ECL peaceman well: fix sign error in gravity correction

This commit is contained in:
Andreas Lauser 2015-01-12 11:55:16 +01:00
parent 1233e12e76
commit e7a80216e8

View File

@ -1178,9 +1178,9 @@ protected:
Valgrind::CheckDefined(refDepth_);
// pressure in the borehole ("hole pressure") at the given location
Scalar ph = pbh + rho*g*(refDepth_ - depth);
Scalar ph = pbh + rho*g*(depth - refDepth_);
// volumetric flux of the phase from the well to the reservoir
// volumetric reservoir rate for the phase
volRates[phaseIdx] = Twj*lambda*(ph - p);
Valgrind::CheckDefined(g);