solvers: Moved fvOption construction into createFields.H for post-processing

This ensures that the fvOptions are constructed for the -postProcessing option so that functionObjects which process fvOption data operate correctly in this mode.
2017-10-24 14:48:43 +01:00
parent d7723455c7
commit 73e24df8e0
83 changed files with 136 additions and 142 deletions
--- a/applications/solvers/basic/laplacianFoam/createFields.H
+++ b/applications/solvers/basic/laplacianFoam/createFields.H
@ -1,37 +1,39 @@
-    Info<< "Reading field T\n" << endl;
+Info<< "Reading field T\n" << endl;

-    volScalarField T
+volScalarField T
+(
+    IOobject
    (
-        IOobject
-        (
-            "T",
-            runTime.timeName(),
-            mesh,
-            IOobject::MUST_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh
-    );
+        "T",
+        runTime.timeName(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
+    ),
+    mesh
+);


-    Info<< "Reading transportProperties\n" << endl;
+Info<< "Reading transportProperties\n" << endl;

-    IOdictionary transportProperties
+IOdictionary transportProperties
+(
+    IOobject
    (
-        IOobject
-        (
-            "transportProperties",
-            runTime.constant(),
-            mesh,
-            IOobject::MUST_READ_IF_MODIFIED,
-            IOobject::NO_WRITE
-        )
-    );
+        "transportProperties",
+        runTime.constant(),
+        mesh,
+        IOobject::MUST_READ_IF_MODIFIED,
+        IOobject::NO_WRITE
+    )
+);


-    Info<< "Reading diffusivity DT\n" << endl;
+Info<< "Reading diffusivity DT\n" << endl;

-    dimensionedScalar DT
-    (
-        transportProperties.lookup("DT")
-    );
+dimensionedScalar DT
+(
+    transportProperties.lookup("DT")
+);
+
+#include "createFvOptions.H"
--- a/applications/solvers/basic/laplacianFoam/laplacianFoam.C
+++ b/applications/solvers/basic/laplacianFoam/laplacianFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -45,7 +45,6 @@ int main(int argc, char *argv[])
    simpleControl simple(mesh);

    #include "createFields.H"
-    #include "createFvOptions.H"

    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

--- a/applications/solvers/basic/scalarTransportFoam/createFields.H
+++ b/applications/solvers/basic/scalarTransportFoam/createFields.H
@ -1,55 +1,57 @@
-    Info<< "Reading field T\n" << endl;
+Info<< "Reading field T\n" << endl;

-    volScalarField T
+volScalarField T
+(
+    IOobject
    (
-        IOobject
-        (
-            "T",
-            runTime.timeName(),
-            mesh,
-            IOobject::MUST_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh
-    );
+        "T",
+        runTime.timeName(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
+    ),
+    mesh
+);


-    Info<< "Reading field U\n" << endl;
+Info<< "Reading field U\n" << endl;

-    volVectorField U
+volVectorField U
+(
+    IOobject
    (
-        IOobject
-        (
-            "U",
-            runTime.timeName(),
-            mesh,
-            IOobject::MUST_READ,
-            IOobject::AUTO_WRITE
-        ),
-        mesh
-    );
+        "U",
+        runTime.timeName(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
+    ),
+    mesh
+);


-    Info<< "Reading transportProperties\n" << endl;
+Info<< "Reading transportProperties\n" << endl;

-    IOdictionary transportProperties
+IOdictionary transportProperties
+(
+    IOobject
    (
-        IOobject
-        (
-            "transportProperties",
-            runTime.constant(),
-            mesh,
-            IOobject::MUST_READ_IF_MODIFIED,
-            IOobject::NO_WRITE
-        )
-    );
+        "transportProperties",
+        runTime.constant(),
+        mesh,
+        IOobject::MUST_READ_IF_MODIFIED,
+        IOobject::NO_WRITE
+    )
+);


-    Info<< "Reading diffusivity DT\n" << endl;
+Info<< "Reading diffusivity DT\n" << endl;

-    dimensionedScalar DT
-    (
-        transportProperties.lookup("DT")
-    );
+dimensionedScalar DT
+(
+    transportProperties.lookup("DT")
+);

-    #include "createPhi.H"
+#include "createPhi.H"
+
+#include "createFvOptions.H"
--- a/applications/solvers/basic/scalarTransportFoam/scalarTransportFoam.C
+++ b/applications/solvers/basic/scalarTransportFoam/scalarTransportFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -44,7 +44,6 @@ int main(int argc, char *argv[])
    simpleControl simple(mesh);

    #include "createFields.H"
