From 94d7ef5b78be00ef5aea1cbef363efe3cfdf9dda Mon Sep 17 00:00:00 2001
From: ekinaci <efekinac@yahoo.com>
Date: Thu, 10 Jan 2019 15:35:51 +0100
Subject: [PATCH] add verbose switch to check reactionHeat per particle

---
 .../subModels/energyModel/reactionHeat/reactionHeat.C    | 9 +++++++++
 .../subModels/energyModel/reactionHeat/reactionHeat.H    | 2 ++
 2 files changed, 11 insertions(+)

diff --git a/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.C b/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.C
index 5ca1c6c6..b91d628d 100644
--- a/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.C
+++ b/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.C
@@ -45,6 +45,7 @@ reactionHeat::reactionHeat
     energyModel(dict,sm),
     propsDict_(dict.subDict(typeName + "Props")),
     interpolation_(propsDict_.lookupOrDefault<bool>("interpolation",false)),
+    verbose_(propsDict_.lookupOrDefault<bool>("verbose",false)),
     mesh_(sm.mesh()),
     maxSource_(1e30),
     reactionHeatName_(propsDict_.lookupOrDefault<word>("reactionHeatName","reactionHeat")),
@@ -97,6 +98,14 @@ void reactionHeat::calcEnergyContribution()
 
     particleCloud_.dataExchangeM().getData(reactionHeatName_,"scalar-atom",reactionHeat_);
 
+    for(int index = 0;index < particleCloud_.numberOfParticles(); ++index)
+    {
+        if (verbose_ && index>=0 && index < 2)
+        {
+            Pout << "reactionHeat = " << reactionHeat_[index][0] << endl;
+        }
+    }
+
     reactionHeatField_.primitiveFieldRef() = 0.0;
     reactionHeatField_.boundaryFieldRef() = 0.0;
 
diff --git a/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.H b/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.H
index 0fa6c805..73054f9a 100644
--- a/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.H
+++ b/src/lagrangian/cfdemParticle/subModels/energyModel/reactionHeat/reactionHeat.H
@@ -47,6 +47,8 @@ protected:
 
     bool interpolation_;
 
+    bool verbose_;
+
     const fvMesh& mesh_;
 
     scalar maxSource_;