er change projy_average to projy_proj
From G.S.
On Fri, Mar 15, 2013 at 2:34 PM, Geoffrey So gsiisg@gmail.com wrote:
I think you can get the number from the projy_average object with projy_average.field_data['CellMass'] and double check that against what you calculate
From G.S.
On Fri, Mar 15, 2013 at 2:26 PM, Elizabeth Tasker < tasker@astro1.sci.hokudai.ac.jp> wrote:
wonderful!
Thank you!
Elizabeth
On Mar 15, 2013, at 5:17 PM, Matthew Turk matthewturk@gmail.com wrote:
Hi Elizabeth,
On Fri, Mar 15, 2013 at 5:15 PM, Elizabeth Tasker tasker@astro1.sci.hokudai.ac.jp wrote:
Hi,
When you create a projection:
projy_proj = pf.h.proj(2, 'Density', source=cloud)
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)?
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.
projy_average = pf.h.proj(2, "CellMass", weight_field="Ones",
source=cloud) will give me an average, but not a total.
You could sum up the cell sizes:
(proj["pdx"] * 2.0)**2.sum()
Note that pdx is half-widths, so we multiply by two. This will need to be unit-converted from code to whatever.
-Matt
Elizabeth _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org