BUGFIX. Store soMax also for the solvent model

This commit is contained in:
Tor Harald Sandve 2017-04-06 14:39:30 +02:00
parent 431abb0012
commit ac36ad5ad1

View File

@ -553,6 +553,8 @@ namespace Opm {
rv = rv.max(zero);
}
sd_.soMax = fluid_.satOilMax();
// Sg is used as primal variable for water only cells.
const double epsilon = std::sqrt(std::numeric_limits<double>::epsilon());
auto watOnly = sw > (1 - epsilon);