First of all, we will keep the time series corresponding to the measured data and to the results of the original model set up in a list. This variable will be used later on to compare them quickly with each other.
| Code Block | ||||
|---|---|---|---|---|
| ||||
# Save measurements and initial model results in a list of time series to compare
timeSeriesToCompare = [measuredTimeSeries, waterlevelResults] |
The model calibration will be done by changing the less reliable parameters from the loaded model. Let’s assume that these factors are the flow boundary condition at the node Boxbergen, and the roughness of the main tributaries leading to Zwolle:
The definitions of boundary conditions object are not known yet in the scripting environment as variables. The complete name of the flow boundary condition at Boxbergen is, as shown in the previous figure, composed by the name of the location and the current value (Boxbergen - Q:25.064 m^3/s). However, the function GetBoundaryDataByName defined in API library SobekWaterFlowFunctions already takes this into account, so you don't have to worry about it. Assign the name of the boundary condition to a new variable bcCalibrationName. Then use that string to assign the Boundary data object to the variable bcCalibration. The property Flow of this object contains the current numeric value of this constant flow boundary. We will save this number in the variable bcInitialValue.
...