-    #include "createFvOptions.H"

    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

--- a/applications/solvers/combustion/PDRFoam/PDRFoam.C
+++ b/applications/solvers/combustion/PDRFoam/PDRFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -93,7 +93,6 @@ int main(int argc, char *argv[])
    #include "readGravitationalAcceleration.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"
    #include "createTimeControls.H"
    #include "compressibleCourantNo.H"
--- a/applications/solvers/combustion/PDRFoam/createFields.H
+++ b/applications/solvers/combustion/PDRFoam/createFields.H
@ -228,3 +228,4 @@ fields.add(thermo.heu());
 flameWrinkling->addXi(fields);

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/combustion/XiFoam/XiEngineFoam/XiEngineFoam.C
+++ b/applications/solvers/combustion/XiFoam/XiEngineFoam/XiEngineFoam.C
@ -76,7 +76,6 @@ int main(int argc, char *argv[])
    #include "readCombustionProperties.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "createRhoUf.H"
    #include "initContinuityErrs.H"
    #include "readEngineTimeControls.H"
--- a/applications/solvers/combustion/XiFoam/XiFoam.C
+++ b/applications/solvers/combustion/XiFoam/XiFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -71,7 +71,6 @@ int main(int argc, char *argv[])
    #include "readCombustionProperties.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"
    #include "createTimeControls.H"
    #include "compressibleCourantNo.H"
--- a/applications/solvers/combustion/XiFoam/createFields.H
+++ b/applications/solvers/combustion/XiFoam/createFields.H
@ -140,3 +140,4 @@ fields.add(thermo.he());
 fields.add(thermo.heu());

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/combustion/coldEngineFoam/coldEngineFoam.C
+++ b/applications/solvers/combustion/coldEngineFoam/coldEngineFoam.C
@ -52,7 +52,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "createRhoUf.H"
    #include "initContinuityErrs.H"
    #include "readEngineTimeControls.H"
--- a/applications/solvers/combustion/coldEngineFoam/createFields.H
+++ b/applications/solvers/combustion/coldEngineFoam/createFields.H
@ -68,3 +68,4 @@ Info<< "Creating field kinetic energy K\n" << endl;
 volScalarField K("K", 0.5*magSqr(U));

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/combustion/fireFoam/createFields.H
+++ b/applications/solvers/combustion/fireFoam/createFields.H
@ -166,3 +166,4 @@ volScalarField K("K", 0.5*magSqr(U));
 #include "createSurfaceFilmModel.H"
 #include "createPyrolysisModel.H"
 #include "createRadiationModel.H"
+#include "createFvOptions.H"
--- a/applications/solvers/combustion/fireFoam/fireFoam.C
+++ b/applications/solvers/combustion/fireFoam/fireFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -54,7 +54,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"
    #include "createTimeControls.H"
    #include "compressibleCourantNo.H"
--- a/applications/solvers/combustion/reactingFoam/createFields.H
+++ b/applications/solvers/combustion/reactingFoam/createFields.H
@ -112,3 +112,4 @@ volScalarField Qdot
 );

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/combustion/reactingFoam/reactingFoam.C
+++ b/applications/solvers/combustion/reactingFoam/reactingFoam.C
@ -53,7 +53,6 @@ int main(int argc, char *argv[])
    #include "initContinuityErrs.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"

    turbulence->validate();

--- a/applications/solvers/combustion/reactingFoam/rhoReactingBuoyantFoam/createFields.H
+++ b/applications/solvers/combustion/reactingFoam/rhoReactingBuoyantFoam/createFields.H
@ -132,3 +132,4 @@ volScalarField Qdot
 );

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/combustion/reactingFoam/rhoReactingBuoyantFoam/rhoReactingBuoyantFoam.C
+++ b/applications/solvers/combustion/reactingFoam/rhoReactingBuoyantFoam/rhoReactingBuoyantFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2012-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2012-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -53,7 +53,6 @@ int main(int argc, char *argv[])
    #include "initContinuityErrs.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"

    turbulence->validate();

