diff --git a/ApplicationCode/ReservoirDataModel/RigGeoMechWellLogExtractor.cpp b/ApplicationCode/ReservoirDataModel/RigGeoMechWellLogExtractor.cpp index 42f26f00a9..bf147f136b 100644 --- a/ApplicationCode/ReservoirDataModel/RigGeoMechWellLogExtractor.cpp +++ b/ApplicationCode/ReservoirDataModel/RigGeoMechWellLogExtractor.cpp @@ -113,7 +113,7 @@ float RigGeoMechWellLogExtractor::calculatePorePressureInSegment(int64_t interse double porePressure = hydroStaticPorePressureBar; // 1: Try pore pressure from the grid - if (porePressure == hydroStaticPorePressureBar && averageSegmentPorePressureBar != std::numeric_limits::infinity()) + if (porePressure == hydroStaticPorePressureBar && averageSegmentPorePressureBar > 0.0) { porePressure = averageSegmentPorePressureBar; } @@ -137,6 +137,8 @@ float RigGeoMechWellLogExtractor::calculatePorePressureInSegment(int64_t interse porePressure = pascalToBar(poreElementPressuresPascal[elmIdx]); } // 4: If no pore-pressure was found, the default value of hydrostatic pore pressure is used. + + CVF_ASSERT(porePressure >= 0.0); return porePressure; }