From ec11bbaac30cdd500538e2d1985c8c3732342edf Mon Sep 17 00:00:00 2001
From: hnil <Halvor.M.Nilsen@sintef.no>
Date: Thu, 29 Feb 2024 11:02:26 +0100
Subject: [PATCH] Avoid deleting derivative i.e. getting "wrong" matrix in case
 of zero flux

---
 opm/simulators/flow/NewTranFluxModule.hpp | 4 ----
 1 file changed, 4 deletions(-)

diff --git a/opm/simulators/flow/NewTranFluxModule.hpp b/opm/simulators/flow/NewTranFluxModule.hpp
index 39f74c58a..08f5d203b 100644
--- a/opm/simulators/flow/NewTranFluxModule.hpp
+++ b/opm/simulators/flow/NewTranFluxModule.hpp
@@ -285,10 +285,6 @@ public:
                                         distZ*g,
                                         thpres,
                                         problem.moduleParams());
-            if (pressureDifferences[phaseIdx] == 0) {
-                volumeFlux[phaseIdx] = 0.0;
-                continue;
-            }
 
             const bool upwindIsInterior = (static_cast<unsigned>(upIdx[phaseIdx]) == interiorDofIdx);
             const IntensiveQuantities& up = upwindIsInterior ? intQuantsIn : intQuantsEx;