itime1=[0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 2.0, 3.0,$
      4.0, 5.0, 7.0, 9.0, 10.0, 15.0, 20.0, 30.0, 40.0, 50.0,$
      70.0, 80.0, 90.0, 100.0, 200.0, 300.0, 400.0, 500.0]

flist=file_search('20050306/lin*',count=numfiles)
;window,0,xs=300,ys=300,retain=2

filetimes=fltarr(numfiles)
counts=fltarr(numfiles)
itime=fltarr(numfiles)
data=fltarr(256,256,numfiles)

for i=0,numfiles-1 do begin
 print,flist(i)
 r=fitshline(file=flist(i),error=e,KEY='ITIME')
 itime(i)=float(strmid(r,11,8))
 temp=readfits(flist(i))
 data(*,*,i)=temp
; tvhist,rebin(temp,256,256)
 counts(i)=mean(temp)
endfor

counts2=counts
for i=0,numfiles-1 do counts(i)=counts(i)-counts(1)


; files 27 and 28 are dark with no light to try to understand stripes
c=poly_fit(itime(10:16),counts(10:16),1,yfit)
predicted=poly(itime(0:26),c)


;plot,itime(0:26),counts(0:26),psym=1,/ynoz
;oplot,itime,predicted

;set_plot,'ps'
;device,/encapsulated,filename='lintest_mar05.eps'
;
;plot,counts2/1.e3,(predicted-counts)/predicted,psym=1,yran=[-0.02,.3],$
;     ystyl=1,chars=1.6,xtit='Measured Total Counts',$
;     ytit='(Predicted-Counts)/Predicted',tit='Linearity Tests Mar 05'
;
;device,/close
;set_plot,'x'

the_end:
end