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 thanks in advance, Cheers Latif
On Tue, May 14, 2013 at 7:21 AM, Latif <latifne@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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi Matt, many thanks for a prompt response. The output looks like this. print np.unique(dd["particle_type"]) [1 5] print np.unique(dd["creation_time"]).size 2586 It seems particle type is 5 which i assumed to be 2. Cheers Latif On Tue, May 14, 2013 at 1:45 PM, Matthew Turk <matthewturk@gmail.com> wrote:
On Tue, May 14, 2013 at 7:21 AM, Latif <latifne@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@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
participants (2)
-
Latif
-
Matthew Turk