filreg = '../FITS/AR9214*'
f = findfile(filreg)
tim = 3600l*strmid(f, 26, 2)+60l*strmid(f, 28, 2)+strmid(f, 30, 2)
tim0 = tim(0)
tim = tim-tim0

  ;;; align the frames, equalize intensity

mkmov, mov, filreg

  ;;; equalize power spectra

enhance, mov

  ;;; remove distortions

rubb, mov, av=5, cell=24, step=16, range=2, /plot, /same

openw, 1, 'mov_642x464x120.raw'
writeu, 1, mov
close, 1

ssx = 16*(sx/16)
ssy = 16*(sy/16)

mm1 = scale_array(mov(0:ssx-1, 0:ssy-1, *))


tt = tim+tim0

hr = tt/3600
mn = (tt MOD 3600)/60
sc = tt MOD 60
tt = nnumber(hr, 2)+":"+nnumber(mn, 2)+":"+nnumber(sc, 2)

show, mm1(*, *, 0), /cent

FOR i=0, num-1 DO BEGIN
    tv, mm1(*, *, i)
    clock, tt(i), /dev, siz=100, pos=[10, 10]
    mm1(*, *, i) = tvrd()
ENDFOR

write_mpeg, '../MPEG/AR9214_02nov00_sp.mpg', mm1

END