surfaceFind: output region/zone number of found face

applications/utilities/surface/surfaceFind/surfaceFind.C

@@ -84,10 +84,10 @@ int main(int argc, char *argv[])
     }
 
     Info<< "Nearest vertex:" << nl
-        << "    index      :" << minIndex << " (in localPoints)" << nl
-        << "    index      :" << surf1.meshPoints()[minIndex]
+        << "    index      : " << minIndex << " (in localPoints)" << nl
+        << "    index      : " << surf1.meshPoints()[minIndex]
         << " (in points)" << nl
-        << "    coordinates:" << localPoints[minIndex] << nl
+        << "    coordinates: " << localPoints[minIndex] << nl
         << endl;
 
     //
@@ -114,14 +114,26 @@ int main(int argc, char *argv[])
     const face& f = surf1[minIndex];
 
     Info<< "Face with nearest centre:" << nl
-        << "    index        :" << minIndex << nl
-        << "    centre       :" << f.centre(points) << nl
-        << "    face         :" << f << nl
-        << "    vertex coords:\n";
+        << "    index        : " << minIndex << nl
+        << "    centre       : " << f.centre(points) << nl
+        << "    face         : " << f << nl
+        << "    vertex coords:" << endl;
+
     forAll(f, fp)
     {
-        Info<< "        " << points[f[fp]] << "\n";
+        Info<< "        " << points[f[fp]] << nl;
+    }
+
+    const List<surfZone>& surfZones = surf1.surfZones();
+    label surfZone = -1;
+    forAll(surfZones, zonei)
+    {
+        if (minIndex >= surfZones[zonei].start())
+        {
+            surfZone = zonei;
+        }
     }
+    Info<< "    zone/region  : " << surfZone << endl;
 
     Info<< nl
         << "End\n" << endl;