// Copyright (C) Stichting Deltares 2016. All rights reserved. // // This file is part of Ringtoets. // // Ringtoets is free software: you can redistribute it and/or modify // it under the terms of the GNU General Public License as published by // the Free Software Foundation, either version 3 of the License, or // (at your option) any later version. // // This program is distributed in the hope that it will be useful, // but WITHOUT ANY WARRANTY; without even the implied warranty of // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the // GNU General Public License for more details. // // You should have received a copy of the GNU General Public License // along with this program. If not, see . // // All names, logos, and references to "Deltares" are registered trademarks of // Stichting Deltares and remain full property of Stichting Deltares at all times. // All rights reserved. using System; using System.Linq; using Core.Common.Base; using Core.Common.Base.Data; using Core.Common.Base.Geometry; using Ringtoets.Common.Data.Calculation; using Ringtoets.Common.Data.DikeProfiles; using Ringtoets.Common.Data.Probabilistics; using Ringtoets.Common.Data.Properties; using Ringtoets.HydraRing.Data; namespace Ringtoets.Common.Data.Structures { /// /// Base class that holds all structures calculation specific input parameters. /// public abstract class StructuresInputBase : Observable, ICalculationInput, IUseBreakWater, IUseForeshore where T : StructureBase { private readonly NormalDistribution modelFactorSuperCriticalFlow; private readonly LogNormalDistribution allowedLevelIncreaseStorage; private readonly VariationCoefficientLogNormalDistribution storageStructureArea; private readonly LogNormalDistribution flowWidthAtBottomProtection; private readonly VariationCoefficientLogNormalDistribution criticalOvertoppingDischarge; private readonly VariationCoefficientNormalDistribution widthFlowApertures; private readonly VariationCoefficientLogNormalDistribution stormDuration; private T structure; private RoundedDouble structureNormalOrientation; private double failureProbabilityStructureWithErosion; private ForeshoreProfile foreshoreProfile; protected StructuresInputBase() { structureNormalOrientation = new RoundedDouble(2, double.NaN); modelFactorSuperCriticalFlow = new NormalDistribution(2) { Mean = (RoundedDouble) 1.1, StandardDeviation = (RoundedDouble) 0.03 }; allowedLevelIncreaseStorage = new LogNormalDistribution(2) { Mean = (RoundedDouble) double.NaN, StandardDeviation = (RoundedDouble) 0.1 }; storageStructureArea = new VariationCoefficientLogNormalDistribution(2) { Mean = (RoundedDouble) double.NaN, CoefficientOfVariation = (RoundedDouble) 0.1 }; flowWidthAtBottomProtection = new LogNormalDistribution(2) { Mean = (RoundedDouble) double.NaN, StandardDeviation = (RoundedDouble) 0.05 }; criticalOvertoppingDischarge = new VariationCoefficientLogNormalDistribution(2) { Mean = (RoundedDouble) double.NaN, CoefficientOfVariation = (RoundedDouble) 0.15 }; failureProbabilityStructureWithErosion = double.NaN; widthFlowApertures = new VariationCoefficientNormalDistribution(2) { Mean = (RoundedDouble) double.NaN, CoefficientOfVariation = (RoundedDouble) 0.05 }; stormDuration = new VariationCoefficientLogNormalDistribution(2) { Mean = (RoundedDouble) 6.0, CoefficientOfVariation = (RoundedDouble) 0.25 }; UpdateForeshoreProperties(); } #region Model factors /// /// Gets or sets the model factor for the super critical flow. /// /// Only sets the mean. public NormalDistribution ModelFactorSuperCriticalFlow { get { return modelFactorSuperCriticalFlow; } set { modelFactorSuperCriticalFlow.Mean = value.Mean; } } #endregion /// /// Synchronizes the input properties with the properties of the structure. /// protected abstract void UpdateStructureProperties(); private static bool ValidProbabilityValue(double probability) { return !double.IsNaN(probability) && probability <= 1 && probability >= 0; } #region Hydraulic data public HydraulicBoundaryLocation HydraulicBoundaryLocation { get; set; } /// /// Gets or sets the storm duration. /// [hrs] /// /// Only sets the mean. public VariationCoefficientLogNormalDistribution StormDuration { get { return stormDuration; } set { stormDuration.Mean = value.Mean; } } #region Schematization /// /// Gets or sets the structure. /// public T Structure { get { return structure; } set { structure = value; UpdateStructureProperties(); } } /// /// Gets or sets the orientation of the normal of the structure. /// [degrees] /// ///Thown when the value for the orientation is not between [0,360] degrees. public RoundedDouble StructureNormalOrientation { get { return structureNormalOrientation; } set { if (double.IsNaN(value)) { structureNormalOrientation = value.ToPrecision(structureNormalOrientation.NumberOfDecimalPlaces); return; } RoundedDouble newOrientationValue = value.ToPrecision(structureNormalOrientation.NumberOfDecimalPlaces); if (newOrientationValue < 0 || newOrientationValue > 360) { throw new ArgumentOutOfRangeException("value", Resources.Orientation_Value_needs_to_be_between_0_and_360); } structureNormalOrientation = newOrientationValue; } } /// /// Gets or sets the allowed increase of level for storage. /// [m] /// public LogNormalDistribution AllowedLevelIncreaseStorage { get { return allowedLevelIncreaseStorage; } set { allowedLevelIncreaseStorage.Mean = value.Mean; allowedLevelIncreaseStorage.StandardDeviation = value.StandardDeviation; } } /// /// Gets or sets the storage structure area. /// [m^2] /// public VariationCoefficientLogNormalDistribution StorageStructureArea { get { return storageStructureArea; } set { storageStructureArea.Mean = value.Mean; storageStructureArea.CoefficientOfVariation = value.CoefficientOfVariation; } } /// /// Gets or sets the flow width at the bottom protection. /// [m] /// public LogNormalDistribution FlowWidthAtBottomProtection { get { return flowWidthAtBottomProtection; } set { flowWidthAtBottomProtection.Mean = value.Mean; flowWidthAtBottomProtection.StandardDeviation = value.StandardDeviation; } } /// /// Gets or sets the critical overtopping discharge. /// [m^3/s/m] /// public VariationCoefficientLogNormalDistribution CriticalOvertoppingDischarge { get { return criticalOvertoppingDischarge; } set { criticalOvertoppingDischarge.Mean = value.Mean; criticalOvertoppingDischarge.CoefficientOfVariation = value.CoefficientOfVariation; } } /// /// Gets or sets the failure probability of structure given erosion. /// [1/year] /// /// Thrown when the value of the probability /// is not between [0, 1]. public double FailureProbabilityStructureWithErosion { get { return failureProbabilityStructureWithErosion; } set { if (!ValidProbabilityValue(value)) { throw new ArgumentOutOfRangeException("value", Resources.FailureProbability_Value_needs_to_be_between_0_and_1); } failureProbabilityStructureWithErosion = value; } } /// /// Gets or sets the width of flow apertures. /// [m] /// public VariationCoefficientNormalDistribution WidthFlowApertures { get { return widthFlowApertures; } set { widthFlowApertures.Mean = value.Mean; widthFlowApertures.CoefficientOfVariation = value.CoefficientOfVariation; } } #endregion #endregion #region Foreshore Profile /// /// Gets or sets the foreshore profile. /// public ForeshoreProfile ForeshoreProfile { get { return foreshoreProfile; } set { foreshoreProfile = value; UpdateForeshoreProperties(); } } /// /// Gets or sets whether the needs to be taken into account. /// public bool UseBreakWater { get; set; } /// /// Gets or sets whether the needs to be taken into account. /// public bool UseForeshore { get; set; } /// /// Gets the geometry of the foreshore. /// public RoundedPoint2DCollection ForeshoreGeometry { get { return foreshoreProfile != null ? foreshoreProfile.Geometry : new RoundedPoint2DCollection(2, Enumerable.Empty()); } } /// /// Gets the . /// public BreakWater BreakWater { get; private set; } private void UpdateForeshoreProperties() { if (foreshoreProfile == null) { UseForeshore = false; UseBreakWater = false; BreakWater = GetDefaultBreakWaterProperties(); } else { UseForeshore = foreshoreProfile.Geometry.Count() > 1; UseBreakWater = foreshoreProfile.HasBreakWater; BreakWater = foreshoreProfile.HasBreakWater ? new BreakWater(foreshoreProfile.BreakWater.Type, foreshoreProfile.BreakWater.Height) : GetDefaultBreakWaterProperties(); } } private static BreakWater GetDefaultBreakWaterProperties() { return new BreakWater(BreakWaterType.Dam, 0.0); } #endregion } }