filreg='../FITS/AR9214_ca_01nov00_10[34]* ../FITS/AR9214_ca_01nov00_105[0-6]*'
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

  ;;; one is missing

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

tmp = mov
n_neu = num+1
mov = intarr(sx, sy, n_neu)
mov(*, *, 0:52) = tmp(*, *, 0:52)
mov(*, *, 54:*) = tmp(*, *, 53:*)
mov(*, *, 53) = total(tmp(*, *, 52:53), 3)/2

t_neu = lonarr(n_neu)
t_neu(0:52) = tim(0:52)
t_neu(54:*) = tim(53:*)
t_neu(53) = (tim(52)+tim(53))/2

  ;;; find rms-values

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

  ;;; fine adjustment

ref = total(mov, 3)
sh = fltarr(2, n_neu)

FOR i=0, n_neu-1 DO BEGIN
    sh(*, i) = shc(ref, mov(*, *, i), /n2, /filt, /int)
    print, i, sh(*, i)
ENDFOR
sh(0, *) = sh(0, *)-min(sh(0, *))
sh(1, *) = sh(1, *)-min(sh(1, *))

s1 = max(sh(0, *))
IF ((sx-s1)/2. NE (sx-s1)/2) THEN s1 = s1+1
s2 = max(sh(1, *))
IF ((sy-s2)/2. NE (sy-s2)/2) THEN s2 = s2+1

FOR i=0, n_neu-1 DO BEGIN
    mov(*, *, i) = f_shift(mov(*, *, i), sh(0, i), sh(1, i))
ENDFOR

mov = mov(s1:*, s2:*, *)
s = size(mov)
sx = s(1)
sy = s(2)


  ;;; equalize power spectra

enhance, mov

  ;;; remove distortions

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

  ;;; Write FITS files

head = fitscard('SIMPLE', 'T', '(logical) file is basic FITS format')
add, head, fitscard('BITPIX', 16, '(integer) number of bits per pixel')
add, head, fitscard('NAXIS', 2)
add, head, fitscard('NAXIS1', 0)
add, head, fitscard('NAXIS2', 0)
add, head, fitscard('TIME', 0)
head = [head, 'END                            / end of fits header data']

FOR i=0, n_neu-1 DO BEGIN
    head(5) = fitscard('TIME', 1000*(t_neu(i)+tim0))
    wrtfits, reform(mov(*, *, i)), head=head, $
      file='file'+nnumber(i+1, 2)+'a.fits', /swap
ENDFOR

END