--- a/applications/solvers/combustion/reactingFoam/rhoReactingFoam/createFields.H
+++ b/applications/solvers/combustion/reactingFoam/rhoReactingFoam/createFields.H
@ -113,3 +113,4 @@ volScalarField Qdot
 );

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/combustion/reactingFoam/rhoReactingFoam/rhoReactingFoam.C
+++ b/applications/solvers/combustion/reactingFoam/rhoReactingFoam/rhoReactingFoam.C
@ -54,7 +54,6 @@ int main(int argc, char *argv[])
    #include "initContinuityErrs.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"

    turbulence->validate();

--- a/applications/solvers/compressible/rhoPimpleFoam/createFields.H
+++ b/applications/solvers/compressible/rhoPimpleFoam/createFields.H
@ -73,3 +73,4 @@ Info<< "Creating field kinetic energy K\n" << endl;
 volScalarField K("K", 0.5*magSqr(U));

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/rhoPimpleDyMFoam.C
+++ b/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleDyMFoam/rhoPimpleDyMFoam.C
@ -58,7 +58,6 @@ int main(int argc, char *argv[])
    #include "initContinuityErrs.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "createRhoUf.H"
    #include "createControls.H"

--- a/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleFoam.C
+++ b/applications/solvers/compressible/rhoPimpleFoam/rhoPimpleFoam.C
@ -57,7 +57,6 @@ int main(int argc, char *argv[])
    #include "initContinuityErrs.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"

    turbulence->validate();

--- a/applications/solvers/compressible/rhoSimpleFoam/createFields.H
+++ b/applications/solvers/compressible/rhoSimpleFoam/createFields.H
@ -57,3 +57,4 @@ autoPtr<compressible::turbulenceModel> turbulence
 dimensionedScalar initialMass = fvc::domainIntegrate(rho);

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/compressible/rhoSimpleFoam/rhoPorousSimpleFoam/rhoPorousSimpleFoam.C
+++ b/applications/solvers/compressible/rhoSimpleFoam/rhoPorousSimpleFoam/rhoPorousSimpleFoam.C
@ -49,7 +49,6 @@ int main(int argc, char *argv[])
    #include "createMesh.H"
    #include "createControl.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "createZones.H"
    #include "initContinuityErrs.H"

--- a/applications/solvers/compressible/rhoSimpleFoam/rhoSimpleFoam.C
+++ b/applications/solvers/compressible/rhoSimpleFoam/rhoSimpleFoam.C
@ -48,7 +48,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    turbulence->validate();
--- a/applications/solvers/compressible/sonicFoam/createFields.H
+++ b/applications/solvers/compressible/sonicFoam/createFields.H
@ -54,3 +54,4 @@ Info<< "Creating field kinetic energy K\n" << endl;
 volScalarField K("K", 0.5*magSqr(U));

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/compressible/sonicFoam/sonicDyMFoam/sonicDyMFoam.C
+++ b/applications/solvers/compressible/sonicFoam/sonicDyMFoam/sonicDyMFoam.C
@ -54,7 +54,6 @@ int main(int argc, char *argv[])
    #include "createControls.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "createRhoUf.H"
    #include "compressibleCourantNo.H"
    #include "setInitialDeltaT.H"
--- a/applications/solvers/compressible/sonicFoam/sonicFoam.C
+++ b/applications/solvers/compressible/sonicFoam/sonicFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -51,7 +51,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    turbulence->validate();
--- a/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/buoyantBoussinesqPimpleFoam.C
+++ b/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/buoyantBoussinesqPimpleFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -63,7 +63,6 @@ int main(int argc, char *argv[])
    #include "createMesh.H"
    #include "createControl.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "createTimeControls.H"
    #include "CourantNo.H"
    #include "setInitialDeltaT.H"
--- a/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/createFields.H
+++ b/applications/solvers/heatTransfer/buoyantBoussinesqPimpleFoam/createFields.H
@ -123,3 +123,4 @@ mesh.setFluxRequired(p_rgh.name());

 #include "createMRF.H"
 #include "createIncompressibleRadiationModel.H"
+#include "createFvOptions.H"
--- a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/buoyantBoussinesqSimpleFoam.C
+++ b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/buoyantBoussinesqSimpleFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -63,7 +63,6 @@ int main(int argc, char *argv[])
    #include "createMesh.H"
    #include "createControl.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    turbulence->validate();
--- a/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/createFields.H
+++ b/applications/solvers/heatTransfer/buoyantBoussinesqSimpleFoam/createFields.H
@ -123,3 +123,4 @@ mesh.setFluxRequired(p_rgh.name());

 #include "createMRF.H"
 #include "createIncompressibleRadiationModel.H"
