From 084e5aad19f1002cc3e191b893045edf0a9c0673 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Tue, 21 Jul 2009 13:29:39 +0100
Subject: [PATCH] rewriting empty boundary conditions

---
 .../isoSurface/isoSurfaceTemplates.C            | 17 ++++++++++++-----
 1 file changed, 12 insertions(+), 5 deletions(-)

diff --git a/src/sampling/sampledSurface/isoSurface/isoSurfaceTemplates.C b/src/sampling/sampledSurface/isoSurface/isoSurfaceTemplates.C
index d7903551c9..bddf6fedf2 100644
--- a/src/sampling/sampledSurface/isoSurface/isoSurfaceTemplates.C
+++ b/src/sampling/sampledSurface/isoSurface/isoSurfaceTemplates.C
@@ -96,11 +96,18 @@ Foam::isoSurface::adaptPatchFields
                     sliceFld
                 )
             );
-            sliceFld.boundaryField()[patchI] ==
-                mesh.boundary()[patchI].patchInternalField
-                (
-                    sliceFld
-                );
+
+            // Note: cannot use patchInternalField since uses emptyFvPatch::size
+            // Do our own internalField instead.
+            const unallocLabelList& faceCells =
+                mesh.boundary()[patchI].patch().faceCells();
+
+            Field<Type>& pfld = sliceFld.boundaryField()[patchI];
+            pfld.setSize(faceCells.size());
+            forAll(faceCells, i)
+            {
+                pfld[i] = sliceFld[faceCells[i]];
+            }
         }
         else if (isA<cyclicPolyPatch>(pp))
         {