diff --git a/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.C b/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.C
index 54173150..61d1a530 100644
--- a/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.C
+++ b/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.C
@@ -47,6 +47,7 @@ reactionHeat::reactionHeat
     propsDict_(dict.subDict(typeName + "Props")),
     interpolation_(propsDict_.lookupOrDefault<bool>("interpolation",false)),
     verbose_(propsDict_.lookupOrDefault<bool>("verbose",false)),
+    execution_(true),
     mesh_(sm.mesh()),
     maxSource_(1e30),
     reactionHeatName_(propsDict_.lookupOrDefault<word>("reactionHeatName","reactionHeat")),
@@ -94,7 +95,8 @@ reactionHeat::~reactionHeat()
 void reactionHeat::calcEnergyContribution()
 {
     loopCounter_++;
-    if (loopCounter_ % Nevery_ != 0)
+    execution_ = (particleCloud_.dataExchangeM().couplingStep() % Nevery_ == 0);
+    if (!execution_)
     {
         return;
     }
diff --git a/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.H b/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.H
index 003c57bf..ef5d9397 100644
--- a/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.H
+++ b/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.H
@@ -50,6 +50,8 @@ protected:
 
     bool verbose_;
 
+    bool execution_;
+
     const fvMesh& mesh_;
 
     scalar maxSource_;