diff --git a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom/system/controlDict b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom/system/controlDict
index cc11f83e1..0d4fed148 100644
--- a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom/system/controlDict
+++ b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom/system/controlDict
@@ -21,7 +21,7 @@ libs
 );
 application     buoyantBoussinesqSimpleFoam;
 
-startFrom       latestTime;
+startFrom       startTime;
 
 startTime       0;
 
@@ -41,7 +41,7 @@ writeFormat     ascii;
 
 writePrecision  6;
 
-writeCompression uncompressed;
+writeCompression compressed;
 
 timeFormat      general;
 
@@ -49,5 +49,54 @@ timePrecision   6;
 
 runTimeModifiable yes;
 
+functions
+(
+    minMaxU
+    {
+        type        minMaxField;
+
+        // Where to load it from (if not already in solver)
+        functionObjectLibs ("libfieldFunctionObjects.so");
+
+        name U;
+    }
+    minMaxP
+    {
+        type        minMaxField;
+
+        // Where to load it from (if not already in solver)
+        functionObjectLibs ("libfieldFunctionObjects.so");
+
+        name p;
+    }
+    minMaxT
+    {
+        type        minMaxField;
+
+        // Where to load it from (if not already in solver)
+        functionObjectLibs ("libfieldFunctionObjects.so");
+
+        name T;
+    }
+    minMaxK
+    {
+        type        minMaxField;
+
+        // Where to load it from (if not already in solver)
+        functionObjectLibs ("libfieldFunctionObjects.so");
+
+        name k;
+    }
+    minMaxEpsilon
+    {
+        type        minMaxField;
+
+        // Where to load it from (if not already in solver)
+        functionObjectLibs ("libfieldFunctionObjects.so");
+
+        name epsilon;
+    }
+);
+
 
 // ************************************************************************* //
diff --git a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom/system/fvSchemes b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom/system/fvSchemes
index 8b0408600..159eeca39 100644
--- a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom/system/fvSchemes
+++ b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom/system/fvSchemes
@@ -32,20 +32,12 @@ divSchemes
     div(phi,T)      Gauss upwind;
     div(phi,k)      Gauss upwind;
     div(phi,epsilon) Gauss upwind;
-    div(phi,R)      Gauss upwind;
-    div(R)          Gauss linear;
     div((nuEff*dev(T(grad(U))))) Gauss linear;
 }
 
 laplacianSchemes
 {
-    default         none;
-    laplacian(nuEff,U) Gauss linear corrected;
-    laplacian((1|A(U)),p) Gauss linear corrected;
-    laplacian(kappaEff,T) Gauss linear corrected;
-    laplacian(DkEff,k) Gauss linear corrected;
-    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
-    laplacian(DREff,R) Gauss linear corrected;
+    default         Gauss linear corrected;
 }
 
 interpolationSchemes
diff --git a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom/system/fvSolution b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom/system/fvSolution
index 5cd951b30..91fb0bf44 100644
--- a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom/system/fvSolution
+++ b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/hotRoom/system/fvSolution
@@ -24,48 +24,37 @@ solvers
         tolerance       1e-08;
         relTol          0.01;
     }
-
     U
     {
         solver          BiCGStab;
         preconditioner  DILU;
         tolerance       1e-05;
-        relTol          0.1;
+        relTol          0;
     }
-
     T
     {
         solver          BiCGStab;
         preconditioner  DILU;
         tolerance       1e-05;
-        relTol          0.1;
+        relTol          0;
     }
-
     k
     {
         solver          BiCGStab;
         preconditioner  DILU;
         tolerance       1e-05;
-        relTol          0.1;
+        relTol          0;
     }
-
     epsilon
     {
         solver          BiCGStab;
         preconditioner  DILU;
         tolerance       1e-05;
-        relTol          0.1;
-    }
-
-    R
-    {
-        solver          BiCGStab;
-        preconditioner  DILU;
-        tolerance       1e-05;
-        relTol          0.1;
+        relTol          0;
     }
 }
 
