I think you can get the number from the projy_average object with<div>projy_average.field_data['CellMass']</div><div>and double check that against what you calculate</div><div><br></div><div>From</div><div>G.S.<br>
<br><div class="gmail_quote">On Fri, Mar 15, 2013 at 2:26 PM, Elizabeth Tasker <span dir="ltr"><<a href="mailto:tasker@astro1.sci.hokudai.ac.jp" target="_blank">tasker@astro1.sci.hokudai.ac.jp</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">wonderful!<br>
<br>
Thank you!<br>
<span class="HOEnZb"><font color="#888888"><br>
Elizabeth<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
On Mar 15, 2013, at 5:17 PM, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>> wrote:<br>
<br>
> Hi Elizabeth,<br>
><br>
> On Fri, Mar 15, 2013 at 5:15 PM, Elizabeth Tasker<br>
> <<a href="mailto:tasker@astro1.sci.hokudai.ac.jp">tasker@astro1.sci.hokudai.ac.jp</a>> wrote:<br>
>> Hi,<br>
>><br>
>> When you create a projection:<br>
>><br>
>> projy_proj = pf.h.proj(2, 'Density', source=cloud)<br>
>><br>
>> is there a way of knowing the sum of the quantity that has been integrated over? (or equally, knowing the integration depth for each cell, so the total quantity can be calculated)?<br>
>><br>
>> I'd like to sum up the mass above a certain surface density threshold in my source, "cloud". From the above, projy_proj will tell me which regions of my cloud are above a threshold, but I can't think of an easy way to know how much mass has been summed over in that cell.<br>

>><br>
>> projy_average = pf.h.proj(2, "CellMass", weight_field="Ones", source=cloud) will give me an average, but not a total.<br>
><br>
> You could sum up the cell sizes:<br>
><br>
> (proj["pdx"] * 2.0)**2.sum()<br>
><br>
> Note that pdx is half-widths, so we multiply by two.  This will need<br>
> to be unit-converted from code to whatever.<br>
><br>
> -Matt<br>
><br>
>><br>
>> Elizabeth<br>
>> _______________________________________________<br>
>> yt-users mailing list<br>
>> <a href="mailto:yt-users@lists.spacepope.org">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>
> _______________________________________________<br>
> yt-users mailing list<br>
> <a href="mailto:yt-users@lists.spacepope.org">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>
<br>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">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>
</div></div></blockquote></div><br></div>