<div dir="ltr">Hi all, <div><br></div><div>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><a href="http://yt-project.org/doc/analyzing/filtering.html">http://yt-project.org/doc/analyzing/filtering.html</a></div><div>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><br></div><div>Here is the quick check that should work on IsoGal:</div><div><br></div><div><div>ds = yt.load("/data001/stonnes/BPr16vs30arcfo/DD"+loop[i]+"/BPr16vs30arcfo"+loop[i])</div><div>    </div><div>    disk3 = ds.disk([0.5,0.5,0.5],[0,0,-1.0],(40,"kpc"),(10,"kpc"))</div><div>    </div><div>    diskn = disk3.cut_region(["obj['Metal_Density']/obj['density'] > 0.5"])</div><div>    diska = diskn.cut_region(["obj['density'] > 1e-25"])</div><div>    diskd = disk3.cut_region(["(obj['Metal_Density']/obj['density'] > 0.5) & (obj['density'] > 1e-25)"])</div><div>    </div><div>    print diska.quantities.total_quantity([("gas","cell_mass")])</div><div>    print disk3.quantities.total_quantity([("gas","cell_mass")])</div><div>    print diskn.quantities.total_quantity([("gas","cell_mass")])</div><div>    print diskd.quantities.total_quantity([("gas","cell_mass")])</div></div><div><br></div><div>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><br></div><div>Best,</div><div>Stephanie<br clear="all"><div><div class="gmail_signature"><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>