readSaltdenTable

This commit is contained in:
goncalvesmachadoc
2022-11-14 16:36:45 +01:00
parent 00d480fcd2
commit c826598759

View File

@@ -344,6 +344,14 @@ public:
return params_.saltsolTable_[pvtnumRegionIdx];
}
static const Scalar saltdenTable(const ElementContext& elemCtx,
unsigned scvIdx,
unsigned timeIdx)
{
unsigned pvtnumRegionIdx = elemCtx.problem().pvtRegionIndex(elemCtx, scvIdx, timeIdx);
return params_.saltdenTable_[pvtnumRegionIdx];
}
static bool hasBDensityTables()
{
return !params_.bdensityTable_.empty();
@@ -417,6 +425,9 @@ public:
const auto& saltsolTable = BrineModule::saltsolTable(elemCtx, dofIdx, timeIdx);
saltSolubility_ = saltsolTable;
const auto& saltdenTable = BrineModule::saltdenTable(elemCtx, dofIdx, timeIdx);
saltDensity_ = saltdenTable;
if (priVars.primaryVarsMeaningBrine() == PrimaryVariables::Sp) {
saltSaturation_ = priVars.makeEvaluation(saltConcentrationIdx, timeIdx);
fs.setSaltConcentration(saltSolubility_);