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