[yt-users] Gravitational Boundedness with Gas and Particles

[Nathan Goldbaum]

Hi Jason,

If you take a look at the validator in the clump finder that is using the
cython `gravitational_binding_energy` function, you'll see that there is
already an option to use particles in the gravitational binding energy
calculation. See the `_gravitationally_bound` function in
`yt/analysis_modules/level_sets/clump_validators.py`.

When you create your master clump, all you need to do is set
`use_particles=True` when you add the `gravitationally_bound` validator to
the master clump. In fact, the example in the yt docs does this:

http://yt-project.org/docs/dev/analyzing/analysis_modules/clump_finding.html#clump-validators

Hope that helps, let us know if you need any additional advice.

-Nathan

On Thu, Nov 17, 2016 at 1:50 PM, Jason Galyardt <jason.galyardt at gmail.com>
wrote:

> Hi all,
>
> I'd like to use the yt clump finder to identify gravitationally bound
> clumps. However, my FLASH simulation includes both gas (on the mesh) and
> massive particles. I would like to use both of these gravitational sources
> in the determination of boundedness, but it's not clear to me that the
> current implementation of gravitational_binding_energy() in
> yt/utilities/lib/misc_utilities.pyx does this. I suspect that it would be
> relatively straightforward to add the ('deposit', 'io_mass') field to the
> gravitational_binding_energy() function, but I'm not entirely sure how to
> access this field in a safe manner in Cython.
>
> Has anyone already solved this problem?
>
> Thanks,
> Jason
>
> -----
> Jason Galyardt
> University of Georgia
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20161117/52dcd869/attachment-0001.html>