diff --git a/applications/solvers/combustion/chemFoam/plasmaChemFoam/createBaseFields.H b/applications/solvers/combustion/chemFoam/plasmaChemFoam/createBaseFields.H
index 0762f270..471377b9 100644
--- a/applications/solvers/combustion/chemFoam/plasmaChemFoam/createBaseFields.H
+++ b/applications/solvers/combustion/chemFoam/plasmaChemFoam/createBaseFields.H
@@ -53,5 +53,39 @@ Info<< "Creating base fields for time " << runTime.timeName() << endl;
     );
 
     T.write();
+
+    volScalarField ne
+    (
+        IOobject
+        (
+            "ne",
+            runTime.timeName(),
+            mesh,
+            IOobject::READ_IF_PRESENT,
+            IOobject::NO_WRITE,
+            false
+        ),
+        mesh,
+        dimensionedScalar("ne", dimTemperature, ne0)
+    );
+
+    ne.write();
+
+    volScalarField Te
+    (
+        IOobject
+        (
+            "Te",
+            runTime.timeName(),
+            mesh,
+            IOobject::READ_IF_PRESENT,
+            IOobject::NO_WRITE,
+            false
+        ),
+        mesh,
+        dimensionedScalar("Te", dimTemperature, Te0)
+    );
+
+    Te.write();
 }
 
diff --git a/applications/solvers/combustion/chemFoam/plasmaChemFoam/createFields.H b/applications/solvers/combustion/chemFoam/plasmaChemFoam/createFields.H
index 8855e450..b9d99ef5 100644
--- a/applications/solvers/combustion/chemFoam/plasmaChemFoam/createFields.H
+++ b/applications/solvers/combustion/chemFoam/plasmaChemFoam/createFields.H
@@ -20,9 +20,40 @@
 
     scalar p0 = readScalar(initialConditions.lookup("p"));
     scalar T0 = readScalar(initialConditions.lookup("T"));
+    scalar ne0 = readScalar(initialConditions.lookup("ne"));
+    scalar Te0 = readScalar(initialConditions.lookup("Te"));
+
+    // read CSV
+    CSV<scalar> EnProfile("EnProfile", initialConditions);
 
     #include "createBaseFields.H"
 
+    volScalarField ne
+    (
+        IOobject
+        (
+            "ne",
+            runTime.timeName(),
+            mesh,
+            IOobject::MUST_READ,
+            IOobject::AUTO_WRITE
+        ),
+        mesh
+    );
+
+    volScalarField Te
+    (
+        IOobject
+        (
+            "Te",
+            runTime.timeName(),
+            mesh,
+            IOobject::MUST_READ,
+            IOobject::AUTO_WRITE
+        ),
+        mesh
+    );
+
     Info<< nl << "Reading thermophysicalProperties" << endl;
     autoPtr<psiChemistryModel> pChemistry(psiChemistryModel::New(mesh));
 
diff --git a/applications/solvers/combustion/chemFoam/plasmaChemFoam/plasmaChemFoam.C b/applications/solvers/combustion/chemFoam/plasmaChemFoam/plasmaChemFoam.C
index 4183045a..e41ba4ea 100644
--- a/applications/solvers/combustion/chemFoam/plasmaChemFoam/plasmaChemFoam.C
+++ b/applications/solvers/combustion/chemFoam/plasmaChemFoam/plasmaChemFoam.C
@@ -43,6 +43,8 @@ Description
 #include "basicMultiComponentMixture.H"
 #include "cellModeller.H"
 
+#include "CSV.H"
+
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 int main(int argc, char *argv[])