;
;
;

dpath='/usr2/mpenn/02sep23/'

dark0100=swap_endian(rfits(dpath+'dark_0100.d.i.001'))
dark0150=swap_endian(rfits(dpath+'dark_0150.d.i.001'))
dark0200=swap_endian(rfits(dpath+'dark_0200.d.i.001'))
dark0800=swap_endian(rfits(dpath+'dark_0800.d.i.001'))
dark1000=swap_endian(rfits(dpath+'dark_1000.d.i.001'))
dark1200=swap_endian(rfits(dpath+'dark_1200.d.i.001'))
dark1400=swap_endian(rfits(dpath+'dark_1400.d.i.001'))
dark1600=swap_endian(rfits(dpath+'dark_1600.d.i.001'))
dark2100=swap_endian(rfits(dpath+'dark_2100.d.i.001'))
dark2300=swap_endian(rfits(dpath+'dark_2300.d.i.001'))
dark2500=swap_endian(rfits(dpath+'dark_2500.d.i.001'))

;!p.multi=[0,1,3]

x000=[100,150,200]
y000=[mean(swap_endian(rfits(dpath+'disk_100.o.i.001'))-dark0100),$
      mean(swap_endian(rfits(dpath+'disk_150.o.i.001'))-dark0150),$
      mean(swap_endian(rfits(dpath+'disk_200.o.i.001'))-dark0200)]
c000=poly_fit(x000,y000,1,yfit000)

x080=[0800,1000,1200]
y080=[mean(swap_endian(rfits(dpath+'sky_080_0800.o.i.001'))-dark0800),$
      mean(swap_endian(rfits(dpath+'sky_080_1000.o.i.001'))-dark1000),$
      mean(swap_endian(rfits(dpath+'sky_080_1200.o.i.001'))-dark1200)]
c080=poly_fit(x080,y080,1,yfit080)

;plot,x080,y080,psym=2,/ynoz
;oplot,x080,yfit080

x160=[0800,1000,1200]
y160=[mean(swap_endian(rfits(dpath+'sky_160_0800.o.i.001'))-dark0800),$
      mean(swap_endian(rfits(dpath+'sky_160_1000.o.i.001'))-dark1000),$
      mean(swap_endian(rfits(dpath+'sky_160_1200.o.i.001'))-dark1200)]
c160=poly_fit(x160,y160,1,yfit160)

;plot,x160,y160,psym=2,/ynoz
;oplot,x160,yfit160

x240=[1200,1400,1600,2100,2300,2500]
y240=[mean(swap_endian(rfits(dpath+'sky_240_1200.o.i.001'))-dark1200),$
      mean(swap_endian(rfits(dpath+'sky_240_1400.o.i.001'))-dark1400),$
      mean(swap_endian(rfits(dpath+'sky_240_1600.o.i.001'))-dark1600),$
      mean(swap_endian(rfits(dpath+'sky240_2100.o.i.001'))-dark2100),$
      mean(swap_endian(rfits(dpath+'sky240_2300.o.i.001'))-dark2300),$
      mean(swap_endian(rfits(dpath+'sky240_2500.o.i.001'))-dark2500)]
c240=poly_fit(x240,y240,1,yfit240)

;plot,x240,y240,psym=2,/ynoz,chars=1.6,xtit='Exposure [msec]',$
    ;ytit='Mean value (frame-dark)'
;oplot,x240,y240,linestyle=1
;oplot,x240,yfit240

;goto,the_end

;!p.multi=0

xsky=[80,160,240]/400.+1.
ysky=[c080(1)/c000(1),c160(1)/c000(1),c240(1)/c000(1)]

plot,xsky,ysky,psym=2,xtit='Radial height',ytit='Normalized Flux',$
     chars=1.6,xran=[1.15,1.65],xstyl=1
oplot,xsky,ysky
oplot,[1.15,1.65],[0,0],linestyl=1


the_end:
end