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=[1000,1075]
  device,set_pixel_depth=24,decomposed=1
endif else begin
  set_plot,'x'
  window,xsize=1000,ysize=1075
  device,decomposed=1
endelse
print,'Windows done'

 restore,/verb,'int_h1_415-1904.save'
 restore,/verb,'byt_h1_415-1904.save'
 restore,/verb,'bht_h1_415-1904.save'
s=size(times)
nt=s[1]
print,nt

!p.thick=2
!x.thick=2
!y.thick=2
!p.charsize=1.5
!p.charthick=2.0
 
for t=1710,nt-1 do begin

  sx=nxt[t]/2
  sy=0
  int1=bytscl(shift(reform(intt[0:nxt[t]-1,0:nxt[t]-1,t]),sx,sy))
  by1=bytscl(shift(reform(by01t[0:nxt[t]-1,0:nxt[t]-1,t]),sx,sy))
  by1p=bytscl(shift(reform(by01t[0:nxt[t]-1,0:nxt[t]-1,t])>0.,sx,sy))
  by1m=bytscl(shift(reform(-by01t[0:nxt[t]-1,0:nxt[t]-1,t])>0.,sx,sy))
  bh1=bytscl(shift(reform(bh01t[0:nxt[t]-1,0:nxt[t]-1,t]),sx,sy))
  dx=48./nxt[t]
  erase
; tv,congrid(int1[.17*nxt[t]:.72*nxt[t],.35*nxt[t]:.9*nxt[t]]*.3,1000,1000),channel=0,0
; tv,congrid(by1[.17*nxt[t]:.72*nxt[t],.35*nxt[t]:.9*nxt[t]],1000,1000),channel=1,0
  tv,congrid(by1p[.17*nxt[t]:.72*nxt[t],.35*nxt[t]:.9*nxt[t]],1000,1000),channel=1,0
  tv,congrid(by1m[.17*nxt[t]:.72*nxt[t],.35*nxt[t]:.9*nxt[t]],1000,1000),channel=3,0
  tv,congrid(bh1[.17*nxt[t]:.72*nxt[t],.35*nxt[t]:.9*nxt[t]],1000,1000),channel=2,0
  plot,[0,.55*nxt[t]*dx],[0,.55*nxt[t]*dx],/noerase,xtit='!5Mm',xmarg=[0,0],ymar=[5.0,0],xst=1,yst=1,/nodata,channel=0
  xyouts,45,15,'Time'+string(times[t]/.6/60.,format='(F6.2)')+' hrs',/device

  im=tvrd(/true)

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

return
end