From 1eb9ab123843ffd24dd27a32946eb93a06f18aa1 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Wed, 19 Feb 2014 15:00:51 +0000
Subject: [PATCH] BUG: snappyHexMesh: put internal faces only in layerFaces

---
 .../autoHexMesh/autoHexMeshDriver/autoLayerDriver.C    | 10 +++++-----
 1 file changed, 5 insertions(+), 5 deletions(-)

diff --git a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C
index 8c10c44001..08c8f7e530 100644
--- a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C
+++ b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C
@@ -2463,8 +2463,8 @@ void Foam::autoLayerDriver::getLayerCellsFaces
         if (layer.size())
         {
             // Layer contains both original boundary face and new boundary
-            // face so is nLayers+1
-            forAll(layer, i)
+            // face so is nLayers+1. Leave out old internal face.
+            for (label i = 1; i < layer.size(); i++)
             {
                 faceRealThickness[layer[i]] = realThickness;
             }
@@ -2608,7 +2608,7 @@ bool Foam::autoLayerDriver::writeLayerData
         }
         {
             label nAdded = 0;
-            forAll(faceRealThickness, faceI)
+            for (label faceI = 0; faceI < mesh.nInternalFaces(); faceI++)
             {
                 if (faceRealThickness[faceI] > 0)
                 {
@@ -2617,7 +2617,7 @@ bool Foam::autoLayerDriver::writeLayerData
             }
 
             faceSet layerFacesSet(mesh, "layerFaces", nAdded);
-            forAll(faceRealThickness, faceI)
+            for (label faceI = 0; faceI < mesh.nInternalFaces(); faceI++)
             {
                 if (faceRealThickness[faceI] > 0)
                 {
@@ -3405,7 +3405,7 @@ void Foam::autoLayerDriver::addLayers
             addedCellSet.write();
 
             faceSet layerFacesSet(newMesh, "layerFaces", newMesh.nFaces()/100);
-            forAll(faceRealThickness, faceI)
+            for (label faceI = 0; faceI < newMesh.nInternalFaces(); faceI++)
             {
                 if (faceRealThickness[faceI] > 0)
                 {