diff --git a/applications/solvers/combustion/plasmaReactingFoam/PhiEqn.H b/applications/solvers/combustion/plasmaReactingFoam/PhiEqn.H
index dc857322..440bb8fa 100644
--- a/applications/solvers/combustion/plasmaReactingFoam/PhiEqn.H
+++ b/applications/solvers/combustion/plasmaReactingFoam/PhiEqn.H
@@ -7,5 +7,8 @@
     );
 
     E = -fvc::grad(Phi);
+
+    ng = p / T * (NA / R);
+
     En = mag(E) / (ng);
 }
diff --git a/applications/solvers/combustion/plasmaReactingFoam/YEqn.H b/applications/solvers/combustion/plasmaReactingFoam/YEqn.H
index 4e2177be..d870f6bd 100644
--- a/applications/solvers/combustion/plasmaReactingFoam/YEqn.H
+++ b/applications/solvers/combustion/plasmaReactingFoam/YEqn.H
@@ -17,10 +17,11 @@ tmp<fv::convectionScheme<scalar> > mvConvection
 
     forAll(rho, celli)
     {
-        const scalar Eni = En[celli];
-        De[celli] = DeN.value(Eni);
-        mue[celli] = mueN.value(Eni);
-        Te[celli] = TeOfEn.value(Eni);
+        const scalar Eni = En[celli] * EnFac;
+
+        De[celli] = DeN.value(Eni) * DeNFac;
+        mue[celli] = mueN.value(Eni) * mueNFac;
+        Te[celli] = TeOfEn.value(Eni) * TeFac;
     }
     De.correctBoundaryConditions();
     mue.correctBoundaryConditions();
diff --git a/applications/solvers/combustion/plasmaReactingFoam/createFields.H b/applications/solvers/combustion/plasmaReactingFoam/createFields.H
index 8fe9826e..994759cc 100644
--- a/applications/solvers/combustion/plasmaReactingFoam/createFields.H
+++ b/applications/solvers/combustion/plasmaReactingFoam/createFields.H
@@ -17,9 +17,30 @@ dimensionedScalar epsilon0
     physicalProperties.lookup("epsilon0")
 );
 
+// Convert E/n in SI unit to table unit.
+// Default V m^2 => Td
+scalar EnFac (
+        physicalProperties.lookupOrDefault("EnFac", 1.0e21)
+        );
+
 CSV<scalar> mueN("mueN", physicalProperties);
+// Convert mu_e * n_g value from table into SI unit.
+scalar mueNFac (
+        physicalProperties.lookupOrDefault("mueNFac", 1.0)
+        );
+
 CSV<scalar> DeN("DeN", physicalProperties);
+// Convert D_e * n_g value from table into SI unit.
+scalar DeNFac (
+        physicalProperties.lookupOrDefault("DeNFac", 1.0)
+        );
+
 CSV<scalar> TeOfEn("TeOfEn", physicalProperties);
+// Convert T_e value from table into SI unit.
+scalar TeFac (
+        physicalProperties.lookupOrDefault("TeFac", 1.0)
+        );
+
 
 Info<< "Reading field Phi\n" << endl;
 volScalarField Phi