Index: DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilProfileHelper.cs =================================================================== diff -u -r4759 -r5163 --- DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilProfileHelper.cs (.../SoilProfileHelper.cs) (revision 4759) +++ DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilProfileHelper.cs (.../SoilProfileHelper.cs) (revision 5163) @@ -53,19 +53,19 @@ } // 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) - // { - // var 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;