pro offband3
;version for Barr filter 1/3/02
!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'

ohtrans=0.45
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
index1a=where(ww1 le 3200)
sum1a=total(r1(index1a))*diff1
print,diff1,sum1,sum1a

index2=where(w2 ge 3300 and w2 le 4600)
len2=n_elements(index2)
r2=ray2(index2)*2.
ww2=w2(index2)
q2=ww2/1.e4
od=6.0
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
for i=0,len3-1 do begin
	if (ww3(i) le 5000) then od=6.0
	if (ww3(i) gt 5000 and ww3(i) le 6000) then od=3.8
	if (ww3(i) gt 6000) then od=3.5
	t3=10^(-od)
	r3(i)=r3(i)*t3*q3(i)
	endfor
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='Barr filter 1/3/02'
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
oob=100.*(sum1-sum1a+sum2+2*sum3)/sum1a
soob=string(format='(f7.4)',oob)
xyouts,5000.,0.02,'Out of band = '+soob+'%',size=1.0

datetag
end