OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-11-27 19:50:16 -06:00
Code Issues Packages Projects Releases Wiki Activity

[bugfix] use DUNE compatible way of computing cell center depth.

Browse Source
This commit is contained in:
Robert Kloefkorn 2016-10-25 16:40:39 +02:00
parent 8784c81a4e
commit 26e5abf699
1 changed files with 20 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
applications/ebos/eclproblem.hh
View File

@ -805,10 +805,25 @@ private:
static bool enableEclOutput_()
{ return EWOMS_GET_PARAM(TypeTag, bool, EnableEclOutput); }
Scalar cellCenterDepth( const Element& element ) const
{
typedef typename Element :: Geometry Geometry;
static constexpr int zCoord = Geometry ::dimension - 1;
Scalar zz = 0.0;
const Geometry geometry = element.geometry();
const int corners = geometry.corners();
for (int i=0; i<corners; ++i)
{
zz += geometry.corner( i )[ zCoord ];
}
return zz/Scalar(corners);
}
void updateElementDepths_()
{
const auto& gridManager = this->simulator().gridManager();
const auto& grid = gridManager.grid();
const auto& gridView = gridManager.gridView();
const auto& elemMapper = this->elementMapper();;
@ -818,13 +833,14 @@ private:
auto elemIt = gridView.template begin</*codim=*/0>();
const auto& elemEndIt = gridView.template end</*codim=*/0>();
for (; elemIt != elemEndIt; ++elemIt) {
const Element& element = *elemIt;
#if DUNE_VERSION_NEWER(DUNE_COMMON, 2,4)
unsigned elemIdx = elemMapper.index(*elemIt);
const unsigned int elemIdx = elemMapper.index(element);
#else
unsigned elemIdx = elemMapper.map(*elemIt);
const unsigned int elemIdx = elemMapper.map(element);
#endif
elementCenterDepth_[elemIdx] = grid.cellCenterDepth(elemIdx);
elementCenterDepth_[elemIdx] = cellCenterDepth( element );
}
}