OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Initialise saturation points for Pc scaling

Browse Source 
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 6f75afc62b
commit 4da2c5862d
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/core/props/satfunc/SaturationPropsFromDeck_impl.hpp
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];