From 6f1102cb5d498459b966607e185b00b8b795e07d Mon Sep 17 00:00:00 2001
From: Martin Beaudoin <beaudoin.martin@ireq.ca>
Date: Wed, 28 May 2014 01:11:17 -0400
Subject: [PATCH] Add options for modifying the global controlDict
 values/switches from the command-line

The following command-line options are available:
 -DebugSwitches        key=value
 -DimensionedConstants key=value
 -InfoSwitches         key=value
 -OptimisationSwitches key=value
 -Tolerances           key=value

In order to modify more than one key/value pair for a given category, you will
need to specify a given option more than once.

If an invalid/unknown key is specified, a complete list of valid keys for the
application will be written at the console.

Martin Beaudoin 05/2014
---
 etc/controlDict                               |   2 +-
 src/Pstream/gamma/Pstream.C                   |   2 +-
 src/Pstream/mpi/Pstream.C                     |   2 +-
 .../fvPatchFields/fvPatchVectorNFields.C      |   2 +-
 .../fvsPatchFields/fvsPatchVectorNFields.C    |   2 +-
 .../autoHexMeshDriver/autoRefineDriver.C      |   6 +-
 .../meshRefinement/meshRefinementBaffles.C    |  14 +-
 .../meshRefinement/meshRefinementRefine.C     |  16 +-
 .../dynamicTopoFvMeshCoupled.C                |  43 ++-
 .../slidingInterface/coupleSlidingInterface.C |  17 +-
 .../slidingInterface/slidingInterface.C       |  10 +-
 .../slidingInterface/slidingInterface.H       |  17 +-
 .../slidingInterfaceProjectPoints.C           |  55 ++-
 .../polyTopoChange/polyTopoChange.C           |   2 +-
 .../mesquiteMotionSolver.C                    |   9 +-
 src/finiteArea/faMesh/faMesh.C                |  10 +-
 src/finiteArea/faMesh/faMesh.H                |   2 +-
 .../constraint/cyclic/cyclicFaPatch.C         |  17 +-
 .../constraint/cyclic/cyclicFaPatch.H         |   3 +-
 .../faPatchField/faPatchFields.C              |   2 +-
 .../faePatchField/faePatchFields.C            |   2 +-
 .../finiteArea/faSchemes/faSchemes.C          |   2 +-
 .../cfdTools/general/adjustPhi/adjustPhi.C    |  10 +-
 .../fvPatchField/fvPatchFields.C              |   2 +-
 .../fvsPatchField/fvsPatchFields.C            |   2 +-
 .../finiteVolume/fvSchemes/fvSchemes.C        |   2 +-
 .../constraint/cyclic/cyclicFvPatch.C         |   6 +-
 .../cellPointWeight/cellPointWeight.C         |   2 +-
 src/foam/algorithms/MeshWave/FaceCellWave.C   |  20 +-
 src/foam/algorithms/MeshWave/FaceCellWave.H   |   8 +-
 .../octree/octree/octreeDataBoundBox.C        |   7 +-
 .../octree/octree/octreeDataBoundBox.H        |   3 +-
 src/foam/db/IOstreams/IOstreams.C             |   7 +-
 src/foam/db/IOstreams/IOstreams/IOstream.H    |  11 +-
 src/foam/db/IOstreams/Pstreams/Pstream.C      |  13 +-
 src/foam/db/IOstreams/Pstreams/Pstream.H      |   5 +-
 .../Pstreams/PstreamCombineReduceOps.H        |   2 +-
 .../db/IOstreams/Pstreams/PstreamReduceOps.H  |   4 +-
 .../IOstreams/Pstreams/combineGatherScatter.C |  12 +-
 .../db/IOstreams/Pstreams/gatherScatter.C     |   4 +-
 .../db/IOstreams/Pstreams/gatherScatterList.C |   4 +-
 src/foam/db/IOstreams/Sstreams/OSstreamI.H    |   2 +-
 src/foam/db/Time/TimeIO.C                     |   3 +
 .../includeEntry/includeEntry.C               |   7 +-
 .../includeIfPresentEntry.C                   |   7 +-
 .../inputModeEntry/inputModeEntry.C           |   7 +-
 .../functionEntries/removeEntry/removeEntry.C |   7 +-
 src/foam/db/error/messageStream.C             |   2 +-
 .../functionObject/functionObject.C           |   2 +-
 src/foam/db/regIOobject/regIOobject.C         |   7 +-
 src/foam/db/regIOobject/regIOobject.H         |   3 +-
 src/foam/db/regIOobject/regIOobjectRead.C     |   4 +-
 src/foam/db/scalarRange/scalarRange.C         |   2 +-
 src/foam/db/typeInfo/className.H              |  15 +-
 src/foam/global/JobInfo/JobInfo.C             |   8 +-
 src/foam/global/JobInfo/JobInfo.H             |   3 +-
 src/foam/global/argList/argList.C             |  72 +++-
 src/foam/global/argList/argList.H             |  24 ++
 .../global/controlSwitches/controlSwitches.C  | 145 +++++++
 .../global/controlSwitches/controlSwitches.H  | 358 ++++++++++++++++++
 .../global/controlSwitches/controlSwitchesI.H |  79 ++++
 src/foam/global/controlSwitches/safe_bool.H   | 100 +++++
 src/foam/global/controlSwitches/safe_boolI.H  |  79 ++++
 src/foam/global/debug/debug.C                 | 326 +++++++++++++++-
 src/foam/global/debug/debug.H                 |  32 +-
 .../GGIInterpolation/GGIInterpolation.H       |  15 +-
 .../GGIInterpolationPolygonIntersection.C     |   8 +-
 .../GGIInterpolationQuickRejectTests.C        |  34 +-
 .../GGIInterpolationWeights.C                 |  17 +-
 .../CalcPatchToPatchWeights.C                 |   6 +-
 .../PatchToPatchInterpolation.C               |  17 +-
 .../PatchToPatchInterpolation.H               |   5 +-
 .../BlockLduMatrix/BlockLduMatrix.C           |   9 +-
 .../BlockLduMatrix/BlockLduMatrix.H           |   3 +-
 src/foam/matrices/solution/solution.C         |   2 +-
 src/foam/meshes/meshShapes/cell/cell.C        |   2 +-
 .../polyMesh/globalMeshData/globalMeshData.C  |  10 +-
 .../polyMesh/globalMeshData/globalMeshData.H  |   3 +-
 .../basic/coupled/coupledPolyPatch.H          |   2 +-
 .../constraint/cyclic/cyclicPolyPatch.C       |  10 +-
 .../constraint/processor/processorPolyPatch.C |   4 +-
 .../polyPatches/polyPatch/polyPatch.C         |  11 +-
 .../polyPatches/polyPatch/polyPatch.H         |   4 +-
 .../PrimitivePatch/PrimitivePatch.C           |   9 +-
 .../PrimitivePatch/PrimitivePatch.H           |   3 +-
 .../PrimitivePatchProjectPoints.C             |   2 +-
 src/foam/meshes/primitiveMesh/primitiveMesh.H |  13 +-
 .../primitiveMeshCheck/primitiveMeshCheck.C   |  75 ++--
 .../primitiveShapes/triangle/intersection.C   |  16 +-
 .../primitiveShapes/triangle/intersection.H   |  17 +-
 .../primitiveShapes/triangle/triangleI.H      |  20 +-
 .../primitives/strings/fileName/fileName.C    |   2 +-
 src/foam/primitives/strings/string/string.C   |   2 +-
 src/foam/primitives/strings/word/word.C       |   2 +-
 src/lagrangian/basic/Particle/Particle.H      |   3 +-
 .../COxidationKineticDiffusionLimitedRate.C   |   4 +-
 .../COxidationMurphyShaddix.C                 |   6 +-
 .../dieselSpray/parcel/setRelaxationTimes.C   |   2 +-
 .../atomizationModel/LISA/LISA.C              |   2 +-
 .../blobsSheetAtomization.C                   |   2 +-
 .../breakupModel/reitzDiwakar/reitzDiwakar.C  |   2 +-
 .../parcels/Templates/DsmcParcel/DsmcParcel.H |   5 +-
 .../KinematicParcel/KinematicParcel.H         |   5 +-
 .../ReactingMultiphaseParcel.C                |   2 +-
 .../Templates/ReactingParcel/ReactingParcel.C |   8 +-
 .../LiquidEvaporation/LiquidEvaporation.C     |   4 +-
 .../SingleKineticRateDevolatilisation.C       |   2 +-
 src/lduSolvers/amg/amgPolicy/pamgPolicy.C     |   8 +-
 src/lduSolvers/amg/amgPolicy/pamgPolicy.H     |   3 +-
 .../surfaceIntersection/surfaceIntersection.C |  12 +-
 .../calculateMeshToMeshAddressing.C           |   8 +-
 .../meshToMesh/calculateMeshToMeshWeights.C   |   2 +-
 .../meshToMesh/meshToMesh.C                   |   6 +-
 .../meshToMesh/meshToMesh.H                   |   5 +-
 .../liquidMixture/liquidMixture.C             |   2 +-
 .../radiationConstants/radiationConstants.C   |  10 +-
 .../radiationConstants/radiationConstants.H   |   1 +
 .../specie/specie/specie.C                    |  23 +-
 .../specie/specie/specie.H                    |   8 +-
 .../specie/specie/specieI.H                   |   2 +-
 .../specie/thermo/eConst/eConstThermoI.H      |   2 +-
 .../thermo/hPolynomial/hPolynomialThermo.C    |   4 +-
 .../specie/thermo/janaf/janafThermoI.H        |  10 +-
 .../specie/thermo/specieThermo/specieThermo.C |  19 +-
 .../specie/thermo/specieThermo/specieThermo.H |   7 +-
 .../thermo/specieThermo/specieThermoI.H       |  22 +-
 .../specie/transport/const/constTransportI.H  |   4 +-
 .../polynomial/polynomialTransportI.H         |   4 +-
 .../sutherland/sutherlandTransportI.H         |   4 +-
 129 files changed, 1779 insertions(+), 412 deletions(-)
 create mode 100644 src/foam/global/controlSwitches/controlSwitches.C
 create mode 100644 src/foam/global/controlSwitches/controlSwitches.H
 create mode 100644 src/foam/global/controlSwitches/controlSwitchesI.H
 create mode 100644 src/foam/global/controlSwitches/safe_bool.H
 create mode 100644 src/foam/global/controlSwitches/safe_boolI.H

diff --git a/etc/controlDict b/etc/controlDict
index b3ed77d82..bbcbea948 100644
--- a/etc/controlDict
+++ b/etc/controlDict
@@ -943,7 +943,7 @@ Tolerances
     slidingEdgeEndCutoffTol         0.0001;
     slidingEdgeCoPlanarTol          0.8;
 
-//     GGIAreaErrorTol                 1e-8;
+     GGIAreaErrorTol                 1e-8;
 }
 
 DimensionedConstants
diff --git a/src/Pstream/gamma/Pstream.C b/src/Pstream/gamma/Pstream.C
index f5f849421..9cee828e1 100644
--- a/src/Pstream/gamma/Pstream.C
+++ b/src/Pstream/gamma/Pstream.C
@@ -320,7 +320,7 @@ bool Pstream::init(int& argc, char**& argv)
 
     Info<< "GAMMA Pstream initialized with:" << nl
         << "    floatTransfer         : " << floatTransfer << nl
-        << "    nProcsSimpleSum       : " << nProcsSimpleSum << nl
+        << "    nProcsSimpleSum       : " << nProcsSimpleSum() << nl
         << "    scheduledTransfer     : " << Pstream::scheduledTransfer << nl
         << Foam::endl;
 
diff --git a/src/Pstream/mpi/Pstream.C b/src/Pstream/mpi/Pstream.C
index b7cddd6ec..1bb8ab404 100644
--- a/src/Pstream/mpi/Pstream.C
+++ b/src/Pstream/mpi/Pstream.C
@@ -154,7 +154,7 @@ void Foam::reduce(scalar& Value, const sumOp<scalar>& bop)
         return;
     }
 
-    if (Pstream::nProcs() <= Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() <= Pstream::nProcsSimpleSum())
     {
         if (Pstream::master())
         {
diff --git a/src/VectorN/finiteVolume/fields/fvPatchFields/fvPatchVectorNFields.C b/src/VectorN/finiteVolume/fields/fvPatchFields/fvPatchVectorNFields.C
index 0bfb01b2c..e50138760 100644
--- a/src/VectorN/finiteVolume/fields/fvPatchFields/fvPatchVectorNFields.C
+++ b/src/VectorN/finiteVolume/fields/fvPatchFields/fvPatchVectorNFields.C
@@ -38,7 +38,7 @@ defineNamedTemplateTypeNameAndDebug(fvPatchTypeField, 0);                     \
 template<>                                                                    \
 int fvPatchTypeField::disallowGenericFvPatchField                             \
 (                                                                             \
-    debug::debugSwitch("disallowGenericFvPatchField", 0)                      \
+    debug::debugSwitchFromDict("disallowGenericFvPatchField", 0)                      \
 );                                                                            \
 defineTemplateRunTimeSelectionTable(fvPatchTypeField, patch);                 \
 defineTemplateRunTimeSelectionTable(fvPatchTypeField, patchMapper);           \
diff --git a/src/VectorN/finiteVolume/fields/fvsPatchFields/fvsPatchVectorNFields.C b/src/VectorN/finiteVolume/fields/fvsPatchFields/fvsPatchVectorNFields.C
index 1c8c36503..108a80ca6 100644
--- a/src/VectorN/finiteVolume/fields/fvsPatchFields/fvsPatchVectorNFields.C
+++ b/src/VectorN/finiteVolume/fields/fvsPatchFields/fvsPatchVectorNFields.C
@@ -38,7 +38,7 @@ defineNamedTemplateTypeNameAndDebug(fvsPatchTypeField, 0);                    \
 template<>                                                                    \
 int fvsPatchTypeField::disallowDefaultFvsPatchField                           \
 (                                                                             \
-    debug::debugSwitch("disallowDefaultFvsPatchField", 0)                     \
+    debug::debugSwitchFromDict("disallowDefaultFvsPatchField", 0)                     \
 );                                                                            \
 defineTemplateRunTimeSelectionTable(fvsPatchTypeField, patch);                \
 defineTemplateRunTimeSelectionTable(fvsPatchTypeField, patchMapper);          \
diff --git a/src/autoMesh/autoHexMesh/autoHexMeshDriver/autoRefineDriver.C b/src/autoMesh/autoHexMesh/autoHexMeshDriver/autoRefineDriver.C
index 4e0e6b5f2..139cb4710 100644
--- a/src/autoMesh/autoHexMesh/autoHexMeshDriver/autoRefineDriver.C
+++ b/src/autoMesh/autoHexMesh/autoHexMeshDriver/autoRefineDriver.C
@@ -393,7 +393,7 @@ void Foam::autoRefineDriver::removeInsideCells
     {
         Pout<< "Writing subsetted mesh to time "
             << meshRefiner_.timeName() << '.' << endl;
-        meshRefiner_.write(debug, mesh.time().path()/meshRefiner_.timeName());
+        meshRefiner_.write(debug(), mesh.time().path()/meshRefiner_.timeName());
         Pout<< "Dumped mesh in = "
             << mesh.time().cpuTimeIncrement() << " s\n" << nl << endl;
     }
@@ -645,7 +645,7 @@ void Foam::autoRefineDriver::zonify
                 << meshRefiner_.timeName() << '.' << endl;
             meshRefiner_.write
             (
-                debug,
+                debug(),
                 mesh.time().path()/meshRefiner_.timeName()
             );
         }
@@ -735,7 +735,7 @@ void Foam::autoRefineDriver::splitAndMergeBaffles
     {
         Pout<< "Writing handleProblemCells mesh to time "
             << meshRefiner_.timeName() << '.' << endl;
-        meshRefiner_.write(debug, mesh.time().path()/meshRefiner_.timeName());
+        meshRefiner_.write(debug(), mesh.time().path()/meshRefiner_.timeName());
     }
 }
 
