diff --git a/src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinementBaffles.C b/src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinementBaffles.C index 9798046ca3c9c5ab91b6c59ea7294fb458fb97c1..0d9cecd039d5c852a972049b07b1091571e18367 100644 --- a/src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinementBaffles.C +++ b/src/mesh/autoMesh/autoHexMesh/meshRefinement/meshRefinementBaffles.C @@ -1234,7 +1234,7 @@ void Foam::meshRefinement::findCellZoneTopo { label surfI = namedSurfaceIndex[faceI]; - if (surfI != -1 && surfaceToCellZone[surfI] != -1) + if (surfI != -1) { // Calculate region to zone from cellRegions on either side // of internal face. @@ -1286,7 +1286,7 @@ void Foam::meshRefinement::findCellZoneTopo label surfI = namedSurfaceIndex[faceI]; - if (surfI != -1 && surfaceToCellZone[surfI] != -1) + if (surfI != -1) { bool changedCell = calcRegionToZone ( @@ -2522,6 +2522,32 @@ Foam::autoPtr<Foam::mapPolyMesh> Foam::meshRefinement::zonify mesh_.setInstance(oldInstance()); } + // Print some stats (note: zones are synchronised) + if (mesh_.cellZones().size() > 0) + { + Info<< "CellZones:" << endl; + forAll(mesh_.cellZones(), zoneI) + { + const cellZone& cz = mesh_.cellZones()[zoneI]; + Info<< " " << cz.name() + << "\tsize:" << returnReduce(cz.size(), sumOp<label>()) + << endl; + } + Info<< endl; + } + if (mesh_.faceZones().size() > 0) + { + Info<< "FaceZones:" << endl; + forAll(mesh_.faceZones(), zoneI) + { + const faceZone& fz = mesh_.faceZones()[zoneI]; + Info<< " " << fz.name() + << "\tsize:" << returnReduce(fz.size(), sumOp<label>()) + << endl; + } + Info<< endl; + } + // None of the faces has changed, only the zones. Still... updateMesh(map, labelList());