diff --git a/applications/test/volPointInterpolation/Test-volPointInterpolation.C b/applications/test/volPointInterpolation/Test-volPointInterpolation.C
index 1bbcd6da5e484edf333b620098f1bccc4362144d..96ffeeac377449a728a529de0f7c92180592cb22 100644
--- a/applications/test/volPointInterpolation/Test-volPointInterpolation.C
+++ b/applications/test/volPointInterpolation/Test-volPointInterpolation.C
@@ -33,11 +33,10 @@ Application
 
 int main(int argc, char *argv[])
 {
+    #include "setRootCase.H"
 
-#   include "setRootCase.H"
-
-#   include "createTime.H"
-#   include "createMesh.H"
+    #include "createTime.H"
+    #include "createMesh.H"
 
     Info<< "Reading field p\n" << endl;
     volScalarField p
@@ -67,9 +66,29 @@ int main(int argc, char *argv[])
         mesh
     );
 
+    const pointMesh& pMesh = pointMesh::New(mesh);
+    const pointBoundaryMesh& pbm = pMesh.boundary();
+
+    Info<< "pointMesh boundary" << nl;
+    forAll(pbm, patchI)
+    {
+        Info<< "patch=" << pbm[patchI].name()
+            << ", type=" << pbm[patchI].type()
+            << ", coupled=" << pbm[patchI].coupled()
+            << endl;
+    }
+
     const volPointInterpolation& pInterp = volPointInterpolation::New(mesh);
 
+
     pointScalarField pp(pInterp.interpolate(p));
+    Info<< pp.name() << " boundary" << endl;
+    forAll(pp.boundaryField(), patchI)
+    {
+        Info<< pbm[patchI].name() << " coupled="
+            << pp.boundaryField()[patchI].coupled()<< endl;
+    }
+
     pp.write();
 
     pointVectorField pU(pInterp.interpolate(U));