+
 SIMPLE
 {
     nNonOrthogonalCorrectors 0;
@@ -81,7 +70,6 @@ relaxationFactors
     T               0.7;
     k               0.7;
     epsilon         0.7;
-    R               0.7;
 }
 
 
diff --git a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/constant/polyMesh/boundary b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/constant/polyMesh/boundary
index e69a9c087..9dea11824 100644
--- a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/constant/polyMesh/boundary
+++ b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/constant/polyMesh/boundary
@@ -21,55 +21,55 @@ FoamFile
     {
         type            empty;
         nFaces          0;
-        startFace       59065;
+        startFace       60436;
     }
     minX
     {
         type            empty;
         nFaces          0;
-        startFace       59065;
+        startFace       60436;
     }
     maxX
     {
         type            empty;
         nFaces          0;
-        startFace       59065;
+        startFace       60436;
     }
     minY
     {
         type            empty;
         nFaces          0;
-        startFace       59065;
+        startFace       60436;
     }
     ground
     {
         type            wall;
         nFaces          590;
-        startFace       59065;
+        startFace       60436;
     }
     maxZ
     {
         type            empty;
         nFaces          0;
-        startFace       59655;
+        startFace       61026;
     }
     igloo_region0
     {
         type            wall;
         nFaces          2260;
-        startFace       59655;
+        startFace       61026;
     }
     twoFridgeFreezers_seal_0
     {
         type            wall;
         nFaces          1344;
-        startFace       61915;
+        startFace       63286;
     }
     twoFridgeFreezers_herring_1
     {
         type            wall;
         nFaces          1116;
-        startFace       63259;
+        startFace       64630;
     }
 )
 
diff --git a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/system/fvSchemes b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/system/fvSchemes
index 8b0408600..159eeca39 100644
--- a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/system/fvSchemes
+++ b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/system/fvSchemes
@@ -32,20 +32,12 @@ divSchemes
     div(phi,T)      Gauss upwind;
     div(phi,k)      Gauss upwind;
     div(phi,epsilon) Gauss upwind;
-    div(phi,R)      Gauss upwind;
-    div(R)          Gauss linear;
     div((nuEff*dev(T(grad(U))))) Gauss linear;
 }
 
 laplacianSchemes
 {
-    default         none;
-    laplacian(nuEff,U) Gauss linear corrected;
-    laplacian((1|A(U)),p) Gauss linear corrected;
-    laplacian(kappaEff,T) Gauss linear corrected;
-    laplacian(DkEff,k) Gauss linear corrected;
-    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
-    laplacian(DREff,R) Gauss linear corrected;
+    default         Gauss linear corrected;
 }
 
 interpolationSchemes
diff --git a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/system/fvSolution b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/system/fvSolution
index 4ac68f9fe..0e79914be 100644
--- a/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/system/fvSolution
+++ b/tutorials/heatTransfer/buoyantBoussinesqSimpleFoam/iglooWithFridges/system/fvSolution
@@ -24,48 +24,37 @@ solvers
         tolerance       1e-08;
         relTol          0.01;
     }
-
     U
     {
         solver          BiCGStab;
         preconditioner  DILU;
         tolerance       1e-05;
-        relTol          0.1;
+        relTol          0;
     }
-
     T
     {
         solver          BiCGStab;
         preconditioner  DILU;
         tolerance       1e-05;
-        relTol          0.1;
+        relTol          0;
     }
-
     k
     {
         solver          BiCGStab;
         preconditioner  DILU;
         tolerance       1e-05;
-        relTol          0.1;
+        relTol          0;
     }
-
     epsilon
     {
         solver          BiCGStab;
         preconditioner  DILU;
         tolerance       1e-05;
-        relTol          0.1;
-    }
-
-    R
-    {
-        solver          BiCGStab;
-        preconditioner  DILU;
-        tolerance       1e-05;
-        relTol          0.1;
+        relTol          0;
     }
 }
 
+
 SIMPLE
 {
     nNonOrthogonalCorrectors 0;
@@ -81,7 +70,6 @@ relaxationFactors
     T               0.7;
     k               0.7;
     epsilon         0.7;
-    R               0.7;
 }