giving the friction pressure loss a direction.

This commit is contained in:
Kai Bao 2017-10-02 14:51:54 +02:00
parent ccc21e9dcd
commit aa8b1ae0f8

View File

@ -1644,7 +1644,9 @@ namespace Opm
const double area = segmentSet()[seg].crossArea();
const double diameter = segmentSet()[seg].internalDiameter();
return mswellhelpers::frictionPressureLoss(length, diameter, area, roughness, density, mass_rate, visc);
const double sign = mass_rate < 0. ? 1.0 : - 1.0;
return sign * mswellhelpers::frictionPressureLoss(length, diameter, area, roughness, density, mass_rate, visc);
}