program testmap

	include		'ftspar.h'

	character	cfile*256
	character	cfits*256

	integer		ndim(2)			/0,0/
	real		zfile(73,37)
	real		zfits(73,37)

	cfile = '/big/oftp/dat/map/wso_hcss/wso2010_345_1.dat'
	cfits = '/home/soft/temp/wso2010_345_1.fts.gz'

	nmax = 73*37
	call ArrR4Bad(nmax,zfile)
	call ArrR4Bad(nmax,zfits)

	i = BRead_WSO(cfile, nmax, zfile, nlng, nlat)
	print * , nLng,nLat
	call whatisr4(nlng*nlat,zfile,'zfile')
	
	nmax = 73*37
	nLng = 0
	nLat = 0
	i = BRead_WSO(cfits, nmax, zfits, nlng, nlat)
	print *, nLng,nlat
	call whatisr4(nlng*nlat,zfits,'zfits')

	call ArrR4MinusArrR4(-nLng*nLat,zfile,zfits,zfile)

	call whatisr4(nlng*nlat,zfile,'zfile')
	
	end