How to get cell volume?

[wuming]

Hi, Dr. Dietiker,

Which keywords can we use to get cell center coordinates?

I noticed that they can be output in the MFIX GUI.

I have searched for them in the source code, but I haven’t found them yet. Two years ago, you answered my similar question, but I am not sure if they can be used to output the cell center coordinates like the option in the GUI does.

How to have access to all cell data?

You can access the East, North and Top face locations of a scalar cell with use cutcell, only: xg_e, yg_n, zg_t at the top of the routine (with other “use” statements). Then if you want the location in a cell that has I,J,K indices, use XG_E(I) , YG_N(J), and ZG_T(K) . If you need the cell center coordinate, you have to offset the face location by half the cell size (dx,dy,dz). These are in the “geometry” module: USE geometry, only: dx,dy,dz so the cell center would be at xc = XG_E(I) - 0.…