Hi Stephen, I guess it has to read in all the particles before it can decide which ones are stars, forgot about that part. I'll give it a try and will let you know how it went, thanks.
sfr = StarFormationRate(pf, star_mass=sm, star_creation_time=ct, volume=dd.volume('mpc'))
You might need to make sure that 'sm' and 'ct' are arrays when you pass them in. I don't remember, so I'd just do it to be careful!
-- Stephen Skory firstname.lastname@example.org http://stephenskory.com/ 510.621.3687 (google voice)
yt-users mailing list email@example.com http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org