Force solvent saturation between [0,1]

This commit is contained in:
Tor Harald Sandve 2017-08-23 09:50:40 +02:00
parent 55fc987af8
commit 6716ca63c9

View File

@ -706,6 +706,7 @@ namespace Opm {
if (has_solvent_) {
double& ss = priVars[Indices::solventSaturationIdx];
ss -= satScaleFactor * dss;
ss = std::min(std::max(ss, 0.0),1.0);
}
if (has_polymer_) {
double& c = priVars[Indices::polymerConcentrationIdx];