ENH: checkMesh: write wall distance. See #71

applications/utilities/mesh/manipulation/checkMesh/checkMesh.C

@@ -160,7 +160,8 @@ int main(int argc, char *argv[])
             "cellVolume",
             "cellVolumeRatio",
             "minTetVolume",
-            "cellRegion"
+            "cellRegion",
+            "wallDistance"
         });
     }
 

applications/utilities/mesh/manipulation/checkMesh/writeFields.C

@@ -6,6 +6,7 @@
 #include "syncTools.H"
 #include "tetPointRef.H"
 #include "regionSplit.H"
+#include "wallDist.H"
 
 using namespace Foam;
 
@@ -281,6 +282,7 @@ void Foam::writeFields
 
     // Aspect ratio
     // ~~~~~~~~~~~~
+
     if (selectedFields.found("aspectRatio"))
     {
         volScalarField aspectRatio
@@ -498,5 +500,19 @@ void Foam::writeFields
         cellRegion.write();
     }
 
+    if (selectedFields.found("wallDistance"))
+    {
+        // Wall distance
+        //const volScalarField& y = wallDist::New(mesh).y();
+        volScalarField y("wallDistance", wallDist::New(mesh).y());
+        Info<< "    Writing wall distance to " << y.name() << endl;
+        y.write();
+
+        // Wall-reflection vectors
+        //const volVectorField& n = wallDist::New(mesh).n();
+        //Info<< "    Writing wall normal to " << n.name() << endl;
+        //n.write();
+    }
+
     Info<< endl;
 }