diff --git a/NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp b/NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp
index de2feaa642ab6fe0a05e708ccc790b32b6ead68e..ef3dfaa923cb46ba117bf088dd89676adea7ca53 100644
--- a/NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp
+++ b/NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp
@@ -577,6 +577,21 @@ void mlawNonLocalPorousCoupledLaw::checkCoalescence(IPNonLocalPorosity* q1, cons
     else
     {
       // In the bulk:
+/*
+      // Only growth model is allowed
+      q1Thom->getRefToCoalescenceOnsetFlag() = false;
+      // Reset onset variables to their default value
+      q1Thom->getRefToPorosityAtCoalescenceOnset() = 0.;
+      q1Thom->getRefToLigamentRatioAtCoalescenceOnset() = 0.;
+      q1Thom->getRefToAspectRatioAtCoalescenceOnset() = 0.;
+      q1Thom->getRefToShapeFactorAtCoalescenceOnset() = 0.;
+
+      q1Thom->getRefToCoalescenceActiveFlag() = false;
+      q1Thom->getRefToCrackOffsetOnCft() = 1.;
+      q1Thom->getRefToAccelerateRate() = 1.;
+      q1Thom->getRefToYieldOffset() = 0.;
+
+*/
       if (!q0Thom->getCoalescenceOnsetFlag()){
         double yieldThomason = _mlawCoales->yieldFunction(kcorEq,pcor,R,yieldfV,q0,q1,T);
 
@@ -630,6 +645,11 @@ void mlawNonLocalPorousCoupledLaw::checkCoalescence(IPNonLocalPorosity* q1, cons
         q1Thom->getRefToAccelerateRate() = q0Thom->getAccelerateRate();
       }
     }
+    
+    
+    
+    // In the bulk: end
+    
   }
   else{
     // If crack transition is not used: no distinction between bulk and interface ipvs