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 ________________________________________________________