From 550a743b41774d93a7eae4abf51ffbcf9f80dde0 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs@hunt.opencfd.co.uk>
Date: Tue, 9 Jun 2009 07:52:25 +0100
Subject: [PATCH] extrapolate intersection distance

---
 src/autoMesh/autoHexMesh/meshRefinement/meshRefinement.C | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/autoMesh/autoHexMesh/meshRefinement/meshRefinement.C b/src/autoMesh/autoHexMesh/meshRefinement/meshRefinement.C
index 9ce1e3db76..eadc22ba3c 100644
--- a/src/autoMesh/autoHexMesh/meshRefinement/meshRefinement.C
+++ b/src/autoMesh/autoHexMesh/meshRefinement/meshRefinement.C
@@ -117,13 +117,13 @@ void Foam::meshRefinement::calcNeighbourData
                 // Extrapolate the face centre.
                 vector fn = faceAreas[i];
                 fn /= mag(fn)+VSMALL;
-                // Normal distance from face centre to cell centre
-                scalar d = ((faceCentres[i] - cellCentres[i]) & fn);
 
                 label own = faceCells[i];
                 label ownLevel = cellLevel[own];
                 label faceLevel = meshCutter_.getAnchorLevel(pp.start()+i);
 
+                // Normal distance from face centre to cell centre
+                scalar d = ((faceCentres[i] - cellCentres[own]) & fn);
                 if (faceLevel > ownLevel)
                 {
                     // Other cell more refined. Adjust normal distance