From 8b52119b97f654be21dc8b5636572c587f3330f6 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Wed, 19 Feb 2014 11:22:54 +0000
Subject: [PATCH] BUG: smappyHexMesh: report correct layerFaces

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

diff --git a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C
index 8c10c44001..8259e9ac22 100644
--- a/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C
+++ b/src/mesh/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C
@@ -2463,8 +2463,14 @@ 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.
+            if (!mesh.isInternalFace(layer[0]))
+            {
+                FatalErrorIn("autoLayerDriver::getLayerCellsFaces()")
+                    << "problem." << layer[0] << exit(FatalError);
+            }
+
+            for (label i = 1; i < layer.size(); i++)
             {
                 faceRealThickness[layer[i]] = realThickness;
             }
@@ -2608,7 +2614,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 +2623,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 +3411,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)
                 {