Bugfix: use soMax consistent with BlackoilIntensiveQuantities::update().

This commit is contained in:
Atgeirr Flø Rasmussen
2023-05-26 13:20:37 +02:00
parent 65739e15d6
commit e43cb70e2c
@@ -690,7 +690,7 @@ public:
if (sg < -eps && s > 0.0 && FluidSystem::enableDissolvedGas()) {
const Scalar& po = (*this)[pressureSwitchIdx];
setPrimaryVarsMeaningGas(GasMeaning::Rs);
Scalar soMax = problem.maxOilSaturation(globalDofIdx);
Scalar soMax = std::max(s, problem.maxOilSaturation(globalDofIdx));
Scalar rsMax = problem.maxGasDissolutionFactor(/*timeIdx=*/0, globalDofIdx);
Scalar rsSat = enableExtbo ? ExtboModule::rs(pvtRegionIndex(),
po,