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 -- 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