; BBTHROUGH.PRO

PRINT,'Enter the name of the passband directory'
dir=''
READ,dir
pbfile='/jasper/cnsr3_data1/passband/'+dir+'/signal.za'
rdcols,pbfile,n,wavelength,reltrans

wavelength=wavelength*10.	; convert to Angstroms
mxw=MAX(wavelength,MIN=mnw)
uwl=FINDGEN(mxw-mnw+1)+mnw	; 1 A steps from min to max
mxr=MAX(reltrans,MIN=mnr)
reltrans=(reltrans-mnr)/(mxr-mnr)	; scale from 0 to 1
srt=SPLINE(wavelength,reltrans,uwl)	; spline to 1 A interval

PRINT,'Enter the black-body temperature and scale factor'
READ,t,sf

planck,uwl,[t,sf],bb
tp=TOTAL(srt*bb)
PRINT,tp,FORMAT='("Total throughput is ",G10.4," kR, for unit peak transmission")'

END