;
;
;

dpath='/usr2/mpenn/02sep20/'

dark=fltarr(256,256,10)
disk=fltarr(256,256,10)
time=fltarr(10)

for i=0,9 do begin
 temp=swap_endian(rfits(dpath+'disk400_'$
                        +string((i+1)*10,format='(i3.3)')+'.d.i.001'))
 dark(*,*,i)=temp
 temp=swap_endian(rfits(dpath+'disk400_'$
                        +string((i+1)*10,format='(i3.3)')+'.o.i.001'))
 disk(*,*,i)=temp
 time(i)=float((i+1)*10.)
endfor

sun=disk-dark

darkmean=fltarr(10)
sunmean=fltarr(10)
for i=0,9 do begin
 darkmean(i)=mean(dark(*,*,i))
 sunmean(i)=mean(sun(*,*,i))
endfor

cdark=poly_fit(time,darkmean,1,darkfit)
;plot,time,darkmean,psym=2,xtit='Exposure [msec]',ytit='Dark counts',$
;     chars=1.6,/ynoz
;oplot,time,darkfit


;csun=poly_fit(time,sunmean,1,sunfit)
;xx=findgen(120)-20
;yy=poly(xx,csun)
;plot,time,sunmean,psym=2,xtit='Exposure [msec]',ytit='Solar counts',$
;     chars=1.6,xran=[-20,105],xstyl=1,yran=[0,750],ystyl=1
;oplot,time,sunfit
;oplot,xx,yy,linestyle=1
;oplot,[0,0],[0,750]

time=time+14.1
csun=poly_fit(time(0:2),sunmean(0:2),1,sunfit)
yy=poly(time,csun)
plot,time,(sunmean-yy)/sunmean,xtit='Exposure [msec]',$
     ytit='(Sig - Fit) / Sig',psym=2,chars=1.6
oplot,time,(sunmean-yy)/sunmean

;time=time+14.1
;
;logsig=alog(sunmean)/alog(10.)
;logtime=alog(time)/alog(10.)
;
;clog=poly_fit(logtime,logsig,1,yfitlog)
;
;plot,logtime,logsig,psym=2



the_end:
end