[yt-users] Halo positions in merger trees

[Jean-Claude Waizmann]

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 at physics.ucsd.edu           o__  Stephen Skory
> http://physics.ucsd.edu/~sskory/ _.>/ _Graduate Student
> ________________________________(_)_\(_)_______________
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org