From 1cb1e8475a7f6306df0d9078980a774ddba4365a Mon Sep 17 00:00:00 2001
From: Liu Ming <miliu@statoil.com>
Date: Fri, 26 Sep 2014 14:05:17 +0800
Subject: [PATCH] add PhasePresence, make muOil() function happy.

---
 opm/polymer/fullyimplicit/utilities.cpp | 8 +++++++-
 1 file changed, 7 insertions(+), 1 deletion(-)

diff --git a/opm/polymer/fullyimplicit/utilities.cpp b/opm/polymer/fullyimplicit/utilities.cpp
index c75185ac0..8b029b255 100644
--- a/opm/polymer/fullyimplicit/utilities.cpp
+++ b/opm/polymer/fullyimplicit/utilities.cpp
@@ -245,12 +245,18 @@ namespace Opm
 				src[cell] = 1.0;
 			}
 		}
+        //Add PhasePresence make muOil() happy.
+        std::vector<PhasePresence> phaseCondition;
+        for (int c = 0; c < num_cells; ++c) {
+            phaseCondition[c].setFreeWater();
+            phaseCondition[c].setFreeOil();
+        }
 		const Selector<double> src_selector(src);
 		const V one = V::Constant(num_cells, 1.0);
 		const V zero = V::Zero(num_cells);
 		const std::vector<V> kr = props.relperm(sw, so, zero, cells);
 		const V muw = props.muWat(p, cells);
-		const V muo = props.muOil(p, zero, cells);
+		const V muo = props.muOil(p, zero, phaseCondition, cells);
         const V krw_eff = polymer_props.effectiveRelPerm(c, cmax, kr[0]);
 		const V inv_muw_eff = polymer_props.effectiveInvWaterVisc(c, muw.data());
 		std::vector<V> mob(np);