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: implement support for the PORV keyword

Browse Source
This commit is contained in:
Andreas Lauser 2015-02-13 14:01:07 +01:00
parent 7a729e73ac
commit df12c481c3
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

@ -710,6 +710,11 @@ private:
////////////////////////////////
// compute the porosity
if (!eclState->hasDoubleGridProperty("PORO") && !eclState->hasDoubleGridProperty("PORV"))
OPM_THROW(std::runtime_error,
"Can't read the porosity from the ECL state object. "
"(The PORO and PORV keywords are missing)");
if (eclState->hasDoubleGridProperty("PORO")) {
const std::vector<double> &poroData =
eclState->getDoubleGridProperty("PORO")->getData();
@ -719,10 +724,21 @@ private:
porosity_[dofIdx] = poroData[cartesianElemIdx];
}
}
else
OPM_THROW(std::runtime_error,
"Can't read the porosity from the ECL state object. "
"(The PORO keyword is missing)");
// overwrite the porosity using the PORV keyword for the elements for which PORV
// is defined...
if (eclState->hasDoubleGridProperty("PORV")) {
const std::vector<double> &porvData =
eclState->getDoubleGridProperty("PORV")->getData();
for (size_t dofIdx = 0; dofIdx < numDof; ++ dofIdx) {
int cartesianElemIdx = gridManager.cartesianCellId(dofIdx);
if (std::isfinite(porvData[cartesianElemIdx])) {
Scalar dofVolume = this->simulator().model().dofTotalVolume(dofIdx);
porosity_[dofIdx] = porvData[cartesianElemIdx]/dofVolume;
}
}
}
// apply the NTG keyword to the porosity
if (eclState->hasDoubleGridProperty("NTG")) {