From 60fa49504fe95837a63ed318a31ab43168c86e75 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Fri, 22 Mar 2019 13:31:02 +0000
Subject: [PATCH] Calculate Co number only on large live faces

---
 .../immersedBoundary/include/immersedBoundaryCourantNo.H        | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/immersedBoundary/immersedBoundary/include/immersedBoundaryCourantNo.H b/src/immersedBoundary/immersedBoundary/include/immersedBoundaryCourantNo.H
index 9608db2c8..41707b63a 100644
--- a/src/immersedBoundary/immersedBoundary/include/immersedBoundaryCourantNo.H
+++ b/src/immersedBoundary/immersedBoundary/include/immersedBoundaryCourantNo.H
@@ -36,7 +36,7 @@ scalar velMag = 0.0;
 
 if (mesh.nInternalFaces())
 {
-    surfaceScalarField magPhi = mag(sGamma*phi);
+    surfaceScalarField magPhi = pos(sGamma - 0.5)*mag(phi);
 
     surfaceScalarField SfUfbyDelta =
         mesh.surfaceInterpolation::deltaCoeffs()*magPhi;