correcting the determination of out-of-table-range velocity

comparing the location to the line instead of comparing the velocity
value directly.
This commit is contained in:
Kai Bao
2015-08-16 19:55:01 +02:00
committed by Kai Bao
parent fc2223605e
commit d9cb0066df

View File

@@ -505,7 +505,8 @@ namespace Opm
return false; // failed in finding the solution.
}
} else {
if (std::abs(water_vel[i]) < maxShearVel) {
// check if the failure in finding the shear multiplier is due to too big water velocity.
if ((logWaterVelO - logShearVRF.back()) < logShearWaterVel.back()) {
std::cout << " the veclocity is " << water_vel[i] << std::endl;
std::cout << " max shear velocity is " << maxShearVel << std::endl;
std::cerr << " something wrong happend in finding segment" << std::endl;