<div dir="ltr"><div>Hi all,<br><br></div><div>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.<br><br>Has anyone already solved this problem?<br><br></div><div>Thanks,<br></div><div>Jason <br><br>-----<br></div><div>Jason Galyardt <br></div><div>University of Georgia<br><br></div></div>