filreg = '../FITS/AR*'

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

  ;;; align the frames, equalize intensity

mkmov, mov, filreg, /clip16

  ;;; equalize power spectra

enhance, mov

  ;;; remove distortions

rubb, mov, av=5, cell=32, step=24, range=3, plot=3, /same

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

FOR i=0, num-1 DO BEGIN
    av = (avg(mov(0:249, 0:299, i)) + avg(mov(840:*, 175:554, i)))/2.
    mov(*, *, i)=1000./av*mov(*, *, i)
ENDFOR

p1 = fix(total(mov(*, *, 44:45), 3)/2.+.5)
p2 = fix(total(mov(*, *, 82:83), 3)/2.+.5)

num = num+2

rawfile = 'mov_'+strtrim(sx, 2)+'x'+strtrim(sy, 2)+'x'+strtrim(num, 2)+'.raw'

mov(*, *, 54) = total(mov(*, *, [53, 55]), 3)/2

openw, 1, rawfile
writeu, 1, mov(*, *, 0:44)
writeu, 1, p1
writeu, 1, mov(*, *, 45:82)
writeu, 1, p2
writeu, 1, mov(*, *, 83:*)
close, 1

tim = [tim(0:44), fix(avg(tim(44:45))), tim(45:82), $
       fix(avg(tim(82:83))), tim(83:*)]

mov = 0b
conefilt_file, rawfile, sx, sy, num, 725*0.071, 30, 4.5, PART=512

mov = intarr(sx, sy, num)
openr, 1, 'sum.dat'
readu, 1, mov
close, 1

mm1 = scale_array(mov)

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/AR9407_01apr01.mpg', mm1, tmp='/data/pit/tmp/mpgtmp/'

mm1 = congrid(mov, sx/2, sy/2, num)
mm1 = scale_array(mm1)

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

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

write_mpeg, '../MPEG/AR9407_01apr01_small.mpg', mm1, $
  tmp='/data/pit/tmp/mpgtmp/'

END