<div dir="ltr"><div>Also, with the idea that "disk" is actually taking a cylinder from -0.385 to 0.385--the reason I wouldn't be getting twice the number of cells is because 1) the tracer fraction cut (specific scalar) and 2) the box only goes to -0.325. <br>
<br></div>-Stephanie<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Aug 4, 2014 at 2:40 PM, Stephanie Tonnesen <span dir="ltr"><<a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>I haven't done the projection, but changed the disk definition to<br><br></div>disk = pf.h.disk([0.0,0.0,0.0][0.0,0.0,1.0],1.2,0.1925)<br>
<br></div>(more or less what Nathan suggested) and end up with basically the same values for the "Ones" count. <br>
<br></div>I thought that in the disk definition the z input was the total height of the cylinder?<span class="HOEnZb"><font color="#888888"><br><br></font></span></div><span class="HOEnZb"><font color="#888888">-Stephanie<br>
</font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Aug 4, 2014 at 2:36 PM, Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">So tracep25 includes about half of the cells that tdp25 does.<div><br></div><div>What if you do this instead:</div>
<div><br></div><div><span style="font-family:arial,sans-serif;font-size:13.333333969116211px">tracerp25 = alld.cut_region(["grid['</span><span style="font-family:arial,sans-serif;font-size:13.333333969116211px">specific_scalar[0]'] > 0.25","grid['z'] < 0.385","grid['z'] > -0.385","grid['cyl_Rcode'] < 1.2"])</span><br>
</div></div><div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Aug 4, 2014 at 11:34 AM, Stephanie Tonnesen <span dir="ltr"><<a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Hi John,<br><br>When I change the tracerp25 selections to all have <= or >= symbols, the value basically doesn't change. With the gt/lt and equal-to "Ones" is definitely different:<br>
<br></div>tracerp25 is 1099285<br></div>tdp25 is 2028955<br><br></div>So the selection criteria is definitely more stringent on <br>tracerp25 = alld.cut_region(["grid['<div>specific_scalar[0]'] > 0.25","grid['z'] <= 0.1925","grid['z'] >= -0.1925","grid['cyl_Rcode'] <= 1.2"])<span><font color="#888888"><br>
<br></font></span></div><span><font color="#888888"><div>-Stephanie<br></div><br></font></span></div><div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Aug 4, 2014 at 2:11 PM, John ZuHone <span dir="ltr"><<a href="mailto:jzuhone@gmail.com" target="_blank">jzuhone@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Stephanie,<br>
<br>
I suspect the problem is that these are not *exactly* the same regions. Have you tried checking the result of:<br>
<br>
tdp25["Ones"].sum()<br>
<br>
tracerp25["Ones"].sum()<br>
<br>
which would be a quick way to determine the number of cells in each region.<br>
<br>
What might be going on (but I have yet to verify) is that the disk selection container object is using >= and <= for the zone selection, whereas you used > and < in the definition of your cut region. You might try swapping those out and seeing if that changes things.<br>
<br>
Best,<br>
<br>
John Z<br>
<div><div><br>
On Aug 4, 2014, at 1:15 PM, Stephanie Tonnesen <<a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a>> wrote:<br>
<br>
> Hi all,<br>
><br>
> I am trying to get the total mass in a disk region of my simulation. I am using athena data and yt 2.7-dev. I thought I was defining the region in the same way using two methods:<br>
><br>
> alld = pf.h.all_data()<br>
> tracerp25 = alld.cut_region(["grid['specific_scalar[0]'] > 0.25","grid['z'] < 0.1925","grid['z'] > -0.1925","grid['cyl_Rcode'] < 1.2"])<br>
><br>
> disk = pf.h.disk([0.0,0.0,0.0],[0.0,0.0,1.0],1.2,0.385)<br>
> tdp25 = disk.cut_region(["grid['specific_scalar[0]'] > 0.25"])<br>
><br>
> but when I print the total mass I get answers that differ by 10%<br>
><br>
> print tracerp25.quantities["TotalQuantity"]("CellMassMsun")<br>
> print tdp25.quantities["TotalQuantity"]("CellMassMsun")<br>
><br>
> >><a href="tel:%5B8019047231" value="+18019047231" target="_blank">[8019047231</a>.5192223]<br>
> >><a href="tel:%5B8886968784" value="+18886968784" target="_blank">[8886968784</a>.4710732]<br>
><br>
> Can someone point out where my difference is?<br>
><br>
> Thanks!<br>
><br>
> Stephanie<br>
</div></div>> _______________________________________________<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>
<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></div>
</div></div><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>
<br></blockquote></div><br></div>
</div></div><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>
<br></blockquote></div><br></div>
</div></div></blockquote></div><br></div>