[yt-users] Halo positions in merger trees
Stephen Skory
stephenskory at yahoo.com
Thu Aug 6 09:18:45 PDT 2009
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
________________________________(_)_\(_)_______________
More information about the yt-users
mailing list