Demo Iteration Project
The Demo Iteration project is available for all users and can be found in the main menu under Edit projects. This project does not count towards your license.
This project is intended for hydrologists and other water experts.
This project showcases the rainfall overlay in a iteration mode.
The demo show the iteration functionality that can be used to adjust the ground infiltration speed using a combo formula.
Calculations using the Iteration option
The demo project simulates a Dutch polder to demonstrate how rainfall infiltration can be modeled using the iteration option. A hypothetical formula is applied to calculate the infiltration rate per time step, based on the saturation level of the unsaturated zone. This iterative process determines the soil's saturation step by step, which in turn influences the potential infiltration.
The purpose of this demo project is to illustrate the capabilities of this iterative calculation method.
Setup of the Demo Project
The demo includes several essential components for performing these calculations. Each component is explained, along with the rationale for its use.
Components:
Combo Overlays Rainfall Overlay Iteration Collector Indicator Number of Iterations
Combo Overlays
Several overlays need to be prepared as inputs for the various iteration steps. Each step's outputs are set up for use in the subsequent iteration. The combo overlays include:
Combo Surface U and Combo Surface V: These represent the initial velocity of surface water in the x and y directions. When water flows downhill, the velocity at the end of an iteration is used as the initial velocity for the next iteration. Formula: IF(EQ(ITERATION, 0), 0, A)
Combo Fraction and Combo Groundwater Depth: These calculate inputs based on previous iterations unless it’s iteration 0, where global values are used:
Saturation percentage: GLOBAL_DEFAULT_FRACTION Groundwater level: GLOBAL_DEFAULT_TABLE Formula examples: IF(EQ(ITERATION, 0), GLOBAL_DEFAULT_FRACTION, IF(EQ(A, NO_DATA), GLOBAL_DEFAULT_FRACTION, A)) IF(EQ(ITERATION, 0), GLOBAL_DEFAULT_TABLE, A)
Combo Surface: Determines input from previous iterations unless it’s iteration 0. Formula: IF(EQ(ITERATION, 0), 0, A)
Combo Infiltration Speed: Calculates infiltration speed based on soil saturation using the Ground BOFEK ID overlay to identify soil type. Formula example: MAX(0, SUB(GLOBAL_DEFAULT_INFIL, SWITCH(B, 0, 1004, MUL(10, SUB(A, GLOBAL_DEFAULT_FRACTION)), 1006, MUL(5, SUB(A, GLOBAL_DEFAULT_FRACTION)))))
Rainfall Overlay
Default settings are used for the rainfall overlay, with some exceptions:
Rain Event Duration: The total duration is specified and later divided into iterations. Water System – Water Level: Set initial water level and UV flow using the Surface combo overlay. Groundwater: Initialize groundwater level using the Groundwater Depth combo overlay. Unsaturated Zone: Initialize saturation and storage percentages per grid cell using the Fraction combo overlay. Generated Output Overlays:
Surface Last Value Rain Ground Last Value Ground Watertable Ground Last Unsaturated Storage Ground Last Unsaturated Fraction Surface Last U Surface Last V The number of timeframes is set to 1.
Iteration Collector
Displays results of each iteration step. The input overlay is the Rain Event.
Modulus: Set to 1 by default, meaning one output is saved per iteration as a timeframe in the collector.
Indicator
Handles partial updates for specific timeframes. Example formulas:
UPDATE_OVERLAY_PARTIAL_END_SEC_WHERE_OVERLAY_IS_10 Queries like SELECT_ITERATIONS_WHERE and SELECT_ITERATION_WHERE retrieve the number of iterations and the current iteration.
Number of Iterations
Set through the interface (default is 1).
Notes and Limitations
Debugging cannot be downloaded per iteration but can be saved for each iteration. Specific updates (e.g., updating attributes for certain building types in iteration 5) can be executed using the API and indicator but not through TQL for non-initial iterations.
Exploring the project
To begin exploring the project, see the outlet located on the left side.