From e43cb70e2c51bf99c3bce14cbcceb2fc15606ef4 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= <atgeirr@sintef.no>
Date: Fri, 26 May 2023 13:20:37 +0200
Subject: [PATCH] Bugfix: use soMax consistent with
 BlackoilIntensiveQuantities::update().

---
 opm/models/blackoil/blackoilprimaryvariables.hh | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/opm/models/blackoil/blackoilprimaryvariables.hh b/opm/models/blackoil/blackoilprimaryvariables.hh
index e1874ee39..88aba097c 100644
--- a/opm/models/blackoil/blackoilprimaryvariables.hh
+++ b/opm/models/blackoil/blackoilprimaryvariables.hh
@@ -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,