Index: DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilProfileHelper.cs =================================================================== diff -u -r3240 -r3241 --- DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilProfileHelper.cs (.../SoilProfileHelper.cs) (revision 3240) +++ DamEngine/trunk/src/Deltares.DamEngine.Data/Geotechnics/SoilProfileHelper.cs (.../SoilProfileHelper.cs) (revision 3241) @@ -25,19 +25,21 @@ }; soilProfile.Layers.Insert(0, topLayer); } - 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); - } + // 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); + // } soilProfile.BottomLevel = originalSoilProfile1D.BottomLevel; return soilProfile;