OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

ebos: use grid.cellCenterDepth() instead of the Z coordinate of the centroid

Browse Source 
grid.cellCenterDepth() is the average of the Z-coordinates of the
element vertices which, for distorted elements, is slightly different
from the depth of the centroid.

IMO this is conceptually a change for the worse, but ECL likes to do a
lot of things inconsistently, so let's budge.
This commit is contained in:
Andreas Lauser 2016-08-05 14:43:56 +02:00
parent 679651bca1
commit 6e9b13ad5a
1 changed files with 8 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
applications/ebos/eclfluxmodule.hh
View File

@ -210,6 +210,7 @@ protected:
Valgrind::SetUndefined(*this);
const auto& problem = elemCtx.problem();
const auto& grid = elemCtx.simulator().gridManager().grid();
const auto& stencil = elemCtx.stencil(timeIdx);
const auto& scvf = stencil.interiorFace(scvfIdx);
@ -231,8 +232,13 @@ protected:
const auto &intQuantsIn = elemCtx.intensiveQuantities(interiorDofIdx_, timeIdx);
const auto &intQuantsEx = elemCtx.intensiveQuantities(exteriorDofIdx_, timeIdx);
Scalar zIn = elemCtx.pos(interiorDofIdx_, timeIdx)[dimWorld - 1];
Scalar zEx = elemCtx.pos(exteriorDofIdx_, timeIdx)[dimWorld - 1];
// this is quite hacky because the dune grid interface does not provide a
// cellCenterDepth() method (so this code is specific to the OPM derived
// grids). The "good" solution would be to take the Z coordinate of the element
// centers, but since ECL seems to like to be inconsistent on that front, it
// needs to be done like here...
Scalar zIn = grid.cellCenterDepth(I);
Scalar zEx = grid.cellCenterDepth(J);
// the distances from the DOF's depths. (i.e., the additional depth of the
// exterior DOF)