PRO movie,device=device

if (n_elements(device) eq 0) then device='x'
if (device eq 'z') then begin
  set_plot,'z'
  device,set_resolution=[1120,1000]
  device,set_pixel_depth=24,decomposed=0
endif else begin
  set_plot,'x'
  window,xsize=1120,ysize=1000
  device,pseudo_color=8,decompose=0
endelse
print,'Windows done'

 restore,/verb,'byt_h1_415-1926b.save'
 restore,/verb,'vxt_h1_415-1926b.save'
 restore,/verb,'vzt_h1_415-1926b.save'
s=size(times)
nt=s[1]
print,nt

!p.background=255
!p.thick=2
!x.thick=2
!y.thick=2
!p.charsize=1.5
!p.charthick=2.0
 
bmin=-1.7
bmax=1.7

for t=1710,nt-1 do begin

  sx=nxt[t]/2
  sy=0
  by1=shift(reform(by01t[0:nxt[t]-1,0:nxt[t]-1,t]),sx,sy)
  vx1=shift(reform(vx01t[0:nxt[t]-1,0:nxt[t]-1,t]),sx,sy)
  vz1=shift(reform(vz01t[0:nxt[t]-1,0:nxt[t]-1,t]),sx,sy)
  dx=48./nxt[t]
  im_vel_scl,by1[.17*nxt[t]:.72*nxt[t],.35*nxt[t]:.9*nxt[t]]<bmax>bmin, $
  vx1[.17*nxt[t]:.72*nxt[t],.35*nxt[t]:.9*nxt[t]],vz1[.17*nxt[t]:.72*nxt[t],.35*nxt[t]:.9*nxt[t]], $
  /horiz,dx=dx,ytitle='!5Y (Mm)',xtitle='!5X (Mm)',colorv=255,colort=2,num=500,len=.5,dim=900
  loadct,0
  colorbar_scl,yr=[bmin,bmax],ytit='!5B!dvert!n (kG)',color=2,ymar=[3,2],x0=0.95,chars=1.5
  xyouts,45,11,'Time'+string(times[t]/.6/60.,format='(F6.2)')+' hrs',/device,color=2 

  im=tvrd(/true)

  if (t le 9) then begin
    write_jpeg,'byvxz_h1.000'+string(t,format='(i1)')+'.jpg',im,/true
  endif else if (t le 99) then begin
    write_jpeg,'byvxz_h1.00'+string(t,format='(i2)')+'.jpg',im,/true
  endif else if (t le 999) then begin
    write_jpeg,'byvxz_h1.0'+string(t,format='(i3)')+'.jpg',im,/true
  endif else begin
    write_jpeg,'byvxz_h1.'+string(t,format='(i4)')+'.jpg',im,/true
  endelse 
  wait,.01
endfor

return
end