diff --git a/src/autoMesh/autoHexMesh/meshRefinement/meshRefinementBaffles.C b/src/autoMesh/autoHexMesh/meshRefinement/meshRefinementBaffles.C
index 65cc227d6..3a8f31ba7 100644
--- a/src/autoMesh/autoHexMesh/meshRefinement/meshRefinementBaffles.C
+++ b/src/autoMesh/autoHexMesh/meshRefinement/meshRefinementBaffles.C
@@ -229,7 +229,7 @@ void Foam::meshRefinement::getBafflePatches
 {
     autoPtr<OFstream> str;
     label vertI = 0;
-    if (debug&OBJINTERSECTIONS)
+    if (debug & OBJINTERSECTIONS)
     {
         str.reset
         (
@@ -1499,13 +1499,13 @@ void Foam::meshRefinement::baffleAndSplitMesh
     Info<< "Created baffles in = "
         << runTime.cpuTimeIncrement() << " s\n" << nl << endl;
 
-    printMeshInfo(debug, "After introducing baffles");
+    printMeshInfo(debug(), "After introducing baffles");
 
     if (debug)
     {
         Pout<< "Writing baffled mesh to time " << timeName()
             << endl;
-        write(debug, runTime.path()/"baffles");
+        write(debug(), runTime.path()/"baffles");
         Pout<< "Dumped debug data in = "
             << runTime.cpuTimeIncrement() << " s\n" << nl << endl;
     }
@@ -1581,13 +1581,13 @@ void Foam::meshRefinement::baffleAndSplitMesh
         Info<< "Created baffles in = "
             << runTime.cpuTimeIncrement() << " s\n" << nl << endl;
 
-        printMeshInfo(debug, "After introducing baffles");
+        printMeshInfo(debug(), "After introducing baffles");
 
         if (debug)
         {
             Pout<< "Writing extra baffled mesh to time "
                 << timeName() << endl;
-            write(debug, runTime.path()/"extraBaffles");
+            write(debug(), runTime.path()/"extraBaffles");
             Pout<< "Dumped debug data in = "
                 << runTime.cpuTimeIncrement() << " s\n" << nl << endl;
         }
@@ -1617,13 +1617,13 @@ void Foam::meshRefinement::baffleAndSplitMesh
     Info<< "Split mesh in = "
         << runTime.cpuTimeIncrement() << " s\n" << nl << endl;
 
-    printMeshInfo(debug, "After subsetting");
+    printMeshInfo(debug(), "After subsetting");
 
     if (debug)
     {
         Pout<< "Writing subsetted mesh to time " << timeName()
             << endl;
-        write(debug, runTime.path()/timeName());
+        write(debug(), runTime.path()/timeName());
         Pout<< "Dumped debug data in = "
             << runTime.cpuTimeIncrement() << " s\n" << nl << endl;
     }
diff --git a/src/autoMesh/autoHexMesh/meshRefinement/meshRefinementRefine.C b/src/autoMesh/autoHexMesh/meshRefinement/meshRefinementRefine.C
index 858020005..85d0996ed 100644
--- a/src/autoMesh/autoHexMesh/meshRefinement/meshRefinementRefine.C
+++ b/src/autoMesh/autoHexMesh/meshRefinement/meshRefinementRefine.C
@@ -1269,7 +1269,7 @@ Foam::meshRefinement::refineAndBalance
             << " mesh to time " << timeName() << endl;
         write
         (
-            debug,
+            debug(),
             mesh_.time().path()
            /timeName()
         );
@@ -1282,7 +1282,7 @@ Foam::meshRefinement::refineAndBalance
 
     Info<< "Refined mesh in = "
         << mesh_.time().cpuTimeIncrement() << " s" << endl;
-    printMeshInfo(debug, "After refinement " + msg);
+    printMeshInfo(debug(), "After refinement " + msg);
 
 
     // Load balancing
@@ -1324,7 +1324,7 @@ Foam::meshRefinement::refineAndBalance
             Info<< "Balanced mesh in = "
                 << mesh_.time().cpuTimeIncrement() << " s" << endl;
 
-            printMeshInfo(debug, "After balancing " + msg);
+            printMeshInfo(debug(), "After balancing " + msg);
 
 
             if (debug)
@@ -1333,7 +1333,7 @@ Foam::meshRefinement::refineAndBalance
                     << " mesh to time " << timeName() << endl;
                 write
                 (
-                    debug,
+                    debug(),
                     mesh_.time().path()/timeName()
                 );
                 Pout<< "Dumped debug data in = "
@@ -1447,7 +1447,7 @@ Foam::meshRefinement::balanceAndRefine
         //    Info<< endl;
         //}
 
-        printMeshInfo(debug, "After balancing " + msg);
+        printMeshInfo(debug(), "After balancing " + msg);
 
         if (debug)
         {
@@ -1455,7 +1455,7 @@ Foam::meshRefinement::balanceAndRefine
                 << " mesh to time " << timeName() << endl;
             write
             (
-                debug,
+                debug(),
                 mesh_.time().path()/timeName()
             );
             Pout<< "Dumped debug data in = "
@@ -1478,7 +1478,7 @@ Foam::meshRefinement::balanceAndRefine
             << " mesh to time " << timeName() << endl;
         write
         (
-            debug,
+            debug(),
             mesh_.time().path()
            /timeName()
         );
@@ -1504,7 +1504,7 @@ Foam::meshRefinement::balanceAndRefine
     //    Info<< endl;
     //}
 
-    printMeshInfo(debug, "After refinement " + msg);
+    printMeshInfo(debug(), "After refinement " + msg);
 
     return distMap;
 }
diff --git a/src/dynamicMesh/dynamicFvMesh/dynamicTopoFvMesh/dynamicTopoFvMeshCoupled.C b/src/dynamicMesh/dynamicFvMesh/dynamicTopoFvMesh/dynamicTopoFvMeshCoupled.C
index d5be28ff9..097b5df7e 100644
--- a/src/dynamicMesh/dynamicFvMesh/dynamicTopoFvMesh/dynamicTopoFvMeshCoupled.C
+++ b/src/dynamicMesh/dynamicFvMesh/dynamicTopoFvMesh/dynamicTopoFvMeshCoupled.C
@@ -56,6 +56,7 @@ Author
 #include "processorPolyPatch.H"
 #include "decompositionMethod.H"
 #include "mapDistributePolyMesh.H"
+#include "controlSwitches.H"
 
 namespace Foam
 {
@@ -4894,7 +4895,11 @@ bool dynamicTopoFvMesh::checkCoupledBoundaries(bool report) const
     bool sizeError = false, misMatchError = false;
 
     // Check if a geometric tolerance has been specified.
-    scalar gTol = debug::tolerances("processorMatchTol", 1e-4);
+    const Foam::debug::tolerancesSwitch gTol
+    (
+	"processorMatchTol",
+	1e-4
+    );
 
     // Maintain a list of master / neighbour anchors
     List<vectorField> mAnchors(boundary.size());
@@ -5053,7 +5058,7 @@ bool dynamicTopoFvMesh::checkCoupledBoundaries(bool report) const
                     scalar rMagSf = mag(half1Areas[faceI]);
                     scalar avSf = 0.5 * (fMagSf + rMagSf);
 
-                    if (mag(fMagSf - rMagSf)/avSf > gTol)
+                    if (mag(fMagSf - rMagSf)/avSf > gTol())
                     {
                         misMatchError = true;
 
@@ -5071,7 +5076,7 @@ bool dynamicTopoFvMesh::checkCoupledBoundaries(bool report) const
                         Foam::max
                         (
                             pTol,
-                            gTol * mag(lP[lF[faceI][0]] - lC[faceI])
+                            gTol() * mag(lP[lF[faceI][0]] - lC[faceI])
                         )
                     );
                 }
@@ -5217,7 +5222,7 @@ bool dynamicTopoFvMesh::checkCoupledBoundaries(bool report) const
             scalar nbrMagSf = mag(fAreas[pI][faceI]);
             scalar avSf = 0.5 * (magSf + nbrMagSf);
 
-            if (mag(magSf - nbrMagSf)/avSf > gTol)
+            if (mag(magSf - nbrMagSf)/avSf > gTol())
             {
                 misMatchError = true;
 
@@ -5236,7 +5241,7 @@ bool dynamicTopoFvMesh::checkCoupledBoundaries(bool report) const
                 Foam::max
                 (
                     pTol,
-                    gTol *
+                    gTol() *
                     mag
                     (
                         myPoints[myFaces[faceI][0]]
@@ -6116,10 +6121,14 @@ void dynamicTopoFvMesh::buildLocalCoupledMaps()
 
     // Check if a geometric tolerance has been specified.
     const boundBox& box = polyMesh::bounds();
-    scalar relTol = debug::tolerances("patchFaceMatchTol", 1e-4);
+    const Foam::debug::tolerancesSwitch relTol
+    (
+	"patchFaceMatchTol",
+	1e-4
+    );
 
     // Compute tolerance
-    scalar tol = relTol * box.mag();
+    scalar tol = relTol() * box.mag();
 
     const polyBoundaryMesh& boundary = boundaryMesh();
 
@@ -6181,7 +6190,7 @@ void dynamicTopoFvMesh::buildLocalCoupledMaps()
                 FatalErrorIn("void dynamicTopoFvMesh::buildLocalCoupledMaps()")
                     << " Failed to match all points"
                     << " within a tolerance of: " << tol << nl
-                    << " relTol: " << relTol << nl
+                    << " relTol: " << relTol() << nl
                     << abort(FatalError);
             }
 
@@ -6843,10 +6852,14 @@ void dynamicTopoFvMesh::buildProcessorCoupledMaps()
             const Map<label>& pMap = cMap.entityMap(coupleMap::POINT);
 
             // Fetch relative tolerance
-            scalar relTol = debug::tolerances("processorMatchTol", 1e-4);
+	    const Foam::debug::tolerancesSwitch relTol
+	    (
+		"processorMatchTol",
+		1e-4
+	    );
 
             // Compute tolerance
-            scalar tol = relTol * box.mag();
+            scalar tol = relTol() * box.mag();
 
             forAllConstIter(Map<label>, pMap, pIter)
             {
@@ -8028,7 +8041,11 @@ bool dynamicTopoFvMesh::syncCoupledBoundaryOrdering
     bool anyChange = false, failedPatchMatch = false;
 
     // Fetch tolerance
-    scalar matchTol = Foam::debug::tolerances("patchFaceMatchTol", 1e-4);
+    const Foam::debug::tolerancesSwitch matchTol
+    (
+	"patchFaceMatchTol",
+	1e-4
+    );
 
     // Calculate centres and tolerances for any slave patches
     List<scalarField> slaveTols(nPatches_);
@@ -8090,7 +8107,7 @@ bool dynamicTopoFvMesh::syncCoupledBoundaryOrdering
                     maxLen = max(maxLen, mag(points_[sFace[fpI]] - sfc));
                 }
 
-                slaveTols[slavePatch][fI] = matchTol*maxLen;
+                slaveTols[slavePatch][fI] = matchTol()*maxLen;
             }
 
             // For cyclics, additionally test for halves,
@@ -8442,7 +8459,7 @@ bool dynamicTopoFvMesh::syncCoupledBoundaryOrdering
                     maxLen = max(maxLen, mag(points_[checkFace[fpI]] - fc));
                 }
 
-                slaveTols[pI][fI] = matchTol*maxLen;
+                slaveTols[pI][fI] = matchTol()*maxLen;
             }
 
             // Write out my centres to disk
diff --git a/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/coupleSlidingInterface.C b/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/coupleSlidingInterface.C
index bf9e069cf..6afda9c70 100644
--- a/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/coupleSlidingInterface.C
+++ b/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/coupleSlidingInterface.C
@@ -47,12 +47,13 @@ Author
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
-const Foam::scalar Foam::slidingInterface::edgeCoPlanarTol_
+const Foam::debug::tolerancesSwitch
+Foam::slidingInterface::edgeCoPlanarTol_
 (
-    debug::tolerances("slidingEdgeCoPlanarTol", 0.8)
+    "slidingEdgeCoPlanarTol",
+    0.8
 );
 
-
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
 // Index of debug signs:
@@ -636,8 +637,8 @@ void Foam::slidingInterface::coupleInterface(polyTopoChange& ref) const
 
                     if
                     (
-                        cutOnMaster > edgeEndCutoffTol_
-                     && cutOnMaster < 1.0 - edgeEndCutoffTol_
+                        cutOnMaster > edgeEndCutoffTol_()
+		     && cutOnMaster < 1.0 - edgeEndCutoffTol_()
                     )
                     {
                         // Master is cut, check the slave
@@ -663,7 +664,7 @@ void Foam::slidingInterface::coupleInterface(polyTopoChange& ref) const
                             // Calculate merge tolerance from the
                             // target edge length
                             scalar mergeTol =
-                                edgeCoPlanarTol_*mag(b - a);
+                                edgeCoPlanarTol_()*mag(b - a);
 //                             Pout<< "cutOnMaster: " << cutOnMaster
 //                                 << " masterCutPoint: " << masterCutPoint
 //                                 << " slaveCutPoint: " << slaveCut.hitPoint()
@@ -676,8 +677,8 @@ void Foam::slidingInterface::coupleInterface(polyTopoChange& ref) const
 //                                 << endl;
                             if
                             (
-                                cutOnSlave > edgeEndCutoffTol_
-                             && cutOnSlave < 1.0 - edgeEndCutoffTol_
+                                cutOnSlave > edgeEndCutoffTol_()
+			     && cutOnSlave < 1.0 - edgeEndCutoffTol_()
                              && slaveCut.distance() < mergeTol
                             )
                             {
diff --git a/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/slidingInterface.C b/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/slidingInterface.C
index cba1ae963..d568e26d2 100644
--- a/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/slidingInterface.C
+++ b/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/slidingInterface.C
@@ -631,8 +631,8 @@ void Foam::slidingInterface::modifyMotionPoints(pointField& motionPoints) const
 
                     if
                     (
-                        cutOnMaster > edgeEndCutoffTol_
-                     && cutOnMaster < 1.0 - edgeEndCutoffTol_
+                        cutOnMaster > edgeEndCutoffTol_()
+		     && cutOnMaster < 1.0 - edgeEndCutoffTol_()
                     )
                     {
                         // Master is cut, check the slave
@@ -658,12 +658,12 @@ void Foam::slidingInterface::modifyMotionPoints(pointField& motionPoints) const
                             // Calculate merge tolerance from the
                             // target edge length
                             scalar mergeTol =
-                                edgeCoPlanarTol_*mag(b - a);
+                                edgeCoPlanarTol_()*mag(b - a);
 
                             if
                             (
-                                cutOnSlave > edgeEndCutoffTol_
-                             && cutOnSlave < 1.0 - edgeEndCutoffTol_
+                                cutOnSlave > edgeEndCutoffTol_()
+			     && cutOnSlave < 1.0 - edgeEndCutoffTol_()
                              && slaveCut.distance() < mergeTol
                             )
                             {
diff --git a/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/slidingInterface.H b/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/slidingInterface.H
index e53f2531d..aeecc54bd 100644
--- a/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/slidingInterface.H
+++ b/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/slidingInterface.H
@@ -63,6 +63,7 @@ SourceFiles
 #include "ZoneIDs.H"
 #include "intersection.H"
 #include "Pair.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -261,28 +262,28 @@ private:
     // Static data members
 
         //- Point merge tolerance
-        static const scalar pointMergeTol_;
+        static const Foam::debug::tolerancesSwitch pointMergeTol_;
 
         //- Edge merge tolerance
-        static const scalar edgeMergeTol_;
+        static const Foam::debug::tolerancesSwitch edgeMergeTol_;
 
         //- Estimated number of faces an edge goes through
-        static const label nFacesPerSlaveEdge_;
+        static const Foam::label nFacesPerSlaveEdge_;
 
         //- Edge-face interaction escape limit
-        static const label edgeFaceEscapeLimit_;
+        static const Foam::label edgeFaceEscapeLimit_;
 
         //- Integral match point adjustment tolerance
-        static const scalar integralAdjTol_;
+        static const Foam::debug::tolerancesSwitch integralAdjTol_;
 
         //- Edge intersection master catch fraction
-        static const scalar edgeMasterCatchFraction_;
+        static const Foam::debug::tolerancesSwitch edgeMasterCatchFraction_;
 
         //- Edge intersection co-planar tolerance
-        static const scalar edgeCoPlanarTol_;
+        static const Foam::debug::tolerancesSwitch edgeCoPlanarTol_;
 
         //- Edge end cut-off tolerance
-        static const scalar edgeEndCutoffTol_;
+        static const Foam::debug::tolerancesSwitch edgeEndCutoffTol_;
 
 
 public:
diff --git a/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/slidingInterfaceProjectPoints.C b/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/slidingInterfaceProjectPoints.C
index a0c455ab1..0d44fdecf 100644
--- a/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/slidingInterfaceProjectPoints.C
+++ b/src/dynamicMesh/dynamicMesh/polyMeshModifiers/slidingInterface/slidingInterfaceProjectPoints.C
@@ -44,31 +44,42 @@ Author
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
-const Foam::scalar Foam::slidingInterface::pointMergeTol_
+const Foam::debug::tolerancesSwitch
+Foam::slidingInterface::pointMergeTol_
 (
-    debug::tolerances("slidingPointMergeTol", 0.2)
+    "slidingPointMergeTol",
+    0.2
 );
 
-const Foam::scalar Foam::slidingInterface::edgeMergeTol_
+const Foam::debug::tolerancesSwitch
+Foam::slidingInterface::edgeMergeTol_
 (
-    debug::tolerances("slidingEdgeMergeTol", 0.05)
+    "slidingEdgeMergeTol",
+    0.05
 );
 
-const Foam::scalar Foam::slidingInterface::integralAdjTol_
+const Foam::debug::tolerancesSwitch
+Foam::slidingInterface::integralAdjTol_
 (
-    debug::tolerances("slidingIntegralAdjTol", 0.15)
+    "slidingIntegralAdjTol",
+    0.15
 );
 
-const Foam::scalar Foam::slidingInterface::edgeMasterCatchFraction_
+const Foam::debug::tolerancesSwitch
+Foam::slidingInterface::edgeMasterCatchFraction_
 (
-    debug::tolerances("slidingEdgeMasterCatchFraction", 0.4)
+    "slidingEdgeMasterCatchFraction",
+    0.4
 );
 
-const Foam::scalar Foam::slidingInterface::edgeEndCutoffTol_
+const Foam::debug::tolerancesSwitch
+Foam::slidingInterface::edgeEndCutoffTol_
 (
-    debug::tolerances("slidingEdgeEndCutoffTol", 0.0001)
+    "slidingEdgeEndCutoffTol",
+    0.0001
 );
 
+
 const Foam::label Foam::slidingInterface::nFacesPerSlaveEdge_ = 5;
 
 // Increased limit for extreme 20-1 cutting.  HJ, 19/Dec/2008
@@ -105,7 +116,13 @@ bool Foam::slidingInterface::projectPoints() const
             << name() << " : "
             << "Projecting slave points onto master surface using ";
 
-        if (debug::optimisationSwitch("nSquaredProjection", 0) > 0)
+	const Foam::debug::optimisationSwitch nSquaredProjection
+	(
+	    "nSquaredProjection",
+	    0
+	);
+
+        if (nSquaredProjection() > 0)
         {
             Pout<< "N-squared point projection" << endl;
         }
@@ -358,7 +375,7 @@ bool Foam::slidingInterface::projectPoints() const
 
                     // Calculate the tolerance
                     const scalar mergeTol =
-                        integralAdjTol_*minSlavePointLength[pointI];
+                        integralAdjTol_()*minSlavePointLength[pointI];
 
                     // Adjust the hit
                     if (mag(nearPoint - missPoint) < mergeTol)
@@ -576,7 +593,7 @@ bool Foam::slidingInterface::projectPoints() const
 
                 // Calculate the tolerance
                 const scalar mergeTol =
-                    pointMergeTol_*
+                    pointMergeTol_()*
                     min
                     (
                         minSlavePointLength[pointI],
@@ -678,7 +695,7 @@ bool Foam::slidingInterface::projectPoints() const
                     minMasterFaceLength[slavePointFaceHits[pointI].hitObject()]
                 );
 
-            const scalar mergeTol = pointMergeTol_*mergeLength;
+            const scalar mergeTol = pointMergeTol_()*mergeLength;
 
             scalar minDistance = GREAT;
 
@@ -739,7 +756,7 @@ bool Foam::slidingInterface::projectPoints() const
 
                     // Calculate the tolerance
                     const scalar mergeTol =
-                        pointMergeTol_*
+                        pointMergeTol_()*
                         min
                         (
                             minSlavePointLength[pointI],
@@ -1051,7 +1068,7 @@ bool Foam::slidingInterface::projectPoints() const
             // Calculated as a combination of travel distance in projection and
             // edge length
             scalar slaveCatchDist =
-                edgeMasterCatchFraction_*edgeMag
+                edgeMasterCatchFraction_()*edgeMag
               + 0.5*
                 (
                     mag
@@ -1144,9 +1161,9 @@ bool Foam::slidingInterface::projectPoints() const
                     // Not a point hit, check for edge
                     if
                     (
-                        cutOnSlave > edgeEndCutoffTol_
-                     && cutOnSlave < 1.0 - edgeEndCutoffTol_ // check edge cut
-                     && distInEdgePlane < edgeMergeTol_*edgeMag // merge plane
+                        cutOnSlave > edgeEndCutoffTol_()
+		     && cutOnSlave < 1.0 - edgeEndCutoffTol_() // check edge cut
+		     && distInEdgePlane < edgeMergeTol_()*edgeMag // merge plane
                      && edgeLineHit.distance()
                       < min
                         (
diff --git a/src/dynamicMesh/dynamicMesh/polyTopoChange/polyTopoChange/polyTopoChange.C b/src/dynamicMesh/dynamicMesh/polyTopoChange/polyTopoChange/polyTopoChange.C
index c64f7c92c..89fb0a4dc 100644
--- a/src/dynamicMesh/dynamicMesh/polyTopoChange/polyTopoChange/polyTopoChange.C
+++ b/src/dynamicMesh/dynamicMesh/polyTopoChange/polyTopoChange/polyTopoChange.C
@@ -36,7 +36,7 @@ const Foam::label Foam::polyTopoChange::cellFraction = 10;
 
 int Foam::polyTopoChange::debug
 (
-    Foam::debug::debugSwitch("polyTopoChange", 0)
+    Foam::debug::debugSwitchFromDict("polyTopoChange", 0)
 );
 
 
diff --git a/src/dynamicMesh/meshMotion/mesquiteMotionSolver/mesquiteMotionSolver.C b/src/dynamicMesh/meshMotion/mesquiteMotionSolver/mesquiteMotionSolver.C
index 6196a4adc..e96cd0b9c 100644
--- a/src/dynamicMesh/meshMotion/mesquiteMotionSolver/mesquiteMotionSolver.C
+++ b/src/dynamicMesh/meshMotion/mesquiteMotionSolver/mesquiteMotionSolver.C
@@ -32,6 +32,7 @@ License
 #include "processorPolyPatch.H"
 #include "addToRunTimeSelectionTable.H"
 #include "polyMesh.H"
+#include "controlSwitches.H"
 
 #include "matchPoints.H"
 #include "DimensionedField.H"
@@ -1872,10 +1873,14 @@ void mesquiteMotionSolver::initParallelSurfaceSmoothing()
     const boundBox& box = mesh().bounds();
 
     // Fetch relative tolerance
-    scalar relTol = debug::tolerances("processorMatchTol", 1e-4);
+    const Foam::debug::tolerancesSwitch relTol
+    (
+	"processorMatchTol",
+	1e-4
+    );
 
     // Compute tolerance
-    scalar tol = relTol * box.mag();
+    scalar tol = relTol() * box.mag();
 
     // Wait for all transfers to complete.
     OPstream::waitRequests();
diff --git a/src/finiteArea/faMesh/faMesh.C b/src/finiteArea/faMesh/faMesh.C
index d46aef71d..2b357aaaa 100644
--- a/src/finiteArea/faMesh/faMesh.C
+++ b/src/finiteArea/faMesh/faMesh.C
@@ -38,6 +38,7 @@ Description
 #include "wedgeFaPatch.H"
 #include "faPatchData.H"
 #include "SortableList.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
@@ -48,11 +49,14 @@ namespace Foam
 
 Foam::word Foam::faMesh::meshSubDir = "faMesh";
 
-const bool Foam::faMesh::quadricsFit_
+const Foam::debug::optimisationSwitch
+Foam::faMesh::quadricsFit_
 (
-    debug::optimisationSwitch("quadricsFit", 0) > 0
+    "quadricsFit",
+    0
 );
 
+
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
 void Foam::faMesh::setPrimitiveMeshData()
@@ -1189,7 +1193,7 @@ const Foam::vectorField& Foam::faMesh::pointAreaNormals() const
     {
         calcPointAreaNormals();
 
-        if (quadricsFit_)
+        if (quadricsFit_() > 0)
         {
             calcPointAreaNormalsByQuadricsFit();
         }
diff --git a/src/finiteArea/faMesh/faMesh.H b/src/finiteArea/faMesh/faMesh.H
index c7fa001db..9ea68058c 100644
--- a/src/finiteArea/faMesh/faMesh.H
+++ b/src/finiteArea/faMesh/faMesh.H
@@ -174,7 +174,7 @@ class faMesh
 
     // Static Private Data
 
-        static const bool quadricsFit_;
+        static const Foam::debug::optimisationSwitch quadricsFit_;
 
 
     // Private Member Functions
diff --git a/src/finiteArea/faMesh/faPatches/constraint/cyclic/cyclicFaPatch.C b/src/finiteArea/faMesh/faPatches/constraint/cyclic/cyclicFaPatch.C
index 32c8177a9..299cd04bb 100644
--- a/src/finiteArea/faMesh/faPatches/constraint/cyclic/cyclicFaPatch.C
+++ b/src/finiteArea/faMesh/faPatches/constraint/cyclic/cyclicFaPatch.C
@@ -38,9 +38,10 @@ namespace Foam
 defineTypeNameAndDebug(cyclicFaPatch, 0);
 addToRunTimeSelectionTable(faPatch, cyclicFaPatch, dictionary);
 
-const scalar cyclicFaPatch::matchTol_
+const Foam::debug::tolerancesSwitch cyclicFaPatch::matchTol_
 (
-    debug::tolerances("patchFaceMatchTol", 1e-3)
+    "patchFaceMatchTol",
+    1e-3
 );
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
@@ -86,7 +87,7 @@ void Foam::cyclicFaPatch::calcTransforms()
             else if
             (
                 mag(magLe - nbrMagLe)/avLe
-              > matchTol_
+	      > matchTol_()
             )
             {
                 // Error in area matching.  Find largest error
@@ -102,7 +103,7 @@ void Foam::cyclicFaPatch::calcTransforms()
         }
 
         // Check for error in edge matching
-        if (maxMatchError > matchTol_)
+        if (maxMatchError > matchTol_())
         {
             label nbrEdgei = errorEdge + size()/2;
             scalar magLe = mag(half0Normals[errorEdge]);
@@ -120,7 +121,7 @@ void Foam::cyclicFaPatch::calcTransforms()
                 << "patch:" << name()
                 << " my area:" << magLe
                 << " neighbour area:" << nbrMagLe
-                << " matching tolerance:" << matchTol_
+                << " matching tolerance:" << matchTol_()
                 << endl
                 << "Mesh edge:" << start() + errorEdge
                 << endl
@@ -175,7 +176,7 @@ void cyclicFaPatch::makeWeights(scalarField& w) const
         if
         (
             mag(magL[edgei] - magL[edgei + sizeby2])/avL
-          > matchTol_
+	  > matchTol_()
         )
         {
             // Found error.  Look for largest matching error
@@ -197,7 +198,7 @@ void cyclicFaPatch::makeWeights(scalarField& w) const
     }
 
     // Check for error in matching
-    if (maxMatchError > polyPatch::matchTol_)
+    if (maxMatchError > polyPatch::matchTol_())
     {
         scalar avL = (magL[errorEdge] + magL[errorEdge + sizeby2])/2.0;
 
@@ -207,7 +208,7 @@ void cyclicFaPatch::makeWeights(scalarField& w) const
             << 100*mag(magL[errorEdge] - magL[errorEdge + sizeby2])/avL
             << "% -- possible edge ordering problem." << nl
             << "Cyclic area match tolerance = "
-            << polyPatch::matchTol_ << " patch: " << name()
+            << polyPatch::matchTol_() << " patch: " << name()
             << abort(FatalError);
     }
 }
diff --git a/src/finiteArea/faMesh/faPatches/constraint/cyclic/cyclicFaPatch.H b/src/finiteArea/faMesh/faPatches/constraint/cyclic/cyclicFaPatch.H
index ff6a7aacd..fbb227c33 100644
--- a/src/finiteArea/faMesh/faPatches/constraint/cyclic/cyclicFaPatch.H
+++ b/src/finiteArea/faMesh/faPatches/constraint/cyclic/cyclicFaPatch.H
@@ -38,6 +38,7 @@ SourceFiles
 #include "coupledFaPatch.H"
 #include "cyclicLduInterface.H"
 #include "cyclicPolyPatch.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -65,7 +66,7 @@ protected:
 
         //- Relative tolerance (for geometric matching). Is factor of
         //  maximum edge length per face.
-        static const scalar matchTol_;
+        static const debug::tolerancesSwitch matchTol_;
 
     // Protected Member functions
 
diff --git a/src/finiteArea/fields/faPatchFields/faPatchField/faPatchFields.C b/src/finiteArea/fields/faPatchFields/faPatchField/faPatchFields.C
index 6a9b05892..e528d1ea9 100644
--- a/src/finiteArea/fields/faPatchFields/faPatchField/faPatchFields.C
+++ b/src/finiteArea/fields/faPatchFields/faPatchField/faPatchFields.C
@@ -38,7 +38,7 @@ defineNamedTemplateTypeNameAndDebug(faPatchTypeField, 0);                     \
 template<>                                                                    \
 int faPatchTypeField::disallowDefaultFaPatchField                             \
 (                                                                             \
-    debug::debugSwitch("disallowDefaultFaPatchField", 0)                      \
+    debug::debugSwitchFromDict("disallowDefaultFaPatchField", 0)                      \
 );                                                                            \
 defineTemplateRunTimeSelectionTable(faPatchTypeField, patch);                 \
 defineTemplateRunTimeSelectionTable(faPatchTypeField, patchMapper);           \
diff --git a/src/finiteArea/fields/faePatchFields/faePatchField/faePatchFields.C b/src/finiteArea/fields/faePatchFields/faePatchField/faePatchFields.C
index 14c3b9ec1..2c744a959 100644
--- a/src/finiteArea/fields/faePatchFields/faePatchField/faePatchFields.C
+++ b/src/finiteArea/fields/faePatchFields/faePatchField/faePatchFields.C
@@ -41,7 +41,7 @@ defineNamedTemplateTypeNameAndDebug(faePatchTypeField, 0);                    \
 template<>                                                                    \
 int faePatchTypeField::disallowDefaultFaePatchField                           \
 (                                                                             \
-    debug::debugSwitch("disallowDefaultFaePatchField", 0)                     \
+    debug::debugSwitchFromDict("disallowDefaultFaePatchField", 0)                     \
 );                                                                            \
 defineTemplateRunTimeSelectionTable(faePatchTypeField, patch);                \
 defineTemplateRunTimeSelectionTable(faePatchTypeField, patchMapper);          \
diff --git a/src/finiteArea/finiteArea/faSchemes/faSchemes.C b/src/finiteArea/finiteArea/faSchemes/faSchemes.C
index f06d5857d..f16953c84 100644
--- a/src/finiteArea/finiteArea/faSchemes/faSchemes.C
+++ b/src/finiteArea/finiteArea/faSchemes/faSchemes.C
@@ -35,7 +35,7 @@ Description
 namespace Foam
 {
 
-int faSchemes::debug(Foam::debug::debugSwitch("faSchemes", false));
+int faSchemes::debug(Foam::debug::debugSwitchFromDict("faSchemes", false));
 
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
diff --git a/src/finiteVolume/cfdTools/general/adjustPhi/adjustPhi.C b/src/finiteVolume/cfdTools/general/adjustPhi/adjustPhi.C
index eb40b8ab7..eb158a7e8 100644
--- a/src/finiteVolume/cfdTools/general/adjustPhi/adjustPhi.C
+++ b/src/finiteVolume/cfdTools/general/adjustPhi/adjustPhi.C
@@ -29,6 +29,7 @@ License
 #include "processorFvsPatchFields.H"
 #include "inletOutletFvPatchFields.H"
 #include "fvc.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * Global Functions  * * * * * * * * * * * * * //
 
@@ -115,8 +116,11 @@ bool Foam::adjustPhi
 
         scalar massCorr = 1.0;
 
-        static const scalar closedDomainTol =
-            debug::tolerances("closedDomainTol", 1e-10);
+        static const Foam::debug::tolerancesSwitch closedDomainTol
+	(
+	    "closedDomainTol",
+	    1e-10
+	);
 
         if (mag(adjustableMassOut) > SMALL)
         {
@@ -125,7 +129,7 @@ bool Foam::adjustPhi
         else if
         (
             mag(fixedMassOut - massIn)
-          > closedDomainTol*Foam::max(1.0, mag(massIn))
+	  > closedDomainTol()*Foam::max(1.0, mag(massIn))
         )
         {
             phi.write();
diff --git a/src/finiteVolume/fields/fvPatchFields/fvPatchField/fvPatchFields.C b/src/finiteVolume/fields/fvPatchFields/fvPatchField/fvPatchFields.C
index 172a6e486..4adfef227 100644
--- a/src/finiteVolume/fields/fvPatchFields/fvPatchField/fvPatchFields.C
+++ b/src/finiteVolume/fields/fvPatchFields/fvPatchField/fvPatchFields.C
@@ -38,7 +38,7 @@ defineNamedTemplateTypeNameAndDebug(fvPatchTypeField, 0);                     \
 template<>                                                                    \
 int fvPatchTypeField::disallowGenericFvPatchField                             \
 (                                                                             \
-    debug::debugSwitch("disallowGenericFvPatchField", 0)                      \
+    debug::debugSwitchFromDict("disallowGenericFvPatchField", 0)                      \
 );                                                                            \
 defineTemplateRunTimeSelectionTable(fvPatchTypeField, patch);                 \
 defineTemplateRunTimeSelectionTable(fvPatchTypeField, patchMapper);           \
diff --git a/src/finiteVolume/fields/fvsPatchFields/fvsPatchField/fvsPatchFields.C b/src/finiteVolume/fields/fvsPatchFields/fvsPatchField/fvsPatchFields.C
index 5c049a8ec..8e89526f9 100644
--- a/src/finiteVolume/fields/fvsPatchFields/fvsPatchField/fvsPatchFields.C
+++ b/src/finiteVolume/fields/fvsPatchFields/fvsPatchField/fvsPatchFields.C
@@ -38,7 +38,7 @@ defineNamedTemplateTypeNameAndDebug(fvsPatchTypeField, 0);                    \
 template<>                                                                    \
 int fvsPatchTypeField::disallowDefaultFvsPatchField                           \
 (                                                                             \
-    debug::debugSwitch("disallowDefaultFvsPatchField", 0)                     \
+    debug::debugSwitchFromDict("disallowDefaultFvsPatchField", 0)                     \
 );                                                                            \
 defineTemplateRunTimeSelectionTable(fvsPatchTypeField, patch);                \
 defineTemplateRunTimeSelectionTable(fvsPatchTypeField, patchMapper);          \
diff --git a/src/finiteVolume/finiteVolume/fvSchemes/fvSchemes.C b/src/finiteVolume/finiteVolume/fvSchemes/fvSchemes.C
index c8b2c6e50..4a9b35359 100644
--- a/src/finiteVolume/finiteVolume/fvSchemes/fvSchemes.C
+++ b/src/finiteVolume/finiteVolume/fvSchemes/fvSchemes.C
@@ -29,7 +29,7 @@ License
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
 
-int Foam::fvSchemes::debug(Foam::debug::debugSwitch("fvSchemes", false));
+int Foam::fvSchemes::debug(Foam::debug::debugSwitchFromDict("fvSchemes", false));
 
 // * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * * //
 
diff --git a/src/finiteVolume/fvMesh/fvPatches/constraint/cyclic/cyclicFvPatch.C b/src/finiteVolume/fvMesh/fvPatches/constraint/cyclic/cyclicFvPatch.C
index e88c0617d..9f3695c90 100644
--- a/src/finiteVolume/fvMesh/fvPatches/constraint/cyclic/cyclicFvPatch.C
+++ b/src/finiteVolume/fvMesh/fvPatches/constraint/cyclic/cyclicFvPatch.C
@@ -63,7 +63,7 @@ void cyclicFvPatch::makeWeights(scalarField& w) const
         if
         (
             mag(magFa[facei] - magFa[facei + sizeby2])/avFa
-          > polyPatch::matchTol_
+          > polyPatch::matchTol_()
         )
         {
             // Found error.  Look for largest matching error
@@ -85,7 +85,7 @@ void cyclicFvPatch::makeWeights(scalarField& w) const
     }
 
     // Check for error in matching
-    if (maxMatchError > polyPatch::matchTol_)
+    if (maxMatchError > polyPatch::matchTol_())
     {
         scalar avFa = (magFa[errorFace] + magFa[errorFace + sizeby2])/2.0;
 
@@ -95,7 +95,7 @@ void cyclicFvPatch::makeWeights(scalarField& w) const
             << 100*mag(magFa[errorFace] - magFa[errorFace + sizeby2])/avFa
             << "% -- possible face ordering problem." << nl
             << "Cyclic area match tolerance = "
-            << polyPatch::matchTol_ << " patch: " << name()
+            << polyPatch::matchTol_() << " patch: " << name()
             << abort(FatalError);
     }
 }
diff --git a/src/finiteVolume/interpolation/interpolation/interpolationCellPoint/cellPointWeight/cellPointWeight.C b/src/finiteVolume/interpolation/interpolation/interpolationCellPoint/cellPointWeight/cellPointWeight.C
index 28207c123..b22f76113 100644
--- a/src/finiteVolume/interpolation/interpolation/interpolationCellPoint/cellPointWeight/cellPointWeight.C
+++ b/src/finiteVolume/interpolation/interpolation/interpolationCellPoint/cellPointWeight/cellPointWeight.C
@@ -28,7 +28,7 @@ License
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
-int Foam::cellPointWeight::debug(debug::debugSwitch("cellPointWeight", 0));
+int Foam::cellPointWeight::debug(debug::debugSwitchFromDict("cellPointWeight", 0));
 Foam::scalar Foam::cellPointWeight::tol(SMALL);
 
 // * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * * //
diff --git a/src/foam/algorithms/MeshWave/FaceCellWave.C b/src/foam/algorithms/MeshWave/FaceCellWave.C
index 96ee1f2cf..d929ed525 100644
--- a/src/foam/algorithms/MeshWave/FaceCellWave.C
+++ b/src/foam/algorithms/MeshWave/FaceCellWave.C
@@ -36,15 +36,17 @@ License
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
 template <class Type>
-Foam::scalar Foam::FaceCellWave<Type>::geomTol_
+Foam::debug::tolerancesSwitch Foam::FaceCellWave<Type>::geomTol_
 (
-    debug::tolerances("FaceCellWaveGeomTol", 1e-6)
+    "FaceCellWaveGeomTol",
+    1e-6
 );
 
 template <class Type>
-Foam::scalar Foam::FaceCellWave<Type>::propagationTol_
+Foam::debug::tolerancesSwitch Foam::FaceCellWave<Type>::propagationTol_
 (
-    debug::tolerances("FaceCellWavePropagationTol", 0.01)
+    "FaceCellWavePropagationTol",
+    0.01
 );
 
 // Write to ostream
@@ -264,7 +266,7 @@ void Foam::FaceCellWave<Type>::checkCyclic(const polyPatch& patch) const
         label i1 = patch.start() + patchFaceI;
         label i2 = i1 + cycOffset;
 
-        if (!allFaceInfo_[i1].sameGeometry(mesh_, allFaceInfo_[i2], geomTol_))
+        if (!allFaceInfo_[i1].sameGeometry(mesh_, allFaceInfo_[i2], geomTol_()))
         {
             FatalErrorIn("FaceCellWave<Type>::checkCyclic(const polyPatch&)")
                 << "problem: i:" << i1 << "  otheri:" << i2
@@ -359,7 +361,7 @@ void Foam::FaceCellWave<Type>::mergeFaceInfo
             (
                 meshFaceI,
                 neighbourWallInfo,
-                propagationTol_,
+                propagationTol_(),
                 currentWallInfo
             );
         }
@@ -934,7 +936,7 @@ Foam::label Foam::FaceCellWave<Type>::faceToCell()
                 cellI,
                 faceI,
                 neighbourWallInfo,
-                propagationTol_,
+                propagationTol_(),
                 currentWallInfo
             );
         }
@@ -952,7 +954,7 @@ Foam::label Foam::FaceCellWave<Type>::faceToCell()
                     cellI,
                     faceI,
                     neighbourWallInfo,
-                    propagationTol_,
+                    propagationTol_(),
                     currentWallInfo2
                 );
             }
@@ -1018,7 +1020,7 @@ Foam::label Foam::FaceCellWave<Type>::cellToFace()
                     faceI,
                     cellI,
                     neighbourWallInfo,
-                    propagationTol_,
+                    propagationTol_(),
                     currentWallInfo
                 );
             }
diff --git a/src/foam/algorithms/MeshWave/FaceCellWave.H b/src/foam/algorithms/MeshWave/FaceCellWave.H
index 6535f338e..c42ff5f9d 100644
--- a/src/foam/algorithms/MeshWave/FaceCellWave.H
+++ b/src/foam/algorithms/MeshWave/FaceCellWave.H
@@ -50,6 +50,7 @@ SourceFiles
 #include "boolList.H"
 #include "labelList.H"
 #include "primitiveFieldsFwd.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -267,9 +268,8 @@ class FaceCellWave
 
       // Private static data
 
-            static scalar geomTol_;
-            static scalar propagationTol_;
-
+            static debug::tolerancesSwitch geomTol_;
+            static debug::tolerancesSwitch propagationTol_;
 
 public:
 
@@ -278,7 +278,7 @@ public:
         //- Access to tolerance
         static scalar propagationTol()
         {
-            return propagationTol_;
+            return propagationTol_();
         }
 
         //- Change tolerance
diff --git a/src/foam/algorithms/octree/octree/octreeDataBoundBox.C b/src/foam/algorithms/octree/octree/octreeDataBoundBox.C
index 45480f99e..66a0b92d8 100644
--- a/src/foam/algorithms/octree/octree/octreeDataBoundBox.C
+++ b/src/foam/algorithms/octree/octree/octreeDataBoundBox.C
@@ -35,12 +35,13 @@ License
 
 defineTypeNameAndDebug(Foam::octreeDataBoundBox, 0);
 
-Foam::scalar Foam::octreeDataBoundBox::tol
+const Foam::debug::tolerancesSwitch
+Foam::octreeDataBoundBox::tol
 (
-    debug::tolerances("octreeDataBoundBoxTol", 1e-6)
+    "octreeDataBoundBoxTol",
+    1e-6
 );
 
-
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
 Foam::octreeDataBoundBox::octreeDataBoundBox
diff --git a/src/foam/algorithms/octree/octree/octreeDataBoundBox.H b/src/foam/algorithms/octree/octree/octreeDataBoundBox.H
index 77ce6ccae..36e489651 100644
--- a/src/foam/algorithms/octree/octree/octreeDataBoundBox.H
+++ b/src/foam/algorithms/octree/octree/octreeDataBoundBox.H
@@ -49,6 +49,7 @@ Author
 #include "point.H"
 #include "className.H"
 #include "linePointRef.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -69,7 +70,7 @@ class octreeDataBoundBox
     // Static data
 
         //- Tolerance on linear dimensions
-        static scalar tol;
+        static const Foam::debug::tolerancesSwitch tol;
 
 
     // Private data
diff --git a/src/foam/db/IOstreams/IOstreams.C b/src/foam/db/IOstreams/IOstreams.C
index aa9523d4c..d539b126f 100644
--- a/src/foam/db/IOstreams/IOstreams.C
+++ b/src/foam/db/IOstreams/IOstreams.C
@@ -36,8 +36,13 @@ namespace Foam
 
 const IOstream::versionNumber IOstream::originalVersion(0.5);
 const IOstream::versionNumber IOstream::currentVersion(2.0);
-unsigned int IOstream::precision_(debug::infoSwitch("writePrecision", 6));
 
+Foam::debug::infoSwitch
+IOstream::precision_
+(
+    "writePrecision",
+    6
+);
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 // Global IO streams
diff --git a/src/foam/db/IOstreams/IOstreams/IOstream.H b/src/foam/db/IOstreams/IOstreams/IOstream.H
index 937e2bd99..f0cd3ae87 100644
--- a/src/foam/db/IOstreams/IOstreams/IOstream.H
+++ b/src/foam/db/IOstreams/IOstreams/IOstream.H
@@ -49,9 +49,13 @@ SourceFiles
 #include "scalar.H"
 #include "fileName.H"
 #include "InfoProxy.H"
+#include "controlSwitches.H"
+
+#include "debug.H"
 
 #include <iostream>
 
+
 #if __GNUC__ < 3
 #   define ios_base ios
 #endif
@@ -211,8 +215,7 @@ public:
         static const versionNumber currentVersion;
 
         //- Default precision
-        static unsigned int precision_;
-
+        static Foam::debug::infoSwitch precision_;
 
 private:
 
@@ -466,13 +469,13 @@ public:
             //- Return the default precision
             static unsigned int defaultPrecision()
             {
-                return precision_;
+                return precision_();
             }
 
             //- Reset the default precision (and return old precision)
             static unsigned int defaultPrecision(unsigned int p)
             {
-                unsigned int precision0 = precision_;
+                unsigned int precision0 = precision_();
                 precision_ = p;
                 return precision0;
             }
diff --git a/src/foam/db/IOstreams/Pstreams/Pstream.C b/src/foam/db/IOstreams/Pstreams/Pstream.C
index a57613aa0..1181f6a48 100644
--- a/src/foam/db/IOstreams/Pstreams/Pstream.C
+++ b/src/foam/db/IOstreams/Pstreams/Pstream.C
@@ -227,16 +227,20 @@ Foam::List<Foam::Pstream::commsStruct> Foam::Pstream::treeCommunication_(0);
 // Should compact transfer be used in which floats replace doubles
 // reducing the bandwidth requirement at the expense of some loss
 // in accuracy
-bool Foam::Pstream::floatTransfer
+const Foam::debug::optimisationSwitch
+Foam::Pstream::floatTransfer
 (
-    debug::optimisationSwitch("floatTransfer", 0)
+    "floatTransfer",
+    0
 );
 
 // Number of processors at which the reduce algorithm changes from linear to
 // tree
-int Foam::Pstream::nProcsSimpleSum
+const Foam::debug::optimisationSwitch
+Foam::Pstream::nProcsSimpleSum
 (
-    debug::optimisationSwitch("nProcsSimpleSum", 16)
+    "nProcsSimpleSum",
+    16
 );
 
 // Default commsType
@@ -245,5 +249,4 @@ Foam::Pstream::commsTypes Foam::Pstream::defaultCommsType
     commsTypeNames.read(debug::optimisationSwitches().lookup("commsType"))
 );
 
-
 // ************************************************************************* //
diff --git a/src/foam/db/IOstreams/Pstreams/Pstream.H b/src/foam/db/IOstreams/Pstreams/Pstream.H
index 523d40ea1..211eff246 100644
--- a/src/foam/db/IOstreams/Pstreams/Pstream.H
+++ b/src/foam/db/IOstreams/Pstreams/Pstream.H
@@ -45,6 +45,7 @@ SourceFiles
 #include "HashTable.H"
 #include "string.H"
 #include "NamedEnum.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -224,11 +225,11 @@ public:
         //- Should compact transfer be used in which floats replace doubles
         //  reducing the bandwidth requirement at the expense of some loss
         //  in accuracy
-        static bool floatTransfer;
+        static const Foam::debug::optimisationSwitch floatTransfer;
 
         //- Number of processors at which the sum algorithm changes from linear
         //  to tree
-        static int nProcsSimpleSum;
+        static const Foam::debug::optimisationSwitch nProcsSimpleSum;
 
         //- Default commsType
         static commsTypes defaultCommsType;
diff --git a/src/foam/db/IOstreams/Pstreams/PstreamCombineReduceOps.H b/src/foam/db/IOstreams/Pstreams/PstreamCombineReduceOps.H
index 5c8a260af..8aea0dd64 100644
--- a/src/foam/db/IOstreams/Pstreams/PstreamCombineReduceOps.H
+++ b/src/foam/db/IOstreams/Pstreams/PstreamCombineReduceOps.H
@@ -62,7 +62,7 @@ void combineReduce
 template <class T, class CombineOp>
 void combineReduce(T& Value, const CombineOp& cop)
 {
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         Pstream::combineGather(Pstream::linearCommunication(), Value, cop);
         Pstream::combineScatter(Pstream::linearCommunication(), Value);
diff --git a/src/foam/db/IOstreams/Pstreams/PstreamReduceOps.H b/src/foam/db/IOstreams/Pstreams/PstreamReduceOps.H
index d46b6cf19..8ea79c8bc 100644
--- a/src/foam/db/IOstreams/Pstreams/PstreamReduceOps.H
+++ b/src/foam/db/IOstreams/Pstreams/PstreamReduceOps.H
@@ -58,7 +58,7 @@ void reduce
     const BinaryOp& bop
 )
 {
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         reduce(Pstream::linearCommunication(), Value, bop);
     }
@@ -79,7 +79,7 @@ T returnReduce
 {
     T WorkValue(Value);
 
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         reduce(Pstream::linearCommunication(), WorkValue, bop);
     }
diff --git a/src/foam/db/IOstreams/Pstreams/combineGatherScatter.C b/src/foam/db/IOstreams/Pstreams/combineGatherScatter.C
index 826b00587..289731a47 100644
--- a/src/foam/db/IOstreams/Pstreams/combineGatherScatter.C
+++ b/src/foam/db/IOstreams/Pstreams/combineGatherScatter.C
@@ -129,7 +129,7 @@ void Pstream::combineGather
 template <class T, class CombineOp>
 void Pstream::combineGather(T& Value, const CombineOp& cop)
 {
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         combineGather(Pstream::linearCommunication(), Value, cop);
     }
@@ -207,7 +207,7 @@ void Pstream::combineScatter(const List<Pstream::commsStruct>& comms, T& Value)
 template <class T>
 void Pstream::combineScatter(T& Value)
 {
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         combineScatter(Pstream::linearCommunication(), Value);
     }
@@ -313,7 +313,7 @@ void Pstream::listCombineGather
 template <class T, class CombineOp>
 void Pstream::listCombineGather(List<T>& Values, const CombineOp& cop)
 {
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         listCombineGather(Pstream::linearCommunication(), Values, cop);
     }
@@ -395,7 +395,7 @@ void Pstream::listCombineScatter
 template <class T>
 void Pstream::listCombineScatter(List<T>& Values)
 {
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         listCombineScatter(Pstream::linearCommunication(), Values);
     }
@@ -480,7 +480,7 @@ void Pstream::mapCombineGather
 template <class Container, class CombineOp>
 void Pstream::mapCombineGather(Container& Values, const CombineOp& cop)
 {
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         mapCombineGather(Pstream::linearCommunication(), Values, cop);
     }
@@ -536,7 +536,7 @@ void Pstream::mapCombineScatter
 template <class Container>
 void Pstream::mapCombineScatter(Container& Values)
 {
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         mapCombineScatter(Pstream::linearCommunication(), Values);
     }
diff --git a/src/foam/db/IOstreams/Pstreams/gatherScatter.C b/src/foam/db/IOstreams/Pstreams/gatherScatter.C
index 73b7a2ec2..5f9e99bf9 100644
--- a/src/foam/db/IOstreams/Pstreams/gatherScatter.C
+++ b/src/foam/db/IOstreams/Pstreams/gatherScatter.C
@@ -103,7 +103,7 @@ void Pstream::gather
 template <class T, class BinaryOp>
 void Pstream::gather(T& Value, const BinaryOp& bop)
 {
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         gather(Pstream::linearCommunication(), Value, bop);
     }
@@ -168,7 +168,7 @@ void Pstream::scatter(const List<Pstream::commsStruct>& comms, T& Value)
 template <class T>
 void Pstream::scatter(T& Value)
 {
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         scatter(Pstream::linearCommunication(), Value);
     }
diff --git a/src/foam/db/IOstreams/Pstreams/gatherScatterList.C b/src/foam/db/IOstreams/Pstreams/gatherScatterList.C
index 1c6a36b78..4b7dbc5df 100644
--- a/src/foam/db/IOstreams/Pstreams/gatherScatterList.C
+++ b/src/foam/db/IOstreams/Pstreams/gatherScatterList.C
@@ -180,7 +180,7 @@ void Pstream::gatherList
 template <class T>
 void Pstream::gatherList(List<T>& Values)
 {
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         gatherList(Pstream::linearCommunication(), Values);
     }
@@ -305,7 +305,7 @@ void Pstream::scatterList
 template <class T>
 void Pstream::scatterList(List<T>& Values)
 {
-    if (Pstream::nProcs() < Pstream::nProcsSimpleSum)
+    if (Pstream::nProcs() < Pstream::nProcsSimpleSum())
     {
         scatterList(Pstream::linearCommunication(), Values);
     }
diff --git a/src/foam/db/IOstreams/Sstreams/OSstreamI.H b/src/foam/db/IOstreams/Sstreams/OSstreamI.H
index 1ac2412ff..018ec331b 100644
--- a/src/foam/db/IOstreams/Sstreams/OSstreamI.H
+++ b/src/foam/db/IOstreams/Sstreams/OSstreamI.H
@@ -44,7 +44,7 @@ inline Foam::OSstream::OSstream
     {
         setOpened();
         setGood();
-        os_.precision(precision_);
+        os_.precision(precision_());
     }
     else
     {
diff --git a/src/foam/db/Time/TimeIO.C b/src/foam/db/Time/TimeIO.C
index 2965665a4..aef237c56 100644
--- a/src/foam/db/Time/TimeIO.C
+++ b/src/foam/db/Time/TimeIO.C
@@ -33,6 +33,9 @@ License
 
 void Foam::Time::readDict()
 {
+    if (debug)
+	Info << "Time::readDict(): reading " << controlDict_.name() << endl;
+
     if (!deltaTchanged_)
     {
         deltaT_ = readScalar(controlDict_.lookup("deltaT"));
diff --git a/src/foam/db/dictionary/functionEntries/includeEntry/includeEntry.C b/src/foam/db/dictionary/functionEntries/includeEntry/includeEntry.C
index fed644cd3..ade17a0db 100644
--- a/src/foam/db/dictionary/functionEntries/includeEntry/includeEntry.C
+++ b/src/foam/db/dictionary/functionEntries/includeEntry/includeEntry.C
@@ -37,7 +37,12 @@ const Foam::word Foam::functionEntries::includeEntry::typeName
 
 // Don't lookup the debug switch here as the debug switch dictionary
 // might include includeEntry
-int Foam::functionEntries::includeEntry::debug(0);
+Foam::debug::debugSwitch
+Foam::functionEntries::includeEntry::debug
+(
+    "includeEntry",
+    0
+);
 
 namespace Foam
 {
diff --git a/src/foam/db/dictionary/functionEntries/includeIfPresentEntry/includeIfPresentEntry.C b/src/foam/db/dictionary/functionEntries/includeIfPresentEntry/includeIfPresentEntry.C
index d680bcb1a..119ad395a 100644
--- a/src/foam/db/dictionary/functionEntries/includeIfPresentEntry/includeIfPresentEntry.C
+++ b/src/foam/db/dictionary/functionEntries/includeIfPresentEntry/includeIfPresentEntry.C
@@ -37,7 +37,12 @@ const Foam::word Foam::functionEntries::includeIfPresentEntry::typeName
 
 // Don't lookup the debug switch here as the debug switch dictionary
 // might include includeIfPresentEntry
-int Foam::functionEntries::includeIfPresentEntry::debug(0);
+Foam::debug::debugSwitch
+Foam::functionEntries::includeIfPresentEntry::debug
+(
+    "includeIfPresentEntry",
+    0
+);
 
 namespace Foam
 {
diff --git a/src/foam/db/dictionary/functionEntries/inputModeEntry/inputModeEntry.C b/src/foam/db/dictionary/functionEntries/inputModeEntry/inputModeEntry.C
index 308523d4f..3d7a562ba 100644
--- a/src/foam/db/dictionary/functionEntries/inputModeEntry/inputModeEntry.C
+++ b/src/foam/db/dictionary/functionEntries/inputModeEntry/inputModeEntry.C
@@ -36,7 +36,12 @@ const Foam::word Foam::functionEntries::inputModeEntry::typeName
 
 // Don't lookup the debug switch here as the debug switch dictionary
 // might include inputModeEntries
-int Foam::functionEntries::inputModeEntry::debug(0);
+Foam::debug::debugSwitch
+Foam::functionEntries::inputModeEntry::debug
+(
+    "inputModeEntry",
+    0
+);
 
 Foam::functionEntries::inputModeEntry::inputMode
     Foam::functionEntries::inputModeEntry::mode_(MERGE);
diff --git a/src/foam/db/dictionary/functionEntries/removeEntry/removeEntry.C b/src/foam/db/dictionary/functionEntries/removeEntry/removeEntry.C
index 85b3aec44..9dadf8fed 100644
--- a/src/foam/db/dictionary/functionEntries/removeEntry/removeEntry.C
+++ b/src/foam/db/dictionary/functionEntries/removeEntry/removeEntry.C
@@ -39,7 +39,12 @@ const Foam::word Foam::functionEntries::removeEntry::typeName
 
 // Don't lookup the debug switch here as the debug switch dictionary
 // might include removeEntry
-int Foam::functionEntries::removeEntry::debug(0);
+Foam::debug::debugSwitch
+Foam::functionEntries::removeEntry::debug
+(
+    "removeEntry",
+    0
+);
 
 namespace Foam
 {
diff --git a/src/foam/db/error/messageStream.C b/src/foam/db/error/messageStream.C
index d08686c93..0056946e2 100644
--- a/src/foam/db/error/messageStream.C
+++ b/src/foam/db/error/messageStream.C
@@ -29,7 +29,7 @@ License
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
-int Foam::messageStream::level(Foam::debug::debugSwitch("level", 2));
+int Foam::messageStream::level(Foam::debug::debugSwitchFromDict("level", 2));
 
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/src/foam/db/functionObjects/functionObject/functionObject.C b/src/foam/db/functionObjects/functionObject/functionObject.C
index c7aa5bad0..0790ecfbe 100644
--- a/src/foam/db/functionObjects/functionObject/functionObject.C
+++ b/src/foam/db/functionObjects/functionObject/functionObject.C
@@ -30,7 +30,7 @@ License
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
 defineRunTimeSelectionTable(Foam::functionObject, dictionary);
-int Foam::functionObject::debug(Foam::debug::debugSwitch("functionObject", 0));
+int Foam::functionObject::debug(Foam::debug::debugSwitchFromDict("functionObject", 0));
 
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
diff --git a/src/foam/db/regIOobject/regIOobject.C b/src/foam/db/regIOobject/regIOobject.C
index 09388e410..1d9167ab0 100644
--- a/src/foam/db/regIOobject/regIOobject.C
+++ b/src/foam/db/regIOobject/regIOobject.C
@@ -31,12 +31,13 @@ License
 
 defineTypeNameAndDebug(Foam::regIOobject, 0);
 
-int Foam::regIOobject::fileModificationSkew
+const Foam::debug::optimisationSwitch
+Foam::regIOobject::fileModificationSkew
 (
-    Foam::debug::optimisationSwitch("fileModificationSkew", 30)
+    "fileModificationSkew",
+    30
 );
 
-
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
 // Construct from IOobject
diff --git a/src/foam/db/regIOobject/regIOobject.H b/src/foam/db/regIOobject/regIOobject.H
index 3832caed5..fddd39645 100644
--- a/src/foam/db/regIOobject/regIOobject.H
+++ b/src/foam/db/regIOobject/regIOobject.H
@@ -41,6 +41,7 @@ SourceFiles
 #include "IOobject.H"
 #include "typeInfo.H"
 #include "OSspecific.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -93,7 +94,7 @@ public:
         //- Runtime type information
         TypeName("regIOobject");
 
-        static int fileModificationSkew;
+        static const Foam::debug::optimisationSwitch fileModificationSkew;
 
 
     // Constructors
diff --git a/src/foam/db/regIOobject/regIOobjectRead.C b/src/foam/db/regIOobject/regIOobjectRead.C
index e82db74af..f69ed0975 100644
--- a/src/foam/db/regIOobject/regIOobjectRead.C
+++ b/src/foam/db/regIOobject/regIOobjectRead.C
@@ -154,7 +154,7 @@ bool Foam::regIOobject::modified() const
     return
     (
         lastModified_
-     && lastModified(filePath()) > (lastModified_ + fileModificationSkew)
+     && lastModified(filePath()) > (lastModified_ + fileModificationSkew())
     );
 }
 
@@ -167,7 +167,7 @@ bool Foam::regIOobject::readIfModified()
 
         bool readFile = false;
 
-        if (newTimeStamp > (lastModified_ + fileModificationSkew))
+        if (newTimeStamp > (lastModified_ + fileModificationSkew()))
         {
             readFile = true;
         }
diff --git a/src/foam/db/scalarRange/scalarRange.C b/src/foam/db/scalarRange/scalarRange.C
index e077e6faa..7d8e2db70 100644
--- a/src/foam/db/scalarRange/scalarRange.C
+++ b/src/foam/db/scalarRange/scalarRange.C
@@ -28,7 +28,7 @@ License
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
-int Foam::scalarRange::debug(::Foam::debug::debugSwitch("scalarRange", 0));
+int Foam::scalarRange::debug(::Foam::debug::debugSwitchFromDict("scalarRange", 0));
 
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
diff --git a/src/foam/db/typeInfo/className.H b/src/foam/db/typeInfo/className.H
index c50f2cd6e..bccc1f403 100644
--- a/src/foam/db/typeInfo/className.H
+++ b/src/foam/db/typeInfo/className.H
@@ -31,6 +31,7 @@ Description
 
 #include "word.H"
 #include "debug.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -72,13 +73,13 @@ public:                                                                       \
 //  Also declares debug information.
 #define ClassName(TypeNameString)                                             \
     ClassNameNoDebug(TypeNameString);                                         \
-    static int debug
+    static Foam::debug::debugSwitch debug;
 
 //- Add typeName information from argument @a TypeNameString to a namespace.
 //  Also declares debug information.
 #define NamespaceName(TypeNameString)                                         \
     NamespaceNameNoDebug(TypeNameString);                                     \
-    extern int debug
+    extern Foam::debug::debugSwitch debug;
 
 //- Add typeName information from argument @a TypeNameString to a template class.
 //  Also declares debug information.
@@ -133,21 +134,23 @@ public:                                                                       \
 
 //- Define the debug information, lookup as @a Name
 #define defineDebugSwitchWithName(Type, Name, DebugSwitch)                    \
-    int Type::debug(::Foam::debug::debugSwitch(Name, DebugSwitch))
+    ::Foam::debug::debugSwitch                                                \
+    Type::debug(std::string(Name), DebugSwitch)
 
 //- Define the debug information
 #define defineDebugSwitch(Type, DebugSwitch)                                  \
-    defineDebugSwitchWithName(Type, Type::typeName_(), DebugSwitch)
+    defineDebugSwitchWithName(Type, Type::typeName_(), DebugSwitch);
 
 #ifdef __INTEL_COMPILER
 //- Define the debug information for templates, lookup as @a Name
 # define defineTemplateDebugSwitchWithName(Type, Name, DebugSwitch)           \
-    defineDebugSwitchWithName(Type, Name, DebugSwitch)
+    defineDebugSwitchWithName(Type, Name, DebugSwitch);
+
 #else
 //- Define the debug information for templates, lookup as @a Name
 # define defineTemplateDebugSwitchWithName(Type, Name, DebugSwitch)           \
     template<>                                                                \
-    defineDebugSwitchWithName(Type, Name, DebugSwitch)
+    defineDebugSwitchWithName(Type, Name, DebugSwitch);
 #endif
 
 //- Define the debug information for templates
diff --git a/src/foam/global/JobInfo/JobInfo.C b/src/foam/global/JobInfo/JobInfo.C
index a1eddf347..410512492 100644
--- a/src/foam/global/JobInfo/JobInfo.C
+++ b/src/foam/global/JobInfo/JobInfo.C
@@ -31,7 +31,13 @@ License
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
-bool Foam::JobInfo::writeJobInfo(Foam::debug::infoSwitch("writeJobInfo", 0));
+Foam::debug::infoSwitch
+Foam::JobInfo::writeJobInfo
+(
+    "writeJobInfo",
+    0
+);
+
 Foam::JobInfo Foam::jobInfo;
 
 
diff --git a/src/foam/global/JobInfo/JobInfo.H b/src/foam/global/JobInfo/JobInfo.H
index f2ffde334..2ec736669 100644
--- a/src/foam/global/JobInfo/JobInfo.H
+++ b/src/foam/global/JobInfo/JobInfo.H
@@ -42,6 +42,7 @@ SourceFiles
 #include "dictionary.H"
 #include "fileName.H"
 #include "cpuTime.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -70,7 +71,7 @@ class JobInfo
 public:
 
     static bool constructed;
-    static bool writeJobInfo;
+    static Foam::debug::infoSwitch writeJobInfo;
 
     // Constructors
 
diff --git a/src/foam/global/argList/argList.C b/src/foam/global/argList/argList.C
index 7cfa7d6eb..445a86d73 100644
--- a/src/foam/global/argList/argList.C
+++ b/src/foam/global/argList/argList.C
@@ -32,6 +32,8 @@ License
 #include "IOobject.H"
 #include "JobInfo.H"
 #include "labelList.H"
+#include "SortableList.H"
+//#include "dimensionedConstants.H"
 
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
@@ -39,8 +41,32 @@ License
 Foam::SLList<Foam::string> Foam::argList::validArgs;
 Foam::HashTable<Foam::string> Foam::argList::validOptions;
 Foam::HashTable<Foam::string> Foam::argList::validParOptions;
+Foam::word Foam::argList::appDictName_("");
 bool Foam::argList::bannerEnabled(true);
 
+template<>
+const char*
+Foam::NamedEnum
+<
+    Foam::debug::globalControlDictSwitchSet,
+    5
+>::names[] =
+{
+    "DebugSwitches",
+    "InfoSwitches",
+    "OptimisationSwitches",
+    "Tolerances",
+    "DimensionedConstants"
+};
+
+
+const Foam::NamedEnum
+<
+    Foam::debug::globalControlDictSwitchSet,
+    5
+>
+Foam::argList::globalControlDictSwitchSetNames_;
+
 
 Foam::argList::initValidTables::initValidTables()
 {
@@ -48,7 +74,11 @@ Foam::argList::initValidTables::initValidTables()
     validOptions.set("parallel", "");
     validParOptions.set("parallel", "");
     validOptions.set("noFunctionObjects", "");
-
+    validOptions.set("DebugSwitches", "key=value");
+    validOptions.set("InfoSwitches", "key=value");
+    validOptions.set("OptimisationSwitches", "key=value");
+    validOptions.set("Tolerances", "key=value");
+    validOptions.set("DimensionedConstants", "key=value");
     Pstream::addValidParOptions(validParOptions);
 }
 
@@ -308,6 +338,7 @@ Foam::argList::argList
         FatalError.exit();
     }
 
+    // From here, we consider the command-line arguments to be valid
 
     string dateString = clock::date();
     string timeString = clock::clockTime();
@@ -328,6 +359,19 @@ Foam::argList::argList
             << endl;
     }
 
+    // Command-line override for central controlDict's variables
+    forAll(globalControlDictSwitchSetNames_, cI)
+    {
+	word switchSetName = globalControlDictSwitchSetNames_.names[cI];
+
+	if( optionFound(switchSetName) )
+	    Foam::debug::updateCentralDictVars
+	    (
+		globalControlDictSwitchSetNames_[switchSetName],
+		option(switchSetName)
+	    );
+    }
+
     jobInfo.add("startDate", dateString);
     jobInfo.add("startTime", timeString);
     jobInfo.add("userName", userName());
@@ -540,7 +584,7 @@ Foam::argList::argList
             Info<< "Slaves : " << slaveProcs << nl
                 << "Pstream initialized with:" << nl
                 << "    floatTransfer     : " << Pstream::floatTransfer << nl
-                << "    nProcsSimpleSum   : " << Pstream::nProcsSimpleSum << nl
+                << "    nProcsSimpleSum   : " << Pstream::nProcsSimpleSum() << nl
                 << "    commsType         : "
                 << Pstream::commsTypeNames[Pstream::defaultCommsType]
                 << endl;
@@ -568,6 +612,13 @@ Foam::argList::argList
         Info<< endl;
         IOobject::writeDivider(Info);
     }
+
+    // If the macro AppSpecificDictionary is used, one can
+    // modify the application-specific dictionnary using the
+    // command-line parameter -appDict
+    if(appDictName_ != "")
+      optionReadIfPresent("appDict", appDictName_);
+
 }
 
 
@@ -608,22 +659,31 @@ void Foam::argList::printUsage() const
         Info<< " <" << iter().c_str() << '>';
     }
 
+    int i=0;
+    SortableList<Foam::string> sortedValidOptions(validOptions.size());
+
     for
     (
         HashTable<string>::iterator iter = validOptions.begin();
         iter != validOptions.end();
-        ++iter
+        ++iter, ++i
     )
     {
-        Info<< " [-" << iter.key();
+	OStringStream keyValuePair;
+        keyValuePair << "[-" << iter.key();
 
         if (iter().size())
         {
-            Info<< ' ' << iter().c_str();
+            keyValuePair<< ' ' << iter().c_str();
         }
 
-        Info<< ']';
+        keyValuePair<< ']';
+	sortedValidOptions[i]= keyValuePair.str();
     }
+    sortedValidOptions.sort();
+
+    forAll(sortedValidOptions, sI)
+	Info << " " << sortedValidOptions[sI].c_str();
 
     // place help/doc/srcDoc options of the way at the end,
     // but with an extra space to separate it a little
diff --git a/src/foam/global/argList/argList.H b/src/foam/global/argList/argList.H
index 3a6855a37..edd32ebd8 100644
--- a/src/foam/global/argList/argList.H
+++ b/src/foam/global/argList/argList.H
@@ -126,6 +126,9 @@ public:
 
     // Static data members
 
+        //- Central dictionary switchSet names
+        static const NamedEnum<Foam::debug::globalControlDictSwitchSet, 5> globalControlDictSwitchSetNames_;
+
         //- A list of valid (mandatory) arguments
         static SLList<string> validArgs;
 
@@ -135,6 +138,9 @@ public:
         //- A list of valid parallel options
         static HashTable<string> validParOptions;
 
+        //- Name of the application-specific dictionary
+        static word appDictName_;
+
         //! @cond ignoreDocumentation
         class initValidTables
         {
@@ -263,6 +269,11 @@ public:
                 return readList<T>(optionLookup(opt)());
             }
 
+            //- Return appDictName_
+            const word appDictName() const
+            {
+                return appDictName_;
+            }
 
         // Edit
 
@@ -290,8 +301,21 @@ public:
 
             //- Check root path and case path
             bool checkRootCase() const;
+
 };
 
+// Allow overriding the application-specific dictionary using the command-line
+// parameter -appDict. The dictionary file must still be located in the case's
+// "system" directory, or accessible using a path relative to the case's
+// "system" directory.
+// Example:
+//    cellSet                                     # Default will be cellSetDict
+//    cellSet -appDict cellSetDict                # Same as above
+//    cellSet -appDict cellSetDict.rotorOnly
+//
+#define  AppSpecificDictionary(DefaultAppDict)                            \
+    argList::appDictName_ = (DefaultAppDict);		           \
+    argList::validOptions.insert("appDict", argList::appDictName_);
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
diff --git a/src/foam/global/controlSwitches/controlSwitches.C b/src/foam/global/controlSwitches/controlSwitches.C
new file mode 100644
index 000000000..f84d845f7
--- /dev/null
+++ b/src/foam/global/controlSwitches/controlSwitches.C
@@ -0,0 +1,145 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | foam-extend: Open Source CFD
+   \\    /   O peration     |
+    \\  /    A nd           | For copyright notice see file Copyright
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of foam-extend.
+
+    foam-extend is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation, either version 3 of the License, or (at your
+    option) any later version.
+
+    foam-extend is distributed in the hope that it will be useful, but
+    WITHOUT ANY WARRANTY; without even the implied warranty of
+    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+    General Public License for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
+
+Author
+    Martin Beaudoin, Hydro-Quebec, 2014.  All rights reserved
+
+\*---------------------------------------------------------------------------*/
+
+#include "controlSwitches.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+namespace debug
+{
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+template<class Type>
+Foam::debug::controlSwitches<Type>::controlSwitches()
+:
+    switchValue_(Type(0))
+{
+}
+
+template<class T>
+Foam::debug::controlSwitches<T>::controlSwitches(const T& switchValue)
+:
+    switchValue_(switchValue)
+{}
+
+
+template<class T>
+Foam::debug::controlSwitches<T>::controlSwitches
+(
+    const std::string& switchName,
+    const T& switchValue,
+    globalControlDictSwitchSet switchSet,
+    std::map<std::string, std::list<controlSwitches<T> *> >** switchValuesTable
+)
+:
+    switchSet_(switchSet),
+    switchName_(switchName),
+    switchValue_(switchValue)
+{
+    // Register the switch in its list
+    if (*switchValuesTable == NULL)
+    {
+	*switchValuesTable = new std::map<std::string, std::list<controlSwitches<T> *> >();
+    }
+
+    switchValuesTable_ = *switchValuesTable;
+    std::map<std::string, std::list<controlSwitches<T> *> >&switchValues = *switchValuesTable_;
+
+    // Memorize this switch object address
+    if (switchValues.find(switchName) != switchValues.end())
+    {
+	switchValues[switchName].push_back(this);
+    }
+    else
+    {
+	std::list<controlSwitches<T>* > pList;
+	pList.push_back(this);
+	switchValues.insert(std::pair<std::string, std::list<controlSwitches<T>* > >(switchName, pList));
+    }
+}
+
+
+template<class T>
+Foam::debug::controlSwitches<T>::controlSwitches(const Foam::debug::controlSwitches<T>& csw)
+:
+    switchName_(csw.switchName_),
+    switchValue_(csw.switchValue_)
+{}
+
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+template<class T>
+Foam::debug::controlSwitches<T>::~controlSwitches()
+{
+    // Unregister the switch from its list
+    std::map<std::string, std::list<controlSwitches<T> *> >&switchValuesTable = *switchValuesTable_;
+
+    // Remove entry or key if pointers list is empty
+    switchValuesTable[switchName_].remove(this);
+
+    // Replace the updated list
+    if(switchValuesTable[switchName_].size() == 0)
+    {
+	switchValuesTable.erase(switchName_);
+    }
+}
+
+
+// * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * * //
+
+template<class T>
+void Foam::debug::controlSwitches<T>::operator=(const Foam::debug::controlSwitches<T>& rhs)
+{
+    // Check for assignment to self
+    if (this == &rhs)
+    {
+	std::cerr << "Foam::debug::controlSwitches<T>::operator=(const Foam::controlSwitches<T>&)"
+		  << "--> FOAM FATAL ERROR: "
+		  << "Attempted assignment to self"
+		  << std::endl;
+	std::abort();
+	exit(-1);
+    }
+    else
+    {
+	switchValue_ = rhs.switchValue_;
+    }
+}
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace debug
+
+} // End namespace Foam
+
+// ************************************************************************* //
diff --git a/src/foam/global/controlSwitches/controlSwitches.H b/src/foam/global/controlSwitches/controlSwitches.H
new file mode 100644
index 000000000..0da85a3f4
--- /dev/null
+++ b/src/foam/global/controlSwitches/controlSwitches.H
@@ -0,0 +1,358 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | foam-extend: Open Source CFD
+   \\    /   O peration     |
+    \\  /    A nd           | For copyright notice see file Copyright
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of foam-extend.
+
+    foam-extend is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation, either version 3 of the License, or (at your
+    option) any later version.
+
+    foam-extend is distributed in the hope that it will be useful, but
+    WITHOUT ANY WARRANTY; without even the implied warranty of
+    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+    General Public License for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::controlSwitches
+
+Author
+    Martin Beaudoin, Hydro-Quebec, 2014.  All rights reserved
+
+Description
+    Handling of runtime controlSwitches
+    - Debug switches
+    - Info switches
+    - Optimisation switches
+    - Tolerances switches
+    - DimensionedConstants values
+
+    This class is a specialization of the safe_bool class, where additional
+    boolean operators are allowed. This class is also offering access to
+    information like the switch name, etc.
+
+Note:
+    We will be using std::map and std::list for the containers for this class.
+    In some sections of the code where controlSwitches are being used, we are at
+    a relatively low level, so we cannot use HashTable and DynamicList.
+    MB 23/05/2014
+
+SourceFiles
+    controlSwitchesI.H
+    controlSwitches.C
+    controlSwitchesIO.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef controlSwitches_H
+#define controlSwitches_H
+
+#include <list>
+#include <map>
+#include <string>
+#include "safe_bool.H"
+
+#include "scalar.H"
+#include "debug.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+namespace debug
+{
+
+// Forward declaration of classes
+template<class T> class controlSwitches;
+
+// Typedefs
+// Hashtables to static class attributes addresses holding the
+// runtime debug/info/optimisation/tolerances values
+typedef std::map<std::string, std::list<controlSwitches<int> *> >         ListDebugControlSwitches;
+typedef std::map<std::string, std::list<controlSwitches<int> *> >         ListInfoControlSwitches;
+typedef std::map<std::string, std::list<controlSwitches<int> *> >         ListOptimisationControlSwitches;
+typedef std::map<std::string, std::list<controlSwitches<Foam::scalar> *> > ListTolerancesControlSwitches;
+typedef std::map<std::string, std::list<controlSwitches<Foam::scalar> *> > ListConstantsControlSwitches;
+
+extern ListDebugControlSwitches* debugSwitchValues_;
+extern ListInfoControlSwitches* infoSwitchValues_;
+extern ListOptimisationControlSwitches* optimisationSwitchValues_;
+extern ListTolerancesControlSwitches* tolerancesSwitchValues_;
+extern ListConstantsControlSwitches* constantsSwitchValues_;
+
+/*---------------------------------------------------------------------------*\
+                          Class controlSwitches Declaration
+\*---------------------------------------------------------------------------*/
+
+template<class T>
+class controlSwitches
+    : public safe_bool<controlSwitches<T> > // using the CRTP idiom
+{
+    // Private data
+
+        //- Switch section
+        globalControlDictSwitchSet switchSet_;
+
+        //- Switch name
+        std::string switchName_;
+
+        //- Switch value
+        T switchValue_;
+
+        //- Handle to runTime switches list
+        std::map<std::string, std::list<controlSwitches<T> *> >* switchValuesTable_ ;
+
+    // Private Member Functions
+
+public:
+
+    // Constructors
+
+        //- Construct null
+        controlSwitches();
+
+        //- Construct from components
+        controlSwitches(const T& data);
+
+        controlSwitches
+	(
+	    const std::string& switchName,
+	    const T& switchValue,
+	    globalControlDictSwitchSet switchSet,
+	    std::map<std::string, std::list<controlSwitches<T> *> >** switchesValues
+	);
+
+        //- Construct as copy
+        controlSwitches(const controlSwitches&);
+
+    //- Destructor
+        ~controlSwitches();
+
+
+    // Member Functions
+        // Access
+
+        // Check
+        bool boolean_test() const
+	{
+	    // Perform Boolean logic here
+	    return switchValue_ != 0;
+	}
+
+    // Member Operators
+
+        // Assignement operator
+        void operator=(const T&);
+
+        //- Assignement operator
+        void operator=(const controlSwitches&);
+
+        // & operator
+        const T operator&(const T&);
+
+        // == operator
+        bool operator==(const T&);
+
+        // >= operator
+        bool operator>=(const T&);
+
+        // > operator
+        bool operator>(const T&);
+
+        // < operator
+        bool operator<(const T&);
+
+        // Return value
+        T operator()() const
+	{
+	    return switchValue_;
+	}
+};
+
+// Utility classes
+class debugSwitch
+    : public controlSwitches<int>
+{
+public:
+    debugSwitch
+    (
+	const std::string& switchName,
+	const int& switchValue
+    )
+	:
+	controlSwitches<int>
+	(
+	    switchName,
+	    debug::debugSwitchFromDict(switchName.c_str(), switchValue),
+	    debug::DEBUGSWITCHES,
+	    &debugSwitchValues_
+	)
+    {}
+
+    virtual ~debugSwitch()
+    {
+    }
+
+    // Assignement operator
+    void operator=(const int& rhs)
+    {
+	controlSwitches<int>::operator=(rhs);
+    }
+};
+
+class infoSwitch
+    : public controlSwitches<int>
+{
+public:
+    infoSwitch
+    (
+	const std::string& switchName,
+	const int& switchValue
+    )
+	:
+	controlSwitches<int>
+	(
+	    switchName,
+	    debug::infoSwitchFromDict(switchName.c_str(), switchValue),
+	    debug::INFOSWITCHES,
+	    &infoSwitchValues_
+	)
+    {}
+
+    virtual ~infoSwitch()
+    {
+    }
+
+    // Assignement operator
+    void operator=(const int& rhs)
+    {
+	controlSwitches<int>::operator=(rhs);
+    }
+};
+
+class optimisationSwitch
+    : public controlSwitches<int>
+{
+public:
+    optimisationSwitch
+    (
+	const std::string& switchName,
+	const int& switchValue
+    )
+	:
+	controlSwitches<int>
+	(
+	    switchName,
+	    debug::optimisationSwitchFromDict(switchName.c_str(), switchValue),
+	    debug::OPTIMISATIONSWITCHES,
+	    &optimisationSwitchValues_
+	)
+    {}
+
+    virtual ~optimisationSwitch()
+    {
+    }
+
+    // Assignement operator
+    void operator=(const int& rhs)
+    {
+	controlSwitches<int>::operator=(rhs);
+    }
+
+
+};
+
+class tolerancesSwitch
+    : public controlSwitches<Foam::scalar>
+{
+public:
+    tolerancesSwitch
+    (
+	const std::string& switchName,
+	const Foam::scalar& switchValue
+    )
+	:
+	controlSwitches<Foam::scalar>
+	(
+	    switchName,
+	    debug::tolerancesFromDict(switchName.c_str(), switchValue),
+	    debug::TOLERANCES,
+	    &tolerancesSwitchValues_
+	)
+    {
+    }
+
+    virtual ~tolerancesSwitch()
+    {
+    }
+
+   // Assignement operator
+    void operator=(const Foam::scalar& rhs)
+    {
+	controlSwitches<Foam::scalar>::operator=(rhs);
+    }
+};
+
+class constantsSwitch
+    : public controlSwitches<Foam::scalar>
+{
+public:
+    constantsSwitch
+    (
+	const std::string& switchName,
+	const Foam::scalar& switchValue
+    )
+	:
+	controlSwitches<Foam::scalar>
+	(
+	    switchName,
+	    debug::constantsFromDict(switchName.c_str(), switchValue),
+	    debug::DIMENSIONEDCONSTANTS,
+	    &constantsSwitchValues_
+	)
+    {
+    }
+
+    virtual ~constantsSwitch()
+    {
+    }
+
+   // Assignement operator
+    void operator=(const Foam::scalar& rhs)
+    {
+	controlSwitches<Foam::scalar>::operator=(rhs);
+    }
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace debug
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#include "controlSwitchesI.H"
+
+#ifdef NoRepository
+#   include "controlSwitches.C"
+#endif
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/foam/global/controlSwitches/controlSwitchesI.H b/src/foam/global/controlSwitches/controlSwitchesI.H
new file mode 100644
index 000000000..751f362cb
--- /dev/null
+++ b/src/foam/global/controlSwitches/controlSwitchesI.H
@@ -0,0 +1,79 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | foam-extend: Open Source CFD
+   \\    /   O peration     |
+    \\  /    A nd           | For copyright notice see file Copyright
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of foam-extend.
+
+    foam-extend is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation, either version 3 of the License, or (at your
+    option) any later version.
+
+    foam-extend is distributed in the hope that it will be useful, but
+    WITHOUT ANY WARRANTY; without even the implied warranty of
+    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+    General Public License for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+
+// * * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * //
+
+template<class T>
+inline void Foam::debug::controlSwitches<T>::operator=(const T& rhs)
+{
+    switchValue_ = rhs;
+}
+
+template<class T>
+inline const T Foam::debug::controlSwitches<T>::operator&(const T& rhs)
+{
+    return switchValue_ & rhs;
+}
+
+template<class T>
+inline bool Foam::debug::controlSwitches<T>::operator==(const T& rhs)
+{
+    return switchValue_ == rhs;
+}
+
+template<class T>
+inline bool Foam::debug::controlSwitches<T>::operator>=(const T& rhs)
+{
+    return switchValue_ >= rhs;
+}
+
+template<class T>
+inline bool Foam::debug::controlSwitches<T>::operator>(const T& rhs)
+{
+    return switchValue_ > rhs;
+}
+
+template<class T>
+inline bool Foam::debug::controlSwitches<T>::operator<(const T& rhs)
+{
+    return switchValue_ < rhs;
+}
+
+
+// * * * * * * * * * * * * * * IOstream Operators  * * * * * * * * * * * * * //
+
+
+// * * * * * * * * * * * * * * * Ostream Operator  * * * * * * * * * * * * * //
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+
+// ************************************************************************* //
diff --git a/src/foam/global/controlSwitches/safe_bool.H b/src/foam/global/controlSwitches/safe_bool.H
new file mode 100644
index 000000000..498335076
--- /dev/null
+++ b/src/foam/global/controlSwitches/safe_bool.H
@@ -0,0 +1,100 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | foam-extend: Open Source CFD
+   \\    /   O peration     |
+    \\  /    A nd           | For copyright notice see file Copyright
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of foam-extend.
+
+    foam-extend is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation, either version 3 of the License, or (at your
+    option) any later version.
+
+    foam-extend is distributed in the hope that it will be useful, but
+    WITHOUT ANY WARRANTY; without even the implied warranty of
+    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+    General Public License for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::safe_bool
+
+Description
+    Implementation of the Safe Bool Idiom.
+    This is the version without virtual functions.
+
+    Author: Bjorn Karlsson
+    http://www.artima.com/cppsource/safebool.html
+
+SourceFiles
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef safe_bool_H
+#define safe_bool_H
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+/*---------------------------------------------------------------------------*\
+                          Class safe_bool_base Declaration
+\*---------------------------------------------------------------------------*/
+class safe_bool_base
+{
+    public:
+        typedef void (safe_bool_base::*bool_type)() const;
+        void this_type_does_not_support_comparisons() const {}
+
+    protected:
+
+        safe_bool_base() {}
+        safe_bool_base(const safe_bool_base&) {}
+        safe_bool_base& operator=(const safe_bool_base&) {return *this;}
+        ~safe_bool_base() {}
+};
+
+/*---------------------------------------------------------------------------*\
+                          Class safe_bool Declaration
+\*---------------------------------------------------------------------------*/
+
+template <typename T=void>
+class safe_bool
+:
+    private safe_bool_base
+{
+  // private or protected inheritance is very important here as it triggers the
+  // access control violation in main.
+
+public:
+    operator bool_type() const
+    {
+	return (static_cast<const T*>(this))->boolean_test()
+	    ? &safe_bool_base::this_type_does_not_support_comparisons : 0;
+    }
+
+protected:
+    //- Destructor
+    ~safe_bool()
+    {}
+};
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#include "safe_boolI.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/foam/global/controlSwitches/safe_boolI.H b/src/foam/global/controlSwitches/safe_boolI.H
new file mode 100644
index 000000000..455a9114b
--- /dev/null
+++ b/src/foam/global/controlSwitches/safe_boolI.H
@@ -0,0 +1,79 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | foam-extend: Open Source CFD
+   \\    /   O peration     |
+    \\  /    A nd           | For copyright notice see file Copyright
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of foam-extend.
+
+    foam-extend is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by the
+    Free Software Foundation, either version 3 of the License, or (at your
+    option) any later version.
+
+    foam-extend is distributed in the hope that it will be useful, but
+    WITHOUT ANY WARRANTY; without even the implied warranty of
+    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+    General Public License for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+// * * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * //
+
+template <typename T>
+inline bool operator==(const Foam::safe_bool<T>& lhs, bool b)
+{
+    if (b)
+    {
+	if (lhs) return true;
+	else return false;
+    }
+    else
+    {
+	if (lhs) return false;
+	else return true;
+    }
+}
+
+template <typename T>
+inline bool operator==(bool b, const Foam::safe_bool<T>& rhs)
+{
+    if (b)
+    {
+	if (rhs) return true;
+	else return false;
+    }
+    else
+    {
+	if (rhs) return false;
+	else return true;
+    }
+}
+
+
+template <typename T, typename U>
+inline bool operator==(const Foam::safe_bool<T>& lhs,const Foam::safe_bool<U>& rhs)
+{
+    lhs.this_type_does_not_support_comparisons();
+    return false;
+}
+
+template <typename T,typename U>
+inline  bool operator!=(const Foam::safe_bool<T>& lhs,const Foam::safe_bool<U>& rhs)
+{
+    lhs.this_type_does_not_support_comparisons();
+    return false;
+}
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+
+// ************************************************************************* //
diff --git a/src/foam/global/debug/debug.C b/src/foam/global/debug/debug.C
index 230e9bbe4..b1548f2ad 100644
--- a/src/foam/global/debug/debug.C
+++ b/src/foam/global/debug/debug.C
@@ -30,6 +30,9 @@ Description
 #include "dictionary.H"
 #include "IFstream.H"
 #include "OSspecific.H"
+#include "dimensionedConstants.H"
+#include "SortableList.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -44,6 +47,18 @@ dictionary* debugSwitchesPtr_(NULL);
 dictionary* infoSwitchesPtr_(NULL);
 dictionary* optimisationSwitchesPtr_(NULL);
 dictionary* tolerancesPtr_(NULL);
+dictionary* constantsPtr_(NULL);
+
+// Hashtables to static class attributes addresses holding the
+// runtime debug/info/optimisation/tolerances values
+// This needs to go on the heap so the destructor will not get
+//  called before the object's destructor it is overseeing
+ListDebugControlSwitches* debugSwitchValues_(NULL);
+ListInfoControlSwitches* infoSwitchValues_(NULL);
+ListOptimisationControlSwitches* optimisationSwitchValues_(NULL);
+ListTolerancesControlSwitches* tolerancesSwitchValues_(NULL);
+ListConstantsControlSwitches* constantsSwitchValues_(NULL);
+
 
 // to ensure controlDictPtr_ is deleted at the end of the run
 class deleteControlDictPtr
@@ -149,8 +164,13 @@ Foam::dictionary& Foam::debug::tolerances()
     return switchSet("Tolerances", tolerancesPtr_);
 }
 
+Foam::dictionary& Foam::debug::constants()
+{
+    return switchSet("DimensionedConstants", constantsPtr_);
+}
 
-int Foam::debug::debugSwitch(const char* name, const int defaultValue)
+
+int Foam::debug::debugSwitchFromDict(const char* name, const int defaultValue)
 {
     return debugSwitches().lookupOrAddDefault
     (
@@ -159,7 +179,7 @@ int Foam::debug::debugSwitch(const char* name, const int defaultValue)
 }
 
 
-int Foam::debug::infoSwitch(const char* name, const int defaultValue)
+int Foam::debug::infoSwitchFromDict(const char* name, const int defaultValue)
 {
     return infoSwitches().lookupOrAddDefault
     (
@@ -168,7 +188,7 @@ int Foam::debug::infoSwitch(const char* name, const int defaultValue)
 }
 
 
-int Foam::debug::optimisationSwitch(const char* name, const int defaultValue)
+int Foam::debug::optimisationSwitchFromDict(const char* name, const int defaultValue)
 {
     return optimisationSwitches().lookupOrAddDefault
     (
@@ -177,11 +197,7 @@ int Foam::debug::optimisationSwitch(const char* name, const int defaultValue)
 }
 
 
-double Foam::debug::tolerances
-(
-    const char* name,
-    const double defaultValue
-)
+double Foam::debug::tolerancesFromDict(const char* name, const double defaultValue)
 {
     return tolerances().lookupOrAddDefault
     (
@@ -189,5 +205,299 @@ double Foam::debug::tolerances
     );
 }
 
+double Foam::debug::constantsFromDict(const char* name, const double defaultValue)
+{
+    return tolerances().lookupOrAddDefault
+    (
+        name, defaultValue, false, false
+    );
+}
+
+void Foam::debug::updateCentralDictVars(Foam::debug::globalControlDictSwitchSet globalControlDictSwitchSetName, Foam::string keyValues)
+{
+    Foam::word token;
+
+    // Extract the multiple keyValues, separated by ','
+    std::istringstream issmkvs(keyValues);
+    while ( getline(issmkvs, token, ',') )
+    {
+	// Cleanup
+	Foam::string::stripInvalid<word>(token);
+
+	// Extract the keyValue pair, separated by '='
+	std::istringstream isskv(token);
+
+	while ( getline(isskv, token, '=') )
+	{
+	    Foam::word key = token;
+
+	    switch(globalControlDictSwitchSetName)
+	    {
+	        case DEBUGSWITCHES:
+		{
+		    if (debugSwitchValues_)
+		    {
+			ListDebugControlSwitches& runTimeDebugSwitchValues = *debugSwitchValues_;
+			label oldDebugValue;
+			label newDebugValue;
+			isskv >> newDebugValue;
+
+			if(Foam::debug::debugSwitches().readIfPresent(key, oldDebugValue))
+			{
+			    Info << endl
+				 << "Warning: Modification of DebugSwitch: " << key << endl
+				 << "    Old value: " << oldDebugValue << endl
+				 << "    New value: " << newDebugValue << endl
+				 << endl;
+
+			    Foam::debug::debugSwitches().set(key, newDebugValue);
+
+			    std::list<controlSwitches<int> *> curList = runTimeDebugSwitchValues[key];
+
+			    // Modify all entries for this key
+			    forAllIter(std::list<controlSwitches<int> *>, curList, iterI)
+				*(*iterI) =  newDebugValue;
+			}
+			else
+			{
+			    //  Usage of non-existent DebugSwitches: best to abort right away
+			    SortableList<Foam::word> sortedValidKeys(runTimeDebugSwitchValues.size());
+
+			    int i=0;
+			    forAllIter(ListDebugControlSwitches, runTimeDebugSwitchValues, iterI)
+			    {
+				sortedValidKeys[i++] = iterI->first;
+			    }
+			    sortedValidKeys.sort();
+
+			    FatalError
+				<< "Usage of non-existent DebugSwitches: " << key
+				<< endl << endl
+				<< "Valid entries for this application are: " << sortedValidKeys
+				<< exit(FatalError);
+			}
+		    }
+		    else
+		    {
+			FatalError
+			    << "No DebugSwitches values are available for this application."
+			    << exit(FatalError);
+		    }
+		}
+		break;
+	        case INFOSWITCHES:
+		{
+		    if (infoSwitchValues_)
+		    {
+			ListInfoControlSwitches& runTimeInfoSwitchValues = *infoSwitchValues_;
+			label oldInfoValue;
+			label newInfoValue;
+			isskv >> newInfoValue;
+
+			if(Foam::debug::infoSwitches().readIfPresent(key, oldInfoValue))
+			{
+			    Info << endl
+				 << "Warning: Modification of InfoSwitch: " << key << endl
+				 << "    Old value: " << oldInfoValue << endl
+				 << "    New value: " << newInfoValue << endl
+				 << endl;
+
+			    Foam::debug::infoSwitches().set(key, newInfoValue);
+
+			    std::list<controlSwitches<int> *> curList = runTimeInfoSwitchValues[key];
+
+			    // Modify all entries for this key
+			    forAllIter(std::list<controlSwitches<int> *>, curList, iterI)
+				*(*iterI) =  newInfoValue;
+			}
+			else
+			{
+			    // Usage of non-existent InfoSwitches: best to abort right away
+			    SortableList<Foam::word> sortedValidKeys(runTimeInfoSwitchValues.size());
+
+			    int i=0;
+			    forAllIter(ListInfoControlSwitches, runTimeInfoSwitchValues, iterI)
+			    {
+				sortedValidKeys[i++] = iterI->first;
+			    }
+			    sortedValidKeys.sort();
+
+			    FatalError
+				<< "Usage of non-existent InfoSwitches: " << key
+				<< endl << endl
+				<< "Valid entries for this application are: " << sortedValidKeys
+				<< exit(FatalError);
+			}
+		    }
+		    else
+		    {
+			FatalError
+			    << "No InfoSwitches values are available for this application."
+			    << exit(FatalError);
+		    }
+		}
+		break;
+	        case OPTIMISATIONSWITCHES:
+		{
+		    if (optimisationSwitchValues_)
+		    {
+			ListOptimisationControlSwitches& runTimeOptimisationSwitchValues = *optimisationSwitchValues_;
+			label oldOptimisationValue;
+			label newOptimisationValue;
+			isskv >> newOptimisationValue;
+
+			if(Foam::debug::optimisationSwitches().readIfPresent(key, oldOptimisationValue))
+			{
+			    Info << endl
+				 << "Warning: Modification of OptimisationSwitch value: " << key << endl
+				 << "    Old value: " << oldOptimisationValue << endl
+				 << "    New value: " << newOptimisationValue << endl
+				 << endl;
+
+			    Foam::debug::optimisationSwitches().set(key, newOptimisationValue);
+
+			    std::list<controlSwitches<int> *> curList = runTimeOptimisationSwitchValues[key];
+
+			    // Modify all entries for this key
+			    forAllIter(std::list<controlSwitches<int> *>, curList, iterI)
+				*(*iterI) =  newOptimisationValue;
+			}
+			else
+		        {
+			    // Usage of non-existent OptimisationSwitches: best to abort right away
+			    SortableList<Foam::word> sortedValidKeys(runTimeOptimisationSwitchValues.size());
+
+			    int i=0;
+			    forAllIter(ListOptimisationControlSwitches, runTimeOptimisationSwitchValues, iterI)
+			    {
+				sortedValidKeys[i++] = iterI->first;
+			    }
+			    sortedValidKeys.sort();
+
+			    FatalError
+				<< "Usage of non-existent OptimisationSwitches value: " << key
+				<< endl << endl
+				<< "Valid entries for this application are: " << sortedValidKeys
+				<< exit(FatalError);
+			}
+		    }
+		    else
+		    {
+			FatalError
+			    << "No OptimisationSwitches values are available for this application."
+			    << exit(FatalError);
+		    }
+		}
+		break;
+	        case TOLERANCES:
+		{
+		    if (tolerancesSwitchValues_)
+		    {
+			ListTolerancesControlSwitches& runTimeTolerancesSwitchValues = *tolerancesSwitchValues_;
+			scalar oldTolerancesValue;
+			scalar newTolerancesValue;
+			isskv >> newTolerancesValue;
+
+			if(Foam::debug::tolerances().readIfPresent(key, oldTolerancesValue))
+			{
+			    Info << endl
+				 << "Warning: Modification of Tolerance value: " << key << endl
+				 << "    Old value: " << oldTolerancesValue << endl
+				 << "    New value: " << newTolerancesValue << endl
+				 << endl;
+
+			    Foam::debug::tolerances().set(key, newTolerancesValue);
+
+			    std::list<controlSwitches<scalar> *> curList = runTimeTolerancesSwitchValues[key];
+
+			    // Modify all entries for this key
+			    forAllIter(std::list<controlSwitches<scalar> *>, curList, iterI)
+				*(*iterI) =  newTolerancesValue;
+			}
+			else
+			{
+			    // Usage of non-existent Tolerances: best to abort right away
+			    SortableList<Foam::word> sortedValidKeys(runTimeTolerancesSwitchValues.size());
+
+			    int i=0;
+			    forAllIter(ListTolerancesControlSwitches, runTimeTolerancesSwitchValues, iterI)
+				{
+				    sortedValidKeys[i++] = iterI->first;
+				}
+			    sortedValidKeys.sort();
+
+			    FatalError
+				<< "Usage of non-existent Tolerances value: " << key
+				<< endl << endl
+				<< "Valid entries for this application are: " << sortedValidKeys
+				<< exit(FatalError);
+			}
+		    }
+		    else
+		    {
+			FatalError
+			    << "No Tolerances values are available for this application."
+			    << exit(FatalError);
+		    }
+		}
+		break;
+	        case DIMENSIONEDCONSTANTS:
+		{
+		    if (constantsSwitchValues_)
+		    {
+			ListConstantsControlSwitches& runTimeConstantsSwitchValues = *constantsSwitchValues_;
+			scalar oldDimensionedConstantsValue;
+			scalar newDimensionedConstantsValue;
+			isskv >> newDimensionedConstantsValue;
+
+			if(Foam::dimensionedConstants().readIfPresent(key, oldDimensionedConstantsValue))
+			{
+			    Info << endl
+				 << "Warning: Modification of DimensionedConstant: " << key << endl
+				 << "    Old value: " << oldDimensionedConstantsValue << endl
+				 << "    New value: " << newDimensionedConstantsValue << endl
+				 << endl;
+
+			    Foam::dimensionedConstants().set(key, newDimensionedConstantsValue);
+
+			    std::list<controlSwitches<scalar> *> curList = runTimeConstantsSwitchValues[key];
+
+			    // Modify all entries for this key
+			    forAllIter(std::list<controlSwitches<scalar> *>, curList, iterI)
+				*(*iterI) =  newDimensionedConstantsValue;
+			}
+			else
+			{
+			    // Usage of non-existent DimensionedConstants: best to abort right away
+			    SortableList<Foam::word> sortedValidKeys(runTimeConstantsSwitchValues.size());
+
+			    int i=0;
+			    forAllIter(ListConstantsControlSwitches, runTimeConstantsSwitchValues, iterI)
+			    {
+				sortedValidKeys[i++] = iterI->first;
+			    }
+			    sortedValidKeys.sort();
+
+			    FatalError
+				<< "Usage of non-existent DimensionedConstants: " << key
+				<< endl << endl
+				<< "Valid entries for this application are: " << sortedValidKeys
+				<< exit(FatalError);
+			}
+		    }
+		    else
+		    {
+			FatalError
+			    << "No DimensionedConstants values are available for this application."
+			    << exit(FatalError);
+		    }
+		}
+		break;
+	        default:
+		break;
+	    }
+	}
+    }
+}
 
 // ************************************************************************* //
diff --git a/src/foam/global/debug/debug.H b/src/foam/global/debug/debug.H
index 02d1b23ad..f52d0c3cf 100644
--- a/src/foam/global/debug/debug.H
+++ b/src/foam/global/debug/debug.H
@@ -47,6 +47,16 @@ class dictionary;
 
 namespace debug
 {
+    //- Define central control dictionary switchSet
+    enum globalControlDictSwitchSet
+    {
+	DEBUGSWITCHES,
+	INFOSWITCHES,
+	OPTIMISATIONSWITCHES,
+	TOLERANCES,
+	DIMENSIONEDCONSTANTS
+    };
+
     //- The central control dictionary.
     //  Located in $WM_PROJECT_DIR/etc
     //  @sa Foam::findEtcFile()
@@ -64,21 +74,35 @@ namespace debug
     //- The Tolerances sub-dictionary in the central controlDict.
     dictionary& tolerances();
 
+    //- The DimensionedConstants sub-dictionary in the central controlDict.
+    dictionary& constants();
+
     //- Lookup debug switch or add default value.
-    int debugSwitch(const char* name, const int defaultValue = 0);
+    int debugSwitchFromDict(const char* name, const int defaultValue = 0);
 
     //- Lookup info switch or add default value.
-    int infoSwitch(const char* name, const int defaultValue = 0);
+    int infoSwitchFromDict(const char* name, const int defaultValue = 0);
 
     //- Lookup optimisation switch or add default value.
-    int optimisationSwitch(const char* name, const int defaultValue = 0);
+    int optimisationSwitchFromDict(const char* name, const int defaultValue = 0);
 
     //- Lookup tolerances switch or add default value.
-    double tolerances(const char* name, const double defaultValue = 0);
+    double tolerancesFromDict(const char* name, const double defaultValue = 0);
+
+    //- Lookup tolerances switch or add default value.
+    double constantsFromDict(const char* name, const double defaultValue = 0);
 
     //- Internal function to lookup a sub-dictionary from controlDict.
     dictionary& switchSet(const char* subDictName, dictionary*& subDictPtr);
 
+    //- Update central controlDict variables
+    //  supported keyValues format: key1=value1,key2=value2
+    void updateCentralDictVars
+    (
+	globalControlDictSwitchSet centralDictSwitchSetName,
+	Foam::string keyValues
+    );
+
 } // End namespace debug
 
 
diff --git a/src/foam/interpolations/GGIInterpolation/GGIInterpolation.H b/src/foam/interpolations/GGIInterpolation/GGIInterpolation.H
index 8ac0e1558..9e5d00d54 100644
--- a/src/foam/interpolations/GGIInterpolation/GGIInterpolation.H
+++ b/src/foam/interpolations/GGIInterpolation/GGIInterpolation.H
@@ -53,6 +53,7 @@ SourceFiles
 #include "intersection.H"
 #include "point2D.H"
 #include "NamedEnum.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -188,13 +189,13 @@ class GGIInterpolation
     // Private static data
 
         //- Facet area error tolerance
-        static const scalar areaErrorTol_;
+        static const debug::tolerancesSwitch areaErrorTol_;
 
         //- Facet normal featureCos criteria
-        static const scalar featureCosTol_;
+        static const debug::tolerancesSwitch featureCosTol_;
 
-       //- Facet bound box extension factor
-       static const scalar faceBoundBoxExtendSpanFraction_;
+        //- Facet bound box extension factor
+        static const debug::tolerancesSwitch faceBoundBoxExtendSpanFraction_;
 
        //- The next 3 attributes are parameters controlling the creation
        //  of an octree search engine for the GGI facets neighbourhood
@@ -218,13 +219,13 @@ class GGIInterpolation
        //
 
        //- Octree search: minNlevel parameter for octree constructor
-       static const label octreeSearchMinNLevel_;
+       static const Foam::debug::optimisationSwitch octreeSearchMinNLevel_;
 
        //- Octree search: maxLeafRatio parameter for octree constructor
-       static const scalar octreeSearchMaxLeafRatio_;
+       static const Foam::debug::optimisationSwitch octreeSearchMaxLeafRatio_;
 
        //- Octree search: maxShapeRatio parameter for octree constructor
-       static const scalar octreeSearchMaxShapeRatio_;
+       static const Foam::debug::optimisationSwitch octreeSearchMaxShapeRatio_;
 
 
     // Private Member Functions
diff --git a/src/foam/interpolations/GGIInterpolation/GGIInterpolationPolygonIntersection.C b/src/foam/interpolations/GGIInterpolation/GGIInterpolationPolygonIntersection.C
index dff08b5d0..0e9c9ee4c 100644
--- a/src/foam/interpolations/GGIInterpolation/GGIInterpolationPolygonIntersection.C
+++ b/src/foam/interpolations/GGIInterpolation/GGIInterpolationPolygonIntersection.C
@@ -184,8 +184,8 @@ GGIInterpolation<MasterPatch, SlavePatch>::polygonIntersection
 
         if
         (
-            mag(intersectionArea/clippingArea) < areaErrorTol_
-         || mag(intersectionArea/subjectArea)  < areaErrorTol_
+            mag(intersectionArea/clippingArea) < areaErrorTol_()
+	 || mag(intersectionArea/subjectArea)  < areaErrorTol_()
         )
         {
             WarningIn
@@ -232,7 +232,7 @@ GGIInterpolation<MasterPatch, SlavePatch>::isVertexInsidePolygon
     // We use distErrorTol_ to evaluate a distance factor called
     // epsilon.  That epsilon factor will be used to detect if a point
     // is on a vertex or an edge.
-    const scalar distErrorTol = sqrt(areaErrorTol_);
+    const scalar distErrorTol = sqrt(areaErrorTol_());
 
     HormannAgathos pip(clippingPolygon, distErrorTol);
 
@@ -303,7 +303,7 @@ GGIInterpolation<MasterPatch, SlavePatch>::clipPolygon2DSutherlandHodgman
     (
         clippingPolygon,
         subjectPolygon,
-        sqrt(areaErrorTol_) // = distErrorTol
+        sqrt(areaErrorTol_()) // = distErrorTol
     ).evaluate();
 }
 
diff --git a/src/foam/interpolations/GGIInterpolation/GGIInterpolationQuickRejectTests.C b/src/foam/interpolations/GGIInterpolation/GGIInterpolationQuickRejectTests.C
index e56b27874..5c5fab3a0 100644
--- a/src/foam/interpolations/GGIInterpolation/GGIInterpolationQuickRejectTests.C
+++ b/src/foam/interpolations/GGIInterpolation/GGIInterpolationQuickRejectTests.C
@@ -44,31 +44,35 @@ namespace Foam
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
 template<class MasterPatch, class SlavePatch>
-const scalar
+const Foam::debug::tolerancesSwitch
 GGIInterpolation<MasterPatch, SlavePatch>::faceBoundBoxExtendSpanFraction_
 (
-    debug::tolerances("GGIFaceBoundBoxExtendSpanFraction", 1.0e-2)
+    "GGIFaceBoundBoxExtendSpanFraction",
+    1.0e-2
 );
 
 template<class MasterPatch, class SlavePatch>
-const label
+const Foam::debug::optimisationSwitch
  GGIInterpolation<MasterPatch, SlavePatch>::octreeSearchMinNLevel_
 (
-    debug::optimisationSwitch("GGIOctreeSearchMinNLevel", 3)
+    "GGIOctreeSearchMinNLevel",
+    3
 );
 
 template<class MasterPatch, class SlavePatch>
-const scalar
+const Foam::debug::optimisationSwitch
 GGIInterpolation<MasterPatch, SlavePatch>::octreeSearchMaxLeafRatio_
 (
-    debug::optimisationSwitch("GGIOctreeSearchMaxLeafRatio", 3)
+    "GGIOctreeSearchMaxLeafRatio",
+    3
 );
 
 template<class MasterPatch, class SlavePatch>
-const scalar
+const Foam::debug::optimisationSwitch
 GGIInterpolation<MasterPatch, SlavePatch>::octreeSearchMaxShapeRatio_
 (
-    debug::optimisationSwitch("GGIOctreeSearchMaxShapeRatio", 1)
+    "GGIOctreeSearchMaxShapeRatio",
+    1
 );
 
 
@@ -397,7 +401,7 @@ void GGIInterpolation<MasterPatch, SlavePatch>::findNeighboursAABB
                 scalar featureCos =
                     masterFaceNormals[faceMi] & slaveNormals[faceSi];
 
-                if (mag(featureCos) > featureCosTol_)
+                if (mag(featureCos) > featureCosTol_())
                 {
                     candidateMasterNeighbors[faceMi].append(faceSi);
                 }
@@ -440,7 +444,7 @@ void GGIInterpolation<MasterPatch, SlavePatch>::findNeighboursBBOctree
         treeBoundBox bbFaceMaster(facePoints);
 
         lmasterFaceBB[faceMi] =
-            bbFaceMaster.extend(faceBoundBoxExtendSpanFraction_);
+            bbFaceMaster.extend(faceBoundBoxExtendSpanFraction_());
     }
 
     // Initialize the list of slave patch faces bounding box
@@ -484,7 +488,7 @@ void GGIInterpolation<MasterPatch, SlavePatch>::findNeighboursBBOctree
         treeBoundBox bbFaceSlave(facePoints);
 
         lslaveFaceBB[faceSi] =
-            bbFaceSlave.extend(faceBoundBoxExtendSpanFraction_);
+            bbFaceSlave.extend(faceBoundBoxExtendSpanFraction_());
     }
 
     // Create the slave octreeData, using the boundBox flavor
@@ -499,9 +503,9 @@ void GGIInterpolation<MasterPatch, SlavePatch>::findNeighboursBBOctree
     (
         slaveOverallBB,              // overall search domain
         slaveDataBB,
-        octreeSearchMinNLevel_,      // min number of levels
-        octreeSearchMaxLeafRatio_,   // max avg. size of leaves
-        octreeSearchMaxShapeRatio_   // max avg. duplicity.
+        octreeSearchMinNLevel_(),      // min number of levels
+        octreeSearchMaxLeafRatio_(),   // max avg. size of leaves
+        octreeSearchMaxShapeRatio_()   // max avg. duplicity.
     );
 
     const vectorField& masterFaceNormals = masterPatch_.faceNormals();
@@ -538,7 +542,7 @@ void GGIInterpolation<MasterPatch, SlavePatch>::findNeighboursBBOctree
             scalar featureCos =
                 masterFaceNormals[faceMi] & slaveNormals[faceSi];
 
-            if (mag(featureCos) > featureCosTol_)
+            if (mag(featureCos) > featureCosTol_())
             {
                 candidateMasterNeighbors[faceMi].append(faceSi);
             }
diff --git a/src/foam/interpolations/GGIInterpolation/GGIInterpolationWeights.C b/src/foam/interpolations/GGIInterpolation/GGIInterpolationWeights.C
index f0aa15fc8..a56a551fc 100644
--- a/src/foam/interpolations/GGIInterpolation/GGIInterpolationWeights.C
+++ b/src/foam/interpolations/GGIInterpolation/GGIInterpolationWeights.C
@@ -37,7 +37,6 @@ Modification by:
 #include "objectHit.H"
 #include "boolList.H"
 #include "DynamicList.H"
-
 #include "dimensionedConstants.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -48,19 +47,21 @@ namespace Foam
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
 template<class MasterPatch, class SlavePatch>
-const scalar GGIInterpolation<MasterPatch, SlavePatch>::areaErrorTol_
+const Foam::debug::tolerancesSwitch
+GGIInterpolation<MasterPatch, SlavePatch>::areaErrorTol_
 (
-    debug::tolerances("GGIAreaErrorTol", 1.0e-8)
+    "GGIAreaErrorTol",
+    1.0e-8
 );
 
-
 template<class MasterPatch, class SlavePatch>
-const scalar GGIInterpolation<MasterPatch, SlavePatch>::featureCosTol_
+const Foam::debug::tolerancesSwitch
+GGIInterpolation<MasterPatch, SlavePatch>::featureCosTol_
 (
-    debug::tolerances("GGIFeatureCosTol", 0.8)
+    "GGIFeatureCosTol",
+    0.8
 );
 
-
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
 template<class MasterPatch, class SlavePatch>
@@ -354,7 +355,7 @@ void GGIInterpolation<MasterPatch, SlavePatch>::calcAddressing() const
                 (
                     masterPointsInUV,
                     neighbPointsInUV,
-                    sqrt(areaErrorTol_) // distErrorTol
+                    sqrt(areaErrorTol_()) // distErrorTol
                 )
             )
             {
diff --git a/src/foam/interpolations/patchToPatchInterpolation/CalcPatchToPatchWeights.C b/src/foam/interpolations/patchToPatchInterpolation/CalcPatchToPatchWeights.C
index 605d1291e..8681c2c5f 100644
--- a/src/foam/interpolations/patchToPatchInterpolation/CalcPatchToPatchWeights.C
+++ b/src/foam/interpolations/patchToPatchInterpolation/CalcPatchToPatchWeights.C
@@ -112,7 +112,7 @@ void PatchToPatchInterpolation<FromPatch, ToPatch>::calcPointAddressing() const
             // Grab hit point
             hitPoint = curHit.hitPoint();
         }
-        else if (projectionTol_ > SMALL)
+        else if (projectionTol_() > SMALL)
         {
             // Check for a near miss
             pointHit ph =
@@ -162,7 +162,7 @@ void PatchToPatchInterpolation<FromPatch, ToPatch>::calcPointAddressing() const
                     );
             }
 
-            if (dist < minEdgeLength*projectionTol_)
+            if (dist < minEdgeLength*projectionTol_())
             {
                 // This point is being corrected
                 doWeights = true;
@@ -304,7 +304,7 @@ void PatchToPatchInterpolation<FromPatch, ToPatch>::calcFaceAddressing() const
 
             if
             (
-                m < directHitTol_                           // Direct hit
+                m < directHitTol_()                           // Direct hit
              || neighbours.empty()
             )
             {
diff --git a/src/foam/interpolations/patchToPatchInterpolation/PatchToPatchInterpolation.C b/src/foam/interpolations/patchToPatchInterpolation/PatchToPatchInterpolation.C
index bd05408bf..e8bd550ff 100644
--- a/src/foam/interpolations/patchToPatchInterpolation/PatchToPatchInterpolation.C
+++ b/src/foam/interpolations/patchToPatchInterpolation/PatchToPatchInterpolation.C
@@ -38,19 +38,21 @@ namespace Foam
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
 template<class FromPatch, class ToPatch>
-scalar PatchToPatchInterpolation<FromPatch, ToPatch>::directHitTol_
+const Foam::debug::tolerancesSwitch
+PatchToPatchInterpolation<FromPatch, ToPatch>::directHitTol_
 (
-    debug::tolerances("patchToPatchDirectHit", 1e-5)
+    "patchToPatchDirectHit",
+    1e-5
 );
 
-
 template<class FromPatch, class ToPatch>
-scalar PatchToPatchInterpolation<FromPatch, ToPatch>::projectionTol_
+const Foam::debug::tolerancesSwitch
+PatchToPatchInterpolation<FromPatch, ToPatch>::projectionTol_
 (
-    debug::tolerances("patchToPatchProjectionTol", 0.05)
+    "patchToPatchProjectionTol",
+    0.05
 );
 
-
 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
 
 template<class FromPatch, class ToPatch>
@@ -109,7 +111,8 @@ PatchToPatchInterpolation<FromPatch, ToPatch>::PatchToPatchInterpolation
     faceAddressingPtr_(NULL),
     faceWeightsPtr_(NULL),
     faceDistancePtr_(NULL)
-{}
+{
+}
 
 
 // Construct as copy
diff --git a/src/foam/interpolations/patchToPatchInterpolation/PatchToPatchInterpolation.H b/src/foam/interpolations/patchToPatchInterpolation/PatchToPatchInterpolation.H
index 993723b60..127711311 100644
--- a/src/foam/interpolations/patchToPatchInterpolation/PatchToPatchInterpolation.H
+++ b/src/foam/interpolations/patchToPatchInterpolation/PatchToPatchInterpolation.H
@@ -45,6 +45,7 @@ SourceFiles
 #include "FieldFields.H"
 #include "faceList.H"
 #include "intersection.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -86,10 +87,10 @@ class PatchToPatchInterpolation
 
         //- Relative merge tolerance for projected points missing the target
         //  Expressed as the fraction of min involved edge size
-        static scalar projectionTol_;
+        static const Foam::debug::tolerancesSwitch projectionTol_;
 
         //- Direct hit tolerance
-        static scalar directHitTol_;
+        static const Foam::debug::tolerancesSwitch directHitTol_;
 
 
         // Point addressing
diff --git a/src/foam/matrices/blockLduMatrix/BlockLduMatrix/BlockLduMatrix.C b/src/foam/matrices/blockLduMatrix/BlockLduMatrix/BlockLduMatrix.C
index 2617656e9..6a8e05a3a 100644
--- a/src/foam/matrices/blockLduMatrix/BlockLduMatrix/BlockLduMatrix.C
+++ b/src/foam/matrices/blockLduMatrix/BlockLduMatrix/BlockLduMatrix.C
@@ -38,9 +38,11 @@ Description
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
 template<class Type>
-const Foam::label Foam::BlockLduMatrix<Type>::fixFillIn
+const Foam::debug::optimisationSwitch
+Foam::BlockLduMatrix<Type>::fixFillIn
 (
-    debug::optimisationSwitch("matrixConstraintFillIn", 4)
+    "matrixConstraintFillIn",
+    4
 );
 
 
@@ -56,7 +58,7 @@ Foam::BlockLduMatrix<Type>::BlockLduMatrix(const lduMesh& ldu)
     interfaces_(ldu.interfaces().size()),
     coupleUpper_(ldu.lduAddr().nPatches()),
     coupleLower_(ldu.lduAddr().nPatches()),
-    fixedEqns_(ldu.lduAddr().size()/fixFillIn)
+    fixedEqns_(ldu.lduAddr().size()/fixFillIn())
 {
     const lduAddressing& addr = ldu.lduAddr();
 
@@ -65,6 +67,7 @@ Foam::BlockLduMatrix<Type>::BlockLduMatrix(const lduMesh& ldu)
         coupleUpper_.set(i, new CoeffField<Type>(addr.patchAddr(i).size()));
         coupleLower_.set(i, new CoeffField<Type>(addr.patchAddr(i).size()));
     }
+
 }
 
 
diff --git a/src/foam/matrices/blockLduMatrix/BlockLduMatrix/BlockLduMatrix.H b/src/foam/matrices/blockLduMatrix/BlockLduMatrix/BlockLduMatrix.H
index a5cd36466..7ed521806 100644
--- a/src/foam/matrices/blockLduMatrix/BlockLduMatrix/BlockLduMatrix.H
+++ b/src/foam/matrices/blockLduMatrix/BlockLduMatrix/BlockLduMatrix.H
@@ -53,6 +53,7 @@ SourceFiles
 #include "HashSet.H"
 #include "BlockLduInterfaceFieldPtrsList.H"
 #include "Map.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -135,7 +136,7 @@ private:
 
         //- Matrix constraint fill-in
         //  Equals to the estimated fraction of fixed nodes in the matrix
-        static const label fixFillIn;
+        static const Foam::debug::optimisationSwitch fixFillIn;
 
 
     // Private member functions
diff --git a/src/foam/matrices/solution/solution.C b/src/foam/matrices/solution/solution.C
index e196d1bb5..69ce7f6e1 100644
--- a/src/foam/matrices/solution/solution.C
+++ b/src/foam/matrices/solution/solution.C
@@ -33,7 +33,7 @@ License
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
-int Foam::solution::debug(::Foam::debug::debugSwitch("solution", 0));
+int Foam::solution::debug(::Foam::debug::debugSwitchFromDict("solution", 0));
 
 // List of sub-dictionaries to rewrite
 //! @cond localScope
diff --git a/src/foam/meshes/meshShapes/cell/cell.C b/src/foam/meshes/meshShapes/cell/cell.C
index 80e78772d..82e95f3ba 100644
--- a/src/foam/meshes/meshShapes/cell/cell.C
+++ b/src/foam/meshes/meshShapes/cell/cell.C
@@ -30,7 +30,7 @@ License
 
 const char* const Foam::cell::typeName = "cell";
 
-int Foam::cell::debug(debug::debugSwitch("cell", 0));
+int Foam::cell::debug(debug::debugSwitchFromDict("cell", 0));
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
diff --git a/src/foam/meshes/polyMesh/globalMeshData/globalMeshData.C b/src/foam/meshes/polyMesh/globalMeshData/globalMeshData.C
index 2cdb247fc..d22a8d3fa 100644
--- a/src/foam/meshes/polyMesh/globalMeshData/globalMeshData.C
+++ b/src/foam/meshes/polyMesh/globalMeshData/globalMeshData.C
@@ -42,9 +42,11 @@ License
 defineTypeNameAndDebug(Foam::globalMeshData, 0);
 
 // Geometric matching tolerance. Factor of mesh bounding box.
-const Foam::scalar Foam::globalMeshData::matchTol_
+const Foam::debug::tolerancesSwitch
+Foam::globalMeshData::matchTol_
 (
-    debug::tolerances("globalMeshDataMatchTol", 1e-8)
+    "globalMeshDataMatchTol",
+    1e-8
 );
 
 
@@ -661,7 +663,7 @@ Foam::pointField Foam::globalMeshData::geometricSharedPoints() const
     combineReduce(sharedPoints, plusEqOp<pointField>());
 
     // Merge tolerance
-    scalar tolDim = matchTol_*bb_.mag();
+    scalar tolDim = matchTol_()*bb_.mag();
 
     // And see how many are unique
     labelList pMap;
@@ -729,7 +731,7 @@ void Foam::globalMeshData::updateMesh()
     // Note: boundBox does reduce
     bb_ = boundBox(mesh_.points());
 
-    scalar tolDim = matchTol_*bb_.mag();
+    scalar tolDim = matchTol_()*bb_.mag();
 
     if (debug)
     {
diff --git a/src/foam/meshes/polyMesh/globalMeshData/globalMeshData.H b/src/foam/meshes/polyMesh/globalMeshData/globalMeshData.H
index e4881fa58..c240fe811 100644
--- a/src/foam/meshes/polyMesh/globalMeshData/globalMeshData.H
+++ b/src/foam/meshes/polyMesh/globalMeshData/globalMeshData.H
@@ -75,6 +75,7 @@ SourceFiles
 #include "Switch.H"
 #include "processorTopology.H"
 #include "labelPair.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -230,7 +231,7 @@ public:
     // Static data members
 
         //- Geomtric tolerance (fraction of bounding box)
-        static const Foam::scalar matchTol_;
+        static const Foam::debug::tolerancesSwitch matchTol_;
 
 
     // Constructors
diff --git a/src/foam/meshes/polyMesh/polyPatches/basic/coupled/coupledPolyPatch.H b/src/foam/meshes/polyMesh/polyPatches/basic/coupled/coupledPolyPatch.H
index ac0f30f22..74ddb3902 100644
--- a/src/foam/meshes/polyMesh/polyPatches/basic/coupled/coupledPolyPatch.H
+++ b/src/foam/meshes/polyMesh/polyPatches/basic/coupled/coupledPolyPatch.H
@@ -82,7 +82,7 @@ protected:
             const vectorField& nf,
             const vectorField& nr,
             const scalarField& smallDist,
-            const scalar absTol = matchTol_
+            const scalar absTol = matchTol_()
         ) const;
 
         //- Initialise the calculation of the patch addressing
diff --git a/src/foam/meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C b/src/foam/meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C
index 0afd406df..37fea5997 100644
--- a/src/foam/meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C
+++ b/src/foam/meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C
@@ -165,7 +165,7 @@ void Foam::cyclicPolyPatch::calcTransforms() const
             else if
             (
                 mag(magSf - nbrMagSf)/avSf
-             > polyPatch::matchTol_
+             > polyPatch::matchTol_()
             )
             {
                 // Error in area matching.  Find largest error
@@ -280,7 +280,7 @@ void Foam::cyclicPolyPatch::calcTransforms() const
         }
 
         // Check for error in face matching
-        if (maxMatchError > polyPatch::matchTol_)
+        if (maxMatchError > polyPatch::matchTol_())
         {
             label nbrFacei = errorFace + size()/2;
             scalar magSf = mag(half0Normals[errorFace]);
@@ -298,7 +298,7 @@ void Foam::cyclicPolyPatch::calcTransforms() const
                 << "patch:" << name()
                 << " my area:" << magSf
                 << " neighbour area:" << nbrMagSf
-                << " matching tolerance:" << polyPatch::matchTol_
+                << " matching tolerance:" << polyPatch::matchTol_()
                 << endl
                 << "Mesh face:" << start() + errorFace
                 << " vertices:"
@@ -384,7 +384,7 @@ void Foam::cyclicPolyPatch::calcTransforms() const
 
             // Check max distance between face centre and
             // transformed face centre
-            if (maxRelDistance > sqrt(polyPatch::matchTol_))
+            if (maxRelDistance > sqrt(polyPatch::matchTol_()))
             {
                 SeriousErrorIn
                 (
@@ -633,7 +633,7 @@ void Foam::cyclicPolyPatch::getCentresAndAnchors
         }
     }
 
-    if (mag(n0 & n1) < 1 - polyPatch::matchTol_)
+    if (mag(n0 & n1) < 1 - polyPatch::matchTol_())
     {
         if (debug)
         {
diff --git a/src/foam/meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C b/src/foam/meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C
index 532c54fac..f677525ba 100644
--- a/src/foam/meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C
+++ b/src/foam/meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C
@@ -209,7 +209,7 @@ void Foam::processorPolyPatch::calcGeometry()
                 faceNormals[facei] = point(1, 0, 0);
                 nbrFaceNormals[facei] = faceNormals[facei];
             }
-            else if (mag(magSf - nbrMagSf)/avSf > polyPatch::matchTol_)
+            else if (mag(magSf - nbrMagSf)/avSf > polyPatch::matchTol_())
             {
                 FatalErrorIn
                 (
@@ -220,7 +220,7 @@ void Foam::processorPolyPatch::calcGeometry()
                     << "patch: " << name()
                     << " my area: " << magSf
                     << " neighbour area: " << nbrMagSf
-                    << " matching tolerance: " << polyPatch::matchTol_
+                    << " matching tolerance: " << polyPatch::matchTol_()
                     << endl
                     << "Mesh face: " << start()+facei
                     << " vertices: "
diff --git a/src/foam/meshes/polyMesh/polyPatches/polyPatch/polyPatch.C b/src/foam/meshes/polyMesh/polyPatches/polyPatch/polyPatch.C
index 9a4df152a..de7781df1 100644
--- a/src/foam/meshes/polyMesh/polyPatches/polyPatch/polyPatch.C
+++ b/src/foam/meshes/polyMesh/polyPatches/polyPatch/polyPatch.C
@@ -41,7 +41,7 @@ namespace Foam
 
     int polyPatch::disallowGenericPolyPatch
     (
-        debug::debugSwitch("disallowGenericPolyPatch", 0)
+        debug::debugSwitchFromDict("disallowGenericPolyPatch", 0)
     );
 
     defineRunTimeSelectionTable(polyPatch, word);
@@ -52,12 +52,13 @@ namespace Foam
 }
 
 
-const Foam::scalar Foam::polyPatch::matchTol_
+const Foam::debug::tolerancesSwitch
+Foam::polyPatch::matchTol_
 (
-    debug::tolerances("patchFaceMatchTol", 1e-3)
+    "patchFaceMatchTol",
+    1e-3
 );
 
-
 // * * * * * * * * * * * * Protected Member Functions  * * * * * * * * * * * //
 
 void Foam::polyPatch::movePoints(const pointField& p)
@@ -253,7 +254,7 @@ Foam::scalarField Foam::polyPatch::calcFaceTol
             maxLen = max(maxLen, mag(points[f[fp]] - cc));
         }
 
-        tols[faceI] = matchTol_*maxLen;
+        tols[faceI] = matchTol_()*maxLen;
     }
 
     return tols;
diff --git a/src/foam/meshes/polyMesh/polyPatches/polyPatch/polyPatch.H b/src/foam/meshes/polyMesh/polyPatches/polyPatch/polyPatch.H
index 8241935a0..eeacba954 100644
--- a/src/foam/meshes/polyMesh/polyPatches/polyPatch/polyPatch.H
+++ b/src/foam/meshes/polyMesh/polyPatches/polyPatch/polyPatch.H
@@ -51,6 +51,7 @@ SourceFiles
 #include "primitivePatch.H"
 #include "typeInfo.H"
 #include "runTimeSelectionTables.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -230,8 +231,7 @@ public:
         //- Relative tolerance (for geometric matching)
         //  Calculated as a factor of maximum edge length per face.
         //  HJ, 21/Dec/2006
-        static const scalar matchTol_;
-
+        static const Foam::debug::tolerancesSwitch matchTol_;
 
 
     // Declare run-time constructor selection tables
diff --git a/src/foam/meshes/primitiveMesh/PrimitivePatch/PrimitivePatch.C b/src/foam/meshes/primitiveMesh/PrimitivePatch/PrimitivePatch.C
index fd7110997..e9bc55a3f 100644
--- a/src/foam/meshes/primitiveMesh/PrimitivePatch/PrimitivePatch.C
+++ b/src/foam/meshes/primitiveMesh/PrimitivePatch/PrimitivePatch.C
@@ -24,6 +24,7 @@ License
 \*---------------------------------------------------------------------------*/
 
 #include "Map.H"
+#include "debug.H"
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
@@ -34,11 +35,12 @@ template
     class PointField,
     class PointType
 >
-const bool
+const Foam::debug::optimisationSwitch
 Foam::PrimitivePatch<Face, FaceList, PointField, PointType>::
 nSquaredProjection_
 (
-    debug::optimisationSwitch("nSquaredProjection", 0) > 0
+    "nSquaredProjection",
+    0
 );
 
 
@@ -77,7 +79,8 @@ PrimitivePatch
     faceCentresPtr_(NULL),
     faceNormalsPtr_(NULL),
     pointNormalsPtr_(NULL)
-{}
+{
+}
 
 
 template
diff --git a/src/foam/meshes/primitiveMesh/PrimitivePatch/PrimitivePatch.H b/src/foam/meshes/primitiveMesh/PrimitivePatch/PrimitivePatch.H
index 2a45f6487..63e89d1b4 100644
--- a/src/foam/meshes/primitiveMesh/PrimitivePatch/PrimitivePatch.H
+++ b/src/foam/meshes/primitiveMesh/PrimitivePatch/PrimitivePatch.H
@@ -57,6 +57,7 @@ SourceFiles
 #include "point.H"
 #include "intersection.H"
 #include "HashSet.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -177,7 +178,7 @@ private:
 
     // Static Private Data
 
-        static const bool nSquaredProjection_;
+        static const Foam::debug::optimisationSwitch nSquaredProjection_;
 
 
     // Private Member Functions
diff --git a/src/foam/meshes/primitiveMesh/PrimitivePatch/PrimitivePatchProjectPoints.C b/src/foam/meshes/primitiveMesh/PrimitivePatch/PrimitivePatchProjectPoints.C
index 753b9231b..471d4e111 100644
--- a/src/foam/meshes/primitiveMesh/PrimitivePatch/PrimitivePatchProjectPoints.C
+++ b/src/foam/meshes/primitiveMesh/PrimitivePatch/PrimitivePatchProjectPoints.C
@@ -399,7 +399,7 @@ projectFaceCentres
 
         // Force the full search for the first point to ensure good
         // starting face
-        if (faceI == 0 || nSquaredProjection_)
+        if (faceI == 0 || nSquaredProjection_() > 0)
         {
             doNSquaredSearch = true;
         }
diff --git a/src/foam/meshes/primitiveMesh/primitiveMesh.H b/src/foam/meshes/primitiveMesh/primitiveMesh.H
index d6bc47a84..fa4693c86 100644
--- a/src/foam/meshes/primitiveMesh/primitiveMesh.H
+++ b/src/foam/meshes/primitiveMesh/primitiveMesh.H
@@ -66,6 +66,7 @@ SourceFiles
 #include "HashSet.H"
 #include "Map.H"
 #include "EdgeMap.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -321,22 +322,22 @@ public:
 
             //- Cell closedness warning threshold
             //  set as the fraction of un-closed area to closed area
-            static scalar closedThreshold_;
+            static const Foam::debug::tolerancesSwitch closedThreshold_;
 
             //- Aspect ratio warning threshold
-            static scalar aspectThreshold_;
+            static const Foam::debug::tolerancesSwitch aspectThreshold_;
 
             //- Non-orthogonality warning threshold in deg
-            static scalar nonOrthThreshold_;
+            static Foam::debug::tolerancesSwitch nonOrthThreshold_;
 
             //- Skewness warning threshold
-            static scalar skewThreshold_;
+            static const Foam::debug::tolerancesSwitch skewThreshold_;
 
             //- Face angle threshold
-            static scalar faceAngleThreshold_;
+            static Foam::debug::tolerancesSwitch faceAngleThreshold_;
 
             //- Face flatness threshold
-            static scalar faceFlatnessThreshold_;
+            static const Foam::debug::tolerancesSwitch faceFlatnessThreshold_;
 
 
     // Constructors
diff --git a/src/foam/meshes/primitiveMesh/primitiveMeshCheck/primitiveMeshCheck.C b/src/foam/meshes/primitiveMesh/primitiveMeshCheck/primitiveMeshCheck.C
index c2aba26fa..5b7f4afb4 100644
--- a/src/foam/meshes/primitiveMesh/primitiveMeshCheck/primitiveMeshCheck.C
+++ b/src/foam/meshes/primitiveMesh/primitiveMeshCheck/primitiveMeshCheck.C
@@ -31,35 +31,46 @@ License
 
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
-
-Foam::scalar Foam::primitiveMesh::closedThreshold_
+const Foam::debug::tolerancesSwitch
+Foam::primitiveMesh::closedThreshold_
 (
-    debug::tolerances("primitiveMeshClosedThreshold", 1e-6)
+    "primitiveMeshClosedThreshold",
+    1e-6
 );
 
-Foam::scalar Foam::primitiveMesh::aspectThreshold_
+const Foam::debug::tolerancesSwitch
+Foam::primitiveMesh::aspectThreshold_
 (
-    debug::tolerances("primitiveMeshAspectThreshold", 1000)
+    "primitiveMeshAspectThreshold",
+    1000
 );
 
-Foam::scalar Foam::primitiveMesh::nonOrthThreshold_ // deg
+Foam::debug::tolerancesSwitch
+Foam::primitiveMesh::nonOrthThreshold_ // deg
 (
-    debug::tolerances("primitiveMeshNonOrthThreshold", 70)
+    "primitiveMeshNonOrthThreshold",
+    70
 );
 
-Foam::scalar Foam::primitiveMesh::skewThreshold_
+const Foam::debug::tolerancesSwitch
+Foam::primitiveMesh::skewThreshold_
 (
-    debug::tolerances("primitiveMeshSkewThreshold", 4)
+    "primitiveMeshSkewThreshold",
+    4
 );
 
-Foam::scalar Foam::primitiveMesh::faceAngleThreshold_
+Foam::debug::tolerancesSwitch
+Foam::primitiveMesh::faceAngleThreshold_
 (
-    debug::tolerances("primitiveMeshFaceAngleThreshold", 10)
+    "primitiveMeshFaceAngleThreshold",
+    10
 );
 
-Foam::scalar Foam::primitiveMesh::faceFlatnessThreshold_
+const Foam::debug::tolerancesSwitch
+Foam::primitiveMesh::faceFlatnessThreshold_
 (
-    debug::tolerances("primitiveMeshFaceFlatnessThreshold", 0.8)
+    "primitiveMeshFaceFlatnessThreshold",
+    0.8
 );
 
 
@@ -93,12 +104,12 @@ bool Foam::primitiveMesh::checkClosedBoundary(const bool report) const
 
     vector openness = sumClosed/(sumMagClosedBoundary + VSMALL);
 
-    if (cmptMax(cmptMag(openness)) > closedThreshold_)
+    if (cmptMax(cmptMag(openness)) > closedThreshold_())
     {
         if (debug || report)
         {
             Info<< " ***Boundary openness " << openness
-                << " Threshold = " << closedThreshold_
+                << " Threshold = " << closedThreshold_()
                 << " possible hole in boundary description."
                 << endl;
         }
@@ -110,7 +121,7 @@ bool Foam::primitiveMesh::checkClosedBoundary(const bool report) const
         if (debug || report)
         {
             Info<< "    Boundary openness " << openness
-                << " Threshold = " << closedThreshold_
+                << " Threshold = " << closedThreshold_()
                 << " OK."
                 << endl;
         }
@@ -245,7 +256,7 @@ bool Foam::primitiveMesh::checkClosedCells
 
         maxOpennessCell = max(maxOpennessCell, maxOpenness);
 
-        if (maxOpenness > closedThreshold_)
+        if (maxOpenness > closedThreshold_())
         {
             if (setPtr)
             {
@@ -267,7 +278,7 @@ bool Foam::primitiveMesh::checkClosedCells
 
         maxAspectRatio = max(maxAspectRatio, aspectRatio);
 
-        if (aspectRatio > aspectThreshold_)
+        if (aspectRatio > aspectThreshold_())
         {
             if (aspectSetPtr)
             {
@@ -291,7 +302,7 @@ bool Foam::primitiveMesh::checkClosedCells
         {
             Info<< " ***Open cells found, max cell openness: "
                 << maxOpennessCell << ", number of open cells " << nOpen
-                << " Threshold = " << closedThreshold_
+                << " Threshold = " << closedThreshold_()
                 << endl;
         }
 
@@ -305,7 +316,7 @@ bool Foam::primitiveMesh::checkClosedCells
             Info<< " ***High aspect ratio cells found, Max aspect ratio: "
                 << maxAspectRatio
                 << ", number of cells " << nAspect
-                << " Threshold = " << aspectThreshold_
+                << " Threshold = " << aspectThreshold_()
                 << endl;
         }
 
@@ -471,7 +482,7 @@ bool Foam::primitiveMesh::checkFaceOrthogonality
 
     // Severe nonorthogonality threshold
     const scalar severeNonorthogonalityThreshold =
-        ::cos(nonOrthThreshold_/180.0*mathematicalConstant::pi);
+        ::cos(nonOrthThreshold_()/180.0*mathematicalConstant::pi);
 
     scalar minDDotS = GREAT;
 
@@ -536,7 +547,7 @@ bool Foam::primitiveMesh::checkFaceOrthogonality
                     << ::acos(minDDotS)/mathematicalConstant::pi*180.0
                     << " average: " <<
                     ::acos(sumDDotS/neiSize)/mathematicalConstant::pi*180.0
-                    << " Threshold = " << nonOrthThreshold_
+                    << " Threshold = " << nonOrthThreshold_()
                     << endl;
             }
         }
@@ -707,7 +718,7 @@ bool Foam::primitiveMesh::checkFaceSkewness
 
         // Check if the skewness vector is greater than the PN vector.
         // This does not cause trouble but is a good indication of a poor mesh.
-        if (skewness > skewThreshold_)
+        if (skewness > skewThreshold_())
         {
             if (setPtr)
             {
@@ -755,7 +766,7 @@ bool Foam::primitiveMesh::checkFaceSkewness
 
         // Check if the skewness vector is greater than the PN vector.
         // This does not cause trouble but is a good indication of a poor mesh.
-        if (skewness > skewThreshold_)
+        if (skewness > skewThreshold_())
         {
             if (setPtr)
             {
@@ -781,7 +792,7 @@ bool Foam::primitiveMesh::checkFaceSkewness
         {
             Info<< " ***Max skewness = " << maxSkew
                 << ", " << nWarnSkew << " highly skew faces detected"
-                << " Threshold = " << skewThreshold_
+                << " Threshold = " << skewThreshold_()
                 << endl;
         }
 
@@ -890,18 +901,18 @@ bool Foam::primitiveMesh::checkFaceAngles
             << "checking face angles" << endl;
     }
 
-    if (faceAngleThreshold_ < -SMALL || faceAngleThreshold_ > 180 + SMALL)
+    if (faceAngleThreshold_() < -SMALL || faceAngleThreshold_() > 180 + SMALL)
     {
         FatalErrorIn
         (
             "primitiveMesh::checkFaceAngles(const bool, labelHashSet*)"
         )   << "faceAngleThreshold_ should be [0..180] but is now "
-            << faceAngleThreshold_
+            << faceAngleThreshold_()
             << exit(FatalError);
     }
 
     const scalar maxSin =
-        Foam::sin(faceAngleThreshold_/180.0*mathematicalConstant::pi);
+        Foam::sin(faceAngleThreshold_()/180.0*mathematicalConstant::pi);
 
     const pointField& p = points();
     const faceList& fcs = faces();
@@ -1018,14 +1029,14 @@ bool Foam::primitiveMesh::checkFaceFlatness
             << "checking face flatness" << endl;
     }
 
-    if (faceFlatnessThreshold_ < 0 || faceFlatnessThreshold_ > 1)
+    if (faceFlatnessThreshold_() < 0 || faceFlatnessThreshold_() > 1)
     {
         FatalErrorIn
         (
             "primitiveMesh::checkFaceFlatness"
             "(const bool, labelHashSet*)"
         )   << "faceFlatnessThreshold_ should be [0..1] but is now "
-            << faceFlatnessThreshold_
+            << faceFlatnessThreshold_()
             << exit(FatalError);
     }
 
@@ -1074,7 +1085,7 @@ bool Foam::primitiveMesh::checkFaceFlatness
 
             minFlatness = min(minFlatness, flatness);
 
-            if (flatness < faceFlatnessThreshold_)
+            if (flatness < faceFlatnessThreshold_())
             {
                 nWarped++;
 
@@ -1109,7 +1120,7 @@ bool Foam::primitiveMesh::checkFaceFlatness
         {
             Info<< "   *There are " << nWarped
                 << " faces with ratio between projected and actual area < "
-                << faceFlatnessThreshold_ << endl;
+                << faceFlatnessThreshold_() << endl;
 
             Info<< "    Minimum ratio (minimum flatness, maximum warpage) = "
                 << minFlatness << endl;
diff --git a/src/foam/meshes/primitiveShapes/triangle/intersection.C b/src/foam/meshes/primitiveShapes/triangle/intersection.C
index 686c1a5c1..d9b0cdd05 100644
--- a/src/foam/meshes/primitiveShapes/triangle/intersection.C
+++ b/src/foam/meshes/primitiveShapes/triangle/intersection.C
@@ -29,18 +29,22 @@ Description
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
-const Foam::scalar Foam::intersection::planarTol_
+//Foam::debug::controlSwitches<Foam::scalar>
+Foam::debug::tolerancesSwitch
+Foam::intersection::planarTol_
 (
-    debug::tolerances("intersectionPlanarTol", 0.2)
+    "intersectionPlanarTol",
+    0.2
 );
 
-
-const Foam::scalar Foam::intersection::missTol_
+//const Foam::debug::controlSwitches<Foam::scalar>
+const Foam::debug::tolerancesSwitch
+Foam::intersection::missTol_
 (
-    debug::tolerances("intersectionMissTol", SMALL)
+    "intersectionMissTol",
+    SMALL
 );
 
-
 template<>
 const char* Foam::NamedEnum<Foam::intersection::direction, 2>::names[] =
 {
diff --git a/src/foam/meshes/primitiveShapes/triangle/intersection.H b/src/foam/meshes/primitiveShapes/triangle/intersection.H
index 3b2e103bf..190f6d3b2 100644
--- a/src/foam/meshes/primitiveShapes/triangle/intersection.H
+++ b/src/foam/meshes/primitiveShapes/triangle/intersection.H
@@ -37,6 +37,7 @@ SourceFiles
 
 #include "scalar.H"
 #include "NamedEnum.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -59,10 +60,11 @@ public:
         //  intersects the triangle or not.
         //  Adjusting the tolerance from the code does not allow it to be const
         //  HJ, 5/Jan/2008
-        static const scalar planarTol_;
+        //static Foam::debug::controlSwitches<scalar> planarTol_;
+        static Foam::debug::tolerancesSwitch planarTol_;
 
         //- Intersection miss tolerance, of the order of 1e-15
-        static const scalar missTol_;
+        static const Foam::debug::tolerancesSwitch missTol_;
 
 
         enum direction
@@ -101,9 +103,14 @@ public:
 
             // Ugly violence to allow tolerance to remain a const scalar
             // HJ, 7/Sep/2008
-            scalar oldTol = planarTol_;
-            scalar* ptPtr = const_cast<scalar*>(&planarTol_);
-            (*ptPtr) = t;
+            scalar oldTol = planarTol_();
+
+	    // Is this still necessary? If not a const, then it is not a const.
+            // MB 18/05_2014
+            //scalar* ptPtr = const_cast<scalar*>(&planarTol_);
+            //(*ptPtr) = t;
+
+	    planarTol_ = t;
 
             return oldTol;
         }
diff --git a/src/foam/meshes/primitiveShapes/triangle/triangleI.H b/src/foam/meshes/primitiveShapes/triangle/triangleI.H
index b417df5e5..74fe03498 100644
--- a/src/foam/meshes/primitiveShapes/triangle/triangleI.H
+++ b/src/foam/meshes/primitiveShapes/triangle/triangleI.H
@@ -100,20 +100,20 @@ bool triangle<Point, PointRef>::intersection
 
     // User-controlled projection tolerance
     // HJ, 23/Oct/2007
-    if (Foam::mag(u1) < intersection::missTol_)
+    if (Foam::mag(u1) < intersection::missTol_())
     {
         beta = u0/u2;
         if
         (
-            (beta > -intersection::missTol_)
-         && (beta < 1 + intersection::missTol_)
+            (beta > -intersection::missTol_())
+         && (beta < 1 + intersection::missTol_())
         )
         {
             alpha = (v0 - beta*v2)/v1;
             hit =
             (
-                (alpha > -intersection::missTol_)
-             && ((alpha + beta) < 1 + intersection::missTol_)
+                (alpha > -intersection::missTol_())
+             && ((alpha + beta) < 1 + intersection::missTol_())
             );
         }
     }
@@ -122,15 +122,15 @@ bool triangle<Point, PointRef>::intersection
         beta = (v0*u1 - u0*v1)/det;
         if
         (
-            (beta > -intersection::missTol_)
-         && (beta < 1 + intersection::missTol_)
+            (beta > -intersection::missTol_())
+         && (beta < 1 + intersection::missTol_())
         )
         {
             alpha = (u0 - beta*u2)/u1;
             hit =
             (
-                (alpha > -intersection::missTol_)
-             && ((alpha + beta) < 1 + intersection::missTol_)
+                (alpha > -intersection::missTol_())
+             && ((alpha + beta) < 1 + intersection::missTol_())
             );
         }
     }
@@ -552,7 +552,7 @@ inline pointHit triangle<Point, PointRef>::ray
                 Foam::mag(E1)
             ),
             Foam::mag(c_ - b_)
-        )*intersection::planarTol_;
+        )*intersection::planarTol_();
 
     bool eligible =
         alg == intersection::FULL_RAY
diff --git a/src/foam/primitives/strings/fileName/fileName.C b/src/foam/primitives/strings/fileName/fileName.C
index 2a1e96771..0355b0496 100644
--- a/src/foam/primitives/strings/fileName/fileName.C
+++ b/src/foam/primitives/strings/fileName/fileName.C
@@ -32,7 +32,7 @@ License
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
 const char* const Foam::fileName::typeName = "fileName";
-int Foam::fileName::debug(debug::debugSwitch(fileName::typeName, 0));
+int Foam::fileName::debug(debug::debugSwitchFromDict(fileName::typeName, 0));
 const Foam::fileName Foam::fileName::null;
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
diff --git a/src/foam/primitives/strings/string/string.C b/src/foam/primitives/strings/string/string.C
index 65986eeef..36e792a3d 100644
--- a/src/foam/primitives/strings/string/string.C
+++ b/src/foam/primitives/strings/string/string.C
@@ -29,7 +29,7 @@ License
 /* * * * * * * * * * * * * * * Static Member Data  * * * * * * * * * * * * * */
 
 const char* const Foam::string::typeName = "string";
-int Foam::string::debug(debug::debugSwitch(string::typeName, 0));
+int Foam::string::debug(debug::debugSwitchFromDict(string::typeName, 0));
 const Foam::string Foam::string::null;
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
diff --git a/src/foam/primitives/strings/word/word.C b/src/foam/primitives/strings/word/word.C
index 6177f55f3..713e9a516 100644
--- a/src/foam/primitives/strings/word/word.C
+++ b/src/foam/primitives/strings/word/word.C
@@ -29,7 +29,7 @@ License
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
 const char* const Foam::word::typeName = "word";
-int Foam::word::debug(Foam::debug::debugSwitch(word::typeName, 0));
+int Foam::word::debug(Foam::debug::debugSwitchFromDict(word::typeName, 0));
 const Foam::word Foam::word::null;
 
 // ************************************************************************* //
diff --git a/src/lagrangian/basic/Particle/Particle.H b/src/lagrangian/basic/Particle/Particle.H
index 04ed89ec6..15c4727fe 100644
--- a/src/lagrangian/basic/Particle/Particle.H
+++ b/src/lagrangian/basic/Particle/Particle.H
@@ -489,7 +489,8 @@ public:
     template<>                                                                \
     const Foam::word Particle<Type>::typeName(#Type);                         \
     template<>                                                                \
-    int Particle<Type>::debug(Foam::debug::debugSwitch(#Type, DebugSwitch));
+    Foam::debug::debugSwitch                                                  \
+    Particle<Type>::debug(std::string(#Type), DebugSwitch);
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
diff --git a/src/lagrangian/coalCombustion/submodels/surfaceReactionModel/COxidationKineticDiffusionLimitedRate/COxidationKineticDiffusionLimitedRate.C b/src/lagrangian/coalCombustion/submodels/surfaceReactionModel/COxidationKineticDiffusionLimitedRate/COxidationKineticDiffusionLimitedRate.C
index 020b1dff1..264a4b2f0 100644
--- a/src/lagrangian/coalCombustion/submodels/surfaceReactionModel/COxidationKineticDiffusionLimitedRate/COxidationKineticDiffusionLimitedRate.C
+++ b/src/lagrangian/coalCombustion/submodels/surfaceReactionModel/COxidationKineticDiffusionLimitedRate/COxidationKineticDiffusionLimitedRate.C
@@ -134,13 +134,13 @@ Foam::scalar Foam::COxidationKineticDiffusionLimitedRate<CloudType>::calculate
     const scalar D0 = C1_/d*pow(0.5*(T + Tc), 0.75);
 
     // Kinetic rate
-    const scalar Rk = C2_*exp(-E_/(specie::RR*Tc));
+    const scalar Rk = C2_*exp(-E_/(specie::RR()*Tc));
 
     // Particle surface area
     const scalar Ap = mathematicalConstant::pi*sqr(d);
 
     // Change in C mass [kg]
-    scalar dmC = Ap*rhoc*specie::RR*Tc*YO2/WO2_*D0*Rk/(D0 + Rk);
+    scalar dmC = Ap*rhoc*specie::RR()*Tc*YO2/WO2_*D0*Rk/(D0 + Rk);
 
     // Limit mass transfer by availability of C
     dmC = min(mass*fComb, dmC);
diff --git a/src/lagrangian/coalCombustion/submodels/surfaceReactionModel/COxidationMurphyShaddix/COxidationMurphyShaddix.C b/src/lagrangian/coalCombustion/submodels/surfaceReactionModel/COxidationMurphyShaddix/COxidationMurphyShaddix.C
index 579758a4c..6660585f1 100644
--- a/src/lagrangian/coalCombustion/submodels/surfaceReactionModel/COxidationMurphyShaddix/COxidationMurphyShaddix.C
+++ b/src/lagrangian/coalCombustion/submodels/surfaceReactionModel/COxidationMurphyShaddix/COxidationMurphyShaddix.C
@@ -140,10 +140,10 @@ Foam::scalar Foam::COxidationMurphyShaddix<CloudType>::calculate
     const scalar D = D0_*(rho0_/rhoc)*pow(Tc/T0_, Dn_);
 
     // Far field partial pressure O2 [Pa]
-    const scalar ppO2 = rhoO2/WO2_*specie::RR*Tc;
+    const scalar ppO2 = rhoO2/WO2_*specie::RR()*Tc;
 
     // Total molar concentration of the carrier phase [kmol/m^3]
-    const scalar C = pc/(specie::RR*Tc);
+    const scalar C = pc/(specie::RR()*Tc);
 
     if (debug)
     {
@@ -171,7 +171,7 @@ Foam::scalar Foam::COxidationMurphyShaddix<CloudType>::calculate
     {
         qCsOld = qCs;
         const scalar PO2Surface = ppO2*exp(-(qCs + N)*d/(2*C*D));
-        qCs = A_*exp(-E_/(specie::RR*T))*pow(PO2Surface, n_);
+        qCs = A_*exp(-E_/(specie::RR()*T))*pow(PO2Surface, n_);
         qCs = (100.0*qCs + iter*qCsOld)/(100.0 + iter);
         qCs = min(qCs, qCsLim);
 
diff --git a/src/lagrangian/dieselSpray/parcel/setRelaxationTimes.C b/src/lagrangian/dieselSpray/parcel/setRelaxationTimes.C
index 1e1166583..09c099f47 100644
--- a/src/lagrangian/dieselSpray/parcel/setRelaxationTimes.C
+++ b/src/lagrangian/dieselSpray/parcel/setRelaxationTimes.C
@@ -114,7 +114,7 @@ void parcel::setRelaxationTimes
     // Assume equilibrium at drop-surface => pressure @ surface
     // = vapour pressure to calculate fuel-vapour density @ surface
     scalar pressureAtSurface = fuels.pv(pressure, T(), X());
-    scalar rhoFuelVap = pressureAtSurface*fuels.W(X())/(specie::RR*Tf);
+    scalar rhoFuelVap = pressureAtSurface*fuels.W(X())/(specie::RR()*Tf);
 
     scalarField Xs(sDB.fuels().Xs(pressure, temperature, T(), Xf, X()));
     scalarField Ys(Nf, 0.0);
diff --git a/src/lagrangian/dieselSpray/spraySubModels/atomizationModel/LISA/LISA.C b/src/lagrangian/dieselSpray/spraySubModels/atomizationModel/LISA/LISA.C
index 745e5465d..c55e678ab 100644
--- a/src/lagrangian/dieselSpray/spraySubModels/atomizationModel/LISA/LISA.C
+++ b/src/lagrangian/dieselSpray/spraySubModels/atomizationModel/LISA/LISA.C
@@ -98,7 +98,7 @@ void LISA::atomizeParcel
         Winv += Y[i][cellI]/spray_.gasProperties()[i].W();
     }
 
-    scalar R = specie::RR*Winv;
+    scalar R = specie::RR()*Winv;
 
     // ideal gas law to evaluate density
     scalar rhoAverage = pressure/R/Taverage;
diff --git a/src/lagrangian/dieselSpray/spraySubModels/atomizationModel/blobsSheetAtomization/blobsSheetAtomization.C b/src/lagrangian/dieselSpray/spraySubModels/atomizationModel/blobsSheetAtomization/blobsSheetAtomization.C
index 609335753..26deae9ec 100644
--- a/src/lagrangian/dieselSpray/spraySubModels/atomizationModel/blobsSheetAtomization/blobsSheetAtomization.C
+++ b/src/lagrangian/dieselSpray/spraySubModels/atomizationModel/blobsSheetAtomization/blobsSheetAtomization.C
@@ -94,7 +94,7 @@ void blobsSheetAtomization::atomizeParcel
     {
         Winv += Y[i][cellI]/spray_.gasProperties()[i].W();
     }
-    scalar R = specie::RR*Winv;
+    scalar R = specie::RR()*Winv;
 
     // ideal gas law to evaluate density
     scalar rhoAverage = pressure/R/Taverage;
diff --git a/src/lagrangian/dieselSpray/spraySubModels/breakupModel/reitzDiwakar/reitzDiwakar.C b/src/lagrangian/dieselSpray/spraySubModels/breakupModel/reitzDiwakar/reitzDiwakar.C
index f69a14798..30382601a 100644
--- a/src/lagrangian/dieselSpray/spraySubModels/breakupModel/reitzDiwakar/reitzDiwakar.C
+++ b/src/lagrangian/dieselSpray/spraySubModels/breakupModel/reitzDiwakar/reitzDiwakar.C
@@ -102,7 +102,7 @@ void reitzDiwakar::breakupParcel
         Winv += Y[i][cellI]/spray_.gasProperties()[i].W();
         muAverage += Y[i][cellI]*spray_.gasProperties()[i].mu(Taverage);
     }
-    scalar R = specie::RR*Winv;
+    scalar R = specie::RR()*Winv;
 
     // ideal gas law to evaluate density
     scalar rhoAverage = pressure/R/Taverage;
diff --git a/src/lagrangian/dsmc/parcels/Templates/DsmcParcel/DsmcParcel.H b/src/lagrangian/dsmc/parcels/Templates/DsmcParcel/DsmcParcel.H
index 28f66d4cf..1c8fcd8e5 100644
--- a/src/lagrangian/dsmc/parcels/Templates/DsmcParcel/DsmcParcel.H
+++ b/src/lagrangian/dsmc/parcels/Templates/DsmcParcel/DsmcParcel.H
@@ -336,9 +336,10 @@ public:
     template<>                                                                \
     const Foam::word DsmcParcel<Type>::typeName(#Type);                       \
     template<>                                                                \
-    int DsmcParcel<Type>::debug                                               \
+    Foam::debug::debugSwitch						      \
+    DsmcParcel<Type>::debug						      \
     (                                                                         \
-        Foam::debug::debugSwitch(#Type, DebugSwitch)                          \
+        std::string(#Type), DebugSwitch       				      \
     );
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H
index 26dc950eb..c788100aa 100644
--- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H
+++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H
@@ -549,9 +549,10 @@ public:
     template<>                                                                \
     const Foam::word KinematicParcel<Type>::typeName(#Type);                  \
     template<>                                                                \
-    int KinematicParcel<Type>::debug                                          \
+    Foam::debug::debugSwitch						      \
+    KinematicParcel<Type>::debug                                              \
     (                                                                         \
-        Foam::debug::debugSwitch(#Type, DebugSwitch)                          \
+        std::string(#Type), DebugSwitch                                       \
     );
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
index ad51557e1..231fb8929 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
@@ -535,7 +535,7 @@ void Foam::ReactingMultiphaseParcel<ParcelType>::calcDevolatilisation
     Sh -= dMassTot*td.constProps().LDevol()/dt;
 
     // Molar average molecular weight of carrier mix
-    const scalar Wc = this->rhoc_*specie::RR*this->Tc_/this->pc_;
+    const scalar Wc = this->rhoc_*specie::RR()*this->Tc_/this->pc_;
 
     // Update molar emissions
     forAll(dMassDV, i)
diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
index d1a53d65b..91199f0a4 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
@@ -130,10 +130,10 @@ void Foam::ReactingParcel<ParcelType>::correctSurfaceValues
     Xinf /= sum(Xinf);
 
     // Molar fraction of far field species at particle surface
-    const scalar Xsff = 1.0 - min(sum(Cs)*specie::RR*this->T_/pc_, 1.0);
+    const scalar Xsff = 1.0 - min(sum(Cs)*specie::RR()*this->T_/pc_, 1.0);
 
     // Surface carrier total molar concentration
-    const scalar CsTot = pc_/(specie::RR*this->T_);
+    const scalar CsTot = pc_/(specie::RR()*this->T_);
 
     // Surface carrier composition (molar fraction)
     scalarField Xs(Xinf.size());
@@ -177,7 +177,7 @@ void Foam::ReactingParcel<ParcelType>::correctSurfaceValues
         sumYiCbrtW += Ys[i]*cbrtW;
     }
 
-    rhos *= pc_/(specie::RR*T);
+    rhos *= pc_/(specie::RR()*T);
     mus /= sumYiSqrtW;
     kappa /= sumYiCbrtW;
     Pr = cps*mus/kappa;
@@ -452,7 +452,7 @@ void Foam::ReactingParcel<ParcelType>::calcPhaseChange
     td.cloud().addToMassPhaseChange(this->nParticle_*dMassTot);
 
     // Average molecular weight of carrier mix - assumes perfect gas
-    scalar Wc = this->rhoc_*specie::RR*this->Tc_/this->pc_;
+    scalar Wc = this->rhoc_*specie::RR()*this->Tc_/this->pc_;
 
     forAll(YComponents, i)
     {
diff --git a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C
index 36ce81ee2..ffb07a5de 100644
--- a/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C
+++ b/src/lagrangian/intermediate/submodels/Reacting/PhaseChangeModel/LiquidEvaporation/LiquidEvaporation.C
@@ -177,10 +177,10 @@ void Foam::LiquidEvaporation<CloudType>::calculate
         scalar kc = Sh*Dab/(d + ROOTVSMALL);
 
         // vapour concentration at droplet surface [kmol/m3] at film temperature
-        scalar Cs = pSat/(specie::RR*Ts);
+        scalar Cs = pSat/(specie::RR()*Ts);
 
         // vapour concentration in bulk gas [kmol/m3] at film temperature
-        scalar Cinf = Xc[gid]*pc/(specie::RR*Ts);
+        scalar Cinf = Xc[gid]*pc/(specie::RR()*Ts);
 
         // molar flux of vapour [kmol/m2/s]
         scalar Ni = max(kc*(Cs - Cinf), 0.0);
diff --git a/src/lagrangian/intermediate/submodels/ReactingMultiphase/DevolatilisationModel/SingleKineticRateDevolatilisation/SingleKineticRateDevolatilisation.C b/src/lagrangian/intermediate/submodels/ReactingMultiphase/DevolatilisationModel/SingleKineticRateDevolatilisation/SingleKineticRateDevolatilisation.C
index 2d8647f9e..aa0fbdef7 100644
--- a/src/lagrangian/intermediate/submodels/ReactingMultiphase/DevolatilisationModel/SingleKineticRateDevolatilisation/SingleKineticRateDevolatilisation.C
+++ b/src/lagrangian/intermediate/submodels/ReactingMultiphase/DevolatilisationModel/SingleKineticRateDevolatilisation/SingleKineticRateDevolatilisation.C
@@ -83,7 +83,7 @@ Foam::scalar Foam::SingleKineticRateDevolatilisation<CloudType>::calculate
     }
 
     // Kinetic rate
-    const scalar kappa = A1_*exp(-E_/(specie::RR*T));
+    const scalar kappa = A1_*exp(-E_/(specie::RR()*T));
 
     // Volatile devolatilisation from particle to carrier gas phase
     const scalar dMass = min(dt*kappa*massVolatile, massVolatile);
diff --git a/src/lduSolvers/amg/amgPolicy/pamgPolicy.C b/src/lduSolvers/amg/amgPolicy/pamgPolicy.C
index df4771f56..8786df9da 100644
--- a/src/lduSolvers/amg/amgPolicy/pamgPolicy.C
+++ b/src/lduSolvers/amg/amgPolicy/pamgPolicy.C
@@ -48,9 +48,11 @@ namespace Foam
 } // End namespace Foam
 
 
-const Foam::scalar Foam::pamgPolicy::diagFactor_
+const Foam::debug::tolerancesSwitch
+Foam::pamgPolicy::diagFactor_
 (
-    debug::tolerances("pamgDiagFactor", 1e-8)
+    "pamgDiagFactor",
+    1e-8
 );
 
 
@@ -160,7 +162,7 @@ void Foam::pamgPolicy::calcChild()
         // Go through all upper and lower coefficients and for the ones
         // larger than threshold mark the equations out of cluster zero
 
-        scalarField magScaledDiag = diagFactor_*mag(diag);
+        scalarField magScaledDiag = diagFactor_()*mag(diag);
 
         boolList zeroCluster(diag.size(), true);
 
diff --git a/src/lduSolvers/amg/amgPolicy/pamgPolicy.H b/src/lduSolvers/amg/amgPolicy/pamgPolicy.H
index 451eae50d..46d1f3e6a 100644
--- a/src/lduSolvers/amg/amgPolicy/pamgPolicy.H
+++ b/src/lduSolvers/amg/amgPolicy/pamgPolicy.H
@@ -40,6 +40,7 @@ SourceFiles
 
 #include "amgPolicy.H"
 #include "lduMatrix.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -84,7 +85,7 @@ class pamgPolicy
     // Private Static Data
 
         //- Diagonal scaling factor
-        static const scalar diagFactor_;
+        static const Foam::debug::tolerancesSwitch diagFactor_;
 
 
 public:
diff --git a/src/meshTools/triSurface/booleanOps/surfaceIntersection/surfaceIntersection.C b/src/meshTools/triSurface/booleanOps/surfaceIntersection/surfaceIntersection.C
index 3c3ea563c..d76442c5c 100644
--- a/src/meshTools/triSurface/booleanOps/surfaceIntersection/surfaceIntersection.C
+++ b/src/meshTools/triSurface/booleanOps/surfaceIntersection/surfaceIntersection.C
@@ -598,7 +598,7 @@ void Foam::surfaceIntersection::doCutEdges
     forAll(surf1PointTol, pointI)
     {
         surf1PointTol[pointI] =
-            intersection::planarTol_
+            intersection::planarTol_()
           * minEdgeLen(surf1, pointI);
     }
 
@@ -611,7 +611,7 @@ void Foam::surfaceIntersection::doCutEdges
         point pStart = surf1Pts[e.start()];
         const point& pEnd = surf1Pts[e.end()];
 
-        const point tolVec = intersection::planarTol_*(pEnd-pStart);
+        const point tolVec = intersection::planarTol_()*(pEnd-pStart);
         const scalar tolDim = mag(tolVec);
 
         bool doTrack = false;
@@ -666,7 +666,7 @@ void Foam::surfaceIntersection::doCutEdges
                                     << "  hit:" << pHit.hitPoint()
                                     << "  tolDim:" << tolDim
                                     << "  planarTol:"
-                                    << intersection::planarTol_
+                                    << intersection::planarTol_()
                                     << endl;
                             }
                             allCutPoints.append(pHit.hitPoint());
@@ -1074,7 +1074,7 @@ Foam::surfaceIntersection::surfaceIntersection
     bool hasMerged = mergePoints
     (
         cutPoints_,
-        minEdgeLen*intersection::planarTol_,
+        minEdgeLen*intersection::planarTol_(),
         false,
         pointMap,
         newPoints
@@ -1085,7 +1085,7 @@ Foam::surfaceIntersection::surfaceIntersection
         if (debug)
         {
             Pout<< "Merged:" << hasMerged
-                << "  mergeDist:" << minEdgeLen*intersection::planarTol_
+                << "  mergeDist:" << minEdgeLen*intersection::planarTol_()
                 << "  cutPoints:" << cutPoints_.size()
                 << "  newPoints:" << newPoints.size()
                 << endl;
@@ -1102,7 +1102,7 @@ Foam::surfaceIntersection::surfaceIntersection
             e.start() = pointMap[e.start()];
             e.end() = pointMap[e.end()];
 
-            if (e.mag(cutPoints_) < minEdgeLen*intersection::planarTol_)
+            if (e.mag(cutPoints_) < minEdgeLen*intersection::planarTol_())
             {
                 if (debug)
                 {
diff --git a/src/sampling/meshToMeshInterpolation/meshToMesh/calculateMeshToMeshAddressing.C b/src/sampling/meshToMeshInterpolation/meshToMesh/calculateMeshToMeshAddressing.C
index 296fb963c..8fef6785b 100644
--- a/src/sampling/meshToMeshInterpolation/meshToMesh/calculateMeshToMeshAddressing.C
+++ b/src/sampling/meshToMeshInterpolation/meshToMesh/calculateMeshToMeshAddressing.C
@@ -40,9 +40,11 @@ Description
 namespace Foam
 {
 
-const scalar meshToMesh::cellCentreDistanceTol
+const Foam::debug::tolerancesSwitch
+meshToMesh::cellCentreDistanceTol
 (
-    debug::tolerances("meshToMeshCellCentreDistanceTol", 1e-3)
+    "meshToMeshCellCentreDistanceTol",
+    1e-3
 );
 
 
@@ -263,7 +265,7 @@ void meshToMesh::cellAddresses
 
     forAll (points, toI)
     {
-        scalar localTol = cellCentreDistanceTol;
+        scalar localTol = cellCentreDistanceTol();
 
         bool isBoundary = false;
 
diff --git a/src/sampling/meshToMeshInterpolation/meshToMesh/calculateMeshToMeshWeights.C b/src/sampling/meshToMeshInterpolation/meshToMesh/calculateMeshToMeshWeights.C
index b293f04b8..979cfabad 100644
--- a/src/sampling/meshToMeshInterpolation/meshToMesh/calculateMeshToMeshWeights.C
+++ b/src/sampling/meshToMeshInterpolation/meshToMesh/calculateMeshToMeshWeights.C
@@ -68,7 +68,7 @@ void meshToMesh::calculateInverseDistanceWeights() const
             // pick up the value
             if
             (
-                m < directHitTol                            // Direct hit
+                m < directHitTol()                            // Direct hit
              || neighbours.empty()
             )
             {
diff --git a/src/sampling/meshToMeshInterpolation/meshToMesh/meshToMesh.C b/src/sampling/meshToMeshInterpolation/meshToMesh/meshToMesh.C
index dfbc32335..e94fa6778 100644
--- a/src/sampling/meshToMeshInterpolation/meshToMesh/meshToMesh.C
+++ b/src/sampling/meshToMeshInterpolation/meshToMesh/meshToMesh.C
@@ -36,9 +36,11 @@ namespace Foam
 
 defineTypeNameAndDebug(meshToMesh, 0);
 
-const scalar meshToMesh::directHitTol
+const Foam::debug::tolerancesSwitch
+meshToMesh::directHitTol
 (
-    debug::tolerances("meshToMeshdirectHitTol", 1e-5)
+    "meshToMeshdirectHitTol",
+    1e-5
 );
 
 
diff --git a/src/sampling/meshToMeshInterpolation/meshToMesh/meshToMesh.H b/src/sampling/meshToMeshInterpolation/meshToMesh/meshToMesh.H
index 56b5cc3e0..fbe097671 100644
--- a/src/sampling/meshToMeshInterpolation/meshToMesh/meshToMesh.H
+++ b/src/sampling/meshToMeshInterpolation/meshToMesh/meshToMesh.H
@@ -43,6 +43,7 @@ SourceFiles
 #include "fvPatchMapper.H"
 #include "scalarList.H"
 #include "className.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -111,10 +112,10 @@ class meshToMesh
     // Private static data members
 
         //- Direct hit tolerance
-        static const scalar directHitTol;
+        static const Foam::debug::tolerancesSwitch directHitTol;
 
         //- Cell centre distance required to establish a hit
-        static const scalar cellCentreDistanceTol;
+        static const Foam::debug::tolerancesSwitch cellCentreDistanceTol;
 
 
 public:
diff --git a/src/thermophysicalModels/liquidMixture/liquidMixture/liquidMixture.C b/src/thermophysicalModels/liquidMixture/liquidMixture/liquidMixture.C
index e02f0b8a2..307ce1d1c 100644
--- a/src/thermophysicalModels/liquidMixture/liquidMixture/liquidMixture.C
+++ b/src/thermophysicalModels/liquidMixture/liquidMixture/liquidMixture.C
@@ -132,7 +132,7 @@ Foam::scalar Foam::liquidMixture::Ppc
         Zc += x[i]*properties_[i].Zc();
     }
 
-    return specie::RR*Zc*Tpc(x)/Vc;
+    return specie::RR()*Zc*Tpc(x)/Vc;
 }
 
 
diff --git a/src/thermophysicalModels/radiation/radiationConstants/radiationConstants.C b/src/thermophysicalModels/radiation/radiationConstants/radiationConstants.C
index ec319b2f8..3003c52b9 100644
--- a/src/thermophysicalModels/radiation/radiationConstants/radiationConstants.C
+++ b/src/thermophysicalModels/radiation/radiationConstants/radiationConstants.C
@@ -29,6 +29,13 @@ License
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 //- Stefan-Boltzmann constant (default in [J/(K4 m2 s)])
+const Foam::debug::constantsSwitch
+Foam::radiation::sigmaSB_
+(
+    "sigmaSB",
+    5.670E-08
+);
+
 const Foam::dimensionedScalar Foam::radiation::sigmaSB
 (
     Foam::dimensionedConstant
@@ -38,10 +45,9 @@ const Foam::dimensionedScalar Foam::radiation::sigmaSB
         (
             "sigmaSB",
             dimensionSet(1, 0, -3, -4, 0, 0, 0),
-            5.670E-08
+            sigmaSB_()
         )
     )
 );
 
-
 // ************************************************************************* //
diff --git a/src/thermophysicalModels/radiation/radiationConstants/radiationConstants.H b/src/thermophysicalModels/radiation/radiationConstants/radiationConstants.H
index 7c2086fb5..7c0e08bae 100644
--- a/src/thermophysicalModels/radiation/radiationConstants/radiationConstants.H
+++ b/src/thermophysicalModels/radiation/radiationConstants/radiationConstants.H
@@ -41,6 +41,7 @@ namespace Foam
 namespace radiation
 {
     //- Stefan-Boltzmann constant [J/(K4 m2 s)]
+    extern const Foam::debug::constantsSwitch sigmaSB_;
     extern const dimensionedScalar sigmaSB;
 }
 }
diff --git a/src/thermophysicalModels/specie/specie/specie.C b/src/thermophysicalModels/specie/specie/specie.C
index 958b5dc6d..4617019c4 100644
--- a/src/thermophysicalModels/specie/specie/specie.C
+++ b/src/thermophysicalModels/specie/specie/specie.C
@@ -33,14 +33,29 @@ Description
 /* * * * * * * * * * * * * public constants  * * * * * * * * * * * * */
 
 //- Universal gas constant (default in [J/(kmol K)])
-const Foam::scalar Foam::specie::RR = dimensionedConstant("R", 8314.51);
+const Foam::debug::constantsSwitch
+Foam::specie::RR
+(
+    "R",
+    8314.51
+);
+
 
 //- Standard pressure (default in [Pa])
-const Foam::scalar Foam::specie::Pstd = dimensionedConstant("Pstd", 1.0e5);
+const Foam::debug::constantsSwitch
+Foam::specie::Pstd
+(
+    "Pstd",
+    1.0e5
+);
 
 //- Standard temperature (default in [K])
-const Foam::scalar Foam::specie::Tstd = dimensionedConstant("Tstd", 298.15);
-
+const Foam::debug::constantsSwitch
+Foam::specie::Tstd
+(
+    "Tstd",
+    298.15
+);
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
diff --git a/src/thermophysicalModels/specie/specie/specie.H b/src/thermophysicalModels/specie/specie/specie.H
index af088fe95..92569cee1 100644
--- a/src/thermophysicalModels/specie/specie/specie.H
+++ b/src/thermophysicalModels/specie/specie/specie.H
@@ -38,6 +38,8 @@ SourceFiles
 
 #include "word.H"
 #include "scalar.H"
+#include "debug.H"
+#include "controlSwitches.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -82,13 +84,13 @@ public:
         // Thermodynamic constants
 
             //- Universal gas constant [J/(kmol K)]
-            static const scalar RR;
+            static const Foam::debug::constantsSwitch RR;
 
             //- Standard pressure [Pa]
-            static const scalar Pstd;
+            static const Foam::debug::constantsSwitch Pstd;
 
             //- Standard temperature [K]
-            static const scalar Tstd;
+            static const Foam::debug::constantsSwitch Tstd;
 
 
     // Constructors
diff --git a/src/thermophysicalModels/specie/specie/specieI.H b/src/thermophysicalModels/specie/specie/specieI.H
index ea0d99fff..3a9ec00ad 100644
--- a/src/thermophysicalModels/specie/specie/specieI.H
+++ b/src/thermophysicalModels/specie/specie/specieI.H
@@ -96,7 +96,7 @@ inline scalar specie::nMoles() const
 //- Gas constant [J/(kg K)]
 inline scalar specie::R() const
 {
-    return RR/molWeight_;
+    return RR()/molWeight_;
 }
 
 
diff --git a/src/thermophysicalModels/specie/thermo/eConst/eConstThermoI.H b/src/thermophysicalModels/specie/thermo/eConst/eConstThermoI.H
index f150bab3a..84167e621 100644
--- a/src/thermophysicalModels/specie/thermo/eConst/eConstThermoI.H
+++ b/src/thermophysicalModels/specie/thermo/eConst/eConstThermoI.H
@@ -84,7 +84,7 @@ inline Foam::scalar Foam::eConstThermo<equationOfState>::cp
     const scalar
 ) const
 {
-    return Cv_*this->W() + specie::RR;
+    return Cv_*this->W() + specie::RR();
 }
 
 
diff --git a/src/thermophysicalModels/specie/thermo/hPolynomial/hPolynomialThermo.C b/src/thermophysicalModels/specie/thermo/hPolynomial/hPolynomialThermo.C
index 73178dc27..b0d4d4bbe 100644
--- a/src/thermophysicalModels/specie/thermo/hPolynomial/hPolynomialThermo.C
+++ b/src/thermophysicalModels/specie/thermo/hPolynomial/hPolynomialThermo.C
@@ -49,10 +49,10 @@ Foam::hPolynomialThermo<EquationOfState, PolySize>::hPolynomialThermo
     sPolynomial_ = cpPolynomial_.integrateMinus1();
 
     // Offset h poly so that it is relative to the enthalpy at Tstd
-    hPolynomial_[0] += Hf_ - hPolynomial_.evaluate(specie::Tstd);
+    hPolynomial_[0] += Hf_ - hPolynomial_.evaluate(specie::Tstd());
 
     // Offset s poly so that it is relative to the entropy at Tstd
-    sPolynomial_[0] += Sf_ - sPolynomial_.evaluate(specie::Tstd);
+    sPolynomial_[0] += Sf_ - sPolynomial_.evaluate(specie::Tstd());
 }
 
 
diff --git a/src/thermophysicalModels/specie/thermo/janaf/janafThermoI.H b/src/thermophysicalModels/specie/thermo/janaf/janafThermoI.H
index b9108868c..d236ec4a3 100644
--- a/src/thermophysicalModels/specie/thermo/janaf/janafThermoI.H
+++ b/src/thermophysicalModels/specie/thermo/janaf/janafThermoI.H
@@ -125,7 +125,7 @@ inline Foam::scalar Foam::janafThermo<equationOfState>::cp
 ) const
 {
     const coeffArray& a = coeffs(T);
-    return this->RR*((((a[4]*T + a[3])*T + a[2])*T + a[1])*T + a[0]);
+    return this->RR()*((((a[4]*T + a[3])*T + a[2])*T + a[1])*T + a[0]);
 }
 
 
@@ -136,7 +136,7 @@ inline Foam::scalar Foam::janafThermo<equationOfState>::h
 ) const
 {
     const coeffArray& a = coeffs(T);
-    return this->RR*
+    return this->RR()*
     (
         ((((a[4]/5.0*T + a[3]/4.0)*T + a[2]/3.0)*T + a[1]/2.0)*T + a[0])*T
       + a[5]
@@ -158,8 +158,8 @@ template<class equationOfState>
 inline Foam::scalar Foam::janafThermo<equationOfState>::hc() const
 {
     const coeffArray& a = lowCpCoeffs_;
-    const scalar& Tstd = specie::Tstd;
-    return this->RR*
+    const scalar& Tstd = specie::Tstd();
+    return this->RR()*
     (
         (
             (((a[4]/5.0*Tstd + a[3]/4.0)*Tstd + a[2]/3.0)*Tstd + a[1]/2.0)*Tstd
@@ -180,7 +180,7 @@ inline Foam::scalar Foam::janafThermo<equationOfState>::s
 
     const coeffArray& a = coeffs(T);
     return
-    this->RR*
+    this->RR()*
     (
         (((a[4]/4.0*T + a[3]/3.0)*T + a[2]/2.0)*T + a[1])*T + a[0]*::log(T)
       + a[6]
diff --git a/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.C b/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.C
index 1218ba06b..d5aadcd4f 100644
--- a/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.C
+++ b/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.C
@@ -34,26 +34,31 @@ Description
 /* * * * * * * * * * * * * * * Private Static Data * * * * * * * * * * * * * */
 
 template<class thermo>
-const Foam::scalar Foam::specieThermo<thermo>::tol_
+const Foam::debug::tolerancesSwitch
+Foam::specieThermo<thermo>::tol_
 (
-    debug::tolerances("speciesThermoTol", 1e-4)
+    "speciesThermoTol",
+    1e-4
 );
 
 
 template<class thermo>
-const Foam::scalar Foam::specieThermo<thermo>::TJump_
+const Foam::debug::tolerancesSwitch
+Foam::specieThermo<thermo>::TJump_
 (
-    debug::tolerances("speciesThermoTJump", 20)
+    "speciesThermoTJump",
+    20
 );
 
 
 template<class thermo>
-const int Foam::specieThermo<thermo>::maxIter_
+const Foam::debug::optimisationSwitch
+Foam::specieThermo<thermo>::maxIter_
 (
-    debug::optimisationSwitch("speciesThermoMaxIter", 100)
+    "speciesThermoMaxIter",
+    100
 );
 
-
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
 template<class thermo>
diff --git a/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.H b/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.H
index 0748c62cd..3a612fb91 100644
--- a/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.H
+++ b/src/thermophysicalModels/specie/thermo/specieThermo/specieThermo.H
@@ -39,6 +39,7 @@ SourceFiles
 #define specieThermo_H
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+#include "controlSwitches.H"
 
 namespace Foam
 {
@@ -95,13 +96,13 @@ class specieThermo
     // Private data
 
         //- Convergence tolerance of energy -> temperature inversion functions
-        static const scalar tol_;
+        static const Foam::debug::tolerancesSwitch tol_;
 
         //- Max temperature jump of energy -> temperature inversion functions
-        static const scalar TJump_;
+        static const Foam::debug::tolerancesSwitch TJump_;
 
         //- Max number of iterations in energy->temperature inversion functions
-        static const int maxIter_;
+        static const Foam::debug::optimisationSwitch maxIter_;
 
 
     // Private member functions
diff --git a/src/thermophysicalModels/specie/thermo/specieThermo/specieThermoI.H b/src/thermophysicalModels/specie/thermo/specieThermo/specieThermoI.H
index f936d66e3..cfa8f8767 100644
--- a/src/thermophysicalModels/specie/thermo/specieThermo/specieThermoI.H
+++ b/src/thermophysicalModels/specie/thermo/specieThermo/specieThermoI.H
@@ -48,7 +48,7 @@ inline Foam::scalar Foam::specieThermo<thermo>::T
 {
     scalar Test = T0;
     scalar Tnew = T0;
-    scalar Ttol = T0*tol_;
+    scalar Ttol = T0*tol_();
     int iter = 0;
 
     do
@@ -58,9 +58,9 @@ inline Foam::scalar Foam::specieThermo<thermo>::T
 
         Test = Tnew;
         Tnew = Test
-            - Foam::min(((this->*F)(Test) - f)/(this->*dFdT)(Test), TJump_);
+            - Foam::min(((this->*F)(Test) - f)/(this->*dFdT)(Test), TJump_());
 
-        if (iter++ > maxIter_)
+        if (iter++ > maxIter_())
         {
             // Improvements: graceful exit with recovery.  HJ, 11/Oct/2010
             InfoIn
@@ -102,7 +102,7 @@ inline Foam::specieThermo<thermo>::specieThermo
 template<class thermo>
 inline Foam::scalar Foam::specieThermo<thermo>::cv(const scalar T) const
 {
-    return this->cp(T) - this->RR;
+    return this->cp(T) - this->RR();
 }
 
 
@@ -110,21 +110,21 @@ template<class thermo>
 inline Foam::scalar Foam::specieThermo<thermo>::gamma(const scalar T) const
 {
     scalar CP = this->cp(T);
-    return CP/(CP - this->RR);
+    return CP/(CP - this->RR());
 }
 
 
 template<class thermo>
 inline Foam::scalar Foam::specieThermo<thermo>::e(const scalar T) const
 {
-    return this->h(T) - this->RR*(T - this->Tstd);
+    return this->h(T) - this->RR()*(T - this->Tstd());
 }
 
 
 template<class thermo>
 inline Foam::scalar Foam::specieThermo<thermo>::es(const scalar T) const
 {
-    return this->hs(T) - this->RR*(T - this->Tstd);
+    return this->hs(T) - this->RR()*(T - this->Tstd());
 }
 
 
@@ -208,7 +208,7 @@ inline Foam::scalar Foam::specieThermo<thermo>::A(const scalar T) const
 template<class thermo>
 inline Foam::scalar Foam::specieThermo<thermo>::K(const scalar T) const
 {
-    scalar arg = -this->nMoles()*this->g(T)/(this->RR*T);
+    scalar arg = -this->nMoles()*this->g(T)/(this->RR()*T);
 
     if (arg < 600.0)
     {
@@ -237,7 +237,7 @@ inline Foam::scalar Foam::specieThermo<thermo>::Kc(const scalar T) const
     }
     else
     {
-        return Kp(T)*pow(this->Pstd/(this->RR*T), this->nMoles());
+        return Kp(T)*pow(this->Pstd()/(this->RR()*T), this->nMoles());
     }
 }
 
@@ -255,7 +255,7 @@ inline Foam::scalar Foam::specieThermo<thermo>::Kx
     }
     else
     {
-        return Kp(T)*pow(this->Pstd/p, this->nMoles());
+        return Kp(T)*pow(this->Pstd()/p, this->nMoles());
     }
 }
 
@@ -274,7 +274,7 @@ inline Foam::scalar Foam::specieThermo<thermo>::Kn
     }
     else
     {
-        return Kp(T)*pow(n*this->Pstd/p, this->nMoles());
+        return Kp(T)*pow(n*this->Pstd()/p, this->nMoles());
     }
 }
 
diff --git a/src/thermophysicalModels/specie/transport/const/constTransportI.H b/src/thermophysicalModels/specie/transport/const/constTransportI.H
index cd83016a7..6cebf5b4b 100644
--- a/src/thermophysicalModels/specie/transport/const/constTransportI.H
+++ b/src/thermophysicalModels/specie/transport/const/constTransportI.H
@@ -109,9 +109,9 @@ inline scalar constTransport<thermo>::alpha(const scalar T) const
 {
     scalar Cp_ = this->Cp(T);
 
-    scalar deltaT = T - specie::Tstd;
+    scalar deltaT = T - specie::Tstd();
     scalar CpBar =
-        (deltaT*(this->H(T) - this->H(specie::Tstd)) + Cp_)/(sqr(deltaT) + 1);
+        (deltaT*(this->H(T) - this->H(specie::Tstd())) + Cp_)/(sqr(deltaT) + 1);
 
     return Cp_*mu(T)*rPr/CpBar;
 }
diff --git a/src/thermophysicalModels/specie/transport/polynomial/polynomialTransportI.H b/src/thermophysicalModels/specie/transport/polynomial/polynomialTransportI.H
index 70b3ae158..c80748e43 100644
--- a/src/thermophysicalModels/specie/transport/polynomial/polynomialTransportI.H
+++ b/src/thermophysicalModels/specie/transport/polynomial/polynomialTransportI.H
@@ -116,9 +116,9 @@ inline Foam::scalar Foam::polynomialTransport<Thermo, PolySize>::alpha
     const scalar T
 ) const
 {
-    scalar deltaT = T - specie::Tstd;
+    scalar deltaT = T - specie::Tstd();
     scalar CpBar =
-        (deltaT*(this->H(T) - this->H(specie::Tstd)) + this->Cp(T))
+        (deltaT*(this->H(T) - this->H(specie::Tstd())) + this->Cp(T))
        /(sqr(deltaT) + 1);
 
     return kappa(T)/CpBar;
diff --git a/src/thermophysicalModels/specie/transport/sutherland/sutherlandTransportI.H b/src/thermophysicalModels/specie/transport/sutherland/sutherlandTransportI.H
index ea50e2ed6..14ae677f4 100644
--- a/src/thermophysicalModels/specie/transport/sutherland/sutherlandTransportI.H
+++ b/src/thermophysicalModels/specie/transport/sutherland/sutherlandTransportI.H
@@ -148,9 +148,9 @@ inline scalar sutherlandTransport<thermo>::alpha(const scalar T) const
     scalar R_ = this->R();
     scalar Cp_ = Cv_ + R_;
 
-    scalar deltaT = T - specie::Tstd;
+    scalar deltaT = T - specie::Tstd();
     scalar CpBar =
-        (deltaT*(this->H(T) - this->H(specie::Tstd)) + Cp_)/(sqr(deltaT) + 1);
+        (deltaT*(this->H(T) - this->H(specie::Tstd())) + Cp_)/(sqr(deltaT) + 1);
 
     return mu(T)*Cv_*(1.32 + 1.77*this->R()/Cv_)/CpBar;
 }