filreg = "FITS/AR*_102[4-9]* FITS/AR*_10[3-5]* FITS/AR*_11* FITS/AR*_12[01]*"

  ;;; align the frames, equalize intensity

mkmov, mov, filreg

  ;;; find rms-values

rms = fltarr((size(mov))(3))
FOR i=0, num-1 DO rms(i) = stdev(mov(*, *, i))

  ;;; Sort out bad ones

ix = where(rms GE 178)

num = n_elements(ix)
FOR i=0, num-1 DO mov(*, *, i) = mov(*, *, ix(i))
mov = (temporary(mov))(*, *, 0:num-1)
f = findfile(filreg)
f = f(ix)


tim = 3600*strmid(f, 20, 2)+60l*strmid(f, 22, 2)+strmid(f, 24, 2)
tim = tim-tim(0)

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

  ;;; equalize power spectra

enhance, mov

  ;;; remove distortions

rubb, mov, av=5, cell=24, step=16, range=4, /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)
    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
openr, 1, 'movie.raw'
writeu, 1, mov
close, 1

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

END