+#include "createFvOptions.H"
--- a/applications/solvers/heatTransfer/buoyantPimpleFoam/buoyantPimpleFoam.C
+++ b/applications/solvers/heatTransfer/buoyantPimpleFoam/buoyantPimpleFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -52,7 +52,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"
    #include "createTimeControls.H"
    #include "compressibleCourantNo.H"
--- a/applications/solvers/heatTransfer/buoyantPimpleFoam/createFields.H
+++ b/applications/solvers/heatTransfer/buoyantPimpleFoam/createFields.H
@ -106,3 +106,4 @@ volScalarField K("K", 0.5*magSqr(U));

 #include "createMRF.H"
 #include "createRadiationModel.H"
+#include "createFvOptions.H"
--- a/applications/solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C
+++ b/applications/solvers/heatTransfer/buoyantSimpleFoam/buoyantSimpleFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -49,7 +49,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    turbulence->validate();
--- a/applications/solvers/heatTransfer/buoyantSimpleFoam/createFields.H
+++ b/applications/solvers/heatTransfer/buoyantSimpleFoam/createFields.H
@ -88,3 +88,4 @@ dimensionedScalar totalVolume = sum(mesh.V());

 #include "createMRF.H"
 #include "createRadiationModel.H"
+#include "createFvOptions.H"
--- a/applications/solvers/heatTransfer/thermoFoam/createFields.H
+++ b/applications/solvers/heatTransfer/thermoFoam/createFields.H
@ -54,3 +54,4 @@ Info<< "Creating field kinetic energy K\n" << endl;
 volScalarField K("K", 0.5*magSqr(U));

 #include "createRadiationModel.H"
+#include "createFvOptions.H"
--- a/applications/solvers/heatTransfer/thermoFoam/thermoFoam.C
+++ b/applications/solvers/heatTransfer/thermoFoam/thermoFoam.C
@ -49,7 +49,6 @@ int main(int argc, char *argv[])
    #include "createTime.H"
    #include "createMesh.H"
    #include "createFields.H"
-    #include "createFvOptions.H"

    const volScalarField& alphaEff = talphaEff();

--- a/applications/solvers/incompressible/adjointShapeOptimizationFoam/adjointShapeOptimizationFoam.C
+++ b/applications/solvers/incompressible/adjointShapeOptimizationFoam/adjointShapeOptimizationFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -76,7 +76,6 @@ int main(int argc, char *argv[])
    #include "createMesh.H"
    #include "createControl.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"
    #include "initAdjointContinuityErrs.H"

--- a/applications/solvers/incompressible/adjointShapeOptimizationFoam/createFields.H
+++ b/applications/solvers/incompressible/adjointShapeOptimizationFoam/createFields.H
@ -109,3 +109,5 @@ volScalarField alpha
 );
 zeroCells(alpha, inletCells);
 //zeroCells(alpha, outletCells);
+
+#include "createFvOptions.H"
--- a/applications/solvers/incompressible/boundaryFoam/boundaryFoam.C
+++ b/applications/solvers/incompressible/boundaryFoam/boundaryFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -53,7 +53,6 @@ int main(int argc, char *argv[])
    #include "createTime.H"
    #include "createMesh.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "interrogateWallPatches.H"

    turbulence->validate();
--- a/applications/solvers/incompressible/boundaryFoam/createFields.H
+++ b/applications/solvers/incompressible/boundaryFoam/createFields.H
@ -47,3 +47,5 @@ dimensionedVector gradP
    dimensionSet(0, 1, -2, 0, 0),
    Zero
 );
+
+#include "createFvOptions.H"
--- a/applications/solvers/incompressible/pimpleFoam/SRFPimpleFoam/SRFPimpleFoam.C
+++ b/applications/solvers/incompressible/pimpleFoam/SRFPimpleFoam/SRFPimpleFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -51,7 +51,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createTimeControls.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    turbulence->validate();
--- a/applications/solvers/incompressible/pimpleFoam/SRFPimpleFoam/createFields.H
+++ b/applications/solvers/incompressible/pimpleFoam/SRFPimpleFoam/createFields.H
@ -71,3 +71,5 @@ autoPtr<incompressible::turbulenceModel> turbulence
 (
    incompressible::turbulenceModel::New(U, phi, laminarTransport)
 );
