From 4250d4eda586ab9bb970a9b4d8cacedb407e1bce Mon Sep 17 00:00:00 2001 From: Liu Ming Date: Wed, 8 Oct 2014 10:23:10 +0800 Subject: [PATCH] use phase pressure to compute fluidRecipeFVF(). --- .../fullyimplicit/FullyImplicitBlackoilPolymerSolver_impl.hpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/opm/polymer/fullyimplicit/FullyImplicitBlackoilPolymerSolver_impl.hpp b/opm/polymer/fullyimplicit/FullyImplicitBlackoilPolymerSolver_impl.hpp index efff9e7ee..00de81243 100644 --- a/opm/polymer/fullyimplicit/FullyImplicitBlackoilPolymerSolver_impl.hpp +++ b/opm/polymer/fullyimplicit/FullyImplicitBlackoilPolymerSolver_impl.hpp @@ -599,7 +599,7 @@ namespace { for (int phase = 0; phase < maxnp; ++phase) { if (active_[ phase ]) { const int pos = pu.phase_pos[ phase ]; - rq_[pos].b = fluidReciprocFVF(phase, pressure, rs, rv, cond, cells_); + rq_[pos].b = fluidReciprocFVF(phase, pressure[phase], rs, rv, cond, cells_); rq_[pos].accum[aix] = pv_mult * rq_[pos].b * sat[pos]; // DUMP(rq_[pos].b); // DUMP(rq_[pos].accum[aix]);