diff --git a/applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H b/applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H
index df180fa7..9b643abd 100644
--- a/applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H
+++ b/applications/solvers/cfdemSolverRhoPimpleChem/YEqn.H
@@ -50,8 +50,11 @@ tmp<fv::convectionScheme<scalar> > mvConvection
         }
     }
 
+    // calculate total molar number in volume
     N = Yt*rhoeps/W;
 
+    Info << "molecular weight of mixture" << W << endl;
+
     Y[inertIndex] = scalar(1) - Yt;
     Y[inertIndex].max(0.0);
 }
diff --git a/applications/solvers/cfdemSolverRhoPimpleChem/createFields.H b/applications/solvers/cfdemSolverRhoPimpleChem/createFields.H
index 89b7d252..e7a19636 100644
--- a/applications/solvers/cfdemSolverRhoPimpleChem/createFields.H
+++ b/applications/solvers/cfdemSolverRhoPimpleChem/createFields.H
@@ -18,9 +18,8 @@
     basicSpecieMixture& composition = thermo.composition();
     PtrList<volScalarField>& Y = composition.Y();
 
-    // The molecular weight [kj/mol]  ====== EKI
-    //const volScalarField W(composition.W());
-    const volScalarField& W = composition.W();
+    // The molecular weight [kj/mol]
+    const volScalarField W(composition.W());
 
     const word inertSpecie(thermo.lookup("inertSpecie"));
 
@@ -42,7 +41,7 @@
     }
     fields.add(thermo.he());
 
-    // define total number of moles in volume ====== EKI
+    // Total Mole in Volume
     volScalarField N
     (
         IOobject
@@ -55,7 +54,6 @@
         ),
         mesh
     );
-    // =============================================
 
     volScalarField rho
     (