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