Hi Suoqing, On 18/06/13 19:00, Suoqing Ji wrote:
I assume the quantity 'CylindricalCellMass' is a derived field defined by yourself? I usually define the 2d cylindrical cell volume as followings: [...] yes, that's pretty much the expression I was using.
And for my data set, the output given by YT (use dd.quantities['TotalQuantity']('CylCellMass')[0]) and FLASH subroutine (IO_writeIntegralQuantities) match down to the machine accuracy.
ahh, thank you for your response! It made me realize two things: (a) I made a stupid mistake while sorting different time steps, and for a normal 2d cylindrical FLASH simulation I do have convergence to machine precision between the two methods. (b) I was post-processing a simulation that has a few coarse blocks removed and used as new boundary conditions in order to get a non-rectangular domain (see source/Simulation/SimulationMain/WindTunnel/Simulation_defineDomain.F90 in FLASH source, or the illustration here: http://perso.crans.org/yurchak/i/illustations_yt.pdf ). This seems to be the origin of the issue I've been experiencing (and that's still there). I looks like yt, treats this configuration partly ok as for instance none of leaf blocks has for parent a block that was removed. On the other side though: * there this issue with d.quantities['TotalQuantity'] * a slice shows that there is some plasma propagating through the region that's supposed to be treated as outside of the simulation (with removed blocks) and there is nothing there in VisIt. So, anyway, thanks for narrowing it down. Would you have any suggestions how to deal with this issue? -- Roman