zhyang-dev/OpenFOAM-5.x
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

surfaceFind: output region/zone number of found face

Browse Source
This commit is contained in:
Chris Greenshields
2017-06-02 21:02:01 +01:00
parent 97585eac40
commit f1240a853c
1 changed files with 19 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

27
applications/utilities/surface/surfaceFind/surfaceFind.C
View File

@ -81,10 +81,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;
//
@ -111,15 +111,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<< endl;
Info<< "End\n" << endl;