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

Lookup for extractPorosity and extractDispersion_

Browse Source
This commit is contained in:
Antonella Ritorto 2023-12-18 15:47:00 +01:00
parent acc08799a4
commit 6d19c1dd53
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
ebos/ecltransmissibility_impl.hh
View File

@ -668,9 +668,10 @@ extractPorosity_()
// provided by eclState are one-per-cell of "uncompressed" grid, whereas the
// simulation grid might remove a few elements. (e.g. because it is distributed
// over several processes.)
unsigned numElem = gridView_.size(/*codim=*/0);
const auto& fp = eclState_.fieldProps();
if (fp.has_double("PORO")) {
porosity_ = fp.get_double("PORO");
porosity_ = this-> lookUpData_.assignFieldPropsDoubleOnLeaf(fp,"PORO", numElem);
}
else
throw std::logic_error("Can't read the porosityfrom the ecl state. "
@ -685,8 +686,9 @@ extractDispersion_()
throw std::runtime_error("Dispersion disabled at compile time, but the deck "
"contains the DISPERC keyword.");
}
unsigned numElem = gridView_.size(/*codim=*/0);
const auto& fp = eclState_.fieldProps();
dispersion_ = fp.get_double("DISPERC");
dispersion_ = this-> lookUpData_.assignFieldPropsDoubleOnLeaf(fp,"DISPERC", numElem);
}
template<class Grid, class GridView, class ElementMapper, class CartesianIndexMapper, class Scalar>