Merge pull request #759 from totto82/fix_horizontalscaling

BUG: use the scaled water saturation
This commit is contained in:
Atgeirr Flø Rasmussen 2015-02-18 13:15:17 +01:00
commit 37bba7e3ae

View File

@ -602,7 +602,7 @@ namespace Opm
if (this->phase_usage.phase_used[BlackoilPhases::Aqua]) {
int pos = this->phase_usage.phase_pos[BlackoilPhases::Aqua];
double _sw = epst->wat.scaleSatPc(s[pos], this->smin_[pos], this->smax_[pos]);
pc[pos] = epst->wat.pcFactor*this->pcow_(s[pos]);
pc[pos] = epst->wat.pcFactor*this->pcow_(_sw);
double _dsdsw = epst->wat.scaleSatDerivPc(s[pos], this->smin_[pos], this->smax_[pos]);
dpcds[np*pos + pos] = epst->wat.pcFactor*_dsdsw*this->pcow_.derivative(_sw);
}