diff --git a/src/immersedBoundary/immersedBoundary/immersedPoly/ImmersedFace.C b/src/immersedBoundary/immersedBoundary/immersedPoly/ImmersedFace.C index 1cb1a01cb..c0c997fab 100644 --- a/src/immersedBoundary/immersedBoundary/immersedPoly/ImmersedFace.C +++ b/src/immersedBoundary/immersedBoundary/immersedPoly/ImmersedFace.C @@ -272,19 +272,19 @@ void Foam::ImmersedFace::createSubfaces { drySubface_.setSize(nDry); - // // Check area: if it is very small, reset the face - // if - // ( - // drySubface_.mag(facePointsAndIntersections_) - // < immersedPoly::liveFactor_()* - // localFace.mag(facePointsAndIntersections_) - // ) - // { - // // The face is practically wet - // isAllWet_ = true; + // Check area: if it is very small, reset the face + if + ( + drySubface_.mag(facePointsAndIntersections_) + < immersedPoly::liveFactor_()* + localFace.mag(facePointsAndIntersections_) + ) + { + // The face is practically wet + isAllWet_ = true; - // drySubface_.clear(); - // } + drySubface_.clear(); + } } if (nWet < 3) @@ -299,21 +299,21 @@ void Foam::ImmersedFace::createSubfaces { wetSubface_.setSize(nWet); - // // Check area: if it is very small, reset the face - // // Note: must use relative tolerance, ie divide with original face - // // area magnitude. HJ, 30/Aug/2018 - // if - // ( - // wetSubface_.mag(facePointsAndIntersections_) - // < immersedPoly::liveFactor_()* - // localFace.mag(facePointsAndIntersections_) - // ) - // { - // // The face is practically dry - // isAllDry_ = true; + // Check area: if it is very small, reset the face + // Note: must use relative tolerance, ie divide with original face + // area magnitude. HJ, 30/Aug/2018 + if + ( + wetSubface_.mag(facePointsAndIntersections_) + < immersedPoly::liveFactor_()* + localFace.mag(facePointsAndIntersections_) + ) + { + // The face is practically dry + isAllDry_ = true; - // wetSubface_.clear(); - // } + wetSubface_.clear(); + } } } }