pro movie

xdi = 200 & ydi = 50

xdi1 = 180 & ydi1 = 30

xdi2 = 2*xdi1 & ydi2 = 2*ydi1

n = intarr(1,/nozero) & k = n 

rho = fltarr(xdi, ydi,/nozero) & velx = rho & vely = rho & p = rho

gamma = fltarr(1)

rrho = fltarr(xdi2, ydi2,/nozero) & rvelx = rrho & rvely = rrho & rp = rrho

; Root of the output files

root = 'hzabu.'

; Number of frames

n = 2

loadct,15
inp_freq=1

set_plot,'x'
device,decomposed=0

; Main loop

FOR k=0,n-1 DO BEGIN
  suffix=strcompress(string(k+1000),/remove_all)
  suffix=strmid(suffix,0,4)
 
  filename=root+suffix

  openr,1,filename
  readu,1,rho,velx,vely,p,gamma
  close,1

  rrho   = congrid(rho(1:xdi1-1,1:ydi1-1),xdi2,ydi2)
  rvelx  = congrid(velx(1:xdi1-1,1:ydi1-1),xdi2,ydi2)
  rvely  = congrid(vely(1:xdi1-1,1:ydi1-1),xdi2,ydi2)
  rp     = congrid(p(0:xdi1-1,0:ydi1-1),xdi2,ydi2)

  rrho = reverse(rrho,2)

  tvscl,rrho

  ext=strcompress(string(k+1001),/remove_all)
  ext=strmid(ext,1,3)

  rooti=root+ext

  print,rooti+'.png'

; To generate a .png file
  tgi=bytarr(xdi2,ydi2) & tgi=tvrd(/true)
  write_png,rooti+'.png',tgi

ENDFOR

end