Page History

We will now add the values for the boundary conditions by importing the corresponding data from a csv file with the function ImportBoundaryConditions function  ImportBoundaryConditions  . The information regarding the lateral sources can be similarly imported using the function function  ImportLateralData  . In both cases, for the sake of simplicity, we have assumed that they have constant (not time dependent) values.

ImportBoundaryConditions(rootPath + r"\boundaryConditions.csv", flowModel)
ImportLateralData(rootPath + r"\laterals.csv", flowModel)

We continue by setting the roughness . Starting with of the main roughness section. The following code sets the default roughness of the , using the SetDefaultRoughness function, sets the main roughness section to StricklerKs the type Strickler k_s with a default value of 30.

SetDefaultRoughness(flowModel, "Main", RoughnessType.StricklerKs, 30)

Now Next, we create a new roughness section and "called FloodPlain" and  and make cross-sections "prof_SW2815-SW2870_Bo" and " and prof_D20060515-DP-295" use  use it. To declare roughness sections on cross-sections, you need to declare specify the start and end distance along the cross-section for which that specific value will be used. Because we want these cross-sections to only use our new "FloodPlain" section we use the min and max FloodPlain section along their entire profile, we need to calculate the minimum and maximum horizontal relative coordinates values of the cross-section profile using the "GetMinYMaxYofCrossSectionProfile" functionthe GetMinYMaxYofCrossSectionProfile  Image Added  function. The function GetItemByName  Image Added  will provide us the objects corresponding to the desired cross sections.

sectionFloodPlain = AddNewRoughnessSection(flowModel,"FloodPlain")

crs1 = GetItemByName(network.CrossSections, "prof_SW2815-SW2870_Bo")
crs2 = GetItemByName(network.CrossSections, "prof_D20060515-DP-295")

minY, maxY = GetMinYMaxYofCrossSectionProfile(crs1)
AddCrossSectionRoughness(flowModel, crs1, minY, maxY, sectionFloodPlain)

minY, maxY = GetMinYMaxYofCrossSectionProfile(crs2)
AddCrossSectionRoughness(flowModel, crs2, minY, maxY, sectionFloodPlain)

Now At this point, all we only need to do is to set the initial depth to a default of 3 using the following lines:conditions for the model. Let's assume a uniform water depth value of 3 m in the entire model domain.

SetInitialConditionType(flowModel, InitialConditionType.Depth)
flowModel.DefaultInitialDepth = 3