pro doit

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

img2 = readfits('20100827.005349.mk4.rpb.vig.fts.gz',hdr,/noscale)
img3 = readfits('20100827.005645.mk4.rpb.vig.fts.gz',hdr,/noscale)
img4 = readfits('20100827.010237.mk4.rpb.vig.fts.gz',hdr,/noscale)
img5 = readfits('20100827.010829.mk4.rpb.vig.fts.gz',hdr,/noscale)
img6 = readfits('20100827.011426.mk4.rpb.vig.fts.gz',hdr,/noscale)
img7 = readfits('20100827.011722.mk4.rpb.vig.fts.gz',hdr,/noscale)
img8 = readfits('20100827.013205.mk4.rpb.vig.fts.gz',hdr,/noscale)
img9 = readfits('20100827.013503.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) + float(img10) + float(img11) + float(img12) + float(img13)

;imgvigavg = imgvigavg/16.
imgvigavg = imgvigavg/13.

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

;  Modify extension header to capture image times used in average
fxaddpar,exthdr,'AVGTIME1','0047 0053 0056 0102'
fxaddpar,exthdr,'AVGTIME2','0108 0114 0117 0132'
fxaddpar,exthdr,'AVGTIME3','0135 0140 0152 0155'
fxaddpar,exthdr,'AVGTIME4','0204               '


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


; display result: 
window,0,xs=960,ys=960,retain=2
test = readfits('20100827.004757.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