Fixed bug for computeMc.

This commit is contained in:
Xavier Raynaud 2012-03-26 13:45:14 +02:00
parent 4e3f486d13
commit 023cbb228e

View File

@ -220,7 +220,6 @@ namespace Opm
} }
void computeMc(const double& c, void computeMc(const double& c,
const double* visc,
double& mc, double& mc,
double& dmcdc) const double& dmcdc) const
{ {
@ -228,7 +227,7 @@ namespace Opm
double omega = mix_param_; double omega = mix_param_;
double r = std::pow(viscMult(c_max_), 1 - omega); // viscMult(c_max_)=mu_p/mu_w double r = std::pow(viscMult(c_max_), 1 - omega); // viscMult(c_max_)=mu_p/mu_w
mc = c/(cbar + (1 - cbar)*r); mc = c/(cbar + (1 - cbar)*r);
dmdc = r/std::pow(cbar + (1 - cbar)*r, 2); dmcdc = r/std::pow(cbar + (1 - cbar)*r, 2);
} }