;
; look at the test.fit.18-393-q1.fit 
; through test.fit.24-419-q1.fit
;

goto,plot1

itime=[0.3,1.0,3.0,10.0,30.0,100.0,300.0]
;img_cnt1=393
;img_cnt2=18
img_cnt1=420
img_cnt2=25

img_mean=fltarr(256,256,7)
img_stdv=fltarr(256,256,5)
img_bffr=fltarr(256,256,5)

indx=0
for i=0,4 do begin
 for j=0,4 do begin
;  fname='test.fit.'+string(img_cnt2+i,format='(i2.2)')+'-'+$
;        string(img_cnt1+indx,format='(i3.3)')+'-q1.fit'
  fname='test2.fit.'+string(img_cnt2+i,format='(i2.2)')+'-'+$
        string(img_cnt1+indx,format='(i3.3)')+'-q1.fit'
  print,'reading file: ',fname
  temp=swap_endian(rfits(fname))
;  tvhist,temp
;  wait,2
  img_bffr(*,*,j)=temp
  indx=indx+1
 endfor
 for ii=0,255 do begin
  for jj=0,255 do begin
   img_mean(ii,jj,i)=mean(img_bffr(ii,jj,*))
   img_stdv(ii,jj,i)=stdev(img_bffr(ii,jj,*))
  endfor
 endfor
endfor
  
; temp=swap_endian(rfits('test.fit.23-418-q1.fit'))
 temp=swap_endian(rfits('test2.fit.30-445-q1.fit'))
 img_mean(*,*,5)=temp
; temp=swap_endian(rfits('test.fit.24-419-q1.fit'))
 temp=swap_endian(rfits('test2.fit.31-446-q1.fit'))
 img_mean(*,*,6)=temp

imgmean=fltarr(7)
imgstdv=fltarr(7)
for i=0,6 do begin
 imgmean(i)=mean(img_mean(*,*,i))
 if(i le 4) then imgstdv(i)=mean(img_stdv(*,*,i))
endfor 

plot1:

;window,0,xs=300,ys=300,retain=2
set_plot,'ps'
;device,/encapsulated,filename='test_mean.ps'
device,/encapsulated,filename='test2_mean.ps'
plot,itime,imgmean,/xlog,xtit='Time [sec]',ytit='Image mean [ADU]',$
;     tit='test data',chars=1.6,psym=2,/ynoz
     tit='test2 data',chars=1.6,psym=2,/ynoz
oplot,itime,imgmean,linestyl=2
device,/close
set_plot,'x'

window,1,xs=300,ys=300,retain=2
;set_plot,'ps'
;device,/encapsulated,filename='test_stdev.ps'
plot,imgmean(0:4),imgstdv(0:4),psym=2,xtit='Image mean [ADU]',$
;     ytit='Image Stdev [ADU]',chars=1.6,tit='test data'
     ytit='Image Stdev [ADU]',chars=1.6,tit='test2 data'
oplot,imgmean(0:4),imgstdv(0:4),linestyl=2
;device,/close
;set_plot,'x'


the_end:
end