program test

	integer		XCvarFormat
	character	cFile*100
	include		'filparts.h'
	integer		BRead_WSO

	real		ZFile(10000)
	
	cFile = '/big/oftp/dat/map/wso_hcss/wso1999_360_1.dat.gz'

	I = BRead_WSO(cFile, 10000, ZFile,nLng,nLat)
	print *, nLng,nLat, I
	if (1 .eq. 1) stop

	cFile = ''
	NCoff = 0
	XC = 0.00456

	print *, NCoff, XC
	J = XCvarFormat(NCoff,XC,cFile)
	
	print *, NCoff, J
	print *, cFile(:itrim(cFile))


	end

c	program test
c	include		'sun.h'
c	include		'phys.h'
c
c	ScSun	= 1.
c	S	= .5
c	Elo	= 30.
c	iChan	= 1
c
c	F1 = THOM_WTF(ScSun, S, Elo, iChan)
c
c	ElSun = ElSunDistance(ScSun,S,Elo, SinChi)
c	call ThomsonUBVConst(iChan,U,APM)
c
c	ElSun = sngl(ElSun/SUN__RAU)
c	F2 = ThomsonElectronS10(ElSun,SinChi,U,APM, P)
c
c	print *, (SUN__AU*1E13)*1E-26
c	print *, F1, F2*(SUN__AU*1E13)*1E-26
c
c	end
c
c	function THOM_WTF(ScSun,ElSc,Elo,iChan)
c	real*8	S, S2, C2, C, G
c	real	U(3)	/.8,.77,.62/,
c    &		LUM(3)	/2.421029,2.728978,4.920395/,
c    &		RSun	/4.652402E-3/
c
c	ScSun2 = ScSun*ScSun
c	D2 = 2.*ScSun*cosd(Elo)
c	U1 = U(iChan)
c	U2 = 1.-U1
c
c	ElSun2	= ScSun2+ElSc*(ElSc-D2)	! El-Sun dist by cosine rule
c	ElSun	= sqrt(ElSun2)
c	S	= RSun/ElSun		! Sine(Omega)
c	S2	= S*S
c	C2	= 1.-S2
c	C	= dsqrt(C2)		! Cosine(Omega)
c	G	= dlog((1.+S)/C)
C-------
C cCoefficients (from Billings, p. 150)
c
c	AA	= C*S2					! Eq. (20)
c	BB	= -(1.-3.*S2-C2*(1.+3.*S2)*G/S)/8.	! Eq. (21)
c	CC	= 4./3.-C-(C*C2/3.)			! Eq. (22)
c	DD	= (5.+S2-C2*(5.-S2)*G/S)/8.		! Eq. (23)
c
c	DIT	= U2*CC+U1*DD				! Eq. (18)
c	DITMDIR	= (U2*AA+U1*BB)*(ScSun/ElSun*sind(Elo))**2 ! Eq. (19)
c	DIR	= DIT-DITMDIR
c
cc	THOM_WTF = DIT+DIR
c
c	THOM_WTF = 835.9088*LUM(iChan)/(1-U(iChan)/3.)*(DIT+DIR)
c
c	return			! Brightness [S10]
c	end