Hi all,

I've been using yt to handle FLASH AMR data stored in HDF5 files. I want to cycle through a series of images generated, plotting each as I go using matplotlib. To do this I have (somewhat simplified):

for f in files:

ds = yt.load(f)

slc = ds.slice(2,0) # its 2D r-z data

frb = slc.to_frb(width=(Rmax,'pc'), height=(Zmax, 'pc'), resolution=(NX,NY))

image = np.array(frb['temp'])

imshow(img)

etc.

This works the way I want it to, but is very slow. Is there a way to speed it up?

Thanks,

Jon

________________________________________________________
Jonathan D. Slavin                 Harvard-Smithsonian CfA
jslavin@cfa.harvard.edu       60 Garden Street, MS 83
phone: (617) 496-7981       Cambridge, MA 02138-1516
cell: (781) 363-0035 USA
________________________________________________________