+
+#include "createFvOptions.H"
--- a/applications/solvers/incompressible/pimpleFoam/createFields.H
+++ b/applications/solvers/incompressible/pimpleFoam/createFields.H
@ -43,3 +43,4 @@ autoPtr<incompressible::turbulenceModel> turbulence
 );

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pimpleDyMFoam.C
+++ b/applications/solvers/incompressible/pimpleFoam/pimpleDyMFoam/pimpleDyMFoam.C
@ -53,7 +53,6 @@ int main(int argc, char *argv[])
    #include "createControls.H"
    #include "createFields.H"
    #include "createUf.H"
-    #include "createFvOptions.H"
    #include "CourantNo.H"
    #include "setInitialDeltaT.H"

--- a/applications/solvers/incompressible/pimpleFoam/pimpleFoam.C
+++ b/applications/solvers/incompressible/pimpleFoam/pimpleFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -52,7 +52,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createTimeControls.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    turbulence->validate();
--- a/applications/solvers/incompressible/pisoFoam/createFields.H
+++ b/applications/solvers/incompressible/pisoFoam/createFields.H
@ -43,3 +43,4 @@ autoPtr<incompressible::turbulenceModel> turbulence
 );

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/incompressible/pisoFoam/pisoFoam.C
+++ b/applications/solvers/incompressible/pisoFoam/pisoFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -51,7 +51,6 @@ int main(int argc, char *argv[])
    #include "createMesh.H"
    #include "createControl.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    turbulence->validate();
--- a/applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/SRFSimpleFoam.C
+++ b/applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/SRFSimpleFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -48,7 +48,6 @@ int main(int argc, char *argv[])
    #include "createMesh.H"
    #include "createControl.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    turbulence->validate();
--- a/applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/createFields.H
+++ b/applications/solvers/incompressible/simpleFoam/SRFSimpleFoam/createFields.H
@ -68,3 +68,5 @@ autoPtr<incompressible::turbulenceModel> turbulence
 (
    incompressible::turbulenceModel::New(U, phi, laminarTransport)
 );
+
+#include "createFvOptions.H"
--- a/applications/solvers/lagrangian/coalChemistryFoam/coalChemistryFoam.C
+++ b/applications/solvers/lagrangian/coalChemistryFoam/coalChemistryFoam.C
@ -56,7 +56,6 @@ int main(int argc, char *argv[])
    #include "createTimeControls.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    turbulence->validate();
--- a/applications/solvers/lagrangian/coalChemistryFoam/createFields.H
+++ b/applications/solvers/lagrangian/coalChemistryFoam/createFields.H
@ -148,3 +148,4 @@ volScalarField Qdot
 #include "createMRF.H"
 #include "createClouds.H"
 #include "createRadiationModel.H"
+#include "createFvOptions.H"
--- a/applications/solvers/lagrangian/reactingParcelFoam/createFields.H
+++ b/applications/solvers/lagrangian/reactingParcelFoam/createFields.H
@ -139,3 +139,4 @@ volScalarField Qdot
 #include "createRadiationModel.H"
 #include "createClouds.H"
 #include "createSurfaceFilmModel.H"
+#include "createFvOptions.H"
--- a/applications/solvers/lagrangian/reactingParcelFoam/reactingParcelFoam.C
+++ b/applications/solvers/lagrangian/reactingParcelFoam/reactingParcelFoam.C
@ -56,7 +56,6 @@ int main(int argc, char *argv[])
    #include "createTimeControls.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    turbulence->validate();
--- a/applications/solvers/lagrangian/reactingParcelFoam/simpleReactingParcelFoam/createFields.H
+++ b/applications/solvers/lagrangian/reactingParcelFoam/simpleReactingParcelFoam/createFields.H
@ -120,3 +120,4 @@ volScalarField Qdot
 #include "createMRF.H"
 #include "createRadiationModel.H"
 #include "createClouds.H"
+#include "createFvOptions.H"
--- a/applications/solvers/lagrangian/reactingParcelFoam/simpleReactingParcelFoam/simpleReactingParcelFoam.C
+++ b/applications/solvers/lagrangian/reactingParcelFoam/simpleReactingParcelFoam/simpleReactingParcelFoam.C
@ -52,7 +52,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    turbulence->validate();
--- a/applications/solvers/lagrangian/sprayFoam/createFields.H
+++ b/applications/solvers/lagrangian/sprayFoam/createFields.H
@ -134,3 +134,4 @@ volScalarField Qdot
 #include "createMRF.H"
 #include "createClouds.H"
 #include "createRadiationModel.H"
