<div dir="ltr">Hi Stephanie,<div><br></div><div>I'm not sure off-hand why this is happening.</div><div><br></div><div>It may help one of us to be able to reproduce and understand what you're seeing if you can phrase your question in the form of an example that one of can run locally. This would make use of one of the test datasets on <a href="http://yt-project.org/data">yt-project.org/data</a>, and a full script that we can run locally pasted to <a href="http://paste.yt-project.org/data">paste.yt-project.org/data</a></div><div><br></div><div>Finally, just to clarify, the confusion is that the volumes of P3 and P4 are *bigger* than P2 and identical?</div><div><br></div><div>It might be the case that the derived quantity infrastructure is somehow selecting the data as if it were selecting from ds.all_data() instead of your disk selector. Can you try computing what the total volume of cells in a data object constructed by applying cut_regions on ds.all_data() instead of ds.disk()?</div><div><br></div><div>-Nathan</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 7, 2015 at 5:21 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><div><div><div><div><div>Dear yt-users,<br><br></div>I am running<br><br>Version = 3.2-dev<br>Changeset = dc2467c4eae7 (yt) tip<br><br></div>I have been fiddling around with cut_region in a disk region and have gotten really confused. I was initially trying to see whether defining a disk region such that the cylinder height was on the x-axis changed how I need to make a cut in the cylinder height. Anyway, as I was fiddling with these cuts I found that if I made a "spatial" cut (rather than a different variable cut) I was getting different volumes. Let me show you:<br><br>small bit of background: the whole box is 130 kpc across, so 0.5 0.5 0.5 is 65 65 65 kpc. Also, I am cutting and pasting bits here so the indentation is off. <br><br>ds = yt.load("/data/stonnes/dwarf_nowind_lowmassgal_bigtracer/DD"+loop[i]+"/DD44sMBhn25now"+loop[i])<br><br>r2 = ds.disk([0.5,0.5,0.5],[-1,0,0],(4,"kpc"),(2,"kpc"))<br> r2c = r2.cut_region(["obj['Metal_Density']/obj['density'] > 0.33"])<br><br> r3 = ds.disk([0.5,0.5,0.5],[-1,0,0],(4,"kpc"),(2,"kpc"))<br> r3r = r3.cut_region(["(obj['x'].in_units('kpc') > 63.0)"])<br> r3r = r3r.cut_region(["(obj['x'].in_units('kpc') < 67.0)"])<br> r3c = r3r.cut_region(["obj['Metal_Density']/obj['density'] > 0.33"])<br> <br> r4 = ds.disk([0.5,0.5,0.5],[-1,0,0],(4,"kpc"),(2,"kpc"))<br> r4r = r4.cut_region(["(obj['x'].in_units('kpc') > 61.)"])<br> r4c = r4r.cut_region(["obj['Metal_Density']/obj['density'] > 0.33"])<br><br> P2 = np.append(P2,r2c.quantities.total_quantity([("index","cell_volume")]))<br> P3 = np.append(P3,r3c.quantities.total_quantity([("index","cell_volume")]))<br> P4 = np.append(P4,r4c.quantities.total_quantity([("index","cell_volume")]))<br><br></div>P2 volume is 9.346e-5 : within 2% of my calculation for the volume I should be getting<br></div>P3 volume is 7.4768e-4<br></div>P4 volume is 7.4768e-4<br><br></div>When I print out the x- y- and z- maxes and mins, they are the same for each of the selected regions. <br>(By doing this:<br>m2 = np.append(m2,r2c.quantities.extrema([("index","z")])[0])<br>rho2 = np.append(rho2,r2c.quantities.extrema([("index","z")])[1])<br><br>)<br><br></div>Can anyone help me figure out/tell me what is going on here?<br><br></div>Thanks!<br></div>Stephanie<br><br clear="all"><div><div><div><div><div><div><div><div><div><div><div><div><div><div><div dir="ltr"><div><div><div><font face="garamond,serif">--<br>Dr. Stephanie Tonnesen<br></font></div><font face="garamond,serif">Alvin E. Nashman Postdoctoral Fellow<br></font></div><font face="garamond,serif">Carnegie Observatories, Pasadena, CA<br></font></div><font face="garamond,serif"><a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a><br></font></div></div></div>
</div></div></div></div></div></div></div></div></div></div></div></div></div>
<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" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>