Hi Stephen, I went through my script again, and found indeed an error in the reading procedure. After correction it seems now to work fine and is giving the correct center positions. So in principle one can use the merger-code with enzo-1.0-type output. Thank you again for your effort, I should have checked my script more carefully ... Cheers, Jean-Claude Am Donnerstag, den 06.08.2009, 09:18 -0700 schrieb 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@physics.ucsd.edu o__ Stephen Skory http://physics.ucsd.edu/~sskory/ _.>/ _Graduate Student ________________________________(_)_\(_)_______________ _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org