sx = 1040
sy = 864

g = assoc(1, intarr(sx, sy))
k = assoc(2, intarr(sx, sy))

close, /all
openr, 1, 'alig_filt_g_'+strtrim(sx, 2)+'x'+strtrim(sy, 2)+"*"
openr, 2, 'alig_filt_k_'+strtrim(sx, 2)+'x'+strtrim(sy, 2)+"*"

  ;;; first a nice movie with both g and k

num = 198

;newx = sx/32*16
;newy = sy/32*16

newx = 384
newy = 320

mov = intarr(2*newx, newy, num)
FOR i=0, num-1 DO BEGIN
    mov(0:newx-1, *, i) = rescale(g(i), newx, newy)
    mov(newx:*, *, i) = rescale(k(i), newx, newy)
ENDFOR

  ;;; bytscale so they look the same
mm1 = scale_array(mov)

diff = intarr(newx, newy, num)
FOR i=0, num-1 DO diff(*, *, i) = (mov(0:newx-1, *, i)-mov(newx:*, *, i)) > 0
dm = scale_array(diff)

realx = 3*newy
realy = newx  ;/16*16
rm = bytarr(realx, realy, num)
FOR i=0, num-1 DO BEGIN
    rm(0:newy-1, *, i) = rotate(mm1(0:realy-1, *, i), 1)
    rm(newy:2*newy-1, *, i) = rotate(dm(0:realy-1, *, i), 1)
    rm(2*newy:3*newy-1, *, i) = rotate(mm1(newx:newx+realy-1, *, i), 1)
ENDFOR
stop
write_mpeg, 'AR9669_19oct01_a_small.mpg', rm, tmp="/data/pit/tmp"

END