From 4cadaeb3ccb56a5fcb8294ffa7418360d0bc866a Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Sun, 8 May 2016 20:57:08 +0100
Subject: [PATCH] applications/solvers/multiphase: Added -postProcess option
 See also commit 30e2f912e54a54a370509685197509e846ef3a9a

---
 .../simpleFoam/SRFSimpleFoam/SRFSimpleFoam.C  |  2 +-
 .../{createControls.H => createControl.H}     |  0
 .../porousSimpleFoam/porousSimpleFoam.C       |  2 +-
 .../incompressible/simpleFoam/simpleFoam.C    |  2 +-
 .../cavitatingDyMFoam/cavitatingDyMFoam.C     | 10 +++---
 .../cavitatingFoam/cavitatingFoam.C           | 10 +++---
 .../createControl.H}                          |  0
 .../multiphase/cavitatingFoam/createFields.H  |  2 ++
 .../compressibleInterDyMFoam.C                | 11 +++++--
 .../compressibleInterFoam.C                   | 11 +++++--
 .../createControl.H}                          |  0
 .../compressibleInterFoam/createFields.H      |  4 ---
 .../Make/options                              |  1 +
 .../compressibleMultiphaseInterFoam.C         |  9 ++++--
 .../createControl.H}                          |  0
 .../createFields.H                            |  3 --
 .../multiphase/driftFluxFoam/createControl.H  |  1 +
 .../multiphase/driftFluxFoam/createFields.H   |  6 ----
 .../multiphase/driftFluxFoam/driftFluxFoam.C  | 12 ++++---
 .../multiphase/interFoam/createControl.H      |  1 +
 .../interFoam/interDyMFoam/interDyMFoam.C     |  2 +-
 .../solvers/multiphase/interFoam/interFoam.C  |  2 +-
 .../interMixingFoam/interMixingFoam.C         |  2 +-
 .../interPhaseChangeFoam/createControl.H      |  1 +
 .../interPhaseChangeDyMFoam.C                 |  2 +-
 .../interPhaseChangeFoam.C                    |  2 +-
 .../multiphaseEulerFoam/createControl.H       |  1 +
 .../multiphaseEulerFoam/createFields.H        | 15 ++-------
 .../multiphaseEulerFoam/multiphaseEulerFoam.C | 20 +++++++++---
 .../multiphaseInterDyMFoam.C                  |  2 +-
 .../multiphaseInterFoam/multiphaseInterFoam.C |  2 +-
 .../potentialFreeSurfaceFoam/createControl.H  |  1 +
 .../potentialFreeSurfaceDyMFoam.C             |  2 +-
 .../potentialFreeSurfaceFoam.C                |  2 +-
 .../createControl.H                           |  1 +
 .../createFieldRefs.H                         |  4 +++
 .../createFields.H                            |  5 ---
 .../reactingMultiphaseEulerFoam.C             |  8 ++---
 .../reactingTwoPhaseEulerFoam/createControl.H |  1 +
 .../createFieldRefs.H                         | 28 ++++++++++++++++
 .../reactingTwoPhaseEulerFoam/createFields.H  | 32 +------------------
 .../reactingTwoPhaseEulerFoam.C               |  8 ++---
 .../twoLiquidMixingFoam/Make/options          |  1 +
 .../twoLiquidMixingFoam/createControl.H       |  1 +
 .../twoLiquidMixingFoam/twoLiquidMixingFoam.C |  6 ++--
 .../twoPhaseEulerFoam/createControl.H         |  1 +
 .../twoPhaseEulerFoam/createFieldRefs.H       | 20 ++++++++++++
 .../twoPhaseEulerFoam/createFields.H          | 24 ++------------
 .../twoPhaseEulerFoam/twoPhaseEulerFoam.C     |  9 +++---
 .../functionObjectList/postProcess.H          | 16 +++++-----
 50 files changed, 160 insertions(+), 148 deletions(-)
 rename applications/solvers/incompressible/simpleFoam/{createControls.H => createControl.H} (100%)
 rename applications/solvers/multiphase/{interFoam/createControls.H => cavitatingFoam/createControl.H} (100%)
 rename applications/solvers/multiphase/{interPhaseChangeFoam/createControls.H => compressibleInterFoam/createControl.H} (100%)
 rename applications/solvers/multiphase/{potentialFreeSurfaceFoam/createControls.H => compressibleMultiphaseInterFoam/createControl.H} (100%)
 create mode 100644 applications/solvers/multiphase/driftFluxFoam/createControl.H
 create mode 100644 applications/solvers/multiphase/interFoam/createControl.H
 create mode 100644 applications/solvers/multiphase/interPhaseChangeFoam/createControl.H
 create mode 100644 applications/solvers/multiphase/multiphaseEulerFoam/createControl.H
 create mode 100644 applications/solvers/multiphase/potentialFreeSurfaceFoam/createControl.H
 create mode 100644 applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/createControl.H
 create mode 100644 applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/createFieldRefs.H
 create mode 100644 applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/createControl.H
 create mode 100644 applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/createFieldRefs.H
 create mode 100644 applications/solvers/multiphase/twoLiquidMixingFoam/createControl.H
 create mode 100644 applications/solvers/multiphase/twoPhaseEulerFoam/createControl.H
 create mode 100644 applications/solvers/multiphase/twoPhaseEulerFoam/createFieldRefs.H

