diff --git a/src/mesh/snappyHexMesh/externalDisplacementMeshMover/medialAxisMeshMover.C b/src/mesh/snappyHexMesh/externalDisplacementMeshMover/medialAxisMeshMover.C
index a70742cbf8341e403c90b3d839b0885acea19190..b5ce64a5f78a1331f7d5ffad8cc56d28e6825f4c 100644
--- a/src/mesh/snappyHexMesh/externalDisplacementMeshMover/medialAxisMeshMover.C
+++ b/src/mesh/snappyHexMesh/externalDisplacementMeshMover/medialAxisMeshMover.C
@@ -36,6 +36,7 @@ License
 #include "OBJstream.H"
 #include "PointData.H"
 #include "zeroFixedValuePointPatchFields.H"
+#include "pointSet.H"
 
 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
 
@@ -568,6 +569,33 @@ void Foam::medialAxisMeshMover::update(const dictionary& coeffDict)
         maxInfo.shrink();
         maxPoints.shrink();
 
+
+        if (debug)
+        {
+            mkDir(mesh().time().timePath());
+            OBJstream str(mesh().time().timePath()/"medialSurfacePoints.obj");
+
+            pointSet seedPoints
+            (
+                mesh(),
+                "medialSurfacePoints",
+                maxPoints
+            );
+
+            Info<< typeName
+                << " : Writing estimated medial surface:" << nl << incrIndent
+                << indent << "locations : " << str.name() << nl
+                << indent << "pointSet  : " << seedPoints.name() << nl
+                << decrIndent << endl;
+
+            for (const auto& info : maxInfo)
+            {
+                str.write(info.origin());
+            }
+            seedPoints.write();
+        }
+
+
         // Do all calculations
         PointEdgeWave<pointEdgePoint> medialDistCalc
         (
@@ -658,18 +686,16 @@ void Foam::medialAxisMeshMover::update(const dictionary& coeffDict)
     if (debug)
     {
         Info<< typeName
-            << " : Writing medial axis fields:" << nl
-            << incrIndent
-            << "ratio of medial distance to wall distance : "
+            << " : Writing medial axis fields:" << nl << incrIndent
+            << indent << "ratio of medial distance to wall distance : "
             << medialRatio_.name() << nl
-            << "distance to nearest medial axis           : "
+            << indent << "distance to nearest medial axis           : "
             << medialDist_.name() << nl
-            << "nearest medial axis location              : "
+            << indent << "nearest medial axis location              : "
             << medialVec_.name() << nl
-            << "normal at nearest wall                    : "
+            << indent << "normal at nearest wall                    : "
             << dispVec_.name() << nl
-            << decrIndent << nl
-            << endl;
+            << decrIndent << endl;
 
         dispVec_.write();
         medialRatio_.write();
diff --git a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C
index 440819950b2243dc9c6a8f999cefd91ff1008675..eb0f8d5f35096ab282b4f9a82f8f9c133d0e350f 100644
--- a/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C
+++ b/src/mesh/snappyHexMesh/meshRefinement/meshRefinementBaffles.C
@@ -4194,7 +4194,7 @@ Foam::autoPtr<Foam::mapPolyMesh> Foam::meshRefinement::removeLimitShells
     if (debug&meshRefinement::MESH)
     {
         const_cast<Time&>(mesh_.time())++;
-        Pout<< "Writing mesh after removing limitShells "
+        Pout<< "Writing mesh after removing limitShells"
             << " to time " << timeName() << endl;
         write
         (