Index: DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilSurfaceProfile.cs
===================================================================
diff -u -r3260 -r3261
--- DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilSurfaceProfile.cs	(.../SoilSurfaceProfile.cs)	(revision 3260)
+++ DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilSurfaceProfile.cs	(.../SoilSurfaceProfile.cs)	(revision 3261)
@@ -142,7 +142,8 @@
         }
 
         /// <summary>
-        /// Updates the layers.
+        /// Function to create and update the <see cref="Geometry"/>, <see cref="SoilProfile2D.Surfaces"/>
+        /// based on a 1D profile and a surface line.
         /// </summary>
         private void Create2DGeometryBasedOn1DProfileAndSurfaceLine()
         {
@@ -152,11 +153,10 @@
             var layerData = CreateLayerData(filteredSoilLayers, surfaceLine2, xCoordinatesToTraverse);
 
             BuildGeometryModel(Geometry, layerData, originalSoilProfile1D, surfaceLine2);
-            BuildSoilLayer2D(Surfaces, layerData, Geometry);
+            AssignSoilLayer2DProperties(Surfaces, layerData, Geometry);
         }
 
-
-        private void BuildGeometryModel(GeometryData geometryData, 
+        private static void BuildGeometryModel(GeometryData geometryData, 
             IEnumerable<LayerData> layers, 
             SoilProfile1D profile, 
             SurfaceLine2 surfaceLine)
@@ -177,7 +177,7 @@
             geometryData.Clear();
             geometryData.Left = minX;
             geometryData.Right = maxX;
-            geometryData.Bottom = Math.Min(profile.BottomLevel, surfaceLine2.Geometry.GetMinZ() - 1);
+            geometryData.Bottom = Math.Min(profile.BottomLevel, surfaceLine.Geometry.GetMinZ() - 1);
 
             var pointsToBeAdded = layers.SelectMany(l => l.Points);
             var curvesToBeAdded = layers.SelectMany(l => l.Curves);
@@ -186,21 +186,23 @@
             geometryData.RegenerateGeometry();
         }
 
-        private void BuildSoilLayer2D(ICollection<SoilLayer2D> soilLayer2Ds, 
+        private static void AssignSoilLayer2DProperties(ICollection<SoilLayer2D> surfaces, 
             IEnumerable<LayerData> layers,
             GeometryData geometryData)
         {
-            soilLayer2Ds.Clear();
+            surfaces.Clear();
 
             foreach (GeometrySurface surface in geometryData.Surfaces)
             {
                 GeometryBounds bounds = surface.GetGeometryBounds();
                 double z = (bounds.Top + bounds.Bottom) * 0.5;
 
+                // Get the corresponding layer data. This assumption only holds
+                // because each z coordinate contains ONE single layer definition.
                 LayerData layerData = layers.FirstOrDefault(l => l.IsZLocatedInLayer(z));
                 if (layerData != null)
                 {
-                    Surfaces.Add(new SoilLayer2D
+                    surfaces.Add(new SoilLayer2D
                     {
                         GeometrySurface = surface,
                         IsAquifer = layerData.IsAquifer,
@@ -322,6 +324,12 @@
             return layerData;
         }
 
+        /// <summary>
+        /// Creates a 2D soil layer with points at <paramref name="xCoordinates"/>.
+        /// </summary>
+        /// <param name="soilLayer">The <see cref="SoilLayer1D"/> to create the 2D soil layer for.</param>
+        /// <param name="xCoordinates">The collection of x coordinates to define the 2D soil layer for.</param>
+        /// <returns>A <see cref="LayerData"/> based on the input arguments.</returns>
         private static LayerData CreateSurfaceLineWideSoilLayer(SoilLayer1D soilLayer, IEnumerable<double> xCoordinates)
         {
             // A single square definition does not suffice for the area definition. The 
@@ -336,6 +344,14 @@
             return new LayerData(pointsToConvert, curves, soilLayer.Soil, soilLayer.IsAquifer, soilLayer.WaterpressureInterpolationModel);
         }
 
+        /// <summary>
+        /// Creates a 2D soil layer with points at <paramref name="xCoordinates"/> based on its input arguments.
+        /// </summary>
+        /// <param name="soilLayer">The <see cref="SoilLayer1D"/> to create the 2D soil layer for.</param>
+        /// <param name="surfaceLine">The surface line to create the 2D soil layer with.</param>
+        /// <param name="xCoordinates">The collection of x coordinates to define the 2D soil layer for.</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, 
             SurfaceLine2 surfaceLine, IEnumerable<double> xCoordinates)
         {
@@ -374,7 +390,7 @@
 
                 // Determine if a new area should be started. This should only happen when:
                 if (previousZCoordinate > layerBottomLevel              // The area already started above or at the soil layer bottom
-                    && zCoordinateSurfaceLine < previousZCoordinate     // The surface line has a decreasing trend wrt the previous z coordinate
+                    && zCoordinateSurfaceLine < previousZCoordinate     // The surface line has a decreasing trend w.r.t. the previous z coordinate
                     && zCoordinateSurfaceLine <= layerBottomLevel)      // The surface line crossed the bottom of the layer
                 {
                     enclosedAreas.Add(currentArea);