[yt-users] Halo positions in merger trees

[Stephen Skory]

Jean-Claude,


> Unfortunately, that is terrible news and explains the difference in the
> number of particles that are found. So we will now see what we can do
> about it.

Before I ruin your day completely, there's a simple test to see if this is a problem. Try this:

>>> from yt.mods import *
>>> pf = EnzoStaticOutput('dataset')
>>> sp = pf.h.sphere([0.5]*3, 1.0)
>>> indices = sp["particle_index"]
>>> indices.size
x
>>> uni_indices = na.unique(indices)
>>> uni_indices.size
y

If x == y, things are OK, at least for this dataset.

As I think back to when I discovered this bug, I can't remember if it affected both DM & stars, or just stars. I'd like you to try this so I don't give you wrong information.

If the test works out OK, can you send me the python script?

Sorry for the confusion!

 _______________________________________________________
sskory at physics.ucsd.edu           o__  Stephen Skory
http://physics.ucsd.edu/~sskory/ _.>/ _Graduate Student
________________________________(_)_\(_)_______________