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=[1150,1000]
  device,set_pixel_depth=24,decomposed=0
endif else begin
  set_plot,'x'
  window,xsize=1150,ysize=1000
  device,pseudo_color=8,decompose=0
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.background=255
!p.thick=2
!x.thick=2
!y.thick=2
!p.charsize=1.5
!p.charthick=2.0
 
t=1000
intm=aver(intt[0:nxt[t]-1,0:nxt[t]-1,t:t+10])
imin=min(intt[0:nxt[t]-1,0:nxt[t]-1,t:t+10]/intm,max=imax)
print,imin,imax

for t=1710,nt-1 do begin

  sx=nxt[t]/2
  sy=0
  int1=shift(reform(intt[0:nxt[t]-1,0:nxt[t]-1,t]),sx,sy)
  by1=shift(reform(by01t[0:nxt[t]-1,0:nxt[t]-1,t]),sx,sy)
  bh1=shift(reform(bh01t[0:nxt[t]-1,0:nxt[t]-1,t]),sx,sy)
  dx=48./nxt[t]
  image_cont3,int1[.17*nxt[t]:.72*nxt[t],.35*nxt[t]:.9*nxt[t]]/intm<imax>imin,bh1[.17*nxt[t]:.72*nxt[t],.35*nxt[t]:.9*nxt[t]],by1[.17*nxt[t]:.72*nxt[t],.35*nxt[t]:.9*nxt[t]],/aspect, $
  xr=[0,.55*nxt[t]*dx],yr=[0,.55*nxt[t]*dx],clevelsb=[.1,.2],colorsb=[90,120],clevelsc=[-.2,-.1,.1,.2],colorsc=[180,200,220,250],ytitle='!5Y (Mm)',xtitle='!5X (Mm)'
  loadct,0
  colorbar_scl,yr=[imin,imax],ytitle='!5I/<I>',chars=1.5,x0=.94,color=2
  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,'ibb_h1.000'+string(t,format='(i1)')+'.jpg',im,/true
  endif else if (t le 99) then begin
    write_jpeg,'ibb_h1.00'+string(t,format='(i2)')+'.jpg',im,/true
  endif else if (t le 999) then begin
    write_jpeg,'ibb_h1.0'+string(t,format='(i3)')+'.jpg',im,/true
  endif else begin
    write_jpeg,'ibb_h1.'+string(t,format='(i4)')+'.jpg',im,/true
  endelse 
  wait,.01
endfor

return
end