+#include "createFvOptions.H"
--- a/applications/solvers/lagrangian/sprayFoam/engineFoam/engineFoam.C
+++ b/applications/solvers/lagrangian/sprayFoam/engineFoam/engineFoam.C
@ -56,7 +56,6 @@ int main(int argc, char *argv[])
    #include "readEngineTimeControls.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "createRhoUf.H"
    #include "compressibleCourantNo.H"
    #include "setInitialDeltaT.H"
--- a/applications/solvers/lagrangian/sprayFoam/sprayDyMFoam/sprayDyMFoam.C
+++ b/applications/solvers/lagrangian/sprayFoam/sprayDyMFoam/sprayDyMFoam.C
@ -55,7 +55,6 @@ int main(int argc, char *argv[])
    #include "createFields.H"
    #include "createFieldRefs.H"
    #include "createRhoUf.H"
-    #include "createFvOptions.H"
    #include "compressibleCourantNo.H"
    #include "setInitialDeltaT.H"
    #include "initContinuityErrs.H"
--- a/applications/solvers/lagrangian/sprayFoam/sprayFoam.C
+++ b/applications/solvers/lagrangian/sprayFoam/sprayFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -52,7 +52,6 @@ int main(int argc, char *argv[])
    #include "createTimeControls.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "compressibleCourantNo.H"
    #include "setInitialDeltaT.H"
    #include "initContinuityErrs.H"
--- a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
+++ b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
@ -65,7 +65,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createFields.H"
    #include "createAlphaFluxes.H"
-    #include "createFvOptions.H"
    #include "createUf.H"
    #include "createControls.H"
    #include "CourantNo.H"
--- a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFilmFoam/compressibleInterFilmFoam.C
+++ b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFilmFoam/compressibleInterFilmFoam.C
@ -65,7 +65,6 @@ int main(int argc, char *argv[])
    #include "createFields.H"
    #include "createAlphaFluxes.H"
    #include "createSurfaceFilmModel.H"
-    #include "createFvOptions.H"

    volScalarField& p = mixture.p();
    volScalarField& T = mixture.T();
--- a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
+++ b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
@ -62,7 +62,6 @@ int main(int argc, char *argv[])
    #include "createTimeControls.H"
    #include "createFields.H"
    #include "createAlphaFluxes.H"
-    #include "createFvOptions.H"

    volScalarField& p = mixture.p();
    volScalarField& T = mixture.T();
--- a/applications/solvers/multiphase/compressibleInterFoam/createFields.H
+++ b/applications/solvers/multiphase/compressibleInterFoam/createFields.H
@ -99,3 +99,4 @@ Info<< "Creating field kinetic energy K\n" << endl;
 volScalarField K("K", 0.5*magSqr(U));

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/multiphase/driftFluxFoam/createFields.H
+++ b/applications/solvers/multiphase/driftFluxFoam/createFields.H
@ -132,3 +132,4 @@ mesh.setFluxRequired(alpha1.name());
 tmp<surfaceScalarField> talphaPhiCorr0;

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/multiphase/driftFluxFoam/driftFluxFoam.C
+++ b/applications/solvers/multiphase/driftFluxFoam/driftFluxFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -57,7 +57,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createTimeControls.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    volScalarField& alpha2(mixture.alpha2());
--- a/applications/solvers/multiphase/interFoam/createFields.H
+++ b/applications/solvers/multiphase/interFoam/createFields.H
@ -122,3 +122,4 @@ mesh.setFluxRequired(p_rgh.name());
 mesh.setFluxRequired(alpha1.name());

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
+++ b/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
@ -61,7 +61,6 @@ int main(int argc, char *argv[])
    #include "createDyMControls.H"
    #include "createFields.H"
    #include "createAlphaFluxes.H"
