diff --git a/applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H b/applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H
index 458f343e..3c82ddc8 100644
--- a/applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H
+++ b/applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H
@@ -33,7 +33,6 @@ tmp<fv::convectionScheme<scalar> > mvConvection
                combustion->R(Yi)
               + particleCloud.chemistryM(0).Smi(i)
               + fvOptions(rho, Yi)
-               
             );
 
             YiEqn.relax();
diff --git a/applications/solvers/cfdemSolverRhoPimpleChem/cfdemSolverRhoPimpleChem.C b/applications/solvers/cfdemSolverRhoPimpleChem/cfdemSolverRhoPimpleChem.C
index 33f9901b..a634ef7a 100644
--- a/applications/solvers/cfdemSolverRhoPimpleChem/cfdemSolverRhoPimpleChem.C
+++ b/applications/solvers/cfdemSolverRhoPimpleChem/cfdemSolverRhoPimpleChem.C
@@ -52,7 +52,6 @@ Description
 
 int main(int argc, char *argv[])
 {
-//    #include "postProcess.H"
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
@@ -61,7 +60,6 @@ int main(int argc, char *argv[])
     #include "createRDeltaT.H"
 
     #include "createFields.H"
-    //#include "createFieldRefs.H"
     #include "createFvOptions.H"
     #include "initContinuityErrs.H"
 
diff --git a/applications/solvers/cfdemSolverRhoPimpleChem/createFields.H b/applications/solvers/cfdemSolverRhoPimpleChem/createFields.H
index 2dbfcaa1..3e10274b 100644
--- a/applications/solvers/cfdemSolverRhoPimpleChem/createFields.H
+++ b/applications/solvers/cfdemSolverRhoPimpleChem/createFields.H
@@ -14,7 +14,6 @@
     PtrList<volScalarField>& Y = composition.Y();
 
     // read molecular weight
-    //const volScalarField& W = composition.W();
     volScalarField W(composition.W());
 
     bool propagateInertSpecie = true;
diff --git a/src/lagrangian/cfdemParticle/subModels/chemistryModel/reactantPerParticle/reactantPerParticle.C b/src/lagrangian/cfdemParticle/subModels/chemistryModel/reactantPerParticle/reactantPerParticle.C
index 097d6eb1..0d2151c7 100644
--- a/src/lagrangian/cfdemParticle/subModels/chemistryModel/reactantPerParticle/reactantPerParticle.C
+++ b/src/lagrangian/cfdemParticle/subModels/chemistryModel/reactantPerParticle/reactantPerParticle.C
@@ -56,7 +56,7 @@ reactantPerParticle::reactantPerParticle
     chemistryModel(dict,sm),
     propsDict_(dict.subDict(typeName + "Props")),
     mesh_(sm.mesh()),
-        verbose_(propsDict_.lookupOrDefault<bool>("verbose",false)),
+    verbose_(propsDict_.lookupOrDefault<bool>("verbose",false)),
     reactantPerParticle_(NULL),
     voidfractionFieldName_(propsDict_.lookupOrDefault<word>("voidfractionFieldName","voidfraction")),
     voidfraction_(sm.mesh().lookupObject<volScalarField>(voidfractionFieldName_)),
diff --git a/src/lagrangian/cfdemParticle/subModels/chemistryModel/species/species.C b/src/lagrangian/cfdemParticle/subModels/chemistryModel/species/species.C
index 23f32464..570a9127 100644
--- a/src/lagrangian/cfdemParticle/subModels/chemistryModel/species/species.C
+++ b/src/lagrangian/cfdemParticle/subModels/chemistryModel/species/species.C
@@ -290,8 +290,6 @@ void species::execute()
         }
     }
 
-    //Info << "partTemp from species.C = " << partTemp_[100][0] << endl;
-
         // give DEM data
         particleCloud_.dataExchangeM().giveData(partTempName_, "scalar-atom", partTemp_);
         particleCloud_.dataExchangeM().giveData(partRhoName_,  "scalar-atom", partRho_);
diff --git a/src/lagrangian/cfdemParticle/subModels/energyModel/heatTransferGunn/heatTransferGunn.C b/src/lagrangian/cfdemParticle/subModels/energyModel/heatTransferGunn/heatTransferGunn.C
index 655626bb..3ec830cc 100644
--- a/src/lagrangian/cfdemParticle/subModels/energyModel/heatTransferGunn/heatTransferGunn.C
+++ b/src/lagrangian/cfdemParticle/subModels/energyModel/heatTransferGunn/heatTransferGunn.C
@@ -359,9 +359,8 @@ void heatTransferGunn::calcEnergyContribution()
                     partNu_[index][0] = Nup;
                 }
 
-                 //if(particleCloud_.verbose() && index >=0 && index <2)
-                 if(verbose_ && index >=0 && index <2)
-                 {
+                if(verbose_ && index >=0 && index <2)
+                {
                     Pout << "partHeatFlux = " << partHeatFlux_[index][0] << endl;
                     Pout << "magUr = " << magUr << endl;
                     Pout << "kf0 = " << kf0_ << endl;
@@ -489,9 +488,6 @@ scalar heatTransferGunn::Nusselt(scalar voidfraction, scalar Rep, scalar Pr) con
 
 void heatTransferGunn::heatFlux(label index, scalar h, scalar As, scalar Tfluid, scalar cg3)
 {
-    //Info << "partTemp from HTGunn.C = " << partTemp_[100][0] << endl;
-    //Info << "partTemp = " << partTemp_[index][0] << endl;
-    //Info << "Tfluid =  " << Tfluid << endl;
     scalar hAs = h * As * cg3;
 
     if (particleCloud_.getParticleEffVolFactors())