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

Use only active elements when updating reference porosity

Browse Source
This commit is contained in:
Joakim Hove 2021-02-28 15:48:02 +01:00
parent 81ec44822b
commit d8830c721b
1 changed files with 2 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
ebos/eclproblem.hh
View File

@ -2669,11 +2669,10 @@ private:
referencePorosity_[/*timeIdx=*/0].resize(numDof); referencePorosity_[/*timeIdx=*/0].resize(numDof);
const auto& fp = eclState.fieldProps(); const auto& fp = eclState.fieldProps();
const std::vector<double> porvData = fp.porv(true); const std::vector<double> porvData = fp.porv(false);
const std::vector<int> actnumData = fp.actnum(); const std::vector<int> actnumData = fp.actnum();
for (size_t dofIdx = 0; dofIdx < numDof; ++ dofIdx) { for (size_t dofIdx = 0; dofIdx < numDof; ++ dofIdx) {
unsigned cartElemIdx = vanguard.cartesianIndex(dofIdx); Scalar poreVolume = porvData[dofIdx];
Scalar poreVolume = porvData[cartElemIdx];
// we define the porosity as the accumulated pore volume divided by the // we define the porosity as the accumulated pore volume divided by the
// geometric volume of the element. Note that -- in pathetic cases -- it can // geometric volume of the element. Note that -- in pathetic cases -- it can