Index: DamEngine/trunk/src/Deltares.DamEngine.Calculators/KernelWrappers/DamPipingSellmeijerVnk/DamPipingSellmeijerVnkKernelWrapper.cs =================================================================== diff -u -r902 -r906 --- DamEngine/trunk/src/Deltares.DamEngine.Calculators/KernelWrappers/DamPipingSellmeijerVnk/DamPipingSellmeijerVnkKernelWrapper.cs (.../DamPipingSellmeijerVnkKernelWrapper.cs) (revision 902) +++ DamEngine/trunk/src/Deltares.DamEngine.Calculators/KernelWrappers/DamPipingSellmeijerVnk/DamPipingSellmeijerVnkKernelWrapper.cs (.../DamPipingSellmeijerVnkKernelWrapper.cs) (revision 906) @@ -133,8 +133,6 @@ var xExit = 0.0; var surfaceLevel = 0.0; var dCoverLayer = 0.0; - var d70 = 0.0; - var aquiferHeight = 0.0; var permeabilityKx = 0.0; double? upliftFactor = null; if (upliftLocationAndResult != null) @@ -159,12 +157,18 @@ var referenceLevel = Math.Max(location.PolderLevel, surfaceLevel); Soil inBetweenAquiferlayerSoil = soilProfile1D.BottomAquiferLayer.Soil; - var inBetweenAquiferlayerHeight = 0.0; + double inBetweenAquiferlayerHeight; + var bottomAquiferlayerHeight = soilProfile1D.GetLayerHeight(soilProfile1D.BottomAquiferLayer); if (soilProfile1D.InBetweenAquiferLayer != null) { inBetweenAquiferlayerSoil = soilProfile1D.InBetweenAquiferLayer.Soil; inBetweenAquiferlayerHeight = soilProfile1D.GetLayerHeight(soilProfile1D.InBetweenAquiferLayer); } + else + { + bottomAquiferlayerHeight = bottomAquiferlayerHeight / 2; + inBetweenAquiferlayerHeight = bottomAquiferlayerHeight; + } kernelDataInput = new DamPipingSellmeijerVnkInput() { @@ -181,7 +185,7 @@ WaterViscosity = Physics.WaterViscosity, SurfaceLevel = surfaceLevel, DInBetweenAquiferlayer = inBetweenAquiferlayerHeight, - DBottomAquiferlayer = soilProfile1D.GetLayerHeight(soilProfile1D.BottomAquiferLayer), + DBottomAquiferlayer = bottomAquiferlayerHeight, PermeabilityInBetweenAquiferlayer = inBetweenAquiferlayerSoil.PermeabKx, PermeabilityBottomAquiferlayer = soilProfile1D.BottomAquiferLayer.Soil.PermeabKx, };