Problems in analyzing dark-matter-only simulations with yt

Hi all,
I have problems in using yt with the dm-only output of enzo-1.5. If I run e.g. the hop recipe, than it fails already at:
full_sphere = pf.h.sphere([0.5,0.5,0.5], 1.0)
Using the same script on dm+hydro output works without any problems, so it seems that yt is confused by the output format. Does anyone know how to use yt with the dm-only output.
Thanks, Jean-Claude

Hi Jean-Claude,
Unfortunately, this is a known bug. I'm working on improving DM-only support, but for now I think it's somewhat usable. I believe your error can be fixed if, during instantiation of the EnzoStaticOutput, you specify the type of output it is:
pf = EnzoStaticOutput(path_to_output, data_style=6)
where data_style is the argument that describes the type of data. 5 is normal HDF5, 6 is the "packed" HDF5. You can tell the difference between the two by how the files in the directory are laid out -- do you have one file per processor (packed HDF5) or do you have one file per grid (normal HDF5)?
-Matt
On Tue, Feb 10, 2009 at 6:43 AM, Jean-Claude Waizmann waizmann@ita.uni-heidelberg.de wrote:
Hi all,
I have problems in using yt with the dm-only output of enzo-1.5. If I run e.g. the hop recipe, than it fails already at:
full_sphere = pf.h.sphere([0.5,0.5,0.5], 1.0)
Using the same script on dm+hydro output works without any problems, so it seems that yt is confused by the output format. Does anyone know how to use yt with the dm-only output.
Thanks, Jean-Claude
--
Jean-Claude Waizmann
Institut für Theoretische Astrophysik (ITA/ZAH) Albert-Ueberle-Strasse 2 69120 Heidelberg, Germany
Phone: +49-6221-54-8987 Fax: +49-6221-54-4221 E-mail: waizmann@ita.uni-heidelberg.de
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

Hi Matt,
since I have the packed HDF5 output (RedshiftOutput00xx.cpu00xx files), I modified the recipe to:
pf=lagos.EnzoStaticOutput("$DIR/RD0011/RedshiftOutput0011",data_style=6)
However, the full_sphere = pf.h.sphere([0.5,0.5,0.5], 1.0) command still fails and I still get errors from the HierarchyType.py module.
Thanks, Jean-Claude
participants (2)
-
Jean-Claude Waizmann
-
Matthew Turk