Index: DamEngine/trunk/src/Deltares.DamEngine.Calculators/PlLinesCreator/PlLinesCreator.cs =================================================================== diff -u -r6851 -r6866 --- DamEngine/trunk/src/Deltares.DamEngine.Calculators/PlLinesCreator/PlLinesCreator.cs (.../PlLinesCreator.cs) (revision 6851) +++ DamEngine/trunk/src/Deltares.DamEngine.Calculators/PlLinesCreator/PlLinesCreator.cs (.../PlLinesCreator.cs) (revision 6866) @@ -438,13 +438,14 @@ private bool IsPenetrationZoneCompletelyWithinInfiltrationZoneFor1DGeometry(double penetrationLength) { SoilProfile1D relevantSoilProfile = GetRelevantSoilProfileForAquiferLayersSearch(); - // The infiltration layers are the aquitards situated between the bottom aquifer and: + double bottomLevelLowestInBetweenAquiferCluster = relevantSoilProfile.GetInBetweenAquiferClusters.Last().Item2.TopLevel; + double topLevelBottomAquifer = relevantSoilProfile.BottomAquiferLayer.TopLevel; + // The infiltration layers are the aquitards situated between the bottom aquifer (topLevelBottomAquifer) and: // - the surface line if no in-between aquifer is present - // - the lowest in-between aquifer if at least one in-between aquifer is present + // - the lowest in-between aquifer (bottomLevelLowestInBetweenAquiferCluster) if at least one in-between aquifer is present IList infiltrationLayers = (from SoilLayer1D layer in relevantSoilProfile.Layers - where (relevantSoilProfile.InBetweenAquiferLayer == null || layer.TopLevel < relevantSoilProfile.GetInBetweenAquiferClusters.Last().Item2.TopLevel) && - layer.TopLevel > relevantSoilProfile.BottomAquiferLayer.TopLevel - select layer).ToList(); + where (relevantSoilProfile.InBetweenAquiferLayer == null || layer.TopLevel < bottomLevelLowestInBetweenAquiferCluster) && + layer.TopLevel > topLevelBottomAquifer select layer).ToList(); if (infiltrationLayers.Count > 0) {