Merge pull request #690 from blattms/fix-spe1-brine

[bugfix] Prevent accessing undefined variables for brine without salt precipitation
This commit is contained in:
Markus Blatt
2022-04-21 12:31:21 +02:00
committed by GitHub

View File

@@ -348,7 +348,9 @@ protected:
const double sign = delta >= 0. ? 1. : -1.;
delta = sign * std::min(std::abs(delta), maxTempChange_);
}
else if (enableBrine && pvIdx == Indices::saltConcentrationIdx && currentValue.primaryVarsMeaningBrine() == PrimaryVariables::Sp) {
else if (enableBrine && pvIdx == Indices::saltConcentrationIdx &&
enableSaltPrecipitation &&
currentValue.primaryVarsMeaningBrine() == PrimaryVariables::Sp) {
const double maxSaltSaturationChange = 0.1;
const double sign = delta >= 0. ? 1. : -1.;
delta = sign * std::min(std::abs(delta), maxSaltSaturationChange);