-    #include "createFvOptions.H"

    volScalarField rAU
    (
--- a/applications/solvers/multiphase/interFoam/interFoam.C
+++ b/applications/solvers/multiphase/interFoam/interFoam.C
@ -64,7 +64,6 @@ int main(int argc, char *argv[])
    #include "initContinuityErrs.H"
    #include "createFields.H"
    #include "createAlphaFluxes.H"
-    #include "createFvOptions.H"
    #include "correctPhi.H"

    turbulence->validate();
--- a/applications/solvers/multiphase/interFoam/interMixingFoam/createFields.H
+++ b/applications/solvers/multiphase/interFoam/interMixingFoam/createFields.H
@ -125,3 +125,4 @@ mesh.setFluxRequired(p_rgh.name());
 mesh.setFluxRequired(alpha2.name());

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/multiphase/interFoam/interMixingFoam/interMixingFoam.C
+++ b/applications/solvers/multiphase/interFoam/interMixingFoam/interMixingFoam.C
@ -54,7 +54,6 @@ int main(int argc, char *argv[])
    #include "createTimeControls.H"
    #include "initContinuityErrs.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "correctPhi.H"

    turbulence->validate();
--- a/applications/solvers/multiphase/interPhaseChangeFoam/createFields.H
+++ b/applications/solvers/multiphase/interPhaseChangeFoam/createFields.H
@ -107,3 +107,5 @@ if (p_rgh.needReference())

 mesh.setFluxRequired(p_rgh.name());
 mesh.setFluxRequired(alpha1.name());
+
+#include "createFvOptions.H"
--- a/applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeDyMFoam/interPhaseChangeDyMFoam.C
+++ b/applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeDyMFoam/interPhaseChangeDyMFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -66,7 +66,6 @@ int main(int argc, char *argv[])
    #include "../interFoam/interDyMFoam/createDyMControls.H"
    #include "initContinuityErrs.H"
    #include "createFields.H"
-    #include "createFvOptions.H"

    volScalarField rAU
    (
--- a/applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeFoam.C
+++ b/applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -60,7 +60,6 @@ int main(int argc, char *argv[])
    #include "createMesh.H"
    #include "createControl.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "createTimeControls.H"
    #include "CourantNo.H"
    #include "setInitialDeltaT.H"
--- a/applications/solvers/multiphase/multiphaseInterFoam/createFields.H
+++ b/applications/solvers/multiphase/multiphaseInterFoam/createFields.H
@ -93,3 +93,4 @@ if (p_rgh.needReference())
 mesh.setFluxRequired(p_rgh.name());

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterDyMFoam/multiphaseInterDyMFoam.C
+++ b/applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterDyMFoam/multiphaseInterDyMFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -55,7 +55,6 @@ int main(int argc, char *argv[])
    #include "createTimeControls.H"
    #include "createDyMControls.H"
    #include "createFields.H"
-    #include "createFvOptions.H"

    volScalarField rAU
    (
--- a/applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterFoam.C
+++ b/applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -51,7 +51,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "initContinuityErrs.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "createTimeControls.H"
    #include "correctPhi.H"
    #include "CourantNo.H"
--- a/applications/solvers/multiphase/potentialFreeSurfaceFoam/createFields.H
+++ b/applications/solvers/multiphase/potentialFreeSurfaceFoam/createFields.H
@ -77,3 +77,4 @@ setRefCell(p_gh, pimple.dict(), p_ghRefCell, p_ghRefValue);
 mesh.setFluxRequired(p_gh.name());

 #include "createMRF.H"
+#include "createFvOptions.H"
--- a/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceDyMFoam/potentialFreeSurfaceDyMFoam.C
+++ b/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceDyMFoam/potentialFreeSurfaceDyMFoam.C
@ -60,7 +60,6 @@ int main(int argc, char *argv[])
    #include "createTimeControls.H"
    #include "createDyMControls.H"
    #include "createFields.H"
-    #include "createFvOptions.H"

    volScalarField rAU
    (
--- a/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceFoam.C
+++ b/applications/solvers/multiphase/potentialFreeSurfaceFoam/potentialFreeSurfaceFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -52,7 +52,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createTimeControls.H"
    #include "createFields.H"
-    #include "createFvOptions.H"
    #include "initContinuityErrs.H"

    turbulence->validate();
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/createFields.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/createFields.H
@ -82,3 +82,4 @@ 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"
+#include "createFvOptions.H"
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseEulerFoam.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseEulerFoam.C
@ -2,7 +2,7 @@
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011-2016 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
     \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@ -49,7 +49,6 @@ int main(int argc, char *argv[])
    #include "createControl.H"
    #include "createFields.H"
    #include "createFieldRefs.H"
-    #include "createFvOptions.H"
    #include "createTimeControls.H"
    #include "CourantNos.H"
    #include "setInitialDeltaT.H"