correct method 3

NonLinearSolver/materialLaw/mlawNonLocalPorousCoalescence.cpp

@@ -129,9 +129,12 @@ double mlawNonLocalPorousThomasonLaw::yieldFunction(const double kcorEq, const d
   // Get values in ipvs
   double Cft = q1->getConstRefToIPCoalescence().getConcentrationFactor();
   double CftOffset = q0->getConstRefToIPCoalescence().getCrackOffsetOnCft();
-  if (_CfTOffsetMethod==2 or _CfTOffsetMethod==3){
+  if (_CfTOffsetMethod==2){
     if (!q0->getConstRefToIPCoalescence().getCoalescenceOnsetFlag()){CftOffset = 1.;};
-  };
+  }
+  else if (_CfTOffsetMethod==3){
+    if (!q0->getConstRefToIPCoalescence().getCoalescenceActiveFlag()){CftOffset = 1.;};
+  }
 
   double R0 = _j2IH->getYield0();
   if (_roundedMethod == 1){

NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp

@@ -501,7 +501,7 @@ void mlawNonLocalPorousCoupledLaw::checkCoalescence(IPNonLocalPorosity* q1, cons
             if (_CfTOffsetMethod == 3){
               q1Thom->getRefToCrackOffsetOnCft() = 1. + yieldThomason*R0/(Cft*R);
               
-              if (q1Thom->getRefToCrackOffsetOnCft() > 2.){Msg::Error("High value of Cft Offset = %e",q1Thom->getRefToCrackOffsetOnCft());};
+              if (q1Thom->getRefToCrackOffsetOnCft() > 2.){Msg::Error("High predicted value of Cft Offset = %e",q1Thom->getRefToCrackOffsetOnCft());};
             }
             
           }