solvers: rationalize the construction of the fvOptions

applications/solvers/multiphase/interFoam/LTSInterFoam/LTSInterFoam.C

@@ -59,6 +59,7 @@ int main(int argc, char *argv[])
 
     #include "initContinuityErrs.H"
     #include "createFields.H"
+    #include "createFvOptions.H"
     #include "createPrghCorrTypes.H"
     #include "correctPhi.H"
     #include "CourantNo.H"

applications/solvers/multiphase/interFoam/createFields.H

@@ -128,8 +128,5 @@
     }
 
 
-    fv::IOoptionList fvOptions(mesh);
-
-
     // MULES Correction
     tmp<surfaceScalarField> tphiAlphaCorr0;

applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C

@@ -54,6 +54,7 @@ int main(int argc, char *argv[])
     pimpleControl pimple(mesh);
 
     #include "createFields.H"
+    #include "createFvOptions.H"
     #include "readTimeControls.H"
     #include "createPrghCorrTypes.H"
 

applications/solvers/multiphase/interFoam/interFoam.C

@@ -58,6 +58,7 @@ int main(int argc, char *argv[])
 
     #include "initContinuityErrs.H"
     #include "createFields.H"
+    #include "createFvOptions.H"
     #include "readTimeControls.H"
     #include "createPrghCorrTypes.H"
     #include "correctPhi.H"

applications/solvers/multiphase/interFoam/interMixingFoam/createFields.H

@@ -118,5 +118,3 @@
         );
         p_rgh = p - rho*gh;
     }
-
-    fv::IOoptionList fvOptions(mesh);

applications/solvers/multiphase/interFoam/interMixingFoam/interMixingFoam.C

@@ -52,6 +52,7 @@ int main(int argc, char *argv[])
     #include "readGravitationalAcceleration.H"
     #include "initContinuityErrs.H"
     #include "createFields.H"
+    #include "createFvOptions.H"
     #include "readTimeControls.H"
     #include "CourantNo.H"
     #include "setInitialDeltaT.H"

applications/solvers/multiphase/interFoam/porousInterFoam/porousInterFoam.C

@@ -60,6 +60,7 @@ int main(int argc, char *argv[])
 
     #include "initContinuityErrs.H"
     #include "createFields.H"
+    #include "createFvOptions.H"
     #include "createPorousZones.H"
     #include "readTimeControls.H"
     #include "createPrghCorrTypes.H"

applications/solvers/multiphase/interPhaseChangeFoam/createFields.H

@@ -105,5 +105,3 @@
         );
         p_rgh = p - rho*gh;
     }
-
-    fv::IOoptionList fvOptions(mesh);

applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeDyMFoam/interPhaseChangeDyMFoam.C

@@ -65,6 +65,7 @@ int main(int argc, char *argv[])
     #include "readGravitationalAcceleration.H"
     #include "initContinuityErrs.H"
     #include "createFields.H"
+    #include "createFvOptions.H"
     #include "readTimeControls.H"
     #include "createPcorrTypes.H"
 

applications/solvers/multiphase/interPhaseChangeFoam/interPhaseChangeFoam.C

@@ -63,6 +63,7 @@ int main(int argc, char *argv[])
     #include "readGravitationalAcceleration.H"
     #include "initContinuityErrs.H"
     #include "createFields.H"
+    #include "createFvOptions.H"
     #include "readTimeControls.H"
     #include "createPrghCorrTypes.H"
     #include "CourantNo.H"

applications/solvers/multiphase/multiphaseInterFoam/createFields.H

@@ -91,5 +91,3 @@
             pRefValue - getRefCellValue(p, pRefCell)
         );
     }
-
-    fv::IOoptionList fvOptions(mesh);

applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterDyMFoam/multiphaseInterDyMFoam.C

@@ -52,6 +52,7 @@ int main(int argc, char *argv[])
     pimpleControl pimple(mesh);
 
     #include "createFields.H"
+    #include "createFvOptions.H"
     #include "readTimeControls.H"
     #include "createPrghCorrTypes.H"
 

applications/solvers/multiphase/multiphaseInterFoam/multiphaseInterFoam.C

@@ -51,6 +51,7 @@ int main(int argc, char *argv[])
 
     #include "initContinuityErrs.H"
     #include "createFields.H"
+    #include "createFvOptions.H"
     #include "readTimeControls.H"
     #include "createPrghCorrTypes.H"
     #include "correctPhi.H"