From 65a6604112f500391d84db2d4f8b25ae62d5bd88 Mon Sep 17 00:00:00 2001
From: Ove Saevareid <osae@iris.no>
Date: Mon, 22 Nov 2010 13:07:24 +0100
Subject: [PATCH] Oil relperm: krocw_ evaluated at connate water, krow at sw.

---
 .../blackoil/fluid/FluidMatrixInteractionBlackoil.hpp    | 9 ++-------
 1 file changed, 2 insertions(+), 7 deletions(-)

diff --git a/dune/porsol/blackoil/fluid/FluidMatrixInteractionBlackoil.hpp b/dune/porsol/blackoil/fluid/FluidMatrixInteractionBlackoil.hpp
index 91cbb9a0..798f1b6e 100644
--- a/dune/porsol/blackoil/fluid/FluidMatrixInteractionBlackoil.hpp
+++ b/dune/porsol/blackoil/fluid/FluidMatrixInteractionBlackoil.hpp
@@ -64,12 +64,7 @@ public:
         buildUniformMonotoneTable(sw, krow, samples, krow_);
         buildUniformMonotoneTable(sg, krg,  samples, krg_);
         buildUniformMonotoneTable(sg, krog, samples, krog_);
-        for (int i = 0; i < int(krw.size()); ++i) {
-            if (krw[i] > 0.0) {
-                krocw_ = krow[i];
-                break;
-            }
-        }
+        krocw_ = krow[0]; // At connate water -> ecl. SWOF
 
         // Create tables for pcow and pcog.
         // We must convert the pressures depending on units.
@@ -174,7 +169,7 @@ public:
         Scalar so = saturations[Liquid];
         Scalar krw = params.krw_(sw);
         Scalar krg = params.krg_(sg);
-        Scalar krow = params.krow_(so);
+        Scalar krow = params.krow_(sw);
         Scalar krog = params.krog_(sg);
         Scalar krocw = params.krocw_;
         kr[Aqua] = krw;