PRO Rubb, mov, AVERAGE=n_av, CELLSIZE=cs, STEP=sw, RANGE=range, $
          PLOT=plot, SAME=same
  ;;; do rubbersheeting of the movie to remove last jitter


IF NOT keyword_set(n_av) THEN n_av = 3
n_s = -(n_av-1)/2
n_e = n_s+n_av-1

IF NOT keyword_set(cs) THEN cs = 32
IF NOT keyword_set(sw) THEN sw = 16

s = size(mov)
num = s(3)

FOR img=0, num-1 DO BEGIN
    
    IF img LT -n_s THEN BEGIN
        ref = total(mov(*, *, 0:n_av-1), 3)
    ENDIF ELSE BEGIN
        IF img GT num-1-n_e THEN $
          ref=total(mov(*, *, num-n_av:num-1), 3) $
        ELSE $
          ref=total(mov(*, *, img+n_s:img+n_e), 3)
    ENDELSE
    
    mov(*, *, img) = map_image(ref, mov(*, *, img), cs, sw, $
                               RANGE=range, PLOT=plot, SAME=same)
    
ENDFOR

END