Merge pull request #137 from GitPaean/correcting_determination_in_range_velocity

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

opm/polymer/PolymerProperties.cpp

@@ -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;