From aa319f188e1ab8355bb1c151f585f066303b5b81 Mon Sep 17 00:00:00 2001
From: ignis <ignis@postech.ac.kr>
Date: Mon, 17 Oct 2016 23:01:38 +0900
Subject: [PATCH] plasmaReactingFoam CSV->DataEntry for flexibility and add
 Boltzmann constant

---
 .../plasmaReactingFoam/createFields.H         | 22 ++++++++++++++++---
 1 file changed, 19 insertions(+), 3 deletions(-)

diff --git a/applications/solvers/combustion/plasmaReactingFoam/createFields.H b/applications/solvers/combustion/plasmaReactingFoam/createFields.H
index 074b0cdd..9c1e98d0 100644
--- a/applications/solvers/combustion/plasmaReactingFoam/createFields.H
+++ b/applications/solvers/combustion/plasmaReactingFoam/createFields.H
@@ -23,24 +23,37 @@ scalar EnFac (
         physicalProperties.lookupOrDefault("EnFac", 1.0e21)
         );
 
-CSV<scalar> mueN("mueN", physicalProperties);
+
+autoPtr< DataEntry< scalar > > pmueN (
+        DataEntry<scalar>::New("mueN", physicalProperties));
+const DataEntry<scalar> &mueN(pmueN());
+
 // Convert mu_e * n_g value from table into SI unit.
 scalar mueNFac (
         physicalProperties.lookupOrDefault("mueNFac", 1.0)
         );
 
-CSV<scalar> DeN("DeN", physicalProperties);
+
+autoPtr< DataEntry< scalar > > pDeN (
+        DataEntry<scalar>::New("DeN", physicalProperties));
+const DataEntry<scalar> &DeN(pDeN());
+
 // Convert D_e * n_g value from table into SI unit.
 scalar DeNFac (
         physicalProperties.lookupOrDefault("DeNFac", 1.0)
         );
 
-CSV<scalar> TeOfEn("TeOfEn", physicalProperties);
+
+autoPtr< DataEntry< scalar > > pTeOfEn(
+        DataEntry<scalar>::New ("TeOfEn", physicalProperties));
+const DataEntry<scalar> &TeOfEn(pTeOfEn());
+
 // Convert T_e value from table into SI unit.
 scalar TeFac (
         physicalProperties.lookupOrDefault("TeFac", 1.0)
         );
 
+
 dictionary wallElectronFlux
 (
     physicalProperties.subDict("wallElectronFlux")
@@ -114,6 +127,9 @@ const dimensionedScalar NA = constant::physicoChemical::NA;
 //- Universal gas constant (default in [J/mol/K])
 const dimensionedScalar R = constant::physicoChemical::R;
 
+//- Boltzmann constant
+const dimensionedScalar kB = constant::physicoChemical::k;
+
 //- Pi
 const scalar pi = constant::mathematical::pi;