diff --git a/applications/solvers/combustion/plasmaReactingFoam/createFields.H b/applications/solvers/combustion/plasmaReactingFoam/createFields.H
index 8eb34ed6..794066fb 100644
--- a/applications/solvers/combustion/plasmaReactingFoam/createFields.H
+++ b/applications/solvers/combustion/plasmaReactingFoam/createFields.H
@@ -265,10 +265,11 @@ volScalarField rhoq
         "rhoq",
         runTime.timeName(),
         mesh,
-        IOobject::NO_READ,
+        IOobject::MUST_READ,
         IOobject::AUTO_WRITE
     ),
-    (rho * qc) - (eCharge * ne)
+    mesh
+    // (rho * qc) - (eCharge * ne)
 );
 
 Info<< "Creating field gas number density\n" << endl;