From 92592672fe6203bbee5a5948c647084ad4ca9aba Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?B=C3=A5rd=20Skaflestad?= <Bard.Skaflestad@sintef.no>
Date: Wed, 8 May 2013 15:06:13 +0200
Subject: [PATCH] Use upwind mobility flux.

---
 ImpesTPFAAD.hpp | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/ImpesTPFAAD.hpp b/ImpesTPFAAD.hpp
index 7b3eaf408..8304aa363 100644
--- a/ImpesTPFAAD.hpp
+++ b/ImpesTPFAAD.hpp
@@ -233,7 +233,8 @@ namespace Opm {
 
                 const V   kr = pdepfdata_.phaseRelPerm(phase);
                 const ADB mu = pdepfdata_.phaseViscosity(phase);
-                const ADB flux = (kr / mu) * nkgradp;
+                const ADB mf = upwind.select(kr / mu);
+                const ADB flux = mf * nkgradp;
 
                 const ADB face_B = upwind.select(cell_B);