<div>Hi Stephanie,</div><div><br></div><div>Any chance you can demonstrate this using one of the public datasets on <a href="http://yt-project.org/data">yt-project.org/data</a>? Alternatively you could share one of your datasets that you're having trouble with. Either way, can you please open an issue about it on bitbucket so we don't lose track?</div><div><br></div><div>Nathan</div><div><br><div class="gmail_quote"><div>On Fri, Apr 7, 2017 at 10:17 PM Stephanie Tonnesen <<a href="mailto:stonnes@gmail.com">stonnes@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="gmail_msg">Hi all, <div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg">So I think I found a bug in cut_region.  Right now the directions say you can chain cut regions in two ways:  either listing several cut_regions or using one call and an "&".</div><div class="gmail_msg"><a href="http://yt-project.org/doc/analyzing/filtering.html" class="gmail_msg" target="_blank">http://yt-project.org/doc/analyzing/filtering.html</a></div><div class="gmail_msg">However, for some extremely mysterious reason, when I do a multi-call to cut_region I get wackadoo numbers, like too high by a couple orders of magnitude.  </div><div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg">Here is the quick check that should work on IsoGal:</div><div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg"><div class="gmail_msg">ds = yt.load("/data001/stonnes/BPr16vs30arcfo/DD"+loop[i]+"/BPr16vs30arcfo"+loop[i])</div><div class="gmail_msg">    </div><div class="gmail_msg">    disk3 = ds.disk([0.5,0.5,0.5],[0,0,-1.0],(40,"kpc"),(10,"kpc"))</div><div class="gmail_msg">    </div><div class="gmail_msg">    diskn = disk3.cut_region(["obj['Metal_Density']/obj['density'] > 0.5"])</div><div class="gmail_msg">    diska = diskn.cut_region(["obj['density'] > 1e-25"])</div><div class="gmail_msg">    diskd = disk3.cut_region(["(obj['Metal_Density']/obj['density'] > 0.5) & (obj['density'] > 1e-25)"])</div><div class="gmail_msg">    </div><div class="gmail_msg">    print diska.quantities.total_quantity([("gas","cell_mass")])</div><div class="gmail_msg">    print disk3.quantities.total_quantity([("gas","cell_mass")])</div><div class="gmail_msg">    print diskn.quantities.total_quantity([("gas","cell_mass")])</div><div class="gmail_msg">    print diskd.quantities.total_quantity([("gas","cell_mass")])</div></div><div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg">I am also attaching my whole code, which is a loop, just in case.  So here the diska mass is two orders of magnitude larger than any of the others, and diskd is what I would expect with the two selection criteria.</div><div class="gmail_msg"><br class="gmail_msg"></div><div class="gmail_msg">Best,</div><div class="gmail_msg">Stephanie<br clear="all" class="gmail_msg"><div class="gmail_msg"><div class="m_2269913028756309132gmail_signature gmail_msg"><div class="gmail_msg"><div class="gmail_msg"><div class="gmail_msg"><div class="gmail_msg"><font face="garamond,serif" class="gmail_msg">--<br class="gmail_msg">Dr. Stephanie Tonnesen<br class="gmail_msg"></font></div><font face="garamond,serif" class="gmail_msg">Alvin E. Nashman Postdoctoral Fellow<br class="gmail_msg"></font></div><font face="garamond,serif" class="gmail_msg">Carnegie Observatories, Pasadena, CA<br class="gmail_msg"></font></div><font face="garamond,serif" class="gmail_msg"><a href="mailto:stonnes@gmail.com" class="gmail_msg" target="_blank">stonnes@gmail.com</a><br class="gmail_msg"></font></div></div></div>
</div></div>
_______________________________________________<br class="gmail_msg">
yt-users mailing list<br class="gmail_msg">
<a href="mailto:yt-users@lists.spacepope.org" class="gmail_msg" target="_blank">yt-users@lists.spacepope.org</a><br class="gmail_msg">
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" class="gmail_msg" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br class="gmail_msg">
</blockquote></div></div>