diff --git a/applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/SRFSimpleFoam.C b/applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/SRFSimpleFoam.C
index f7c723a5ce..e3940b7f0f 100644
--- a/applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/SRFSimpleFoam.C
+++ b/applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/SRFSimpleFoam.C
@@ -46,7 +46,7 @@ int main(int argc, char *argv[])
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-    #include "createControls.H"
+    #include "createControl.H"
     #include "createFields.H"
     #include "createFvOptions.H"
     #include "initContinuityErrs.H"
diff --git a/applications/solvers/incompressible/simpleFoam/createControls.H b/applications/solvers/incompressible/simpleFoam/createControl.H
similarity index 100%
rename from applications/solvers/incompressible/simpleFoam/createControls.H
rename to applications/solvers/incompressible/simpleFoam/createControl.H
diff --git a/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/porousSimpleFoam.C b/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/porousSimpleFoam.C
index 4bf481946c..9c30d039b2 100644
--- a/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/porousSimpleFoam.C
+++ b/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/porousSimpleFoam.C
@@ -48,7 +48,7 @@ int main(int argc, char *argv[])
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-    #include "createControls.H"
+    #include "createControl.H"
     #include "createFields.H"
     #include "createPorousZones.H"
     #include "createFvOptions.H"
diff --git a/applications/solvers/incompressible/simpleFoam/simpleFoam.C b/applications/solvers/incompressible/simpleFoam/simpleFoam.C
index f97f8a79b5..28af7112b4 100644
--- a/applications/solvers/incompressible/simpleFoam/simpleFoam.C
+++ b/applications/solvers/incompressible/simpleFoam/simpleFoam.C
@@ -44,7 +44,7 @@ int main(int argc, char *argv[])
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-    #include "createControls.H"
+    #include "createControl.H"
     #include "createFields.H"
     #include "createFvOptions.H"
     #include "initContinuityErrs.H"
