<div dir="ltr"><div><div><div><div><div><div><div><div><div>Hi all,<br><br></div>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:<br>
<br></div><div>alld = pf.h.all_data()<br></div>tracerp25 = alld.cut_region(["grid['specific_scalar[0]'] > 0.25","grid['z'] < 0.1925","grid['z'] > -0.1925","grid['cyl_Rcode'] < 1.2"])<br>
<br></div>disk = pf.h.disk([0.0,0.0,0.0],[0.0,0.0,1.0],1.2,0.385)<br></div>tdp25 = disk.cut_region(["grid['specific_scalar[0]'] > 0.25"])<br><br></div>but when I print the total mass I get answers that differ by 10%<br>
<br></div>print tracerp25.quantities["TotalQuantity"]("CellMassMsun")<br></div>print tdp25.quantities["TotalQuantity"]("CellMassMsun")<br><br>>>[8019047231.5192223]<br>>>[8886968784.4710732]<br>
<br></div>Can someone point out where my difference is?<br><br></div>Thanks!<br><br>Stephanie<br></div>