Hi JC,<div><br></div><div>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.</div><div><br>
</div><div>Britton<br><br><div class="gmail_quote">On Thu, Mar 8, 2012 at 10:25 AM, Jean-Claude Passy <span dir="ltr"><<a href="mailto:jcpassy@gmail.com">jcpassy@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

  

    
  
  <div bgcolor="#FFFFFF" text="#000000">
    Hi all,<br>
    <br>
    As far as I understand, a:<br>
    <br>
    <address>------------------------------------</address>
    <address>region = pf.h.region(...)</address>
    <address>density = region["Density"]</address>
    <address>
      ------------------------------------</address>
    <br>
    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 <br>
    <br>
    <address>
      ------------------------------------</address>
    <address>grid = pf.h.grids[i]</address>
    <address>rhogrid = grid.get_data("Density")</address>
    <address>
      ------------------------------------</address>
    <br>
    but I would like to get a 1D array.<br>
    Thanks for you help,<br>
    <br>
    JC<br>
  </div>

<br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>