;
;
;
goto,point1

dpath='/usr2/mpenn/02sep19/'

dark1=fltarr(256,256,7)
dark2=fltarr(256,256,7)
dark3=fltarr(256,256,7)
dark4=fltarr(256,256,7)

expos=[50,100,150,200,250,300,350]

for i=0,6 do begin
 temp=swap_endian(rfits(dpath+'dark400_'$
                        +string(expos(i),format='(i3.3)')+'.d.i.001'))
 dark1(*,*,i)=temp
 temp=swap_endian(rfits(dpath+'dark400_'$
                        +string(expos(i),format='(i3.3)')+'.o.i.001'))
 dark2(*,*,i)=temp
 temp=swap_endian(rfits(dpath+'dark400b_'$
                        +string(expos(i),format='(i3.3)')+'.d.i.001'))
 dark3(*,*,i)=temp
 temp=swap_endian(rfits(dpath+'dark400b_'$
                        +string(expos(i),format='(i3.3)')+'.o.i.001'))
 dark4(*,*,i)=temp
endfor

darkmean1=fltarr(7)
darkmean2=fltarr(7)
darkmean3=fltarr(7)
darkmean4=fltarr(7)

for i=0,6 do begin
 darkmean1(i)=mean(dark1(*,*,i))
 darkmean2(i)=mean(dark2(*,*,i))
 darkmean3(i)=mean(dark3(*,*,i))
 darkmean4(i)=mean(dark4(*,*,i))
endfor


;plot,expos,darkmean2-darkmean1,psym=1,yran=[-0.5,0.5],$
     ;xtit='Exposure [msec]',ytit='Intensity Diff [ADU]',chars=1.6
;oplot,expos,darkmean2-darkmean1
;oplot,expos,darkmean3-darkmean1,psym=2
;oplot,expos,darkmean3-darkmean1
;oplot,expos,darkmean4-darkmean1,psym=4
;oplot,expos,darkmean4-darkmean1

;cdark=poly_fit(expos,darkmean1,1,ydark)
;plot,expos,darkmean1,psym=2,$
     ;xtit='Exposure [msec]',ytit='Dark counts [ADU]',chars=1.6
;oplot,expos,ydark
;print,cdark

; there is a fit coefficient mistake, not a non-linearity in darks
;cdark=poly_fit(expos(0:2),darkmean1(0:2),1,ydark)
;yy=poly(expos,cdark)
;plot,expos,(darkmean1-yy)/darkmean1

point1:

dpath='/usr2/mpenn/02sep19/'
dark1=fltarr(256,256,7)
dark2=fltarr(256,256,7)
expos=[50,100,150,200,250,300,350]

for i=0,6 do begin
 temp=swap_endian(rfits(dpath+'disk400_'$
                        +string(expos(i),format='(i3.3)')+'.d.i.001'))
 dark1(*,*,i)=temp
 temp=swap_endian(rfits(dpath+'disk400_'$
                        +string(expos(i),format='(i3.3)')+'.o.i.001'))
 dark2(*,*,i)=temp
endfor

darkmean1=fltarr(7)
darkmean2=fltarr(7)
for i=0,6 do begin
 darkmean1(i)=mean(dark1(*,*,i))
 darkmean2(i)=mean(dark2(*,*,i)-dark1(*,*,i))
endfor

;cdark1=poly_fit(expos,darkmean1,1)
;print,cdark1
;cdark2=poly_fit(expos,darkmean2,1,darkfit2)
;print,cdark2
;plot,expos,darkmean2,psym=2,xran=[-20,370],xstyl=1,yran=[0,2400],ystyl=1,$
;     xtit='Exposure [msec]',ytit='Intensity',chars=1.6
;oplot,expos,darkfit2
;xx=findgen(390)-20
;yy=poly(xx,cdark2)
;oplot,xx,yy,linestyle=1
;oplot,[0,0],[0,2400]

cdark2=poly_fit(expos(0:2),darkmean2(0:2),1,darkfit2)
yy=poly(expos,cdark2)
plot,expos,(darkmean2-yy)/darkmean2,psym=2,$
     xtit='Exposure [msec]',ytit='(Sig - Fit) / Sig',chars=1.6
oplot,expos,(darkmean2-yy)/darkmean2,linestyl=1



the_end:
end