avoid Gurson in CBM

NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp

@@ -513,7 +513,7 @@ void mlawNonLocalPorousCoupledLaw::checkCoalescence(IPNonLocalPorosity* q1, cons
             q1Thom->getRefToYieldOffset() = yieldGurson;
             if(yieldGurson > _tol){
               Msg::Info("Change yield surface to Gurson fG =%e",yieldGurson);
-              q1Thom->getRefToCoalescenceActiveFlag() = false; // Gurson will be used at the next time step
+              q1Thom->getRefToCoalescenceActiveFlag() = true; // Gurson will be used at the next time step
               q1Thom->getRefToAccelerateRate() = 1.0;
             }
             else{