diff --git a/applications/solvers/compressible/rhoCentralFoam/BCs/U/maxwellSlipUFvPatchVectorField.C b/applications/solvers/compressible/rhoCentralFoam/BCs/U/maxwellSlipUFvPatchVectorField.C
index 1d5c6713efb44abbde0b27dc7872632723b1e800..4962a338a19c13ef6a518aacbf4e776ac75cf481 100644
--- a/applications/solvers/compressible/rhoCentralFoam/BCs/U/maxwellSlipUFvPatchVectorField.C
+++ b/applications/solvers/compressible/rhoCentralFoam/BCs/U/maxwellSlipUFvPatchVectorField.C
@@ -28,8 +28,7 @@ License
 #include "mathematicalConstants.H"
 #include "fvPatchFieldMapper.H"
 #include "volFields.H"
-#include "surfaceFields.H"
-#include "fvCFD.H"
+#include "fvcGrad.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -129,7 +128,6 @@ maxwellSlipUFvPatchVectorField::maxwellSlipUFvPatchVectorField
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
-// Update the coefficients associated with the patch field
 void maxwellSlipUFvPatchVectorField::updateCoeffs()
 {
     if (updated())
@@ -177,7 +175,6 @@ void maxwellSlipUFvPatchVectorField::updateCoeffs()
 }
 
 
-// Write
 void maxwellSlipUFvPatchVectorField::write(Ostream& os) const
 {
     fvPatchVectorField::write(os);
diff --git a/applications/solvers/compressible/rhoCentralFoam/BCs/U/maxwellSlipUFvPatchVectorField.H b/applications/solvers/compressible/rhoCentralFoam/BCs/U/maxwellSlipUFvPatchVectorField.H
index 1e640891bfc88f89ed70bc8185401efb0698ab78..258e06f01ccd928aa6b2891f5ece629739b8764c 100644
--- a/applications/solvers/compressible/rhoCentralFoam/BCs/U/maxwellSlipUFvPatchVectorField.H
+++ b/applications/solvers/compressible/rhoCentralFoam/BCs/U/maxwellSlipUFvPatchVectorField.H
@@ -45,14 +45,13 @@ namespace Foam
 {
 
 /*---------------------------------------------------------------------------*\
-                Class maxwellSlipUFvPatch Declaration
+                     Class maxwellSlipUFvPatch Declaration
 \*---------------------------------------------------------------------------*/
 
 class maxwellSlipUFvPatchVectorField
 :
     public mixedFixedValueSlipFvPatchVectorField
 {
-
     // Private data
 
         // Accommodation coefficient
@@ -67,6 +66,7 @@ class maxwellSlipUFvPatchVectorField
         // Include boundary curvature term (default on)
         Switch curvature_;
 
+
 public:
 
     //- Runtime type information
@@ -128,6 +128,7 @@ public:
             );
         }
 
+
     // Member functions
 
         //- Update the coefficients associated with the patch field
diff --git a/applications/utilities/mesh/manipulation/checkMesh/checkTopology.C b/applications/utilities/mesh/manipulation/checkMesh/checkTopology.C
index 9e96eec03f776b8335c18d7b33e4d660a4887521..66bbaa23fc3b42eddb7f18eba462372647b6760d 100644
--- a/applications/utilities/mesh/manipulation/checkMesh/checkTopology.C
+++ b/applications/utilities/mesh/manipulation/checkMesh/checkTopology.C
@@ -338,20 +338,18 @@ Foam::label Foam::checkTopology
                 const pointField& pts = pp.points();
                 const labelList& mp = pp.meshPoints();
 
-                boundBox bb;   // zero-sized
                 if (returnReduce(mp.size(), sumOp<label>()) > 0)
                 {
-                    bb.min() = pts[mp[0]];
-                    bb.max() = pts[mp[0]];
-                    for (label i = 1; i < mp.size(); i++)
+                    boundBox bb(point::max, point::min);
+                    forAll (mp, i)
                     {
                         bb.min() = min(bb.min(), pts[mp[i]]);
                         bb.max() = max(bb.max(), pts[mp[i]]);
                     }
                     reduce(bb.min(), minOp<vector>());
                     reduce(bb.max(), maxOp<vector>());
+                    Pout<< ' ' << bb;
                 }
-                Pout<< ' ' << bb;
             }
             Pout<< endl;
         }
diff --git a/applications/utilities/mesh/manipulation/refineMesh/refineMesh.C b/applications/utilities/mesh/manipulation/refineMesh/refineMesh.C
index f890c845466467ee53a7f8e41fa3a74c212e2528..eb33c26427a7291f74cb1093d8e5bb31717f4a0b 100644
--- a/applications/utilities/mesh/manipulation/refineMesh/refineMesh.C
+++ b/applications/utilities/mesh/manipulation/refineMesh/refineMesh.C
@@ -296,6 +296,7 @@ int main(int argc, char *argv[])
     );
 
     #include "addOverwriteOption.H"
+    #include "addRegionOption.H"
     argList::addBoolOption
     (
         "dict",
@@ -305,7 +306,7 @@ int main(int argc, char *argv[])
     #include "setRootCase.H"
     #include "createTime.H"
     runTime.functionObjects().off();
-    #include "createPolyMesh.H"
+    #include "createNamedPolyMesh.H"
     const word oldInstance = mesh.pointsInstance();
 
     printEdgeStats(mesh);
diff --git a/tutorials/compressible/rhoCentralFoam/biconic25-55Run35/0/U b/tutorials/compressible/rhoCentralFoam/biconic25-55Run35/0/U
index d2c86a48344cb2f9a1de664755b3632bfe197b47..dbcd94b421b16ebbd6374ce08ed5578a915402eb 100644
--- a/tutorials/compressible/rhoCentralFoam/biconic25-55Run35/0/U
+++ b/tutorials/compressible/rhoCentralFoam/biconic25-55Run35/0/U
@@ -32416,7 +32416,7 @@ boundaryField
         Uwall           uniform (0 0 0);
         thermalCreep    on;
         curvature       on;
-        refValue        
+        refValue        nonuniform List<vector>
 255
 (
 (3.37115877925426 1.5719971556227 1.63068197702874e-10)
@@ -32676,7 +32676,7 @@ boundaryField
 (0.103666873607679 0 -6.64100738660639e-16)
 )
 ;
-        valueFraction   
+        valueFraction  nonuniform List<scalar> 
 255
 (
 0.00388739013936794