[yt-users] star particles count

[Matthew Turk]

On Tue, May 14, 2013 at 7:21 AM, Latif <latifne at gmail.com> wrote:
> Hi all,
> I have star particles in enzo simulations.  Counting the particles by
> creation time and type gave different results. Is it a problem with Yt or
> comes from simulation? Any ideas?
>
> star_particles = dd["creation_time"] > 0.0
> print "Number of star particles by creation time", star_particles.sum()
> star_particles = dd["particle_type"] == 2
> print "Number of star particles by type", star_particles.sum()
>
>
> Number of star particles by creation time 3818
> Number of star particles by type 0

What's the output of

print np.unique(dd["particle_type"])
print np.unique(dd["creation_time"]).size

?

-Matt

>
> thanks in advance,
> Cheers
> Latif
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>