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


itime=[0.3,1.0,3.0,10.0,30.0,100.0,300.0]
img_cnt1=243
img_cnt2=5

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

indx=0
for i=0,4 do begin
 for j=0,9 do begin
  fname='noise.'+string(img_cnt2+i,format='(i1.1)')+'-'+$
        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
  

imgmean=fltarr(5)
imgstdv=fltarr(5)
for i=0,4 do begin
 imgmean(i)=mean(img_mean(*,*,i))
 imgstdv(i)=mean(img_stdv(*,*,i))
endfor 

plot1:

;window,1,xs=300,ys=300,retain=2
set_plot,'ps'
device,/encapsulated,filename='noise_plot.eps'
plot,imgmean/1.e4,imgstdv,psym=2,xtit='Image mean [1.e4 ADU]',$
     ytit='Image Stdev [ADU]',chars=1.6,tit='noise files',/ynoz
device,/close
set_plot,'x'


the_end:
end