pro offband,ohtrans,od
!p.multi=[0,1,2]
!p.charsize=0.8

!p.psym=10
!x.style=1
!y.range=[0,2000]
!x.range=[3000,7000]
!x.title='Wavelength (A)'
!y.title='Brightness (Rayleighs/A)'
!p.title='Hartley 2 FOS composite'

oh=string(format='(f5.2)',ohtrans)
odep=string(format='(f5.2)',od)

file1='bq02t'
groups_np,file1,w1,ray1
plot,w1,ray1
file2='b303t'
groups_np,file2,w2,ray2
oplot,w2,ray2*2.
file3='b502t'
groups_np,file3,w3,ray3
oplot,w3,ray3

index1=where(w1 ge 3000)
len1=n_elements(index1)
r1=ray1(index1)
ww1=w1(index1)
q1=0.2
t1=ohtrans*ww1/ww1
for i=0,len1-1 do begin
	if(ww1(i) lt 3065 or ww1(i) gt 3125) then t1(i)=t1(i)/1000.
	r1(i)=r1(i)*q1*t1(i)
	endfor
diff1=0.1*(ww1(10)-ww1(0))
sum1=total(r1)*diff1
print,diff1,sum1

index2=where(w2 ge 3300 and w2 le 4600)
len2=n_elements(index2)
r2=ray2(index2)*2.
ww2=w2(index2)
q2=ww2/1.e4
t2=10^(-od)
for i=0,len2-1 do r2(i)=r2(i)*t2*q2(i)
diff2=-0.1*(ww2(10)-ww2(0))
sum2=total(r2)*diff2
print,diff2,sum2

index3=where(w3 ge 4600)
len3=n_elements(index3)
r3=ray3(index3)
ww3=w3(index3)
q3=ww3/1.e4
index6=where(ww3 gt 6000)
q3(index6)=0.6
t3=10^(-od)
for i=0,len3-1 do r3(i)=r3(i)*t3*q3(i)
diff3=-0.1*(ww3(10)-ww3(0))
sum3=total(r3)*diff3
print,diff3,sum3

rmax=max(r1)
!y.range=[1.e-6,1.]
!y.title='Relative signal'
!y.style=1
!p.title='CONTOUR filter estimate'
plot_io,ww1,r1/rmax
;plot,ww1,r1,/ylog
oplot,ww2,r2/rmax
oplot,ww3,r3/rmax
xyouts,5000.,0.1,'T(309 nm) ='+oh
xyouts,5000.,0.02,'Out of band optical depth ='+odep
datetag
end