diff --git a/libsrc/meshing/meshclass.cpp b/libsrc/meshing/meshclass.cpp
index 87017b04..892a5275 100644
--- a/libsrc/meshing/meshclass.cpp
+++ b/libsrc/meshing/meshclass.cpp
@@ -4099,6 +4099,8 @@ namespace netgen
         UpdateTopology();
         timestamp = NextTimeStamp();
       }
+
+    RebuildSurfaceElementLists();
   }
 
   void Mesh :: BuildElementSearchTree ()
@@ -5024,6 +5026,7 @@ namespace netgen
     static int timer = NgProfiler::CreateTimer ("GetSurfaceElementsOfFace");
     NgProfiler::RegionTimer reg (timer);
 
+     /*
      sei.SetSize (0);
      for (SurfaceElementIndex i = 0; i < GetNSE(); i++)
      {
@@ -5033,6 +5036,7 @@ namespace netgen
            sei.Append (i);
         }
      }
+     */
 
      /* Philippose - 01/10/2009
      Commented out the following lines, and activated the originally 
@@ -5041,7 +5045,6 @@ namespace netgen
      */
 
      //      int size1 = sei.Size();
-     /*
      sei.SetSize(0);
 
      SurfaceElementIndex si = facedecoding[facenr-1].firstelement;
@@ -5055,7 +5058,6 @@ namespace netgen
 
         si = (*this)[si].next;
      }
-    */
      
      /*
      // *testout << "with list = " << endl << sei << endl;