Initialise saturation points for Pc scaling

This gives predictable failure modes although the actual behaviour
is (probably) unchanged.
This commit is contained in:
Bård Skaflestad 2014-07-29 11:34:57 +02:00
parent d067ef07a2
commit 59eedcd82b

View File

@ -431,9 +431,9 @@ namespace Opm
} else {
const int wpos = phase_usage_.phase_pos[BlackoilPhases::Aqua];
const int np = BlackoilPhases::MaxNumPhases;
double s[np];
double s[np] = { 0.0 };
s[wpos] = swat;
double pc[np];
double pc[np] = { 0.0 };
funcForCell(cell).evalPc(s, pc, &(eps_transf_[cell]));
if (pc[wpos] > 1.0e-8) {
eps_transf_[cell].wat.pcFactor *= pcow/pc[wpos];