Bugfix: follow opm-material convention for capillary pressure

This commit is contained in:
Atgeirr Flø Rasmussen 2017-08-04 11:39:51 +02:00
parent e9e4416428
commit c39341006d

View File

@ -499,8 +499,8 @@ namespace Opm {
// Compute phase pressures. // Compute phase pressures.
cstate.p[Oil] = constant(poval); cstate.p[Oil] = constant(poval);
cstate.p[Water] = cstate.p[Oil] - cstate.pc[Water]; // pcow = po - pw cstate.p[Water] = cstate.p[Oil] + cstate.pc[Water]; // pcow = pw - po (!) [different from old convention]
cstate.p[Gas] = cstate.p[Oil] + cstate.pc[Gas]; // pcog = pg - po (!) cstate.p[Gas] = cstate.p[Oil] + cstate.pc[Gas]; // pcog = pg - po
// Compute PVT properties. // Compute PVT properties.
cstate.temperature = constant(0.0); // Temperature is not used. cstate.temperature = constant(0.0); // Temperature is not used.