[yt-users] Cutting region confusion

[Stephanie Tonnesen]

Hi all,

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:

alld = pf.h.all_data()
tracerp25 = alld.cut_region(["grid['specific_scalar[0]'] > 0.25","grid['z']
< 0.1925","grid['z'] > -0.1925","grid['cyl_Rcode'] < 1.2"])

disk = pf.h.disk([0.0,0.0,0.0],[0.0,0.0,1.0],1.2,0.385)
tdp25 = disk.cut_region(["grid['specific_scalar[0]'] > 0.25"])

but when I print the total mass I get answers that differ by 10%

print tracerp25.quantities["TotalQuantity"]("CellMassMsun")
print tdp25.quantities["TotalQuantity"]("CellMassMsun")

>>[8019047231.5192223]
>>[8886968784.4710732]

Can someone point out where my difference is?

Thanks!

Stephanie
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140804/39ce573b/attachment.htm>