Index: DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilSurfaceProfile.cs
===================================================================
diff -u -r3268 -r3269
--- DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilSurfaceProfile.cs	(.../SoilSurfaceProfile.cs)	(revision 3268)
+++ DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilSurfaceProfile.cs	(.../SoilSurfaceProfile.cs)	(revision 3269)
@@ -151,9 +151,10 @@
 
             var localSurfaceLinePoints = GetLocalSurfaceLinePoints(filteredSoilLayers, surfaceLine2);
             var layerData = CreateLayerData(filteredSoilLayers, localSurfaceLinePoints);
+            var filteredLayerData = FilterLayerDataWithClosedLoops(layerData);
 
-            BuildGeometryModel(Geometry, layerData, originalSoilProfile1D, surfaceLine2, localSurfaceLinePoints);
-            AssignSoilLayer2DProperties(Surfaces, layerData, Geometry);
+            BuildGeometryModel(Geometry, filteredLayerData, originalSoilProfile1D, surfaceLine2, localSurfaceLinePoints);
+            AssignSoilLayer2DProperties(Surfaces, filteredLayerData, Geometry);
         }
 
         private static void BuildGeometryModel(GeometryData geometryData, 
@@ -333,12 +334,22 @@
 
             return layerData;
         }
+        
+        /// <summary>
+        /// Filters layer data that consist of closed loops.
+        /// </summary>
+        /// <param name="layerData">The collection of <see cref="LayerData"/> to filter.</param>
+        /// <returns>A collection of <see cref="LayerData"/> that define a closed loop.</returns>
+        private static IEnumerable<LayerData> FilterLayerDataWithClosedLoops(IEnumerable<LayerData> layerData)
+        {
+            return layerData.Where(l => l.Curves.Count() > 2);
+        }
 
         /// <summary>
         /// Creates a 2D soil layer with points at <paramref name="surfaceLinePoints"/>.
         /// </summary>
         /// <param name="soilLayer">The <see cref="SoilLayer1D"/> to create the 2D soil layer for.</param>
-        /// <param name="surfaceLinePoints">The collection of x coordinates to define the 2D soil layer for.</param>
+        /// <param name="surfaceLinePoints">The collection of <see cref="Point2D"/> that defines the surface line.</param>
         /// <returns>A <see cref="LayerData"/> based on the input arguments.</returns>
         private static LayerData CreateSurfaceLineWideSoilLayer(SoilLayer1D soilLayer, IEnumerable<Point2D> surfaceLinePoints)
         {
@@ -358,9 +369,8 @@
         /// Creates a 2D soil layer with points at <paramref name="surfaceLinePoints"/> based on its input arguments.
         /// </summary>
         /// <param name="soilLayer">The <see cref="SoilLayer1D"/> to create the 2D soil layer for.</param>
-        /// <param name="surfaceLinePoints">The collection of x coordinates to define the 2D soil layer for.</param>
+        /// <param name="surfaceLinePoints">The collection of <see cref="Point2D"/> that defines the surface line.</param>
         /// <returns>A collection of <see cref="LayerData"/> based on the input arguments.</returns>
-        /// <remarks>The collection only contains <see cref="LayerData"/> with non-empty coordinates.</remarks>
         private static IEnumerable<LayerData> CreateSoilLayerWithSurfaceLineIntersection(SoilLayer1D soilLayer, IEnumerable<Point2D> surfaceLinePoints)
         {
             // Surface cannot be determined if there are not at least 2 coordinates
@@ -383,9 +393,9 @@
             {
                 Point2D currentPoint = pointsArray[i];
                 double xCoordinate = currentPoint.X;
-
-                double nextSurfaceLineZCoordinate = pointsArray[j].Z;
                 double zCoordinateSurfaceLine = currentPoint.Z;
+                double nextSurfaceLineZCoordinate = pointsArray[j].Z;
+
                 double currentZCoordinate = zCoordinateSurfaceLine;
                 if (zCoordinateSurfaceLine > layerTopLevel)
                 {
@@ -442,7 +452,7 @@
             var pointsToConvert = new List<Point2D>(area.Coordinates);
 
             // Determine the additional points need to be added to close the area
-            // The area ends at the top right point. Go clockwise to close the area if necessary by adding 
+            // The area ends at the top right point. Go clockwise to close the area by adding 
             // points corresponding to the layer bottom level.
             var pointsToTraverse = area.Coordinates.Reverse();
             foreach (Point2D point in pointsToTraverse)
@@ -463,7 +473,7 @@
         /// </summary>
         /// <param name="pointsToConvert">The collection points to generate a curve collection for.</param>
         /// <returns>A collection of <see cref="GeometryCurve"/>.</returns>
-        private static IEnumerable<GeometryCurve> GenerateCurves(List<Point2D> pointsToConvert)
+        private static IEnumerable<GeometryCurve> GenerateCurves(IReadOnlyList<Point2D> pointsToConvert)
         {
             Point2D startPoint = pointsToConvert[0];
             var curves = new List<GeometryCurve>();
@@ -482,7 +492,7 @@
         /// </summary>
         /// <param name="pointsToConvert">The collection points to generate a curve collection for.</param>
         /// <returns>A collection of <see cref="GeometryCurve"/>.</returns>
-        private static IEnumerable<GeometryCurve> GenerateClosedLoopCurves(List<Point2D> pointsToConvert)
+        private static IEnumerable<GeometryCurve> GenerateClosedLoopCurves(IReadOnlyList<Point2D> pointsToConvert)
         {
             var curves = GenerateCurves(pointsToConvert).ToList();