restore,'8210_working.sav' &  step = 0 
;restore,'/data/welsch/Saves/ar8210_leka_ilct_r6.sav' &  running=6
;restore,'/data/welsch/Saves/ar8210_george_samp7.sav'

nx = n_elements(data.bz(*,0,0))
ny = n_elements(data.bz(0,*,0))

if (n_elements(step) eq 0) then step = 0
if not(keyword_set(running)) then running = 7

bx = data.bx(*,*,running+step)
by = data.by(*,*,running+step)
bz = data.bz(*,*,running+step)
dbzdt = [data.bz(*,*,step+2*running)-data.bz(*,*,step)]

if (max(abs(vel.u)) ne 0) then vx = vel.u(*,*,step) else $
  vx = vel.vx(*,*,step)
if (max(abs(vel.v)) ne 0) then vy = vel.v(*,*,step) else $
  vy = vel.vy(*,*,step)

;for i=0,50 do begin
;    bx = pad_array(bx)
;    by = pad_array(by)
;    bz = pad_array(bz)
;    dbzdt = pad_array(dbzdt)
;    vx = pad_array(vx)
;    vy = pad_array(vy)
;endfor

nx = n_elements(bz(*,0))
ny = n_elements(bz(0,*))

mag  = fltarr(nx,ny,4)
mag(*,*,0) = bx
mag(*,*,1) = by
mag(*,*,2) = bz
mag(*,*,3) = dbzdt

ivel = fltarr(nx,ny,2)
ivel(*,*,0) = vx
ivel(*,*,1) = vy

if not(keyword_set(thr)) then thr = 100
nzdat = where(abs(data.bz(*,*,step+running)) gt thr)

ilct,ivel*1e+5,mag,data.dx*1e+5,data.t(step+2*running)-data.t(step),$
   check=1, induct_check=induct_check, /verbose

end