Add external forcings, roughness and initial conditions

Create a new region in your script for this section (you can call it Add external forcings, roughness and initial conditions).

We will now add the values for the boundary conditions by importing the corresponding data from a csv file with the function  ImportBoundaryConditions  . The information regarding the lateral sources can be similarly imported using the 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 of the main roughness section. The following code, using the SetDefaultRoughness function, sets the main roughness section to the type Strickler k_s with a default value of 30.

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

Next, we create a new roughness section called FloodPlain and make cross-sections prof_SW2815-SW2870_Bo and prof_D20060515-DP-295 use it. To declare roughness sections on cross-sections, you need to 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 use our new 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    function. The function GetItemByName    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)

At this point, all we need to do is to set the initial 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

Note: The icon  indicates functions that you might find interesting to see how they have been built. You can do this by opening the corresponding library (another normal script python file on its own) in the toolbox and checking the code which defines the function or method you are interested in.