From 09a15c6e583ea4f63b85451d559112670a3eb1ff Mon Sep 17 00:00:00 2001 From: Tor Harald Sandve Date: Wed, 18 Feb 2015 08:11:50 +0100 Subject: [PATCH] BUG: use the scaled water saturation --- opm/core/props/satfunc/SatFuncGwseg.hpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/opm/core/props/satfunc/SatFuncGwseg.hpp b/opm/core/props/satfunc/SatFuncGwseg.hpp index 4a6d0fed..c6acc9ee 100644 --- a/opm/core/props/satfunc/SatFuncGwseg.hpp +++ b/opm/core/props/satfunc/SatFuncGwseg.hpp @@ -602,7 +602,7 @@ namespace Opm if (this->phase_usage.phase_used[BlackoilPhases::Aqua]) { int pos = this->phase_usage.phase_pos[BlackoilPhases::Aqua]; double _sw = epst->wat.scaleSatPc(s[pos], this->smin_[pos], this->smax_[pos]); - pc[pos] = epst->wat.pcFactor*this->pcow_(s[pos]); + pc[pos] = epst->wat.pcFactor*this->pcow_(_sw); double _dsdsw = epst->wat.scaleSatDerivPc(s[pos], this->smin_[pos], this->smax_[pos]); dpcds[np*pos + pos] = epst->wat.pcFactor*_dsdsw*this->pcow_.derivative(_sw); }