PRO movie,data,times,nxt,device=device,datamin=datamin,datamax=datamax,tmin=tmin

if (n_elements(datamin) eq 0) then datamin=min(data)
if (n_elements(datamax) eq 0) then datamax=max(data)
if (n_elements(tmin) eq 0) then tmin=0

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

s=size(times)
nt=s[1]
print,nt
datab=bytscl(data,min=datamin,max=datamax)

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

  sx=nxt[t]/3
; sx=nxt[t]/2
  sy=nxt[t]/4
; sy=0
  databs=shift(reform(datab[0:nxt[t]-1,0:nxt[t]-1,t]),sx,sy)
  ix=0
  iy=0  
  panel,ix,iy,nx=1,ny=1,xdist=1,ydist=1,xleft=7,ybottom=5,ytop=1,xright=17
; Plot entire data set
  plot_image,databs,[0,48],[0,48],ytitle='!5Y (Mm)',xtitle='!5X (Mm)'
; Plot only subset of data
; plot_image,databs[.17*nxt[t]:.72*nxt[t],.35*nxt[t]:.9*nxt[t]], $
; [0,.55*48],[0,.55*48],ytitle='!5Y (Mm)',xtitle='!5X (Mm)'
  colorbar_scl,yr=[datamin,datamax],ytitle='!5B!dvert!n (kG)',chars=1.5,x0=.94,color=2
; ix=1
; iy=0
; panel,ix,iy,nx=2,ny=1,xdist=1,ydist=1,xleft=7,ybottom=5,ytop=1,xright=10
; plot_image,shift(reform(databs[0:630,2016-630:*,t]),sx,sy),[0,15],[0,15],xtitle='!5X (Mm)'
; colorbar_scl,yr=[by01n,by01x],ymar=[5,1],color=1,x0=.96
  xyouts,45,11,'Time'+string(times[t]/.6/60.,format='(F6.2)')+' hrs',/device,color=2 

; ix=0
; iy=1
; panel,ix,iy,nx=2,ny=2,xdist=1,ydist=1,xleft=7,ybottom=6
; plot_image,by1tb[*,*,t],[0,48],[0,48],xtitle='!5X (Mm)'
; ix=1
; iy=1
; panel,ix,iy,nx=2,ny=2,xdist=1,ydist=1,xleft=7,ybottom=6
; plot_image,imb20[*,*,t],[0,48],[0,48],xtitle='!5X (Mm)'

  im=tvrd()

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

return
end