Merge pull request #1125 from totto82/BUGFIXsoMaxSolvent

BUGFIX. Store soMax also for the solvent model
This commit is contained in:
Atgeirr Flø Rasmussen
2017-04-06 21:19:28 +02:00
committed by GitHub

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