Hi Jean-Claude,<br><br>You can get the total mass in the grid in the following way:<br>pf = load(dataset)<br>whole_box = pf.h.all_data()<br>total_mass = whole_box.quantities['TotalQuantity']('CellMassMsun')<br>
<br>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:<br><a href="http://yt.enzotools.org/doc/extensions/analyzing_an_entire_simulation.html">http://yt.enzotools.org/doc/extensions/analyzing_an_entire_simulation.html</a><br>
<br>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:<br>
<a href="http://yt.enzotools.org/doc/cookbook/recipes.html#find-clumps">http://yt.enzotools.org/doc/cookbook/recipes.html#find-clumps</a><br><br>Britton<br><br><div class="gmail_quote">On Mon, May 10, 2010 at 3:01 PM, Jean-Claude Passy <span dir="ltr"><<a href="mailto:jcpassy@gmail.com">jcpassy@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear all,<br>
<br>
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:<br>
   - the total mass in the grid<br>
   - the mass lost by the grid<br>
   - the bound mass in the grid<br>
   - the unbound mass in the grid<br>
<br>
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  :-)  ?<br>
<br>
Thanks for your help,<br>
<br>
<br>
Jean-Claude<br>
<br>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</blockquote></div><br>