diff --git a/ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelCalculator.cpp b/ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelCalculator.cpp
index 89e24b4193..7b8ef268af 100644
--- a/ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelCalculator.cpp
+++ b/ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelCalculator.cpp
@@ -436,6 +436,11 @@ bool RimStimPlanModelCalculator::calculateStressWithGradients( std::vector<doubl
         double Sv_diff               = 0.0;
         double deltaHorizontalStress = poissonsRatio / ( 1.0 - poissonsRatio ) * ( Sv_diff - biot * pressureDiff ) +
                                        ( biot * pressureDiff );
+        if ( std::isnan( deltaHorizontalStress ) || std::isinf( deltaHorizontalStress ) )
+        {
+            RiaLogging::warning( "Invalid horizontal stress delta calculated. Setting to zero." );
+            deltaHorizontalStress = 0.0;
+        }
 
         double depletionStress = Sh_init + deltaHorizontalStress;
         stress.push_back( RiaEclipseUnitTools::barToPsi( depletionStress ) );