[yt-users] tracking the masses

[Britton Smith]

Hi Jean-Claude,

You can get the total mass in the grid in the following way:
pf = load(dataset)
whole_box = pf.h.all_data()
total_mass = whole_box.quantities['TotalQuantity']('CellMassMsun')

You can use the EnzoSimulation class to perform this analysis on all of the
data for one simulation.  Check out a description of that here:
http://yt.enzotools.org/doc/extensions/analyzing_an_entire_simulation.html

That will allow you to track the total mass in the box over time.  As far as
your final two points, you might try using the clump finder to locate bound
objects within your simulation.  You can get the total mass in bound clumps
over time and that should take care of those last points.  Information on
the clump finder can be found here:
http://yt.enzotools.org/doc/cookbook/recipes.html#find-clumps

Britton

On Mon, May 10, 2010 at 3:01 PM, Jean-Claude Passy <jcpassy at gmail.com>wrote:

> Dear all,
>
> what I would like to do is tracking bound/lost/unbound quantities. For
> instance, I would like to be able to plot as a function of time:
>   - the total mass in the grid
>   - the mass lost by the grid
>   - the bound mass in the grid
>   - the unbound mass in the grid
>
> Is there a way to get these quantities with the output files only, or do I
> have to track them during the simulation - which will be bad news  :-)  ?
>
> Thanks for your help,
>
>
> Jean-Claude
>
> _______________________________________________
> 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/20100510/6dc47eab/attachment.html>