filreg = "../FITS/AR*_09085* ../FITS/AR*_0909* ../FITS/AR*_09[1-5]* " + $
  "../FITS/AR*_10[01]* ../FITS/AR*_102[0-5]*"

  ;;; align the frames, equalize intensity

mkmov, mov, filreg

  ;;; find rms-values

num = (size(mov))(3)

f = findfile(filreg)

tim = 3600l*strmid(f, 23, 2)+60l*strmid(f, 25, 2)+strmid(f, 27, 2)
t_null = tim(0)
tim = tim-t_null

n_neu = tim(num-1)/30 + 1
t_neu = indgen(n_neu)*30

  ;;; equalize power spectra

enhance, mov

  ;;; remove distortions

rubb, mov, av=7, cell=24, step=16, range=4, /align, /plot, /same

  ;;; now interpolate to equidistant timesteps

s = size(mov)
sx = s(1)
sy = s(2)

openw, 1, 'mov_interp.raw'
FOR i=0, sx-1 DO BEGIN
    FOR j=0, sy-1 DO BEGIN
        tmp = spline(tim, mov(i, j, *), t_neu)
        writeu, 1, fix(tmp)
    ENDFOR
    writeu, -1, string(i)
ENDFOR
close, 1

delvar, mov
mov = intarr(sx, sy, n_neu)
lin = intarr(n_neu)

openr, 1, 'mov_interp.raw'

FOR i=0, sx-1 DO BEGIN
    FOR j=0, sy-1 DO BEGIN
        readu, 1, lin
        mov(i, j, *) = lin
    ENDFOR
ENDFOR

close, 1
openw, 1, 'mov_interp.raw'
writeu, 1, mov
close, 1
delvar, mov

conefilt, 'mov_interp.raw', sx, sy, n_neu, 420

spawn, 'touch Size_is_'+strtrim(sx, 2)+"x"+strtrim(sy, 2)+"x"+strtrim(n_neu, 2)
mov = intarr(sx, sy, n_neu)

openr, 1, 'sum.dat'
readu, 1, mov
close, 1

mm1 = scale_array(mov, top=2700, bot=200)
 ;;; for mpeg: multiples of 16
mm1 = mm1(4:4+623, 0:431, *)

tt = t_neu + t_null
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)

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

gamma = '1.6'
g = byte(255*(findgen(256)/255)^(1./gamma))

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

mm1 = g(mm1)
write_mpeg, '../MPEG/AR8739_23oct99_sp_g'+gamma+'.mpg', mm1

END