[yt-users] faster method for getting frbs?

[Slavin, Jonathan]

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 at cfa.harvard.edu       60 Garden Street, MS 83
phone: (617) 496-7981       Cambridge, MA 02138-1516
cell: (781) 363-0035             USA
________________________________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20160203/4537583b/attachment.htm>