-> arr = myGroup.create_dataset(name,data=array) (Pdb) name 'DataFields' (Pdb) array 
Interesting. Any ideas on how to fix it? :(
None now, sorry. I'm focusing on the hg version.
As a sanity check, can you do this with the hierarchy-opt branch? To make sure I didn't screw anything up during the merge.
I get exactly the same issues. Which is a kind of good.
self.max_level is being calculated wrong, and my suspicion is that it's getting the grid count, not the maximum level. Can you manually examine the .harrays file to make sure this is the case? Looking at the short output you have below, looks like it might actually be putting the ParentIDs into the level array.
The maximum value in .harrays["/Level"] (I used h5dump to check, to be sure) is the same as self.max_level, which is as it should be, right? I have self.max_level = 2867, max(["/Level"]) = 2867, the number of grids is 5403, MAXLEVEL is 48. Should self.level_stats have 2867 entries, 5403 or 48?
_______________________________________________________ firstname.lastname@example.org o__ Stephen Skory http://physics.ucsd.edu/%7Esskory/ _.>/ _Graduate Student ________________________________(_)_\(_)_______________