The code bellow will consider each individual star so that's a good way of doing it.

Corentin

I've been trying to use a filter with "particle_mass" <  mass resolution cut in between stars and dark matter:

def old(pfilter, data):

        old_only = data[data, "particle_mass"].in_units('Msun')

        filter = old_only < 4e6

        return filter

However, I am uncertain if this is going to consider individual mass particles or combined mass particles in a cell.  

_______________________________________________
yt-users mailing list -- yt-users@python.org
To unsubscribe send an email to yt-users-leave@python.org