OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-28 02:00:59 -06:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #598 from akva2/use_tlmixpar_from_state

Browse Source 
use TLMIXPAR from state
This commit is contained in:
Arne Morten Kvarving 2020-03-04 09:58:25 +01:00 committed by GitHub
commit 913439ec28
1 changed files with 8 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
opm/models/blackoil/blackoilsolventmodules.hh
View File

@ -286,27 +286,21 @@ public:
}
if (deck.hasKeyword("TLMIXPAR")) {
if (!eclState.getTableManager().getTlmixparTable().empty()) {
// resize the attributes of the object
tlMixParamViscosity_.resize(numMiscRegions);
tlMixParamDensity_.resize(numMiscRegions);
assert(numMiscRegions == deck.getKeyword("TLMIXPAR").size());
const auto& tlmixparTable = eclState.getTableManager().getTlmixparTable();
assert(numMiscRegions == tlmixparTable.size());
for (unsigned regionIdx = 0; regionIdx < numMiscRegions; ++regionIdx) {
const auto& tlmixparRecord = deck.getKeyword("TLMIXPAR").getRecord(regionIdx);
const auto& mixParamsViscosityItem = tlmixparRecord.getItem("TL_VISCOSITY_PARAMETER");
const auto& mixParamsDensityItem = tlmixparRecord.getItem("TL_DENSITY_PARAMETER");
if (mixParamsViscosityItem.hasValue(0))
tlMixParamViscosity_[regionIdx] = mixParamsViscosityItem.getSIDouble(0);
else
throw std::invalid_argument("The TL_VISCOSITY parameter can not be defaulted");
if (mixParamsDensityItem.hasValue(0))
tlMixParamDensity_[regionIdx] = mixParamsDensityItem.getSIDouble(0);
else
// Copy data
tlMixParamViscosity_[regionIdx] = tlmixparTable[regionIdx].viscosity;
if (tlmixparTable[regionIdx].density == 0.0)
tlMixParamDensity_[regionIdx] = tlMixParamViscosity_[regionIdx];
else
tlMixParamDensity_[regionIdx] = tlmixparTable[regionIdx].density;
}
}
else