fg = findfile('../FITS/GBand/AR*')
fb = findfile('../FITS/BlueCont/AR*')
;fc = findfile('../FITS/CaIIH/AR*')

fg = fg(8:*)
fb = fb(8:*)
num = n_elements(fg)

;;; (for i in AR*; do fits $i |grep NAXIS2|cut -c 27-30;done)|sort|m
sx = 1238 & sy = 991
gb = intarr(sx, sy, num)
bc = intarr(sx, sy, num)
;ca = intarr(sx, sy, num)
sall = intarr(2, num)
nsum = 3
sh = intarr(2, num)
scl = 5000

par_bc = [-0.3954, 1.0025]
par_ca = [-0.3111, 1.0031]

window, 1, xs=sx, ys=sy

FOR i=0, num-1 DO BEGIN
      ;;; read G band
    rdfits, p, fg(i)
    p1 = reform(p(*, *, 0))
    sgb = size(p1, /dim)
    gb(0, 0, i) = fix(scl*p1)
      ;;; read blue cont
    rdfits, p, fb(i)
    p1 = rot(reform(p(*, *, 0)), par_bc(0), par_bc(1), cub=-.5)
    sbc = size(p1, /dim)
    bc(0, 0, i) = fix(scl*p1)
;      ;;; read Ca II H
;    rdfits, p, fc(i)
;    p1 = rot(reform(p(*, *, 0)), par_ca(0), par_ca(1), cub=-.5)
;    sca = size(p1, /dim)
;    ca(0, 0, i) = fix(scl*p1)
    
      ;;; now shift all triples to match and cut the size to be the same
    
    ssh = intarr(2, 3)
    ssh(*, 1) = shc(gb(*, *, i), bc(*, *, i), /N2, /filt)
;    ssh(*, 2) = shc(gb(*, *, i), ca(*, *, i), /filt)
    ssh(0, *) = ssh(0, *)-max(ssh(0, *))
    ssh(1, *) = ssh(1, *)-max(ssh(1, *))
    gb(*, *, i) = shift(gb(*, *, i), ssh(0, 0), ssh(1, 0))
    bc(*, *, i) = shift(bc(*, *, i), ssh(0, 1), ssh(1, 1))
;    ca(*, *, i) = shift(ca(*, *, i), ssh(0, 2), ssh(1, 2))
    sall(*, i) = (sgb+ssh(*, 0)) < (sbc+ssh(*, 1))  ;; < (sca+ssh(*, 2))
    
    tvscl, gb(*, *, i) & wait, .5
    tvscl, bc(*, *, i) & wait, .5
;    tvscl, ca(*, *, i) > 0 < 10000 & empty
    
    IF i EQ 0 THEN BEGIN
        pp = fltarr(nsum, sall(0, 0), sall(1, 0))
        pp(0, *, *) = float(gb(0:sall(0, 0)-1, 0:sall(1, 0)-1, 0))+ $
          bc(0:sall(0, 0)-1, 0:sall(1, 0)-1, 0)
    ENDIF ELSE BEGIN
        p1=float(gb(0:sall(0, i)-1, 0:sall(1, i)-1, i))+ $
          bc(0:sall(0, i)-1, 0:sall(1, i)-1, i)
        sh(*, i) = shc(total(pp, 1), p1, /filt, /N2)
        sx = (sall(0, i) < sall(0, 0))
        sy = (sall(1, i) < sall(1, 0))
        p1 = p1(0:sx-1, 0:sy-1)
        pp(i MOD nsum, 0:sx-1, 0:sy-1) = shift(p1, sh(0, i), sh(1, i))
        print, i, sh(*, i)
    ENDELSE

ENDFOR
sx = max(sh(0, *))-sh(0, *)
sy = max(sh(1, *))-sh(1, *)

resx = min(sall(0, *)-sx)
resy = min(sall(1, *)-sy)

  ;;; mave even dimensions
resx = resx/2*2
resy = resy/2*2

stop

FOR i=0, num-1 DO BEGIN
    gb(*, *, i) = shift(gb(*, *, i), -sx(i), -sy(i))
    bc(*, *, i) = shift(bc(*, *, i), -sx(i), -sy(i))
;    ca(*, *, i) = shift(ca(*, *, i), -sx(i), -sy(i))
ENDFOR

stop

gb = (temporary(gb))(0:resx-1, 0:resy-1, *)
bc = (temporary(bc))(0:resx-1, 0:resy-1, *)
;ca = (temporary(ca))(0:resx-1, 0:resy-1, *)

stop

enhance, gb
enhance, bc
;enhance, ca

ang = [-0.3111, 1.0031] # transpose(replicate(1., 194))
fc = findfile('../FITS/CaIIH/AR*')
fc = fc(8:*)
mkmov, ca, fc, roi=[256, 0, 1150, 850], sc=5000., ang=ang
enhance, ca


END