Index: DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilProfileHelper.cs =================================================================== diff -u -r5166 -r5239 --- DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilProfileHelper.cs (.../SoilProfileHelper.cs) (revision 5166) +++ DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilProfileHelper.cs (.../SoilProfileHelper.cs) (revision 5239) @@ -51,21 +51,19 @@ }; soilProfile.Layers.Insert(0, topLayer); } - // Todo #Bka: I believe this code should be here to ensure you do not take soil above surface line into account. - // Todo However, it causes a series of errors in the tests which I do not get. Have to look at that together with Tom - // else if (soilProfile.TopLevel > top) - // { - // int bottom = RemoveLayersAboveTop(originalSoilProfile1D.LayerCount, soilProfile, top); - // var topLayer = new SoilLayer1D - // { - // TopLevel = top, - // Soil = originalSoilProfile1D.Layers[bottom].Soil, - // SoilProfile = soilProfile, - // IsAquifer = originalSoilProfile1D.Layers[bottom].IsAquifer, - // BottomLevel = originalSoilProfile1D.Layers[bottom].BottomLevel - // }; - // soilProfile.Layers.Insert(0, topLayer); - // } + else if (soilProfile.TopLevel > top) + { + int bottom = RemoveLayersAboveTop(originalSoilProfile1D.LayerCount, soilProfile, top); + var topLayer = new SoilLayer1D + { + TopLevel = top, + Soil = originalSoilProfile1D.Layers[bottom].Soil, + SoilProfile = soilProfile, + IsAquifer = originalSoilProfile1D.Layers[bottom].IsAquifer, + BottomLevel = originalSoilProfile1D.Layers[bottom].BottomLevel + }; + soilProfile.Layers.Insert(0, topLayer); + } soilProfile.BottomLevel = originalSoilProfile1D.BottomLevel; return soilProfile;