diff --git a/src/OpenFOAM/meshes/polyMesh/polyMeshCheck/polyMeshCheck.C b/src/OpenFOAM/meshes/polyMesh/polyMeshCheck/polyMeshCheck.C index 2cee8c6914d49597167134ad218c6b2dc170db56..5136e354ee8bdebe1ecbc9ad2a4e34d4e027b04c 100644 --- a/src/OpenFOAM/meshes/polyMesh/polyMeshCheck/polyMeshCheck.C +++ b/src/OpenFOAM/meshes/polyMesh/polyMeshCheck/polyMeshCheck.C @@ -123,13 +123,13 @@ bool Foam::polyMesh::checkFaceOrthogonality reduce(severeNonOrth, sumOp<label>()); reduce(errorNonOrth, sumOp<label>()); - const scalar maxNonOrth = radToDeg(::acos(clamp(minDDotS, -1, 1))); - const scalar aveNonOrth = radToDeg(::acos(clamp(sumDDotS/nSummed, -1, 1))); - dictionary& meshDict = const_cast<dictionary&>(data().meshDict()); if (nSummed > 0) { + scalar maxNonOrth = radToDeg(::acos(clamp(minDDotS, -1, 1))); + scalar aveNonOrth = radToDeg(::acos(clamp(sumDDotS/nSummed, -1, 1))); + meshDict.set("maxNonOrth", maxNonOrth); meshDict.set("aveNonOrth", aveNonOrth);