pro wso_testcase

@compile_opt.pro		; On error, return to caller

;+
; NAME:
;	wso_testcase
; PURPOSE:
;	Runs a test case on XuePu brss_csss_scbf program
; CATEGORY:
;	smei/user/xuepu
; CALLING SEQUENCE:
;	www_testcase
; INPUTS:
;	(from file)
; OUTPUTS:
;	(to fits file)
; PROCEDURE:
;	The files wso_1955_145_1_xuepu.fits was sent to us by
;	Xuepu and is his result from running the program on input
;	file wso1955_145_1.dat (using the last 72 of 73 longitudes
;	in the file, i.e. the first longitude was ignored.
;
;	File wso1955_145_1.fits is the file produced by running our
;	version of the program. Differences between the two should
;	be small (of order 10^-6 or better).
; MODIFICATION HISTORY:
;	NOV-2002, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
;-

root = filepath(root=getenv('SMEI'), subdir=['user','xuepu'],'testcase')

f_in  = filepath(root=root, 'wso1955_145_1.dat')
f_out = filepath(root=root, 'wso1955_145_1.fits')
f_ghf = filepath(root=root, 'wso1955_145_1.txt')


brss_csss_scbf, f_in, 'wso', f_ghf, f_out, rcp=1.6

f_xuepu = filepath(root=root, 'wso1955_145_1_xuepu.fits')

a = readfits(f_xuepu)
b = readfits(f_out  )

whatis, a, b, a-b

return  &  end