;;rdfits, mos, '../Mosaic/AR8731_Mosaic.fits'

f=findfile('Mosaic_18oct99*')

rdfits, p, f(0), head=info, /noread
h0 = info([0,1,2,3,4,7,8,9,14,15])
add, h0, fitscard('DATE_OBS', "'Mon Oct 18 1999'", '(character) date of data aquisition'), 5
add, h0, fitscard('TIME_BEG', "'08:32:40'", '(char) start of mosaic'), 6
add, h0, fitscard('TIME_END', "'08:41:27'", '(char) end of mosaic'), 7

add, h0, fitscard('HEL_LAT', '13', '(float) Heliographic latitude (N pos.)'), 8
add, h0, fitscard('HEL_LON', '24', '(float) Heliographic longitude (W pos.)'), 9
add, h0, fitscard('SOL_MU', '0.908', '(float) cosine theta'), 10

alp = float(strmid(info(11), 10, 20))
for i=1, n_elements(f)-1 do begin
    rdfits, p, f(i), head=info, /noread
    alp += float(strmid(info(11), 10, 20))
endfor
alp /= n_elements(f)

add, h0, fitscard('ALPHA', string(alp, form="(f4.1)"), '(float) avg. normalized Fried parameter'), 14
file = 'AR8731-gb-0000-19991018-mosaic.fits'
add, h0, fitscard('FILENAME', "'"+file+"'")

wrtfits, mos, head=h0, scale=1e4

end