diff --git a/applications/solvers/multiphase/cavitatingFoam/cavitatingDyMFoam/cavitatingDyMFoam.C b/applications/solvers/multiphase/cavitatingFoam/cavitatingDyMFoam/cavitatingDyMFoam.C
index 58c1483b24..3220508117 100644
--- a/applications/solvers/multiphase/cavitatingFoam/cavitatingDyMFoam/cavitatingDyMFoam.C
+++ b/applications/solvers/multiphase/cavitatingFoam/cavitatingDyMFoam/cavitatingDyMFoam.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2012-2015 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2012-2016 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -44,14 +44,12 @@ Description
 
 int main(int argc, char *argv[])
 {
-    #include "setRootCase.H"
+    #include "postProcess.H"
 
+    #include "setRootCase.H"
     #include "createTime.H"
     #include "createDynamicFvMesh.H"
-
-    pimpleControl pimple(mesh);
-
-    #include "readThermodynamicProperties.H"
+    #include "createControl.H"
     #include "createControls.H"
     #include "createFields.H"
     #include "createUf.H"
diff --git a/applications/solvers/multiphase/cavitatingFoam/cavitatingFoam.C b/applications/solvers/multiphase/cavitatingFoam/cavitatingFoam.C
index e40105aab7..4f20f53313 100644
--- a/applications/solvers/multiphase/cavitatingFoam/cavitatingFoam.C
+++ b/applications/solvers/multiphase/cavitatingFoam/cavitatingFoam.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2015 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -42,14 +42,12 @@ Description
 
 int main(int argc, char *argv[])
 {
-    #include "setRootCase.H"
+    #include "postProcess.H"
 
+    #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
-    pimpleControl pimple(mesh);
-
-    #include "readThermodynamicProperties.H"
+    #include "createControl.H"
     #include "createControls.H"
     #include "createFields.H"
     #include "CourantNo.H"
diff --git a/applications/solvers/multiphase/interFoam/createControls.H b/applications/solvers/multiphase/cavitatingFoam/createControl.H
similarity index 100%
rename from applications/solvers/multiphase/interFoam/createControls.H
rename to applications/solvers/multiphase/cavitatingFoam/createControl.H
diff --git a/applications/solvers/multiphase/cavitatingFoam/createFields.H b/applications/solvers/multiphase/cavitatingFoam/createFields.H
index f124d8972c..0108b8c56a 100644
--- a/applications/solvers/multiphase/cavitatingFoam/createFields.H
+++ b/applications/solvers/multiphase/cavitatingFoam/createFields.H
@@ -1,3 +1,5 @@
+#include "readThermodynamicProperties.H"
+
 Info<< "Reading field p\n" << endl;
 volScalarField p
 (
diff --git a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
index 521628fa10..6f7ee7ffce 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
+++ b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
@@ -52,19 +52,24 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createDynamicFvMesh.H"
     #include "initContinuityErrs.H"
-
-    pimpleControl pimple(mesh);
-
+    #include "createControl.H"
     #include "createFields.H"
     #include "createUf.H"
     #include "createControls.H"
     #include "CourantNo.H"
     #include "setInitialDeltaT.H"
 
+    volScalarField& p = mixture.p();
+    volScalarField& T = mixture.T();
+    const volScalarField& psi1 = mixture.thermo1().psi();
+    const volScalarField& psi2 = mixture.thermo2().psi();
+
     turbulence->validate();
 
     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
index 3a135b735e..150cc8bdfd 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
+++ b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
@@ -49,17 +49,22 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
-    pimpleControl pimple(mesh);
-
+    #include "createControl.H"
     #include "createTimeControls.H"
     #include "createFields.H"
     #include "CourantNo.H"
     #include "setInitialDeltaT.H"
 
+    volScalarField& p = mixture.p();
+    volScalarField& T = mixture.T();
+    const volScalarField& psi1 = mixture.thermo1().psi();
+    const volScalarField& psi2 = mixture.thermo2().psi();
+
     turbulence->validate();
 
     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/applications/solvers/multiphase/interPhaseChangeFoam/createControls.H b/applications/solvers/multiphase/compressibleInterFoam/createControl.H
similarity index 100%
rename from applications/solvers/multiphase/interPhaseChangeFoam/createControls.H
rename to applications/solvers/multiphase/compressibleInterFoam/createControl.H
diff --git a/applications/solvers/multiphase/compressibleInterFoam/createFields.H b/applications/solvers/multiphase/compressibleInterFoam/createFields.H
index 4862876828..9a49e504c6 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/createFields.H
+++ b/applications/solvers/multiphase/compressibleInterFoam/createFields.H
@@ -36,12 +36,8 @@ volScalarField& alpha2(mixture.alpha2());
 
 Info<< "Reading thermophysical properties\n" << endl;
 
-volScalarField& p = mixture.p();
-volScalarField& T = mixture.T();
 volScalarField& rho1 = mixture.thermo1().rho();
-const volScalarField& psi1 = mixture.thermo1().psi();
 volScalarField& rho2 = mixture.thermo2().rho();
-const volScalarField& psi2 = mixture.thermo2().psi();
 
 volScalarField rho
 (
diff --git a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/Make/options b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/Make/options
index b237353985..f3010e0e57 100644
--- a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/Make/options
+++ b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/Make/options
@@ -1,4 +1,5 @@
 EXE_INC = \
+    -I. \
     -I../interFoam \
     -ImultiphaseMixtureThermo/lnInclude \
     -I$(LIB_SRC)/transportModels/compressible/lnInclude \
diff --git a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/compressibleMultiphaseInterFoam.C b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/compressibleMultiphaseInterFoam.C
index f3d9bf8314..3dc33d6a64 100644
--- a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/compressibleMultiphaseInterFoam.C
+++ b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/compressibleMultiphaseInterFoam.C
@@ -44,17 +44,20 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
-    pimpleControl pimple(mesh);
-
+    #include "createControl.H"
     #include "createTimeControls.H"
     #include "createFields.H"
     #include "CourantNo.H"
     #include "setInitialDeltaT.H"
 
+    volScalarField& p = mixture.p();
+    volScalarField& T = mixture.T();
+
     turbulence->validate();
 
     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/applications/solvers/multiphase/potentialFreeSurfaceFoam/createControls.H b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/createControl.H
similarity index 100%
rename from applications/solvers/multiphase/potentialFreeSurfaceFoam/createControls.H
rename to applications/solvers/multiphase/compressibleMultiphaseInterFoam/createControl.H
diff --git a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/createFields.H b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/createFields.H
index 8ac81a5195..dcb89c35b9 100644
--- a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/createFields.H
+++ b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/createFields.H
@@ -31,9 +31,6 @@ volVectorField U
 Info<< "Constructing multiphaseMixtureThermo\n" << endl;
 multiphaseMixtureThermo mixture(U, phi);
 
-volScalarField& p = mixture.p();
-volScalarField& T = mixture.T();
-
 volScalarField rho
 (
     IOobject
diff --git a/applications/solvers/multiphase/driftFluxFoam/createControl.H b/applications/solvers/multiphase/driftFluxFoam/createControl.H
new file mode 100644
index 0000000000..cf82a59b35
--- /dev/null
+++ b/applications/solvers/multiphase/driftFluxFoam/createControl.H
@@ -0,0 +1 @@
+pimpleControl pimple(mesh);
diff --git a/applications/solvers/multiphase/driftFluxFoam/createFields.H b/applications/solvers/multiphase/driftFluxFoam/createFields.H
index 65d3d38eae..f1b4e78645 100644
--- a/applications/solvers/multiphase/driftFluxFoam/createFields.H
+++ b/applications/solvers/multiphase/driftFluxFoam/createFields.H
@@ -33,10 +33,6 @@ Info<< "Reading incompressibleTwoPhaseInteractingMixture\n" << endl;
 incompressibleTwoPhaseInteractingMixture mixture(U, phi);
 
 volScalarField& alpha1(mixture.alpha1());
-volScalarField& alpha2(mixture.alpha2());
-
-const dimensionedScalar& rho1 = mixture.rhod();
-const dimensionedScalar& rho2 = mixture.rhoc();
 
 // Mixture density
 volScalarField rho
@@ -77,8 +73,6 @@ autoPtr<relativeVelocityModel> UdmModelPtr
     )
 );
 
-relativeVelocityModel& UdmModel(UdmModelPtr());
-
 
 // Construct compressible turbulence model
 autoPtr
diff --git a/applications/solvers/multiphase/driftFluxFoam/driftFluxFoam.C b/applications/solvers/multiphase/driftFluxFoam/driftFluxFoam.C
index c0a2be76a7..990c11b8b0 100644
--- a/applications/solvers/multiphase/driftFluxFoam/driftFluxFoam.C
+++ b/applications/solvers/multiphase/driftFluxFoam/driftFluxFoam.C
@@ -49,19 +49,23 @@ Description
 
 int main(int argc, char *argv[])
 {
-    #include "setRootCase.H"
+    #include "postProcess.H"
 
+    #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
-    pimpleControl pimple(mesh);
-
+    #include "createControl.H"
     #include "createTimeControls.H"
     #include "createFields.H"
     #include "createMRF.H"
     #include "createFvOptions.H"
     #include "initContinuityErrs.H"
 
+    volScalarField& alpha2(mixture.alpha2());
+    const dimensionedScalar& rho1 = mixture.rhod();
+    const dimensionedScalar& rho2 = mixture.rhoc();
+    relativeVelocityModel& UdmModel(UdmModelPtr());
+
     turbulence->validate();
 
     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/applications/solvers/multiphase/interFoam/createControl.H b/applications/solvers/multiphase/interFoam/createControl.H
new file mode 100644
index 0000000000..cf82a59b35
--- /dev/null
+++ b/applications/solvers/multiphase/interFoam/createControl.H
@@ -0,0 +1 @@
+pimpleControl pimple(mesh);
diff --git a/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C b/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
index 15707219a4..07206d3c52 100644
--- a/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
+++ b/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
@@ -57,7 +57,7 @@ int main(int argc, char *argv[])
     #include "createTime.H"
     #include "createDynamicFvMesh.H"
     #include "initContinuityErrs.H"
-    #include "createControls.H"
+    #include "createControl.H"
     #include "createTimeControls.H"
     #include "createDyMControls.H"
     #include "createRDeltaT.H"
diff --git a/applications/solvers/multiphase/interFoam/interFoam.C b/applications/solvers/multiphase/interFoam/interFoam.C
index 8e8c78fff4..c5c4be315b 100644
--- a/applications/solvers/multiphase/interFoam/interFoam.C
+++ b/applications/solvers/multiphase/interFoam/interFoam.C
@@ -60,7 +60,7 @@ int main(int argc, char *argv[])
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-    #include "createControls.H"
+    #include "createControl.H"
     #include "createTimeControls.H"
     #include "createRDeltaT.H"
     #include "initContinuityErrs.H"
diff --git a/applications/solvers/multiphase/interFoam/interMixingFoam/interMixingFoam.C b/applications/solvers/multiphase/interFoam/interMixingFoam/interMixingFoam.C
index 8b59f1ad4f..28d8d6d4bf 100644
--- a/applications/solvers/multiphase/interFoam/interMixingFoam/interMixingFoam.C
+++ b/applications/solvers/multiphase/interFoam/interMixingFoam/interMixingFoam.C
@@ -50,7 +50,7 @@ int main(int argc, char *argv[])
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-    #include "createControls.H"
+    #include "createControl.H"
     #include "createTimeControls.H"
     #include "createRDeltaT.H"
     #include "initContinuityErrs.H"
diff --git a/applications/solvers/multiphase/interPhaseChangeFoam/createControl.H b/applications/solvers/multiphase/interPhaseChangeFoam/createControl.H
new file mode 100644
index 0000000000..cf82a59b35
--- /dev/null
+++ b/applications/solvers/multiphase/interPhaseChangeFoam/createControl.H
@@ -0,0 +1 @@
+pimpleControl pimple(mesh);
diff --git a/applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeDyMFoam/interPhaseChangeDyMFoam.C b/applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeDyMFoam/interPhaseChangeDyMFoam.C
index d12aaad263..f0d9f09917 100644
--- a/applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeDyMFoam/interPhaseChangeDyMFoam.C
+++ b/applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeDyMFoam/interPhaseChangeDyMFoam.C
@@ -61,7 +61,7 @@ int main(int argc, char *argv[])
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createDynamicFvMesh.H"
-    #include "createControls.H"
+    #include "createControl.H"
     #include "createTimeControls.H"
     #include "../interFoam/interDyMFoam/createDyMControls.H"
     #include "initContinuityErrs.H"
diff --git a/applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeFoam.C b/applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeFoam.C
index 2a2d5cbc38..d33efa3fa3 100644
--- a/applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeFoam.C
+++ b/applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeFoam.C
@@ -58,7 +58,7 @@ int main(int argc, char *argv[])
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-    #include "createControls.H"
+    #include "createControl.H"
     #include "createFields.H"
     #include "createFvOptions.H"
     #include "createTimeControls.H"
diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/createControl.H b/applications/solvers/multiphase/multiphaseEulerFoam/createControl.H
new file mode 100644
index 0000000000..cf82a59b35
--- /dev/null
+++ b/applications/solvers/multiphase/multiphaseEulerFoam/createControl.H
@@ -0,0 +1 @@
+pimpleControl pimple(mesh);
diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/createFields.H b/applications/solvers/multiphase/multiphaseEulerFoam/createFields.H
index 212160b83e..a4b9838335 100644
--- a/applications/solvers/multiphase/multiphaseEulerFoam/createFields.H
+++ b/applications/solvers/multiphase/multiphaseEulerFoam/createFields.H
@@ -51,18 +51,6 @@ forAllIter(PtrDictionary<phaseModel>, fluid.phases(), iter)
     phi += fvc::interpolate(alpha)*phase.phi();
 }
 
-scalar slamDampCoeff
-(
-    fluid.lookupOrDefault<scalar>("slamDampCoeff", 1)
-);
-
-dimensionedScalar maxSlamVelocity
-(
-    "maxSlamVelocity",
-    dimVelocity,
-    fluid.lookupOrDefault<scalar>("maxSlamVelocity", GREAT)
-);
-
 
 volScalarField rho
 (
@@ -114,3 +102,6 @@ setRefCell
     pRefValue
 );
 mesh.setFluxRequired(p_rgh.name());
+
+
+#include "createMRFZones.H"
diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam.C b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam.C
index 168c00be7e..c27bf597ed 100644
--- a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam.C
+++ b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam.C
@@ -45,21 +45,31 @@ Description
 
 int main(int argc, char *argv[])
 {
-    #include "setRootCase.H"
+    #include "postProcess.H"
 
+    #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
-    pimpleControl pimple(mesh);
-
+    #include "createControl.H"
     #include "createFields.H"
-    #include "createMRFZones.H"
     #include "initContinuityErrs.H"
     #include "createTimeControls.H"
     #include "correctPhi.H"
     #include "CourantNo.H"
     #include "setInitialDeltaT.H"
 
+    scalar slamDampCoeff
+    (
+        fluid.lookupOrDefault<scalar>("slamDampCoeff", 1)
+    );
+
+    dimensionedScalar maxSlamVelocity
+    (
+        "maxSlamVelocity",
+        dimVelocity,
+        fluid.lookupOrDefault<scalar>("maxSlamVelocity", GREAT)
+    );
+
     turbulence->validate();
 
     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterDyMFoam/multiphaseInterDyMFoam.C b/applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterDyMFoam/multiphaseInterDyMFoam.C
index fa3969d855..9cff08c31f 100644
--- a/applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterDyMFoam/multiphaseInterDyMFoam.C
+++ b/applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterDyMFoam/multiphaseInterDyMFoam.C
@@ -50,7 +50,7 @@ int main(int argc, char *argv[])
     #include "createTime.H"
     #include "createDynamicFvMesh.H"
     #include "initContinuityErrs.H"
-    #include "createControls.H"
+    #include "createControl.H"
     #include "createTimeControls.H"
     #include "createDyMControls.H"
     #include "createFields.H"
diff --git a/applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterFoam.C b/applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterFoam.C
index 933a396fe1..68ea578ff9 100644
--- a/applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterFoam.C
+++ b/applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterFoam.C
@@ -48,7 +48,7 @@ int main(int argc, char *argv[])
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-    #include "createControls.H"
+    #include "createControl.H"
     #include "initContinuityErrs.H"
     #include "createFields.H"
     #include "createMRF.H"
diff --git a/applications/solvers/multiphase/potentialFreeSurfaceFoam/createControl.H b/applications/solvers/multiphase/potentialFreeSurfaceFoam/createControl.H
new file mode 100644
index 0000000000..cf82a59b35
--- /dev/null
+++ b/applications/solvers/multiphase/potentialFreeSurfaceFoam/createControl.H
@@ -0,0 +1 @@
+pimpleControl pimple(mesh);
diff --git a/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceDyMFoam/potentialFreeSurfaceDyMFoam.C b/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceDyMFoam/potentialFreeSurfaceDyMFoam.C
index dea1083919..9b5b6a44ed 100644
--- a/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceDyMFoam/potentialFreeSurfaceDyMFoam.C
+++ b/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceDyMFoam/potentialFreeSurfaceDyMFoam.C
@@ -55,7 +55,7 @@ int main(int argc, char *argv[])
     #include "createTime.H"
     #include "createDynamicFvMesh.H"
     #include "initContinuityErrs.H"
-    #include "createControls.H"
+    #include "createControl.H"
     #include "createTimeControls.H"
     #include "createDyMControls.H"
     #include "createFields.H"
diff --git a/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceFoam.C b/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceFoam.C
index b23fb73cf0..615c76bd52 100644
--- a/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceFoam.C
+++ b/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceFoam.C
@@ -49,7 +49,7 @@ int main(int argc, char *argv[])
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-    #include "createControls.H"
+    #include "createControl.H"
     #include "createTimeControls.H"
     #include "createFields.H"
     #include "createMRF.H"
diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/createControl.H b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/createControl.H
new file mode 100644
index 0000000000..cf82a59b35
--- /dev/null
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/createControl.H
@@ -0,0 +1 @@
+pimpleControl pimple(mesh);
diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/createFieldRefs.H b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/createFieldRefs.H
new file mode 100644
index 0000000000..8a629b017d
--- /dev/null
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/createFieldRefs.H
@@ -0,0 +1,4 @@
+surfaceScalarField& phi = fluid.phi();
+
+const IOMRFZoneList& MRF = fluid.MRF();
+fv::options& fvOptions = fluid.fvOptions();
diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/createFields.H b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/createFields.H
index 6e7b6aacc8..8eafafb2fd 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/createFields.H
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/createFields.H
@@ -10,8 +10,6 @@ autoPtr<multiphaseSystem> fluidPtr
 multiphaseSystem& fluid = fluidPtr();
 multiphaseSystem::phaseModelList& phases = fluid.phases();
 
-surfaceScalarField& phi = fluid.phi();
-
 dimensionedScalar pMin
 (
     "pMin",
@@ -48,6 +46,3 @@ setRefCell
     pRefValue
 );
 mesh.setFluxRequired(p_rgh.name());
-
-const IOMRFZoneList& MRF = fluid.MRF();
-fv::options& fvOptions = fluid.fvOptions();
diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/reactingMultiphaseEulerFoam.C b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/reactingMultiphaseEulerFoam.C
index a98129cc7e..d1b1619567 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/reactingMultiphaseEulerFoam.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/reactingMultiphaseEulerFoam.C
@@ -43,16 +43,16 @@ Description
 
 int main(int argc, char *argv[])
 {
-    #include "setRootCase.H"
+    #include "postProcess.H"
 
+    #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
-    pimpleControl pimple(mesh);
-
+    #include "createControl.H"
     #include "createTimeControls.H"
     #include "createRDeltaT.H"
     #include "createFields.H"
+    #include "createFieldRefs.H"
 
     if (!LTS)
     {
diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/createControl.H b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/createControl.H
new file mode 100644
index 0000000000..cf82a59b35
--- /dev/null
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/createControl.H
@@ -0,0 +1 @@
+pimpleControl pimple(mesh);
diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/createFieldRefs.H b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/createFieldRefs.H
new file mode 100644
index 0000000000..96cbabfe8e
--- /dev/null
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/createFieldRefs.H
@@ -0,0 +1,28 @@
+phaseModel& phase1 = fluid.phase1();
+phaseModel& phase2 = fluid.phase2();
+
+volScalarField& alpha1 = phase1;
+volScalarField& alpha2 = phase2;
+
+volVectorField& U1 = phase1.U();
+surfaceScalarField& phi1 = phase1.phi();
+surfaceScalarField& alphaPhi1 = phase1.alphaPhi();
+surfaceScalarField& alphaRhoPhi1 = phase1.alphaRhoPhi();
+
+volVectorField& U2 = phase2.U();
+surfaceScalarField& phi2 = phase2.phi();
+surfaceScalarField& alphaPhi2 = phase2.alphaPhi();
+surfaceScalarField& alphaRhoPhi2 = phase2.alphaRhoPhi();
+surfaceScalarField& phi = fluid.phi();
+
+rhoThermo& thermo1 = phase1.thermo();
+rhoThermo& thermo2 = phase2.thermo();
+
+volScalarField& rho1 = thermo1.rho();
+const volScalarField& psi1 = thermo1.psi();
+
+volScalarField& rho2 = thermo2.rho();
+const volScalarField& psi2 = thermo2.psi();
+
+const IOMRFZoneList& MRF = fluid.MRF();
+fv::options& fvOptions = fluid.fvOptions();
diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/createFields.H b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/createFields.H
index 75d7de420e..932c05cb0d 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/createFields.H
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/createFields.H
@@ -9,24 +9,6 @@ autoPtr<twoPhaseSystem> fluidPtr
 );
 twoPhaseSystem& fluid = fluidPtr();
 
-phaseModel& phase1 = fluid.phase1();
-phaseModel& phase2 = fluid.phase2();
-
-volScalarField& alpha1 = phase1;
-volScalarField& alpha2 = phase2;
-
-volVectorField& U1 = phase1.U();
-surfaceScalarField& phi1 = phase1.phi();
-surfaceScalarField& alphaPhi1 = phase1.alphaPhi();
-surfaceScalarField& alphaRhoPhi1 = phase1.alphaRhoPhi();
-
-volVectorField& U2 = phase2.U();
-surfaceScalarField& phi2 = phase2.phi();
-surfaceScalarField& alphaPhi2 = phase2.alphaPhi();
-surfaceScalarField& alphaRhoPhi2 = phase2.alphaRhoPhi();
-
-surfaceScalarField& phi = fluid.phi();
-
 dimensionedScalar pMin
 (
     "pMin",
@@ -36,16 +18,7 @@ dimensionedScalar pMin
 
 #include "gh.H"
 
-rhoThermo& thermo1 = phase1.thermo();
-rhoThermo& thermo2 = phase2.thermo();
-
-volScalarField& p = thermo1.p();
-
-volScalarField& rho1 = thermo1.rho();
-const volScalarField& psi1 = thermo1.psi();
-
-volScalarField& rho2 = thermo2.rho();
-const volScalarField& psi2 = thermo2.psi();
+volScalarField& p = fluid.phase1().thermo().p();
 
 Info<< "Reading field p_rgh\n" << endl;
 volScalarField p_rgh
@@ -72,6 +45,3 @@ setRefCell
     pRefValue
 );
 mesh.setFluxRequired(p_rgh.name());
-
-const IOMRFZoneList& MRF = fluid.MRF();
-fv::options& fvOptions = fluid.fvOptions();
diff --git a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/reactingTwoPhaseEulerFoam.C b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/reactingTwoPhaseEulerFoam.C
index 0e5542672d..032b9bce43 100644
--- a/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/reactingTwoPhaseEulerFoam.C
+++ b/applications/solvers/multiphase/reactingEulerFoam/reactingTwoPhaseEulerFoam/reactingTwoPhaseEulerFoam.C
@@ -44,16 +44,16 @@ Description
 
 int main(int argc, char *argv[])
 {
-    #include "setRootCase.H"
+    #include "postProcess.H"
 
+    #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
-    pimpleControl pimple(mesh);
-
+    #include "createControl.H"
     #include "createTimeControls.H"
     #include "createRDeltaT.H"
     #include "createFields.H"
+    #include "createFieldRefs.H"
 
     if (!LTS)
     {
diff --git a/applications/solvers/multiphase/twoLiquidMixingFoam/Make/options b/applications/solvers/multiphase/twoLiquidMixingFoam/Make/options
index 04e5d53a91..4dfc63419d 100644
--- a/applications/solvers/multiphase/twoLiquidMixingFoam/Make/options
+++ b/applications/solvers/multiphase/twoLiquidMixingFoam/Make/options
@@ -1,4 +1,5 @@
 EXE_INC = \
+    -I. \
     -I../interFoam \
     -I$(LIB_SRC)/transportModels/twoPhaseMixture/lnInclude \
     -I$(LIB_SRC)/transportModels \
diff --git a/applications/solvers/multiphase/twoLiquidMixingFoam/createControl.H b/applications/solvers/multiphase/twoLiquidMixingFoam/createControl.H
new file mode 100644
index 0000000000..cf82a59b35
--- /dev/null
+++ b/applications/solvers/multiphase/twoLiquidMixingFoam/createControl.H
@@ -0,0 +1 @@
+pimpleControl pimple(mesh);
diff --git a/applications/solvers/multiphase/twoLiquidMixingFoam/twoLiquidMixingFoam.C b/applications/solvers/multiphase/twoLiquidMixingFoam/twoLiquidMixingFoam.C
index 901be88800..45a6e71469 100644
--- a/applications/solvers/multiphase/twoLiquidMixingFoam/twoLiquidMixingFoam.C
+++ b/applications/solvers/multiphase/twoLiquidMixingFoam/twoLiquidMixingFoam.C
@@ -42,12 +42,12 @@ Description
 
 int main(int argc, char *argv[])
 {
+    #include "postProcess.H"
+
     #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
-    pimpleControl pimple(mesh);
-
+    #include "createControl.H"
     #include "initContinuityErrs.H"
     #include "createFields.H"
     #include "createTimeControls.H"
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/createControl.H b/applications/solvers/multiphase/twoPhaseEulerFoam/createControl.H
new file mode 100644
index 0000000000..cf82a59b35
--- /dev/null
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/createControl.H
@@ -0,0 +1 @@
+pimpleControl pimple(mesh);
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/createFieldRefs.H b/applications/solvers/multiphase/twoPhaseEulerFoam/createFieldRefs.H
new file mode 100644
index 0000000000..71c3e5e39d
--- /dev/null
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/createFieldRefs.H
@@ -0,0 +1,20 @@
+volScalarField& alpha2 = phase2;
+
+surfaceScalarField& phi1 = phase1.phi();
+surfaceScalarField& alphaPhi1 = phase1.alphaPhi();
+surfaceScalarField& alphaRhoPhi1 = phase1.alphaRhoPhi();
+
+surfaceScalarField& phi2 = phase2.phi();
+surfaceScalarField& alphaPhi2 = phase2.alphaPhi();
+surfaceScalarField& alphaRhoPhi2 = phase2.alphaRhoPhi();
+
+surfaceScalarField& phi = fluid.phi();
+
+rhoThermo& thermo1 = phase1.thermo();
+rhoThermo& thermo2 = phase2.thermo();
+
+volScalarField& rho1 = thermo1.rho();
+const volScalarField& psi1 = thermo1.psi();
+
+volScalarField& rho2 = thermo2.rho();
+const volScalarField& psi2 = thermo2.psi();
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/createFields.H b/applications/solvers/multiphase/twoPhaseEulerFoam/createFields.H
index 570d7fd699..b8192fa2ba 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/createFields.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/createFields.H
@@ -9,19 +9,10 @@ phaseModel& phase1 = fluid.phase1();
 phaseModel& phase2 = fluid.phase2();
 
 volScalarField& alpha1 = phase1;
-volScalarField& alpha2 = phase2;
-
 volVectorField& U1 = phase1.U();
-surfaceScalarField& phi1 = phase1.phi();
-surfaceScalarField& alphaPhi1 = phase1.alphaPhi();
-surfaceScalarField& alphaRhoPhi1 = phase1.alphaRhoPhi();
-
 volVectorField& U2 = phase2.U();
-surfaceScalarField& phi2 = phase2.phi();
-surfaceScalarField& alphaPhi2 = phase2.alphaPhi();
-surfaceScalarField& alphaRhoPhi2 = phase2.alphaRhoPhi();
 
-surfaceScalarField& phi = fluid.phi();
+volScalarField& p = phase1.thermo().p();
 
 dimensionedScalar pMin
 (
@@ -32,17 +23,6 @@ dimensionedScalar pMin
 
 #include "gh.H"
 
-rhoThermo& thermo1 = phase1.thermo();
-rhoThermo& thermo2 = phase2.thermo();
-
-volScalarField& p = thermo1.p();
-
-volScalarField& rho1 = thermo1.rho();
-const volScalarField& psi1 = thermo1.psi();
-
-volScalarField& rho2 = thermo2.rho();
-const volScalarField& psi2 = thermo2.psi();
-
 Info<< "Reading field p_rgh\n" << endl;
 volScalarField p_rgh
 (
@@ -100,3 +80,5 @@ volScalarField dpdt
 Info<< "Creating field kinetic energy K\n" << endl;
 volScalarField K1(IOobject::groupName("K", phase1.name()), 0.5*magSqr(U1));
 volScalarField K2(IOobject::groupName("K", phase2.name()), 0.5*magSqr(U2));
+
+#include "createMRF.H"
diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseEulerFoam.C b/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseEulerFoam.C
index 1e6ac19301..d86461087c 100644
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseEulerFoam.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseEulerFoam.C
@@ -41,15 +41,14 @@ Description
 
 int main(int argc, char *argv[])
 {
-    #include "setRootCase.H"
+    #include "postProcess.H"
 
+    #include "setRootCase.H"
     #include "createTime.H"
     #include "createMesh.H"
-
-    pimpleControl pimple(mesh);
-
+    #include "createControl.H"
     #include "createFields.H"
-    #include "createMRF.H"
+    #include "createFieldRefs.H"
     #include "createFvOptions.H"
     #include "createTimeControls.H"
     #include "CourantNos.H"
diff --git a/src/OpenFOAM/db/functionObjects/functionObjectList/postProcess.H b/src/OpenFOAM/db/functionObjects/functionObjectList/postProcess.H
index 3b99c5fa50..2a16ae115a 100644
--- a/src/OpenFOAM/db/functionObjects/functionObjectList/postProcess.H
+++ b/src/OpenFOAM/db/functionObjects/functionObjectList/postProcess.H
@@ -42,12 +42,12 @@ Description
 #define CREATE_MESH createMesh.H
 #endif
 
-#ifndef CREATE_FIELDS_1
-#define CREATE_FIELDS_1 createFields.H
+#ifndef CREATE_FIELDS
+#define CREATE_FIELDS createFields.H
 #endif
 
-#ifndef CREATE_CONTROLS
-#define CREATE_CONTROLS createControls.H
+#ifndef CREATE_CONTROL
+#define CREATE_CONTROL createControl.H
 #endif
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -71,7 +71,7 @@ if (argList::postProcess(argc, argv))
     #include "createTime.H"
     Foam::instantList timeDirs = Foam::timeSelector::select0(runTime, args);
     #include INCLUDE_FILE(CREATE_MESH)
-    #include INCLUDE_FILE(CREATE_CONTROLS)
+    #include INCLUDE_FILE(CREATE_CONTROL)
 
     // Externally stored dictionary for functionObjectList
     // if not constructed from runTime
@@ -96,7 +96,7 @@ if (argList::postProcess(argc, argv))
                 functionObjectList::New(args, runTime, functionObjectsDict);
         }
 
-        #include INCLUDE_FILE(CREATE_FIELDS_1)
+        #include INCLUDE_FILE(CREATE_FIELDS)
 
         #ifdef CREATE_FIELDS_2
         #include INCLUDE_FILE(CREATE_FIELDS_2)
@@ -132,7 +132,7 @@ if (argList::postProcess(argc, argv))
 #undef INCLUDE_FILE2
 
 #undef CREATE_MESH
-#undef CREATE_FIELDS_1
-#undef CREATE_CONTROLS
+#undef CREATE_FIELDS
+#undef CREATE_CONTROL
 
 // ************************************************************************* //
-- 
GitLab