Thanks Britton, this is exactly what I needed. Cheers, JC On 08/03/12 12:36, Britton Smith wrote:

Hi JC,

You can always reshape the 3D data array to a 1D array. You can try something like numpy.reshape. Also, you should be able to just do grid['Density'] without doing get_data.

Britton

On Thu, Mar 8, 2012 at 10:25 AM, Jean-Claude Passy <jcpassy@gmail.com <mailto:jcpassy@gmail.com>> wrote:

Hi all,

As far as I understand, a:

------------------------------------

region = pf.h.region(...)

density = region["Density"]

------------------------------------

will return the density field without considering cells that have subcells. How can I access ALL the data stored in the hdf5 file? I have tried to use

------------------------------------

grid = pf.h.grids[i]

rhogrid = grid.get_data("Density")

------------------------------------

but I would like to get a 1D array. Thanks for you help,

JC

_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org <mailto:yt-users@lists.spacepope.org> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org