use phase pressure to compute fluidRecipeFVF().

This commit is contained in:
Liu Ming 2014-10-08 10:23:10 +08:00
parent 9beaf4e03b
commit 4250d4eda5

View File

@ -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]);