pro doit

img1 = readfits('20100919.200608.mk4.rpb.vig.fts.gz',hdr1,/noscale)
data = readfits('20100919.200608.mk4.rpb.vig.fts.gz',exthdr,/exten,/noscale)

img2 = readfits('20100919.201200.mk4.rpb.vig.fts.gz',hdr,/noscale)
img3 = readfits('20100919.201752.mk4.rpb.vig.fts.gz',hdr,/noscale)
img4 = readfits('20100919.202344.mk4.rpb.vig.fts.gz',hdr,/noscale)
img5 = readfits('20100919.202936.mk4.rpb.vig.fts.gz',hdr,/noscale)
img6 = readfits('20100919.203532.mk4.rpb.vig.fts.gz',hdr,/noscale)
img7 = readfits('20100919.204123.mk4.rpb.vig.fts.gz',hdr,/noscale)
img8 = readfits('20100919.204715.mk4.rpb.vig.fts.gz',hdr,/noscale)
img9 = readfits('20100919.205307.mk4.rpb.vig.fts.gz',hdr,/noscale)
;img10 = readfits('20100827.014055.mk4.rpb.vig.fts.gz',hdr,/noscale)
;img11 = readfits('20100827.015239.mk4.rpb.vig.fts.gz',hdr,/noscale)
;img12 = readfits('20100827.015538.mk4.rpb.vig.fts.gz',hdr,/noscale)
;img13 = readfits('20100827.020427.mk4.rpb.vig.fts.gz',hdr,/noscale)
;img14 = readfits('20100827.211939.mk4.rpb.vig.fts.gz',hdr,/noscale)
;img15 = readfits('20100827.212235.mk4.rpb.vig.fts.gz',hdr,/noscale)
;img16 = readfits('20100827.212530.mk4.rpb.vig.fts.gz',hdr,/noscale)

;************************
;  FOR rpb.vig images:
;************************
;imgvigavg = fltarr(720,384)

;************************
;  FOR rpb.vig images:
;************************
imgvigavg = fltarr(960,960)

;imgvigavg =  float(img1)+float(img2)+float(img3)+float(img4)+float(img5)+float(img6)+float(img7)+float(img8) + float(img9) + float(img10) + float(img11) + float(img12) + float(img13) + float(img14) + float(img15) + float(img16)

imgvigavg =  float(img1)+float(img2)+float(img3)+float(img4)+float(img5)+float(img6)+float(img7)+float(img8)+float(img9)

;imgvigavg = imgvigavg/16.
imgvigavg = imgvigavg/9.

;************************
;  FOR VIG images:
;************************
imgvigavg = fix(imgvigavg)

;  Modify extension header to capture image times used in average
fxaddpar,exthdr,'AVGTIME1','2006 2012 2017 2023'
fxaddpar,exthdr,'AVGTIME2','2029 2035 2041 2047'
fxaddpar,exthdr,'AVGTIME3','2053               '
fxaddpar,exthdr,'AVGTIME4','                   '


writefits,'20100919.200608.mk4.rpb.avg.vig.fts',imgvigavg,hdr1
writefits,'20100919.200608.mk4.rpb.avg.vig.fts',data,exthdr,/app


; display result: 
window,0,xs=960,ys=960,retain=2
test = readfits('20100919.200608.mk4.rpb.avg.vig.fts',hdr)
tv,bytscl(test,min=-120,max=300)
;tv,bytscl(test,min=-2.e-07,max=6.e-07)

end