finish implementation of cft offset method 0 and 1

NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp

@@ -123,7 +123,12 @@ void mlawNonLocalPorousCoupledLaw::setOnsetTriaxialityForRoundedYieldSurface(con
 
 void mlawNonLocalPorousCoupledLaw::setCfTOffsetMethod(const int MethodNumber){
   _CfTOffsetMethod = MethodNumber;
+  if (_CfTOffsetMethod == 0 or _CfTOffsetMethod == 1 or _CfTOffsetMethod == 2){
     Msg::Info("mlawNonLocalPorousCoupledLaw::useCfTOffset: Offset method %d is used",_CfTOffsetMethod);
+  }
+  else{
+    Msg::Fatal("mlawNonLocalPorousCoupledLaw::useCfTOffset: Offset method %d is not existing",_CfTOffsetMethod);
+  }
 };
 
 void mlawNonLocalPorousCoupledLaw::setCrackTransition(const bool fl){
@@ -388,6 +393,10 @@ void mlawNonLocalPorousCoupledLaw::checkCoalescence(IPNonLocalPorosity* q1, cons
         // Check Thomason criterion:
         double yieldThomason = _mlawCoales->yieldFunction(kcorEq,pcor,R,yieldfV,q0,q1,T);
         /* check this value */
+        if(_CfTOffsetMethod == 2)
+        {
+          q1Thom->getRefToCrackOffsetOnCft() = 1. + yieldThomason*R0/(Cft*R);
+        }
         if ((sigmax - Cft*R > _tol) and (yieldThomason > _tol)){
           // If coalescence occurs, update Cft offset value and others variables
           //