diff --git a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/CourantNos.H b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/CourantNos.H
index 248f4c0602f62d6d273b62df33c41b322e90889d..7ab578dd3f7fa72e84d667898238345e090c8af3 100644
--- a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/CourantNos.H
+++ b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/CourantNos.H
@@ -1,4 +1,4 @@
-    #include "CourantNo.H"
+#include "CourantNo.H"
 
 {
     scalar UrCoNum = 0.5*gMax
diff --git a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/createFields.H b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/createFields.H
index 41eb6e0e746630f887b3405a05d500f1f04de1a7..85bb8be45e15ac1dee46e903b06a723fa4860940 100644
--- a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/createFields.H
+++ b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/createFields.H
@@ -71,3 +71,6 @@
         pRefCell,
         pRefValue
     );
+
+    const IOMRFZoneList& MRF = fluid.MRF();
+    fv::IOoptionList& fvOptions = fluid.fvOptions();
diff --git a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/phaseSystems/phaseSystem/phaseSystem.H b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/phaseSystems/phaseSystem/phaseSystem.H
index 691a663ffbe4d604c1d438780153e40b09237755..81bfd708f8ac9dbe652b250bfd6203ff61d9a514 100644
--- a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/phaseSystems/phaseSystem/phaseSystem.H
+++ b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/phaseSystems/phaseSystem/phaseSystem.H
@@ -41,6 +41,8 @@ SourceFiles
 #include "phaseModel.H"
 #include "phasePair.H"
 #include "orderedPhasePair.H"
+#include "HashPtrTable.H"
+
 #include "IOMRFZoneList.H"
 #include "fvIOoptionList.H"
 
diff --git a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/reactingTwoPhaseEulerFoam.C b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/reactingTwoPhaseEulerFoam.C
index f74e3e93f081252a48dbf30da6697af8120980c7..27826f12a1934a8d4c5745c09e0b0e52bc0b938c 100644
--- a/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/reactingTwoPhaseEulerFoam.C
+++ b/applications/solvers/multiphase/reactingTwoPhaseEulerFoam/reactingTwoPhaseEulerFoam.C
@@ -22,7 +22,7 @@ License
     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
 Application
-    twoPhaseEulerFoam
+    reactingTwoPhaseEulerFoam
 
 Description
     Solver for a system of 2 compressible fluid phases with a common pressure,
@@ -36,9 +36,8 @@ Description
 #include "fvCFD.H"
 #include "twoPhaseSystem.H"
 #include "PhaseCompressibleTurbulenceModel.H"
-#include "pimpleControl.H"
 #include "fixedFluxPressureFvPatchScalarField.H"
-#include "HashPtrTable.H"
+#include "pimpleControl.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -72,9 +71,6 @@ int main(int argc, char *argv[])
 
     #include "pUf/createDDtU.H"
 
-    const IOMRFZoneList& MRF = fluid.MRF();
-    fv::IOoptionList& fvOptions = fluid.fvOptions();
-
     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
     Info<< "\nStarting time loop\n" << endl;