nimages = 101
xsize = 81
ysize = 40
mag = 10
openr, lun, file_which('jet.dat'), /get_lun
jet = assoc(lun, bytarr(xsize, ysize))
xinteranimate, set=[xsize * mag, ysize * mag, nimages], /showload
for i = 0L, nimages - 1L do $
xinteranimate, frame=i, image=congrid(jet[i], xsize * mag, ysize * mag)
xinteranimate, /keep_pixmaps
free_lun, lun
end