;
;
;

dpath='/usr2/mpenn/02sep17/'

sunflux100=fltarr(256,256)
sundark100=fltarr(256,256)
sunflux200=fltarr(256,256)
sundark200=fltarr(256,256)

for i=0,63 do begin
 temp=swap_endian(rfits(dpath+'suncent02.d.i.'$
                        +string(i+1,format='(i3.3)')))
 sundark200=sundark200+float(temp)/64.

 temp=swap_endian(rfits(dpath+'suncent02.o.i.'$
                        +string(i+1,format='(i3.3)')))
 sunflux200=sunflux200+float(temp)/64.

 temp=swap_endian(rfits(dpath+'suncent02b.d.i.'$
                        +string(i+1,format='(i3.3)')))
 sundark100=sundark100+float(temp)/64.
 temp=swap_endian(rfits(dpath+'suncent02b.o.i.'$
                        +string(i+1,format='(i3.3)')))
 sunflux100=sunflux100+float(temp)/64.
endfor

sflux100=mean(sunflux100-sundark100)/1250.
sflux200=mean(sunflux200-sundark200)/750.

sky1=fltarr(64)
sky2=fltarr(64)
sky3=fltarr(64)

; must estimate dark value since no dark taken at 2500 milliseconds
;darkval=185.+2500.
; but, dark images not taken so bias effects dominate, unless we
; take differences of the images first 

for i=0,63 do begin
 fname1=dpath+'sky1_1rb.o.i.'+string(i+1,format='(i3.3)')
 fname2=dpath+'sky1_1r.o.i.'+string(i+1,format='(i3.3)')
 sky1(i)=mean(swap_endian(rfits(fname2))-swap_endian(rfits(fname1)))

 fname1=dpath+'sky2rb.o.i.'+string(i+1,format='(i3.3)')
 fname2=dpath+'sky2r.o.i.'+string(i+1,format='(i3.3)')
 sky2(i)=mean(swap_endian(rfits(fname2))-swap_endian(rfits(fname1)))

 fname1=dpath+'sky3rb.o.i.'+string(i+1,format='(i3.3)')
 fname2=dpath+'sky3r.o.i.'+string(i+1,format='(i3.3)')
 sky3(i)=mean(swap_endian(rfits(fname2))-swap_endian(rfits(fname1)))

endfor

rad=[1.1,2.0,3.0]
y=[mean(sky1),mean(sky2),mean(sky3)]/sflux100/2500.*1.e6
e=[stdev(sky1),stdev(sky2),stdev(sky3)]/sflux100/2500.*1.e6

plot,rad,y,psym=2,xtit='Radial Height',ytit='Intensity [millionths]',$
     chars=1.6,xran=[1.0,3.2],xstyl=1,yran=[-300,1000],ystyl=1

oplot,rad,y,linestyl=1
errplot,rad,y-e,y+e



the_end:
end