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@physics.ucsd.edu o__ Stephen Skory http://physics.ucsd.edu/~sskory/ _.>/ _Graduate Student ________________________________(_)_\(_)_______________