MIKE11_UserManual
MIKE11_UserManual
User Guide
Limited Liability
The liability of DHI is limited as specified in Section III of your 'DHI
Software Licence Agreement':
'IN NO EVENT SHALL DHI OR ITS REPRESENTATIVES (AGENTS
AND SUPPLIERS) BE LIABLE FOR ANY DAMAGES WHATSO-
EVER INCLUDING, WITHOUT LIMITATION, SPECIAL, INDIRECT,
INCIDENTAL OR CONSEQUENTIAL DAMAGES OR DAMAGES
FOR LOSS OF BUSINESS PROFITS OR SAVINGS, BUSINESS
INTERRUPTION, LOSS OF BUSINESS INFORMATION OR OTHER
PECUNIARY LOSS ARISING OUT OF THE USE OF OR THE INA-
BILITY TO USE THIS DHI SOFTWARE PRODUCT, EVEN IF DHI
HAS BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGES.
THIS LIMITATION SHALL APPLY TO CLAIMS OF PERSONAL
INJURY TO THE EXTENT PERMITTED BY LAW. SOME COUN-
TRIES OR STATES DO NOT ALLOW THE EXCLUSION OR LIMITA-
TION OF LIABILITY FOR CONSEQUENTIAL, SPECIAL, INDIRECT,
INCIDENTAL DAMAGES AND, ACCORDINGLY, SOME PORTIONS
OF THESE LIMITATIONS MAY NOT APPLY TO YOU. BY YOUR
OPENING OF THIS SEALED PACKAGE OR INSTALLING OR
USING THE SOFTWARE, YOU HAVE ACCEPTED THAT THE
ABOVE LIMITATIONS OR THE MAXIMUM LEGALLY APPLICA-
BLE SUBSET OF THESE LIMITATIONS APPLY TO YOUR PUR-
CHASE OF THIS SOFTWARE.'
Printing History
June 2004
June 2005
December 2006
November 2007
January 2009
December 2011
3
4 MIKE 11
CONTENTS
5
Simulation Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1 SIMULATION EDITOR . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.1 Models . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.1.1 Models . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
1.1.2 Simulation Mode . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.2 Input . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
1.3 Simulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
1.3.1 Simulation Period . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
1.3.2 Initial Conditions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
1.4 Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.5 Start . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
6 MIKE 11
2.3.12 Dambreak Str. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 143
2.3.13 Dambreak Erosion . . . . . . . . . . . . . . . . . . . . . . . . . . 149
2.3.14 User Defined Structure . . . . . . . . . . . . . . . . . . . . . . . 151
2.3.15 Tabulated Structure . . . . . . . . . . . . . . . . . . . . . . . . . 152
2.3.16 Energy Loss . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 154
2.4 Tabular view: Routing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 155
2.4.1 Channel routing . . . . . . . . . . . . . . . . . . . . . . . . . . . 157
2.4.2 Flood control Q and Q-rate . . . . . . . . . . . . . . . . . . . . . 158
2.4.3 Flood control H-Q / H-V curve . . . . . . . . . . . . . . . . . . . 158
2.4.4 Flood control by orifice . . . . . . . . . . . . . . . . . . . . . . . 160
2.4.5 Diversions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 160
2.4.6 Kinematic Routing Method . . . . . . . . . . . . . . . . . . . . . 162
2.5 Tabular view: Runoff / Groundwater Links . . . . . . . . . . . . . . . . . 164
2.5.1 MIKE SHE Links . . . . . . . . . . . . . . . . . . . . . . . . . . . 165
2.5.2 Rainfall-runoff links . . . . . . . . . . . . . . . . . . . . . . . . . 171
2.6 Tabular View: Grid Points . . . . . . . . . . . . . . . . . . . . . . . . . . 172
7
4 BOUNDARY EDITOR . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 229
4.1 Overview of the Boundary Editor . . . . . . . . . . . . . . . . . . . . . . . 229
4.1.1 The Boundary Table - Upper Split Window . . . . . . . . . . . . 230
4.1.2 Specifying the Boundary Description . . . . . . . . . . . . . . . . 231
4.1.3 Specifying the Boundary Type, Data Type and File/Values . . . 236
4.2 Boundary Editor Tools . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 254
4.2.1 Quick set up of Graded Sediment Boundaries . . . . . . . . . . 254
4.2.2 Quick set up of AD Boundaries . . . . . . . . . . . . . . . . . . . 255
4.2.3 Copying Point Source Boundaries . . . . . . . . . . . . . . . . . 256
4.2.4 Scale factor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 257
4.3 Users Upgrading from MIKE 11 Version 2002 or Previous Versions . . . 259
8 MIKE 11
5.7.3 Rainfall . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 305
5.8 Time Series . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 310
5.9 Parameters menu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 314
5.9.1 Enlargement ratio . . . . . . . . . . . . . . . . . . . . . . . . . . 314
5.9.2 Loss Parameters . . . . . . . . . . . . . . . . . . . . . . . . . . . 314
5.9.3 Land use definitions for QLSF method . . . . . . . . . . . . . . 314
5.9.4 Default values for specific method . . . . . . . . . . . . . . . . . 315
5.9.5 Time-fixed combinations . . . . . . . . . . . . . . . . . . . . . . 315
5.9.6 MAW merged output file. . . . . . . . . . . . . . . . . . . . . . . 315
5.10 Basin View . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 315
5.10.1 Activating the Basin View . . . . . . . . . . . . . . . . . . . . . . 315
5.10.2 Importing Layers . . . . . . . . . . . . . . . . . . . . . . . . . . . 316
5.10.3 Basin Work Area . . . . . . . . . . . . . . . . . . . . . . . . . . . 316
5.10.4 Preparing Catchments . . . . . . . . . . . . . . . . . . . . . . . 321
5.10.5 Inserting Rainfall Stations . . . . . . . . . . . . . . . . . . . . . 322
5.10.6 Preparing Thiessen weights . . . . . . . . . . . . . . . . . . . . 322
5.11 Result Presentation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 323
5.12 A Step-by-step procedure for using the RR-Editor . . . . . . . . . . . . 326
9
6.9 Stratification . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 349
6.10 Time Series Output . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 352
6.10.1 Generating Time Series Output Files . . . . . . . . . . . . . . . . 353
6.10.2 Text File Settings . . . . . . . . . . . . . . . . . . . . . . . . . . . 353
6.10.3 Time Period in Output . . . . . . . . . . . . . . . . . . . . . . . . 353
6.10.4 Items in Output . . . . . . . . . . . . . . . . . . . . . . . . . . . . 354
6.10.5 Grid Points in Output . . . . . . . . . . . . . . . . . . . . . . . . . 354
6.11 Maps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 354
6.11.1 A step by step guide to generating two-dimensional maps . . . 361
6.11.2 A step by step guide to generating Digital Elevation Models (DEM)
363
6.12 Groundwater Leakage . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 364
6.13 Reach Lengths . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 365
6.14 Add. Output . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 367
6.14.1 Additional output for QSS with vegetation . . . . . . . . . . . . . 370
6.15 Flood Plain Resistance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 371
6.16 User Def. Marks . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 372
6.16.1 Activation of Bed resistance Triple Zone Approach . . . . . . . . 372
6.17 Encroachment . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 373
6.17.1 Iteration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 374
6.17.2 Location . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 374
6.17.3 Encroachment method . . . . . . . . . . . . . . . . . . . . . . . . 375
6.17.4 Encroachment positions . . . . . . . . . . . . . . . . . . . . . . . 375
6.17.5 Reduction parameters (only encroachment methods 3 to 5) . . 376
6.17.6 Target Values . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 376
6.17.7 Encroachment simulation overview . . . . . . . . . . . . . . . . . 377
6.17.8 Encroachment station overview . . . . . . . . . . . . . . . . . . . 378
6.17.9 General guide lines for carrying out encroachment simulations . 378
6.18 MIKE 12 Parameters . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 378
6.19 MIKE 12 Initial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 379
6.20 Mixing Coefficients . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 381
6.20.1 Water & Water . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 381
6.20.2 Location . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 382
6.20.3 Water & Vegetation . . . . . . . . . . . . . . . . . . . . . . . . . . 382
6.21 W. L. Incr.- Curve . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 382
6.21.1 General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 383
6.21.2 System Definition . . . . . . . . . . . . . . . . . . . . . . . . . . . 384
6.21.3 Tabular view . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 384
6.22 W. L. Incr.- Sand Bars . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 384
6.22.1 General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 385
6.22.2 System Definition . . . . . . . . . . . . . . . . . . . . . . . . . . . 385
6.22.3 Tabular view . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 386
10 MIKE 11
Advection-Dispersion Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . 387
7 ADVECTION-DISPERSION EDITOR . . . . . . . . . . . . . . . . . . . . . . . 389
7.1 Module Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 389
7.1.1 Advection-Dispersion module (AD) . . . . . . . . . . . . . . . . 389
7.1.2 The Advection-Dispersion Equation . . . . . . . . . . . . . . . . 390
7.1.3 Cohesive Sediment Transport module (CST) . . . . . . . . . . 391
7.1.4 Advanced Cohesive Sediment Transport module (ACST) . . . 391
7.2 Components . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 391
7.3 Dispersion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 394
7.4 Init. cond. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 396
7.5 Decay . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 399
7.6 Cohesive ST . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 401
7.6.1 Single Layer Cohesive Model . . . . . . . . . . . . . . . . . . . 401
7.6.2 Multi Layer Cohesive Model . . . . . . . . . . . . . . . . . . . . 403
7.6.3 Cohesive Sediment transport dynamics; Description . . . . . . 405
7.7 Sediment layers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 409
7.7.1 Initial layer thicnkness values from external file . . . . . . . . . 411
7.8 Non-cohesive ST . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 414
7.9 Additional output . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 416
7.10 Dispersion, MIKE 12 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 417
7.11 Additional Output, MIKE 12 . . . . . . . . . . . . . . . . . . . . . . . . . 418
7.12 AD Simulations, Additional features . . . . . . . . . . . . . . . . . . . . . 419
7.12.1 Alternative Solution scheme; The ULTIMATE method . . . . . 419
7.12.2 AD simulation features activated through the MIKE11.Ini file . 420
11
9.0.2 The transport models . . . . . . . . . . . . . . . . . . . . . . . . . 436
9.1 Sediment grain diameter . . . . . . . . . . . . . . . . . . . . . . . . . . . . 437
9.2 Transport model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 438
9.2.1 Model Parameters . . . . . . . . . . . . . . . . . . . . . . . . . . . 440
9.2.2 Special features for specific transport models . . . . . . . . . . . 443
9.2.3 Bottom level update methods . . . . . . . . . . . . . . . . . . . . 444
9.3 Calibration factors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 446
9.4 Data for graded ST . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 447
9.5 Preset distribution of sediment in nodes . . . . . . . . . . . . . . . . . . . 448
9.6 Passive branches . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 450
9.7 Non-Scouring Bed Level . . . . . . . . . . . . . . . . . . . . . . . . . . . . 451
9.8 Initial dune dimensions (Engelund Fredsoe model) . . . . . . . . . . . . 453
9.9 ST Simulations; Additional features . . . . . . . . . . . . . . . . . . . . . 453
9.9.1 Transport layer depth options through ASCII file; TrsDepth.Inp . 454
9.9.2 Shielding option defined through ASCII file; Shield.inp . . . . . . 454
9.9.3 ST simulation features activated through the MIKE11.Ini file . . 455
12 MIKE 11
11 DATA ASSIMILATION EDITOR . . . . . . . . . . . . . . . . . . . . . . . . . . . 475
11.1 General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 476
11.1.1 Module selection . . . . . . . . . . . . . . . . . . . . . . . . . . . 476
11.1.2 Basic parameters . . . . . . . . . . . . . . . . . . . . . . . . . . 477
11.1.3 Forecast . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 478
11.2 Measurements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 478
11.2.1 Measurement location . . . . . . . . . . . . . . . . . . . . . . . 479
11.2.2 Standard deviation . . . . . . . . . . . . . . . . . . . . . . . . . . 479
11.2.3 Weighting function . . . . . . . . . . . . . . . . . . . . . . . . . . 480
11.2.4 Error forecast model . . . . . . . . . . . . . . . . . . . . . . . . . 481
11.3 Equation Editor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 482
11.3.1 General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 482
11.3.2 Parameter definition . . . . . . . . . . . . . . . . . . . . . . . . . 485
11.4 Boundary Statistics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 486
11.4.1 Details . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 487
11.5 Standard deviation editor . . . . . . . . . . . . . . . . . . . . . . . . . . . 488
11.6 Output . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 491
11.7 A step by step guide to uncertainty assessment . . . . . . . . . . . . . 492
11.8 A step by step guide to updating using the Kalman filter method . . . . 494
11.9 A step by step guide to updating using the Weighting function method 495
11.10 Examples . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 496
11.10.1 Uncertainty assessment on hydrodynamic simulation . . . . . 497
11.10.2 Kalman filter updating on hydrodynamic set-up . . . . . . . . . 497
11.10.3 Uncertainty assessment on advection dispersion simulation . 498
11.10.4 Kalman filter updating on advection dispersion set-up . . . . . 499
13
13.8 Detailed Result page . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 517
13.9 Outlined Output page . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 518
Appendix A . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 521
Appendix B . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 529
Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 539
14 MIKE 11
SIMULATION EDITOR
15
16 MIKE 11
1 SIMULATION EDITOR
The simulation editor serves three purposes:
Figure 1.1
The editing of cross sections is a typical example of this link, where the
graphical view of the network editor e.g. is used to present and select cross
sections from the cross section editor.
The linkage requires a file name to be specified for each of the required
editors. The file names are input on the Input Property Page of the simula-
tion editor. An alternative is to select a file from the File Menu which will
recall the appropriate editor. The edit menu can then be used to edit the
objects
Simulation Editor 17
Simulation Editor
1.1 Models
Figure 1.2 The Models tab. Note that the Simulation Mode Box may differ if a
quasi two dimensional steady state solver with vegetation is not
installed.
This page is used to define the simulation models to execute and the simu-
lation mode (unsteady or quasi steady).
1.1.1 Models
The following abbreviations of module names are used:
HD: Hydrodynamic
AD: Advection-Dispersion
18 MIKE 11
Models
RR: Rainfall-Runoff
Some of the model that can be selected are dependent on other modules in
a simulation and it is therefore required to have more modules selected.
This rule for model-dependency is implemented such that once a model is
selected there will be an automatic selection of eventual dependent models
(e.g. Selection of FF-model selects HD-model also, Selection of ECOLab
selects AD-model also etc.).
Water Quality modelling takes place through the ECO Lab model entry
where a variety of Water Quality models can be selected from so-called
ECO Lab templates.
The River Ice Model as listed above is the result of a development project
and the module is not yet completed and publically available as a standard
module. Please consult DHI for the current state if you have interest in
modelling river ice formation etc.
Unsteady
The HD calculations are based on hydrodynamic flow conditions.
Quasi steady
At every time step the calculations are based on steady flow conditions.
Simulation Editor 19
Simulation Editor
20 MIKE 11
Input
1.2 Input
Based on the model selection from the Models Property Page a number of
filename fields becomes active, and the user is required to specify a range
of input file names.
Note that files required for the simulation are indicated by the active fields
on the page. Two exceptions exists though; the fields for the result files.
It is important to notice, that files specified in the Models page are all
input files for a simulation, so the result-file name fields are Not to be
interpretated as being result filenames for the simulation. Instead, these
result-files indicate an input to the simulation as described below.
Simulation Editor 21
Simulation Editor
1.3 Simulation
The simulation property page contains details of simulation time, time
stepping specifications and initial conditions for each of the chosen types
of models.
22 MIKE 11
Simulation
In case the time stepping type is specified as ‘Tabulated’, the time steps
are specified by activating the settings button and selecting a timeserie
from a time series file (dsf0) with Item-type defined as ‘TimeStep’ in the
dialog presented in Figure 1.5.
Simulation Editor 23
Simulation Editor
Figure 1.5 Dialog for specification of the tabulated time step time series file.
The settings for time step adaptation is similarly specified in a menu acti-
vated by the Settings button. The menu is presented as Figure 1.6.
Minimum/maximum/unit
Defines the limits for the adaptation of the time step.
Change ratio
The time step is successively lowered with change ratio until the criterias
specified in this menu are met. The starting value for the adaptation within
the time step is change ratio times the previous time step.
Criterias
The time step adaptation model offers seven criteras which may each be
enabled or disablet and given threshold values.
24 MIKE 11
Simulation
Simulation Editor 25
Simulation Editor
Period
The date and time for the start and end of the simulation period. The
standard windows date time format is used.
26 MIKE 11
Simulation
Type of condition
– Steady State: (HD only).
The initial conditions will be calculated automatically assuming a
steady state condition with discharges and water levels at the
boundaries corresponding to the start time of the simulation.
– Parameter File:
The initial conditions will be taken from the parameter file relevant
to the module in question.
– Hotstart:
The initial conditions will be loaded from an existing result file.
– Steady+Parameter: (HD only).
The initial conditions will be established using both the steady state
and parameter file method. In those grid points where data are spec-
ified in the Initial (p. 332) Property Page of the Hydrodynamic
parameters Editor (p. 331) the initial conditions will be taken from
the parameter file, other grid points will be calculated using the
steady state option.
Hotstart Filename
The name of an existing result file from which initial conditions for the
simulation should be loaded. For hot-start to be valid using control struc-
tures an additional hydrodynamic result file (*addhd.res11) file must be
present in the same location as the standard hydrodynamic result file. Fur-
ther this additional result file should hold structure information as selected
in the additional result tab in the hydrodynamic editor.
Add to File
The results of the current simulation will be added to the end of the hot-
start file. Any information (in the hotstart file) after the simulation start
date will be lost. This part of the file will be replaced by the new simula-
tion results.
Simulation Editor 27
Simulation Editor
1.4 Results
For each of the modules selected on the Models Property Page the user
should specify a filename for saving of the simulation results.
The filename can not be edited if the flag “Add to File” has been selected
on the Simulation Property Page. In this case the selected hotstart file will
become the result file as well.
28 MIKE 11
Start
1.5 Start
If all specified input files exist, the “Start” button can be pressed and the
simulation will commence. The simulation will take place as a separate
process (running the program; ‘Mike11.exe’), and the progress of the sim-
ulation will be reported in the info bar of the editor.
Simulation Editor 29
Simulation Editor
Any error or warning message from the simulation will be presented in the
log-part of the editor and additionally, saved in a file with the same name
as the simulation file and a .log extension.
Upon completion the simulation results can be viewed using MIKE View.
30 MIKE 11
RIVER NETWORK EDITOR
31
32 MIKE 11
2 RIVER NETWORK EDITOR
The River Network Editor is a key editor of MIKE 11 as it supplies the
user with an overview of the current model setup components defined in
the network file it self but also by presenting data which is stored in the
other MIKE 11 editors.
Figure 2.1 The graphical overview display of the MIKE 11 Network Editor (=
the, Graphical View)
The Tabular View of the network editor contains a large group of sub-
pages in which the numerical values behind most of the objects being pre-
sented in the Graphical view are stored for editing. The tabular view con-
tains individual pages with edit fields and tables for digitization points,
river branches definitions and connections, hydraulic structures details
and other data.
The Tabular view can be opened from the Graphical view through the
‘View’->’Tabular view’ option in the main menu bar or simply by using
the hotkey combination; ‘<CTRL>+T’. The Tabular view, branch page is
presented below in Figure 2.2.
Details on features and usage of the Network Editor including both the
Graphical View and the Tabular View is described in the following.
Some of the features available in the Network Editor have been developed
in cooperation with CTI Engineering, CO., Ltd., Japan. Amongst these
are; Tabulated structures, Honma's weir formula, bridges (D’Aubuisson
and submerged bridge), Routing along channels, Outflow from
Dams/retarding basins and the Steady flow with vegetation.
34 MIKE 11
Graphical View
Figure 2.3 Geographical Area coordinates and map projection definition dialog
appearing when creating a new network file (.nwk11)
Note, that additional graphical view sessions can be opened using the New
Window item under the Window Menu.
Editing of the river network (i.e. the points and branches) is undertaken
using the Graphical Editing Toolbar. Editing tools are also found using the
Pop-Up Menu (right mouse button) where features like insert, edit and
delete functions are included. Typically the Pop-Up Menu is used for edit-
ing of cross section geometry, parameters, hydraulic structures and data
stored in other MIKE 11 editors.
Note that to access information from another editor other than the Net-
work Editor, an editor file name must be specified using the Simulation
File Editor.
Figure 2.4 Illustration of right mouse pop up menu from where all data editors
can be accessed.
Select object. This icon activates the selection mode which is also
the default mode. Points, layers and other objects can be selected
by pointing and clicking with the left mouse button. Multiple
36 MIKE 11
Graphical View
Add new points. New points can be added by a point and click
operation using the left mouse button. Multiple points can be
added by pressing the left mouse button and holding it down
while moving the mouse along the desired path. New points will be cre-
ated with a spacing determined by the “minimum digitize distance” speci-
fied in the Mouse Settings property page of Network Settings. Points
added with this tool will always be added as free points, i.e. not connected
to a river branch. Alternatively you can add new points and define a river
branch in one operation using the following tool.
Add points and define branch. This tool creates points and
branches in a single operation. Point and click at successive loca-
tions along a desired path. Points can also be added by pressing
the left mouse button and holding it down while moving. Double click on
the last point to end the branch.
Delete points. This tool deletes both free points and points con-
nected by a branch. Move the cursor over the point (the cursor
will change style to indicate that a point has been detected) and
press the left mouse button to delete. Multiple points can be deleted by
holding the left mouse button down while moving the cursor over the
points.
Move points. This tool moves both free points and points con-
nected by branching. Select the point using the left mouse button
and then drag to the desired location.
Cut branch. This tool divides a single branch into two separate
branches. Move the cursor to the required segment where the
break is required. When the cursor changes style press the left
mouse button once to cut the branch.
Merge branch. This tool merges two separate branches into one.
Move the cursor to the beginning or the end of a branch, click at
this point with the left mouse button and drag to the connection
point on another branch.
Insert point. This tool will insert free points into an existing
branch. Move the cursor to a point on an existing branch, click
with the left mouse button and drag the cursor to the free point for
inclusion into the branch path.
38 MIKE 11
Graphical View
Repeat insert. The repeat insert tool will add a copy of the latest
object (weir, cross section, boundary condition, initial condition
etc.) created using the Insert facility in the Pop-Up Menu. The
repeat insert button is a fast and convenient way of inserting multiple
objects to the river network. The current object type is shown in the status
bar when the repeat insert button is activated. After activating the repeat
insert tool you should point and click once at the desired location of the
new object.
Select & edit. This tool is similar to the Edit facility found in the
Pop-Up Menu. The tool is a fast convenient way of accessing the
various editors required for objects on the river network. To con-
trol the number of editor windows activated use the Select and Edit Set-
tings Property Page of the Network Settings property sheet.
mouse button to delete. Multiple points can be deleted by holding the left
mouse button down while moving the cursor over the points
Insert points to alignment line. This tool will insert free points
into an existing alignment line. Move the cursor to a point on an
existing branch, click with the left mouse button and drag the cur-
sor to the free point for inclusion into the branch path. The free point to be
inserted must be added using the tool “Add new point” in the available in
the toolbar for river network editing.
Add points to alignment line. Using this tool you can add points
to an existing alignment line. Point are added at the upstream or
downstream end. Click once at the point to which you want to add
new points. Then point and click at successive locations along the desired
path.
40 MIKE 11
Graphical View
line, then click once at the upstream end of the second line, and the lines
are merged.
Additional features comprise the printing and the key issue in relation to
the Network Editor is the possibility of importing network data from ext-
renal text-file containing network data for either points, branches or both.
Figure 2.9 The File menu under the Main menu bar
Import
Example:
42 MIKE 11
Graphical View
The data in the file must be of the form shown in Figure 2.10. The file
should contain a section of the type “[AlignmentLine] ... EndSect // Align-
mentLine” for each alignment line.
Figure 2.10 The format used for importing alignment line data.
Figure 2.11 The Edit menu of the Network editor Main menu bar
The following features are applicable only when located in one of the
tables in the Network Editor, Tabular View property pages.
Remove rows
Delete a selected row in a table.
Insert rows
Insert additional row at highlighted location in table.
Find...
Search for a specific string in an active table
Replace...
Search and replace string values in tables.
Plot XY
Plot of two selected column in a XY-like plot.
Duplicate
Duplicate the currently selected item in an overview grid.
44 MIKE 11
Graphical View
Figure 2.12 The View menu of the Network Editor Main menu bar
Project Explorer
General MIKE Zero switch for showing the project explorer as a visible
component in the MIKE Zero interface (Primarily useful if the MIKE 11
editor files are part of a MIKE Zero project)
Project Map
General feature for the MIKE Zero project map feature. Switch for mak-
ing the project map visible within the MIKE Zero GUI.
Start Page
Switch on or off the MIKE Zero Start Page.
Tabular view
Used to open the tabular view of the network file.
Figure 2.13 Dialog for seleting the required longitudinal profile amongst all pos-
sible options available.
Network
The Network entry offers an easy way of turning on and off a number of
the most common used graphical network objects in the graphical view.
More detailed layout of the graphical view see Graphics (p. 55).
Boundary
Here the presentations of the different boundary types can quickly be
turned on or off. For a more detailed layout of the graphical view see
Graphics (p. 55).
Hydrodynamic Parameters
Here the presentations of the different hydrodynamic parameters can
quickly be turned on or off. For a more detailed layout of the graphical
view see Graphics (p. 55).
46 MIKE 11
Graphical View
Swap Profile
The Swap Profile option relates to the plotting of longitudinal profiles,
where the profile - once plotted - can be swapped in orientation (plotted
either Upstream towards Downstream locations or swapped; Downstrem
towards Upstream locations)
Draw Grid
A switch for drawing the grid or not in the graphical view.
Zoom features
Standard zoom features for zooming in and out, previous zoom, panning
or refreshing the view if required.
Export Graphics
The graphical view can be exported in the following ways:
– Copy to Clipboard.
– Save to metafile.
– Save to bitmap.
– Export layer graphics to file.
Switching off the Fixed Aspect Ration allows for distorting the content in
the graphical view such that the view always fills the entire window but
the width-height relation is not locked when adjusting the window size.
Toolbars...
Selection tool for MIKE Zero toolbars .
Status bar
Switch for showing or hiding the status bar in the bottom of the view.
Figure 2.15 Menu for resizing and changine map projection settings for the area
of the graphical view.
48 MIKE 11
Graphical View
The graphical view can be resized by entering the minimum and maxi-
mum coordinates for both the x-axis and the y-axis.
Note: Prior to using the network data in a Shape file for generating (parts
of) the river network it is required to load the same shape file as a back-
ground layer. Otherwise, it is not possible to select the appropriate shape-
file from the file selection drop-down boxes as shown in below.
Shapes file with point information can be used for generating points, and
Shape files with polylines can be used to generate points and branches or
only points. When generating points and branches at together, the
attributes defind in the Shape files can be used to define the river name
and eventually the Topo-ID (Topographical Identifier) for each branch
defined from the Shape file.
50 MIKE 11
Graphical View
Note that the last four options (Alignment lines, Vegetation zones, Dead
water zones due to E/C and Dead water zones behind vegetation) are con-
nected to the special Japanese version of MIKE 11 which requires special
modules in the license to activate.
Note, that in order to use this feature, it is required to have both the net-
work file and a boundary file ‘connected’, which is done by having a sim-
ulation file open in which both the network and the boundary files are
selected in the ‘Input’ page of the simulation editor
Update Chainages
This option is only meaningful if the Auto Update Chainages option is not
selected. The Update Chainages option could be used after having moved
one or several points.
Figure 2.18 The menu in which the number of the first point can be entered.
When joining two network files (see B.1.1 Merging .pfs files (p. 531)) it is
necessary that the number of the points in the two files do not overlap. To
avoid this it is possible to renumber the points internally in one of the net-
work files.
Note that the renumbering will not re-arrange the location of points in the
Tabular view, Points page but only change the internal numbering of
points in the network file which amongst other things are used to deter-
mine which points is used to define each river branch in the network.
Figure 2.19 The Layers menu of the Network editor Main menu bar
Add/remove
The Add/remove dialog as presented in Figure 2.20 offers the features for
managing layers in a network file.
1 ‘Add/Remove Layers’:
The Add/Remove property page includes options for adding or remov-
ing layer files in the network file. Use the button to add an addi-
tioal line to the Layers table and select the layer type to import in the
‘File Type’ drop-down selection box.
Browse for a layer file using the file browsing button; to search
for the layer file to import.
52 MIKE 11
Graphical View
2 ‘Overlay Manager’
If muliple layers have been imported it is possible to arrange the order
of the different layers in the ‘Overlay Manager’ dialog (Figure 2.21) by
selecting a specific layer and using the arrows; to either move
the specific layer up or down in the list of layers.
Layers are displayed in the Graphical View following the order defined
in this Overlay manager page such that layer in line 1 is plotted first,
then layer 2 on top of layer 1, then Layer 3 on top of Layer 1 and 2 etc.
Properties...
The properties dialog (Figure 2.22) defines the appearance and character-
istics for the individual layers imported to the Network editor file .
54 MIKE 11
Graphical View
Figure 2.23 The Settings menu of the Network Editor Main menu bar
Network...
The ‘Network...’ option opens the network settings dialog containing the
following four property pages :
z Graphics
z Mouse
z Network data
z Select and edit
Details on the property pages content follows below.
Graphics
The Graphics property page (Figure 2.24) controls the appearance and
partly the functionality of the graphical view.
The left hand side of the page presents the graphical object possible to
present on the graphical view - organized in a tree structure. Each object in
the tree-view has an entry for defining individual settings for points, lines,
labels etc. Open settings for editing a specific object by opening the tree-
view and highlight the object you wish to edit. Selecting an object allows
for changes in the right hand side of the dialog.
The ‘Display’ tick-box determinates wheter the specific object will be pre-
sented in the Graphical view (in case the object is present in the network
file presently edited).
Mouse
This property page sets the properties for the cursor (mouse) when digitiz-
ing and editing points on the graphical view.
56 MIKE 11
Graphical View
The Digitize distance; ‘Minimum distance between new points’ sets this
minimum distance at which a new point can be defined following a previ-
ous digitized point. Consequently, a smaller values allows for a fine detail
digitization (many points) whereas a large value reduces the number of
points during a digitization process.
Network data
The Network data property page (Figure 2.26) contains global, default set-
tings for misc. parameters related to specific graphical view features .
Connect to:
The automatic connection can either be made to the nearest point or
to the nearest branch segment.
Snap to grid
This facility may be used for snapping points to a user defined grid.
The spacing of the grid may be defined here as well.
58 MIKE 11
Graphical View
Drawing Title
Defines a title, or header string, which will be presented on top of
the Graphical View.
Font...
Feature for defining Font settings of the Graphical View (Figure 2.29).
60 MIKE 11
Graphical View
Figure 2.30 The Window menu of the Network editor, Main menu bar
Addiitonally, all open windows are listed in the lower part of the menu
allowing for an easy location of and access to a specific editor file.
Figure 2.31 The Help menu of the Network editor Main Menu bar
The Help menu follows the standards of Windows programs and offers
two options:
Help Topics...
The Help Topics option opens the MIKE Zero on-line help framework
through which any help item of a specific MIKE Zero product can be
accessed.
Please note, that for information on modules available in the actual license
and other license-details it is required to open the MIKE By DHI License
Administration system dialog. The License administration is found at:
‘Start’->‘Programs’->‘MIKE By DHI 2009’->‘License Administration’
Use hotkey; ‘<CTRL>+T’ to open the tabular view from the graphical
view of the network editor.
2.2.1 Points
The geographical coordinates (x,y) of the digitization points in the river
network may edited here.
Note, that the Points in this dialog is only digitizing points used to define
the appearance of the river network in the graphical view. The digitization
points are not used in the calculation at all, so they must not be mixed up
with calculation points.
Figure 2.32 The points property page, Network editor, Tabular view.
62 MIKE 11
Tabular view: Network
Definitions
The X- and Y-coordinate of the selected point may be edited here.
Attributes
Different attributes are available for editing.
Branch
This displays the river branch to which the present point belongs and is
only for verification purposes.
Chainage type
The chainage may either be chosen as user defined or system defined,
where system defined chainage type means that the chainage of a point
will be calculated as a function of the location in the plan coordinate sys-
tem and the distance to up- and downstream points.
User-defined points are opposite, where the user determines exactly which
chainage a point at a certain location must take (e.g. at a bridge where the
surveyer knows the exact chainage of the bridge measured from an
upstream reference point). The points at the upstream and downstream
ends of branches are best set to ‘User-defined’, and unless greater accu-
racy is essential the chainage is best rounded to the nearest metre.
Chainage
If the chainage type is set to user defined the chainage may be edited using
this box.
Type
The type of the point may be set here. Three types are available:
Overview table
An overview of the points is given in the overview table and all points
data can also be edited in this overview table instead of using the above
mentioned group-boxes for points editing.
2.2.2 Branches
The Branches property page (Figure 2.33) contains parameters defining
the river branches and their connections in the river model.
Please note that it is not possible to create new branches within this page.
If new river branches are to be included in the model then it is required to
shift to the Graphical view and digitize a new river branch in the graphical
View. The creation of a river branch in the Graphical View automatically
inserts a new line in the present Branches page.
Also note, that it is not possible to adjust the chainages for the upstream
and downstream point of the river branch in this page! These chainages
are being ported to the Branches page from the Points page where chain-
ages for digitization points belonging to specific river branches are
defined. Therefore, it is strongly recommended in order for the user to
control the exact length of a river branch to define - in the Points property
page - at least the first and the last points in a river branch as ‘User-
Defined’ points such that an exact chainage can be defined at the upstream
and downstream points of a river branch respectively. With this, you will
be certain that MIKE 11 uses exactly a river branch length as you defined
it should be!
64 MIKE 11
Tabular view: Network
Definitions
Branch name
Name of the actual river branch.
Topo ID
Topo ID is an abreviation of ‘Topographical Identifier’, which refers to the
data defined in the Cross section file (xns11). The cross section file offers
the possibility to store several measurement campaigns (e.g. from differ-
ent years or before and after the construction of a bridge structure). The
main purpose of this is therefore, that it is extremely easy to make an alter-
native scenatio simulation once the river model has been properly defined.
As it is normally a sound assumption to say, that the river digitization and
length etc. does not change then it is possible to change the topographical
data used in a simulation just by changing name of the Topo ID for one or
more river branches (and eventually save the nwk11 with a different name
than the original file)
Note that the String defined for the Topo ID must be identical to a Topo ID
present in the cross section file for the actual river branch. If not, the cal-
culation can not be performed.
Upstr. Ch.
The chainage of the first point in the branch.
In case the chainage value of this point should be adjusted it must be
changed in the ‘Points’ property page.
Downstr. Ch.
The chainage of the last point in the branch.
Again, in case the chainage value of this point should be adjusted it must
be changed in the ‘Points’ property page.
Flow Direction
Indicator for how MIKE 11 should interprete the chainage definition with
respect to the natural flow direction. The default assumption of MIKE 11
is to define flow direction positive with increasing chainages looking
downstream of the river.
Hence, a river survey with river chainages starting with 0 at the river
mouth and chainages increasing upstream should be defined as a ‘Nega-
tive Flow’ branch as illustrated in .
Figure 2.34 Illustration of flow direction definiton. Left figure illustrates situation
where increasing chainages follow positive flow direction and
branch should therefore be defined as ‘Positive Flowdir’ whereas
the right part of the figure illustrates the situation where a survey has
increasing chainages from the river outlet and upstream. This situa-
tion requires the river branch to be defiend as ‘Negative Flowdir’.
66 MIKE 11
Tabular view: Network
Maximum dx
Maximum dx is the maximum grid-spacing (or distance) allowed by the
calculation engine between two adjacent h-points (water level calculation
points). At locations where Cross sections are present in the cross section
file, the calculation engine will always create h-points at these location
and in case the Maximum dx is defined with a value smaller than the dis-
tance between cross sections, the calculation engine will automatically
insert a number of additional calculation points by interpolation in
between existing cross sections such, that the maximum distance between
h-points in the calculation will less than or equal to the defined value of
Maximum dx.
Branch Type
– Regular:
A ‘normal’ branch composed by a number of calculation points
defined from cross sections. A minimum of one cross section is
required
– Link Channel:
No cross sections are required. Instead the parameters given in the
Link channel dialog must be specified using the The ‘Edit Link
Channel Parameters...’ button button.
Note that LINK channels are ONLY to be used when conducting
Hydrodynamic simulations.
– Routing:
A special routing feature is available in the network branch defini-
tion in case a special module is present in the MIKE 11 license. This
feature has been developed as part of a dedicated Japanese version
of MIKE 11 and hence only applicable if a special M11-Japanese
module is purchased and added to the MIKE 11 license.
The Routing option in the branch type column of the Branches dia-
log as well as the ‘Routing’ option in the tree view as described in
section 2.4, is only visible to users holding the required special
module in their license.
No cross sections are required. Only the flow is calculated, no water
levels. See section 2.4 Tabular view: Routing (p. 155).
– Kinematic Routing:
Kinematic Routing can be used to model the hydraulics of upstream
tributaries and secondary river branches, where the main concern is
to route water to the main river system. The Kinematic Routing
method does not facilitate the use of structures at Kinematic Rout-
ing branches. Moreover, the method does not account for backwater
effects.
At Kinematic Routing branches, it is possible to run the model
without information on cross-sections. In turn, this indicates
that Kinematic Routing branches can not be used to model a
looped part of a river network. Employment of Kinematic
Routing branches requires that all branches located upstream
of a Kinematic Routing branch are defined in the same way.
– Stratified:
If stratified flow is to be included in the simulation. The branches
for which this vertical resolution is required are to be specified as
stratified and additional information on stratification details (no of
layers etc.) are defined in the HD Parameter file (see Stratification
(p. 349)).
– MIKE 12:
An option to define a branch as a 2-layer flow branch. MIKE 12 is
though no longer further developed and supported and hence, it
may be beneficial to select a Stratified branch instead if layered
flow are to be included in a MIKE 11 simulation.
Connections
The connection point of one branch to another can be specified here. How-
ever, it is recommended that branch connections be defined using the Con-
nect Branch tool in the Graphical Editing Toolbar.
68 MIKE 11
Tabular view: Network
an ‘open culvert’ type of structure. That is, an overflow weir with a length
to include the friction loss component.
Note: Following the internal formulation of the link channel where no true
cross sections are defined in a link channel branch, the usage is restricted
to Hydro-Dynamic simulations ONLY. Thus, for Advection-Dispersion,
Water-Quality and Sediment-Transport calculations the set-up should be
void of link channels and instead connecting branches between main riv-
ers should be defined out of ‘regular’ branches with cross sections defined
and eventual a weir included to describe embankment levels if this is what
should be represented with the connecting branch.
The link channel dialog (see Figure 2.35) is used for specifying all param-
eters appropriate for the link channel e.g. geometry, head los coefficients
etc.
Geometry
The link-channel geometry comprises typically the definition of a longitu-
dinal geometry of the embankment along the river. The geometry is
defined from the following parameters:
Additional Storage:
Link channels do not contain cross sections and do not contribute to the
storage capacity at nodal points where the link connects to a main
branch. The Additional Storage parameter can be used to avoid zero
storage at nodal points to which only link channels and no regular
channels are connected.
Bed resistance
The bed resistance along the length of a link channel can be described
using Manning's M or Manning's n.
70 MIKE 11
Tabular view: Network
Q/h - relations
Free outflow Q/h-relation must be calculated prior to calculation. To cal-
culate the Q/h relationship, specify the number of relationships required
and press the Calculate button. The result of the calculation will appear in
the table. If any of the parameters defining the link channel geometry, or
loss coefficients are changed the Q/h relations must be re-calculated.
Table values in the Q/H relations can be changed if required but general
recommendation is to leave the table as is from the automatic calculation
when pressing the Calculated button. Only the columns ‘US Type’ and
‘DS Type’ can not be changed. These are only indicators of flow type
derived from the calculated flow parameters in the link channel and the
values of these columns are not included in the calculations.
Purpose
The purpose of using alignment lines is to save geo-referenced informa-
tion in the network editor, and to utilize this information to update infor-
mation in the cross section editor. Alignment lines information in the
network file will influence the simulation results only when transferred to
the cross section editor, and such transfer is requested in the cross section
editor. The information in the cross section editor which is subject to be
updated as the result of transferring alignment line information is:
Definition.
An alignment line is similar to a branch in the sense that it is a line going
through an ordered list of points with x- and y- coordinates. The following
list of types of alignment lines are available:
User Interface
Figure 2.37 shows the property page for alignment lines. Each alignment
line is shown as a row in the overview in the bottom of the dialog, and the
x- and y-coordinates of the points along the actual line (the line in the row
being high lighted in the overview) is shown in the details in the top of the
dialog.
72 MIKE 11
Tabular view: Network
The x- and y-coordinates for the points along the alignment lines can be
edited in three ways: 1) Using the tools available in alignment lines tool
bar in the graphical view (see Tool Bar for Alignment Lines (p. 39)). 2)
Editing the numbers in the tabular view. 3) Using the File menu to import
the coordinates from a text file.
Once the alignment data are added, the information is ready to be trans-
ferred to the cross section editor.
74 MIKE 11
Tabular view: Network
2.2.4 Junctions
The junctions feature is part of the quasi two dimensional steady state with
vegetation module.
Details
Name, Name2 and Name3: The river name of the three rivers meeting at
the junction.
Angle 1 and Angle 2: The direction angle of channel 1 and 2 with respect
to channel 3.
Distance along channel 3 (D): The distance along channel 3 at which the
local water depth should be used for the determination of the water level
in the downstream points of channel 1 and 2.
76 MIKE 11
Tabular view: Structures
ual structures in the composite structure is defined at the exact same loca-
tion (same River Name and Chainage). This example is illustrated in
Figure 2.42.
Regular
Regular structures are internal structures that specify a flow in the river
branch typically based on up- and down- stream flow conditions.
Side Structure
Side structures are special type of a regular structure. Side structures takes
water out of the river network at the location where the side structure is
defined.
– Head losses are excluded from the flow calculation for side struc-
tures by setting internally the flow areas of the cross sections in the
side structure branch as extremely large (1e20 m2) and hence, the
velocity component in these sections (V=Q/A) is a neglible number.
– For culverts and weirs, an assumption of free out- and over-flow
respectively is implemented through an internally defined water
level boundary condition for the outfflow point of the artificial side
structure branch. The boundary condition is defined as a water level
boundary with a low water level.
– Naming convention used for the artificial side branch are:
"SS_"<original branch name>"_"<original chainage>
– In case of multiple side structures defined at the same location,
these will all be included in the same side branch as a ‘composite
side structure’.
– A side structure branch is a simple h-Q-h branch with two cross sec-
tion h-points in each end of the branch and a structure in the centre
Q-point.
Up- and downstream cross sections in the Side Structure Branch are
automatically assigned with a shape identical to the main river cross
section located immediately downstream of the side structure loca-
tion. Upstream cross section of the Side Structure branch is identi-
cal to the section in the main river, but the downstream section has
been shifted 10 m down in order to obtain a free outflow condition
in the Side Structure. See Figure 2.43.
Figure 2.43 Cross sections as defined in the artificial Side Structure branch. Left
part illustrates the cross section in main river downstream of the
chainage defining the Side Structure. Middle part illustrates the
upstream cross section at chainage 0 of the side structure branch
(=copy of main branch section) and right part illustrates the down-
stream cross section of the Side Structure Branch, which is identical
shape as the other sections but only shifted -10 m.
78 MIKE 11
Tabular view: Structures
The downstream cross section of the side structure branch has a user-
defined, additional storage area assigned, which assigns the reservoir data
specified as a level dependent additional storage function (ref. section
3.2.2 Processed data, Tabular View) defined in the downstream (outflow)
h-point of the side structure branch. The significiant difference between a
structure defined as ‘Side Structure’ and ‘Side Structure with Reservoir’ is
therefore, that in the latter case, the water conveyed through the side struc-
ture is stored in the system and can be analyzed, whereas for the structure
type without the Reservoir, the water simply leaves the river model and
disappears completely.
Figure 2.45 Dialog for specifying Additional storage for Side Structures with
Reservoir
– Elevation-area
The Elevation-Area relation specified in the table are used directly
as additional storage area as a function of the water level in the res-
ervoir.
– Elevation-volume
Additional storage area is derived from the Elevation-Storage rela-
tion following the relation:
∆V Table
FlA i + 1 = ------------------
- – FlA i ⋅ 2 (2.1)
∆H Table
where:
80 MIKE 11
Tabular view: Structures
– Head losses are excluded from the flow calculation for side struc-
tures with Reservoir by setting internally the flow areas of the cross
sections in the side structure branch as extremely large (1e20 m2)
and hence, the velocity component in these sections (V=Q/A) is a
neglible number.
– The boundary introduced through the inclusion of the side weir is
internally specified as a no flow boundary. Flow may be included
through specification of point source boundaries at the boundary.
– More side structures with reservoir in the same location will result
in one side branch for each.
– The side branches will be named with the naming convention:
"SSPR"<original branch name>"_"<original chainage>"_"<Struc-
ture ID>.
– ‘Apply Coordinates’. The Apply coordinates offers a possibility for
defining coordinates for the reservoir - which in the case of the arti-
ficially generated side structure branch is identical to the end point
of the side structure branch. Hence, the coordinates defines the
appearance of the side structure branch and a length of the branch
will be calculated according to the specifiec coordinates.
– The cross sections applied in the artificially generated Side Struc-
ture Branch with Reservoir are defined such that the upstream cross
section is a copy of the section in the main river downstream of the
chainage defining the side structure and the downstream cross sec-
tion is a defined as a triangular cross section shape as illustrated in
Figure 2.46.
Figure 2.46 Cross sections as defined in the artificial Side Structure with
Reservoir branch. Left section illustrates cross section in main river
downstream of chainage defined for Side Structure. Middle section
illustrates the upstream section of the side structure branch (copy of
main branch section) and right section illustrates triangular section
applied at downstream h-point of Side Structure Branch
Figure 2.47 Illustration of Structure Geometry definition. The blue area which
potentially conveys water over/through the structure is the geometry
which must be defined.
82 MIKE 11
Tabular view: Structures
For structures which can be plotted the structure dialog contains a section
named “Graphic” with a “Plot...” button. Clicking this will open a separate
graphical window showing two plots. The upper plot being the upstream
cross section and the structure, and the lower plot being the downstream
cross section and the structure. Initially, the centre of the structure is hori-
zontally aligned to marker 2 of the cross section.
Clicking the right-mouse button in the graphical view the user will be
given the following options in addition to the regular zoom and pan
options:
z Settings...:
This will open a separate dialog where the graphical properties of cross
section, structure etc. can be changed.
z Use same axis scaling U/S and D/S:
Selecting this the two plots will be shown with the same scaling of ver-
tical and horizontal axis
z Show next:
The next structure (according to the list in the tabular view) neighbour-
ing cross sections will be shown instead of the current.
z Show previous:
The previous structure (according to the list in the tabular view) neigh-
bouring cross sections will be shown instead of the current.
z Pan structure left:
The structure will be moved to the left in the graphical view.
z Pan structure right:
The structure will be moved to the right in the graphical view
z Increase pan step:
The distance by which the structure moves left/right is increased by a
factor 2.
2.3.6 Weirs
The Weirs page in MIKE 11 includes different types and formulas for
overflow weirs. The weirs property page is presented in Figure 2.49.
84 MIKE 11
Tabular view: Structures
Location
– River Name:
Name of the river branch in which the weir is located.
– Chainage:
Chainage at which the weir is located.
– ID:
String identification of the weir. ID will be used in structure results
and hence makes it easier to identify a specific structure if setup
contains many structures, or if multiple structures are defined at the
same location. It is recommended always to give the weir an ID.
– Type:
The lcation type may be Regular, Side Structure or Side Structure +
Reservoir. See 2.3.3 Structure Types definition (p. 77) for details.
Attributes
– Type:
Type of Overflow Weir formulation:
Special Weir:
Special Weir is chosen in the situation where a Broad Crested weir
formulation is not vaild. For this weir type, the Q/h relationship
table for free overflow must be specified manually by the user.
Weir Formula 1:
A standard weir expression reduced according to the Villemonte
formula is applied. See the Reference Manual, under weir formula.
– Valve:
None:
No valve regulation applies.
Geometry
Geometry definitions applies only for broad crested weir and special weir
types.
Type of Geometry:
– Level-Width:
The weir geometry is specified as a level/width table relative to the
datum.
– Cross Section DB:
The weir geometry is specified in the cross section editor. A cross
section with a matching branch name, Topo ID and chainage must
exist in the applied cross section file. The Topo ID is assumed to be
the same as that specified in the Branches Property page, see Topo
ID (p. 65).
Datum:
Offset which is added to the level column in the level/width table.
Level/Width table:
Weir shape defined as levels and corresponding flow widths. Values in
the levels column must be increasing. See also section 2.3.4 Structure
Geometry definition (p. 82).
86 MIKE 11
Tabular view: Structures
88 MIKE 11
Tabular view: Structures
– Special weir:
Unlike for the broad crested weir, where Q/h relations are calcu-
lated automatically, the user must manually specify Q/h relations
corresponding to free overflow conditions for a Special weir. These
must be specified for both positive and negative flows.
Note that Q/h relations must be re-calculated if any changes are made to
the weir or the cross sections up- or downstream have been altered. Fur-
ther, since a weir in MIKE 11 is defined as a structure causing a contrac-
tion loss and subsequently an expansion loss, some constraints are placed
on the geometry of a broad crested weir. The geometry of the weir must be
such that the cross sectional area at the weir is less than the cross sectional
area at both the upstream and the downstream cross section for all water
levels!
2.3.7 Culverts
Branch name
River Name:
Name of the river branch in which the weir is located.
90 MIKE 11
Tabular view: Structures
Chainage:
Chainage at which the weir is located.
ID:
String identification of the culvert. It is used to identify the culvert if
there are multiple structures at the same location. It is recommended
always to give the culvert an ID.
Type:
The lcation type may be Regular, Side Structure or Side Structure +
Reservoir. See 2.3.3 Structure Types definition (p. 77) for details.
Attributes
Upstream Invert:
Invert level upstream of the culvert.
Downstr. Invert:
Invert level downstream of the culvert.
Length:
Length of the culvert.
Manning’s n:
Manning’s bed resistance number along the culvert (for calculation of
the Friction loss present in a structure with a certain length).
No. of Culverts:
Number of culvert cells.
Valve Regulation:
– None:
No valve regulation applies.
– Only Positive Flow:
Only positive flow is allowed, i.e. whenever the water level down-
stream is higher than upstream the flow through the structure will
be zero.
– Only Negative Flow:
Only negative flow is allowed, i.e. whenever the water level
upstream is higher than downstream the flow through the structure
will be zero.
Section Type: Closed or Open.
Geometry
The cross sectional geometry of a culvert can be specified as:
– Rectangular:
The width and height specify the geometry.
– Circular:
The geometry is specified by the diameter.
– Irregular Level-Width Table:
The geometry is specified using a level/width table. Values in the
level column must be increasing.
– Irregular Depth-Width Table:
The geometry is specified using a depth/width table. Values in the
depth column must be increasing.
– Section DB:
The geometry is specified by a cross section. A cross section with
the same branch name, Topo ID and chainage must exist in the
cross section file. The Topo ID is assumed to be the same as speci-
fied in Topo ID (p. 65).
For further details on structure geometry definition, see section 2.3.4
Structure Geometry definition (p. 82).
Flow Conditions
Once the above parameters and the desired number of Q/h relations have
been filled in the button Calculate Q/h relations can be pressed. The result
of the calculation will appear in the table. If any of the parameters defining
the culvert is changed the user should remember to re-calculate the Q/h
relations. In order to compute the Q/h relation, the nearest upstream and
downstream cross section are used. The cross sections must be within the
distance maximum dx (Maximum dx (p. 67)) defined for the branch in
question. The Q/h relation can not be calculated unless the cross sections
are defined. It is also necessary that the Simulation File is open in order to
load the cross section data from a cross section file.
The Q/h relations are given as Q/y relations (where y is depth above
invert).
The Q/y relations table also shows the type of flow occurring for the actual
levels. Please note that the flow-types are only indicators derived from
92 MIKE 11
Tabular view: Structures
MIKE 11. It will have no effect on the calculation results to change the
flow types.
The possible types are:
2.3.8 Pumps
Location
– Branch Name: Name of the river branch in which the pump is
located.
– Chainage: Chainage at which the pump is located.
– ID: String identification of the pump. Used for identification of the
pump in case of multiple structures at the same location. Specifica-
tion of pump ID is recommended.
– Type:
The lcation type may be Regular, Side Structure or Side Structure +
Reservoir. See 2.3.3 Structure Types definition (p. 77) for details.
Control Parameters
– Start Level: Water level at the inflow that activates the pump.Note
that for pumps with internal outlet the inflow is situated at the previ-
ous h-point (previous with regard to chainage) in case of positive
discharge and at the next h-point with regard to chainage) in case of
negative discharge. The sign of the discharge follows from the spec-
ifications made under Pump Data.
– Stop Level: Water level at which the pump starts closing down.
94 MIKE 11
Tabular view: Structures
Pump Data
– Specification Type:
Fixed Discharge:
Pump rate independent of the local water head expect for the
start/stop control.
Tabulated Characteristic:
Pump rate controlled by specified characteristic (Q-dH-curve) and
the water level difference between upstream water level and outlet
level/downstream water level.
– Discharge:
Pump rate when applying "Fixed Discharge".
– Outlet Level:
The outlet level is used when defining tabulated (dh-Q) characteris-
tics for the pump. It then represents the outlet level on the down-
stream side of the pump.
– Q-dH-curve:
Q-dH-characteristic of the pump. The discharge is determined
through interpolated look up in the table specified. The dH used for
the look-up is given as the difference between Downstream and
Upstream water level in case of submerged outlet and as the differ-
ence between upstream level and outlet level in case of free outlet.
The shift between the two is fully dynamic allowing an outlet to
change from being free to submerged and vice versa.
2.3.9 Bridges
MIKE 11 offers a number of approaches when modelling flow through
bridges. The approach to choose should be based on the assumptions for
the different methods and the requirements of the modelling.
The bridge modelling approaches can be divided into pure free flow meth-
ods and methods which may be combined with submergence/overflow
methods. The pure free flow methods can be further sub-divided into
methods for piers and methods for arches.
96 MIKE 11
Tabular view: Structures
For all bridges there are some common settings consisting of:
98 MIKE 11
Tabular view: Structures
– Name: The river name. Used along with the chainage to identify the
location of the bridge in the network. River name and branch name
are identical,
– Chainage: The location of the bridge in the river. The river cannot
be located at a location where a cross section exists.
– Bridge ID: String identification of the bridge. The string identifica-
tion for the bridge is used in the output, but also for the preprocess-
ing of the bridges. It is important that every bridge in the set-up has
a unique non-empty Bridge ID. The latter ensures that bridges can
be uniquely identified by the calculation kernel of MIKE 11.
– Method: Pull down menu for selecting the free flow method used in
calculating the flow through the bridge. The methods are briefly
described above. Please refer to the reference manual for a more
thorough decsription of the underlying assumptions and governing
equations.
– Options: Some of the bridge methods offer a number of additional
options. These options may include e.g. piers, skewness etc. To
include an option simply select the appropriate box.
– Geometry and Loss factors: Thorugh the use of the Edit button the
user can access the interface for for entering geometric and loss fac-
tor parameters. The Detail button gives access to the edit and or
review of the tables used for determining loss coefficients.
– Graphic: Not available in present release.
– Submergence: (Pressure flow) Available if the Submergence
checkbox is marked (See Options). The user must select one of the
methods: Energy Equation, FHWA WSPRO or MIKE 11 culvert.
For details on FHWA WSPRO see Submergence, FHWA WSPRO.
For a description of culverts please refer to the appropriate section
(2.3.7). The energy equation method is equivalent to the one used
for free surface flow. Note that only some methods allow the sub-
mergence option.
– Culvert no: When choosing a MIKE 11 culvert details of the cul-
vert structure are in the culvert menu (See section 2.3.7). Culvert no
is chosen as the number marked in the overview box in the culvert
menu.
– Bridge level (bottom): The bottom level of the bridge. The geome-
try of the lower part of the bridge consists of the data in bridge cross
sections combined with the bridge level bottom.
Figure 2.57 Overview, Branch2 has three bridge openings - 2, 3 and 4. Marked
in the right part of the overview window.
The overview is split in two. The left hand side gives an overview of the
physical bridges in the set-up and the right the different openings. If a
bridge consists of one opening there will be one entry on the left and one
on the right for that particular bridge. If a bridge has multiple openings
there will only be one entry on the left for that particular bridge and a
number of entries on the right corresponding to the number of openings in
the bridge.
To add a bridge to a set-up simply place the cursor in the filed of the last
bridge and hit the TAB key until a new line appears. Similarly for addi-
100 MIKE 11
Tabular view: Structures
tional openings in the same bridge. As an alternative the insert key on the
keyboard may be used. To sumarize the layout:
– Overview: Left part shows River name, Chainage and Bridge ID.
Right part show methods for the bridge openings.
– Multiple waterway openings: If working with multiple waterway
openings all multiple waterway openings are marked when the
bridge is activated. (See Figure 2.57). In order to add additional
openings, mark a row in the right part of the overview window and
press insert on the keyboard.
The loss / adjustment factor tables are viewed by pressing the Details but-
ton. The default loss factor tables are generated by pressing the Edit but-
ton.
When having default unmarked for a loss factor changes in the loss factor
table will be saved. If default is marked changes will not be saved after
pressing edit.
In the loss factor tables the user can create more columns and rows. Plac-
ing the cursor in the last column (right end) and pressing the arrow button
on the keyboard will create a new column. Pressing the tab button on the
keyboard when having the cursor in last bottom cell creates an additional
line.
The location of the cross sections outside the bridge should be so that any
potential contraction or expansion loss is taken into account. In other
words the optimal location is where the stream lines are parallel prior to a
contraction and post a possible expansion. As a rule of thumb the distance
between the bridge and the cross sections should be of the order one open-
ing width (see Figure 2.58).
Figure 2.58 Location of up- and downstream cross section. 1: Upstream river
cross section. Defined in the cross section editor. 2: Upstream
bridge cross section. Defined in the network editor, bridge geometry.
3: Downstream bridge cross section. Defined in the network editor,
bridge geometry. 4: Downstream river cross section. Defined in the
cross section editor.
– Submergence
– Overflow
– Skewness, Used when the embankments is not perpendicular to the
approaching flow.
– Eccentricity, Used when the bridge opening is eccentrically located
in the river.
– Multiple waterway opening
– Asymmetric opening, Used for individual definition of left and
right abutments.
– Spur dykes
– Piers/piles
Available options for USBPR Bridge:
– Submergence
– Overflow
102 MIKE 11
Tabular view: Structures
– At level z: Only for FHWA WSPRO, opening type II, III, and IV.
Enter the level for witch the Waterway length is measured.
Figure 2.60 Definition sketch of type I opening, vertical embankments and verti-
cal abutments, with or without wingwalls (after Matthai).
Figure 2.62 Definition sketch of type III opening, sloping embankments and
sloping abutments (spillthrough) (after Matthai).
104 MIKE 11
Tabular view: Structures
Figure 2.63 Definition sketch of type IV opening, sloping embankments and ver-
tical abutments with wingwalls (after Matthai).
– Slope: If the slope check box is marked the only the upstream
bridge cross section must be defined. The downstream cross section
is generated be copying the upstream cross section and adding the
slope defined in the slope edit box. Upstream bridge cross section
correspond to section 2 and downstream bridge cross section corre-
spond to section 3. (See Figure 2.58).
– Datum: The water level datum is added to the Z values in the Cross
section table.
– X: Horizontal values for the cross section. Note that the x-values
are evaluated with the up and downstream cross section. As a result
it is important that the four cross sections (See Figure 2.58) are
placed correct in respect to the x-values.
– Z: Vertical level of the cross section point.
– Resistance: Additional resistance in the cross section point. 1 is
resistance corresponding to the manning number.
– Marker: Define cross section embankments and abutments (See
Figure 2.64).
– Cross section definitions: Cross sections 2 and 3 must be defined as
illustrated by the blue arrows in Figure 2.64 - i.e. along the flood-
plain then to the underside (soffit) of the bridge then down the abut-
ment to the channel and similar to the other side. Markers 1 and 3
correspond to floodplain limits and 4 and 5 are bottom of left and
right abutments respectively
Geometry and loss factors are defined for each opening when working
with multiple waterway openings (see Figure 2.57). The position of each
opening and the corresponding stagnation points are defined from the
stagnation point value (if not default) and from the horizontal values in the
bridge cross sections. (See Figure 2.65)
106 MIKE 11
Tabular view: Structures
Loss factor:
108 MIKE 11
Tabular view: Structures
Average depth k y II m or M ( Y a + Y b ) ⁄ ( 2b )
Wingwall kθ IV m or M
In the Loss Factor menu the user can choose to use m or M as axis in the
tables.
Where:
e Eccentricity.
Skewness ∆k φ* M or m
Piers ∆k p* M or m
110 MIKE 11
Tabular view: Structures
In the Loss Factor menu the user can choose to use m or M as axis in the
tables.
Where:
e Degree of eccentricity.
The details of the bridge geometry and location are inserted in the appro-
priate boxes:
BRIDGE_FROUDE_CRITERIA.
Since the equations used for arch bridges are only valid for free flow sub-
mergence and and overflow cannot be selected. Geometry and loss factors
are viewed by pressing the Edit button under Geometry and Loss factors.
112 MIKE 11
Tabular view: Structures
Please note, that Invert level is not specified in the geometrical definition
of the Arch Bridge. Arch Bridge flow calculation depends on the water
levels and depths at up- and downstream side of the bridge, but these are
calculated from the cross sections up- and downstream respectively and
hence is not necessary to define in the parameters page.
The details of the bridge geometry and location are inserted in the appro-
priate boxes:
114 MIKE 11
Tabular view: Structures
Note: If the Froude number downstream of the piers is greater than the cri-
teria (default 0.6) the effect of bridge piers using D’Aubuisson’s formula
is ignored. The criteria value may be changed in the Mike11.ini file by set-
ting the variable:
BRIDGE_FROUDE_CRITERIA.
– Submergence
– Overflow
Geometry and loss factors are viewed by pressing the Edit button under
Geometry and Loss factors.
Figure 2.72 Nagler and Yarnell Bridge piers geometry property page.
Figure 2.73 Nagler and Yarnell Bridge piers Loss factor property page.
116 MIKE 11
Tabular view: Structures
2.3.10 Regulating
118 MIKE 11
Tabular view: Structures
Type:
The lcation type may be Regular, Side Structure or Side Structure +
Reservoir. See 2.3.3 Structure Types definition (p. 77) for details.
Regulation Function
This section is only available when the regulation is specified as a h/Q
function. The function f(J2) is specified in the Regulation Function table
as a series of factors for corresponding values of J2.
Location
Branch name: Name of the river branch in which the structure is located.
Type:
The lcation type may be Regular, Side Structure or Side Structure + Reser-
voir. See 2.3.3 Structure Types definition (p. 77) for details.
Attributes
Gate Type
– Overflow: This gate type corresponds to a variable crested weir.
– Underflow: This gate type corresponds to a vertical sluice gate.
– Discharge: This gate type corresponds to a pump.
– Radial gate: This gate type corresponds to a Tainter gate. In con-
trast to the other gate types a radial gate does not need any informa-
tion about head loss factors. Instead a number of radial gate
parameters must be entered, see Radial Gate Parameters (p. 122).
120 MIKE 11
Tabular view: Structures
Number of gates
The number of identical gates is entered here. This variable is used when a
series of identical gates are simulated.
Underflow CC
This is the contraction coefficient used for underflow gates only. Default
value is 0.63.
Gate width
The width of the gate. Not applicable for gates of the Discharge type.
Sill level
The level of the sill just upstream of the gate. Not applicable for gates of
the Discharge type.
Max. Speed
This variable defines the maximum allowable change in gate level per.
time. (If a discharge gate is chosen the variable defines the maximum
allowable change in discharge per. time). This variable is introduced
because the control strategy defining the variation of the gate level can
result in very rapid changes in gate level. This is probably not realistic,
further it can create instabilities in the computation.
Initial Value
If the Initial Value checkbox is checked the value specified will be used as
initial value
Max Value
If the Max Value checkbox is checked the value specified will interpreted
as the highest possible gate level or in case of a discharge structure the
highest possible pump discharge.
Figure 2.78 The control structure property page when a radial gate has been
selected.
Tune Factor
Discharge calibration factor. The Tune Factor is applied as a multiplica-
tion factor solely on the part of the discharge that flows below the radial
gate (not applied to eventual overtopping gate flow).
122 MIKE 11
Tabular view: Structures
Height
Height above sill level of the overflow gate crest when the gate is closed,
see Figure 2.79.
Radius
Radius of gate, see Figure 2.79.
Trunnion
Height above sill levle of the centre of gate circle, see Figure 2.79.
Tran. Bottom
The Tran. Bottom parameter is used to define the level where a free flow
regime for the gate changes to a transition zone between free flow and
submerged flow regime. Corresponds to yTran,Bottom as defined in Hydrau-
lic Aspects - Radial Gates in the reference manual.
Tran. Depth
The Tran. Depth arameter is used to define the depth (or height) of the
transition zone between free flow and submerged flow regimes and conse-
quently is used to define the level where the flow regime changes to sub-
merged flow. Corresponds to yTran,Depth as defined in Hydraulic Aspects -
Radial Gates in the reference manual.
124 MIKE 11
Tabular view: Structures
With this gate type the flow under the gate is divided into four
flowregimes the choice of which depends on the upstream and down-
stream water level. The four flow regimes are:
Additionally flow over the top of the sluice gate is taken into account
when the water level up stream and/or downstream exceeds the gate level
plus the gate height.
In short the high limit and low limit parameters are used to smoothen the
transitions between fow regimes and coefficient a and exponent b are
pararameters for the flow equations. Please refer to the reference manual
for further details.
Control Definitions
The way the gate operation is calculated is determined from a control
strategy. A control strategy describes how the gate level depends on the
value of a Control Point. For a specific gate it is possible to choose
between an arbitrary number of control strategies by using a list of ‘if’
statements. For each of these statements it is possible to define an arbitrary
number of conditions that all must be evaluated to TRUE if the ‘if’-state-
ment is to be evaluated to TRUE. It is hereby made possible to use differ-
ent operating policies depending on the actual flow regime, time etc.
1 The conditions that must be fulfilled for the strategy to be executed and
2 the control strategy itself.
The control strategy itself is a relationship between an independent varia-
ble (the value of the Control Point) and a dependent variable (the value of
the Target Point).
The reason for using the concept ‘Target Point’ and not just call it gate
level is as follows:
As seen in Figure 2.77 or Figure 2.78 the control definitions section con-
sists of a table. Each line in this table represents the main parameters of an
‘if’-statement.
Priority
As mentioned under Control Definitions (p. 125) it is possible to make
Mike11 choose between an arbitrary number of control strategies. These
control strategies are organised using a list of ‘if’ statements. The control
strategy belonging to the first of these statements that are evaluated to
TRUE will be executed. It is thus of importance for the user to define
which ‘if’-statement that are evaluated first, second, third and so on. This
is enabled by the priority field. In this the user defines the priority of the
‘if’-statement by writing an integer number. By default the first line in the
table will have priority equal to one, the second line will have priority
equal to two and so on. Note that the ‘if’-statement with the lowest prior-
ity always will be evaluated to TRUE. This is because this statement is
connected to the default control strategy that will be executed when all
other ‘if’-statements are evaluated to FALSE.
Calculation Mode
– Tabulated: This is the default calculation mode, which determines
the value of the gate level directly (gate discharge in case of a dis-
charge gate).
126 MIKE 11
Tabular view: Structures
n T T n n
u = α 1 K 1 + -----s + ----d- { y ref – y } (2.2)
T i T s
The variables K, Td, Ti, α1, α2 and α3 are entered by the user, see
Iteration / PID (p. 141). The rest is calculated by Mike11.
– Fully open: If this calculation mode is chosen the gate is fully open.
For an overflow gate this means that the gate level will equal the
Sill Level specified by the user. For underflow gates and radial
gates if will correspond to the Max Level defined by the user. For a
pump this corresponds to the maximum pump capacity specified as
Max Value.
– Close: If this calculation mode is chosen the gate is closed. For an
overflow gate this means that the gate level will equal the Max
Level specified by the user. For underflow gates and radial gates if
will correspond to the Sill Level defined by the user. For a pump
this corresponds to a pump discharge equal to zero.
– Unchanged: If this calculation mode is chosen the gate level or the
pump discharge will remain unchanged.
– Change with: If this calculation mode is chosen the gate level or
the pumps discharge will change with the amount specified in the
Value column in the Control Definitions data section.
– Set equal to: If this calculation mode is chosen the gate level or the
pumps discharge will equal the value specified in the Value column
in the Control Definitions data section.
Control Type
Here the type of control point is chosen.
128 MIKE 11
Tabular view: Structures
– Day of week: Integer expressing the day of the week at the time of
calculation. Monday corresponds to one, tuesday to two and so on.
– Day of month: Integer expressing the day of the month at the time
of calculation.
– Month of year: Integer expressing the month of the year. January
corresponds to one, February to two and so on.
– Year: The year given as an integer value.
– Time after start: This control type is used in control strategies with
a gate operation that can not be interrupted. An example could be a
gate that closes from fully open to fully closed during half an hour
when the water level downstream reaches a certain level. Because it
is not known when the closing procedure is initiated it is not possi-
ble to describe it using a time series. Instead the gate level is
described as a function of time measured relative to the time at
which the procedure was initiated, i.e. the first value of the Control
Type ‘Time after start’ MUST always equal zero. If it is decided to
operate a gate using this control type no other operating policies can
be invoked before the actual gate operation has finished. In the
example this means that no other operating policies can be used
during the half hour it takes to close the gate.
– Concentration: A concentration of any compound.
– Hups: Water level just upstream of the structure.
– Hdws: Water level just downstream of the structure.
– Qups: Discharge just upstream of the structure.
– Qdws: Discharge just downstream of the structure.
– Vol: The volume of water in a point.
– Volups: The volume of water just upstream of the structure.
– Voldws: The volume of water just downstream of the structure.
– BranchVol: The volume of water in a certain part of a river branch
– Depth: The depth in a point.
– Area: The area in a point.
– DepthUps: The depth just upstream of the structure.
– DepthDws: The depth just downstream of the structure.
– AreaUps: The area just upstream of the structure.
– AreaDws: The area just downstream of the structure.
– ThisGate dh: The water level difference across the structure.
Target Type
Here the type of the target point is chosen.
130 MIKE 11
Tabular view: Structures
Type of scaling:
– None: This is the default value. When this is chosen no scaling of
the value of the target point will take place.
– Scaling with internal variable: When this is chosen the value of
the target point will be scaled with the value of a specified internal
variable. See Control Strategy (p. 139) for a list of the internal vari-
ables that can be used as scaling factors.
– Scaling with time series. When this is chosen the value of the tar-
get point is scaled with a factor taken from a time series.
Value:
Specific value of a structure change or setting. This column should be
edited only if the calculation mode is chosen as either Change With or Set
Equal to.
Details...
When pressing the details button a new dialog pops up. This is used to
enter the necessary details in defining the operating rules for control struc-
tures in Mike11.
There are four property pages: ‘Logical Operands’, ‘Control- and Target
point’, ‘Iteration/PID’ and ‘Control Strategy’.
Logical Operands
Note that it is not necessary to enter any logical operands for the ‘if’-state-
ment with the lowest priority. The control strategy belonging to this ‘if’-
statement is the default strategy and will always be executed when all
other ‘if’-statements with higher priority are evaluated to FALSE. As an
example think of a gate where the gate level is a known function of time.
In this case only one control strategy is needed. The control type will be
‘Time’ and the target type will be ‘Gate Level’. Calculation mode is cho-
sen as ‘Direct Gate Operation’. It is not necessary to enter any logical
operands because when only one control strategy is specified this strategy
will have the lowest priority.
132 MIKE 11
Tabular view: Structures
Chainage LO1: This field contains the chainage of the Logical Operand.
Name LO1: This field is used only when LO Type equals ‘Gate Level’,
‘Q_Structure’ or ‘TSLGLC’. Then this field holds the structure ID of the
relevant structure.
134 MIKE 11
Tabular view: Structures
Comp. No.: This field is used only when LO Type equals ‘Concentration’.
The field holds the number of the relevant component.
Branch Name LO2: This field is only used if the LO Type equals ‘dH’
(H1 - H2) or ‘dQ’ (Q1 - Q2). The field holds the name of the branch in
which the H2 or Q2 should be found.
Chainage LO2: This field is only used if the LO Type equals ‘dH’ (H1 -
H2) or ‘dQ’ (Q1 - Q2). The field holds the name of the chainage of the H2
or Q2 point.
Sign: Here the operator used in the logical expression is used. The user
can choose between {<, <=, >, =>, =, <>}.
Use TS-value:
Time Series File: This field holds information about the relevant time
series file in case that the Use TS-value is chosen as ‘Yes’ or in the situa-
tion where the LO Type is chosen to be ‘TS-Scalar’.
Time Series Item: This field holds the name of the item chosen from the
time series file that was selected in the Time Series File field.
Sum of Discharges
i=n
Branch: The name of the branch with the grid point / structure from
which the discharge should be taken.
Struc. Name: In case Structure Discharge has been chosen as the type the
structure ID must be given here.
136 MIKE 11
Tabular view: Structures
Control Type: Here the type of Control Point is chosen. This field is
linked to the Control Type field described in Control Type (p. 128).
Branch, Control Point 1: This field contains the name of the branch with
the control point.
Chainage, Control Point 1: This field contains the chainage of the con-
trol point.
Name, Control Point 1: This field is used only when Control Type equals
‘Gate Level’ or ‘Q_Structure’. The field holds the structure ID of the rele-
vant structure.
Comp. No., Control Point 1: This field is used only when Control Type
equals ‘Concentration’. The field holds the number of the relevant compo-
nent.
Branch, Control Point 2: This field is only used if the Control Type
equals ‘dH ‘(H1 - H2) or ‘dQ’ (Q1 - Q2). The field holds the name of the
branch in which the H2 or Q2 should be found.
Chainage, Control Point 2: This field is only used if the Control Type
equals ‘dH’ (H1 - H2) or ‘dQ’ (Q1 - Q2). The field holds the name of the
chainage of the H2 or Q2 point.
Target Point Type: Here the type of target point is chosen. This field is
linked to the Target Type field described in Target Type (p. 130).
Branch, Target Point 1: This field contains the name of the branch with
the Target point.
Chainage, Target Point 1: This field contains the chainage of the target
point.
Name, Target Point 1: This field is used only when Target Type equals
‘Gate Level’ or ‘Q_Structure’. Then this field holds the structure ID of the
relevant structure.
Comp. No., Target Point 1: This field is used only when Target Type
equals ‘Concentration’. Then this field holds the number of the relevant
component.
Branch, Target Point 2: This field is only used if the Target Type equals
‘dH’ (H1 - H2) or ‘dQ’ (Q1 - Q2). Then this field holds the name of the
branch in which the H2 or Q2 should be found.
Chainage, Target Point 2: This field is only used if the Target Type
equals ‘dH’ (H1 - H2) or ‘dQ’ (Q1 - Q2). Then this field holds the name of
the chainage of the H2 or Q2 point.
Time Series File: This field holds information about the relevant time
series file in case that the Control Type is chosen as ‘Time’. If the button
to the right of this field is pressed it is possible to browse for the file. At
the same time the relevant item in the time series file can be selected.
Time Series Item: This fields hold the name of the item chosen in the
time series file that are selected in the Time Series File field.
138 MIKE 11
Tabular view: Structures
Control Strategy
Here the relationship between the value of the Control Point and the value
of the Target Point are entered. This is done in the table on the left side of
the property page.
Also the information about scaling of the target point are entered here. The
Type of Scaling field is linked to the Type of Scaling field described in
Target Type (p. 130). Below this field there are two sections: A ‘Scaling,
Internal Variables’ section and a ‘Scaling, Time series’ section. Both of
these will be greyed out if ‘None’ is chosen as scaling type.
If Type of Scaling is chosen as ‘Scale with time series’ a dfs0 file contain-
ing the relevant time series can be allocated by pressing the button to the
right of the ‘Time Series File’. At the same time the relevant item in the
dfs0 file can be selected.
Variable Type: The type of internal variables that can be used are:
– V: Velocity in a point.
– Gate level: The level of a gate.
– Concentration: A concentration of any compound.
– Hups: Water level just upstream of the structure.
– Hdws: Water level just downstream of the structure.
– Qups: Discharge just upstream of the structure.
– Qdws: Discharge just downstream of the structure.
– Vol: The volume of water in a point.
– Volups: The volume of water just upstream of the structure.
– Voldws: The volume of water just downstream of the structure.
– BranchVol: The volume of water in a certain part of a river branch
– Depth: The depth in a point.
– Area: The area in a point.
– DepthUps: The depth just upstream of the structure.
– DepthDws: The depth just downstream of the structure.
– AreaUps: The area just upstream of the structure.
– AreaDws: The area just downstream of the structure.
– ThisGate dh: The water level difference across the structure.
– ThisGate Q_Structure: The flow through the structure.
– ThisGate Gate Level: The gate level of the gate.
The use of references to the points up- and downstream of the structure
and reference to the ThisGate enables a faster editing of the strategies
because no location data must be entered
Branch, Scale Point 1: This field contains the name of the branch with
the scaling point.
Chainage, Scale Point 1: This field contains the chainage of the scaling
point.
Name, Scale Point 1: This field is used only when Variable Type equals
Gate Level or Q_Structure. The field holds the structure ID of the relevant
structure.
140 MIKE 11
Tabular view: Structures
Comp No, Scale Point 1: This field is used only when Variable Type
equals Concentration. The field holds the number of the relevant compo-
nent.
Branch, Scale Point 2: This field is only used if the Variable Type equals
dH (H1 - H2) or dQ (Q1 - Q2). The field holds the name of the branch in
which the H2 or Q2 should be found.
Chainage, Scale Point 2: This field is only used if the Variable Type
equals dH (H1 - H2) or dQ (Q1 - Q2). The field holds the name of the
chainage of the H2 or Q2 point.
Iteration / PID
PID-Section: Here the necessary data is entered if the calculation mode is
chosen as PID-operation.
Figure 2.85 The Iteration / PID property page when calculation mode is chosen
as PID-Operation.
where TPRequired is the required value of the target point, TPAct is the actual
value of the target point, LimitLow is the amount that the actual value of the
target point can be smaller than the required target point and LimitHigh is
the amount that the actual value of the target point can be larger that the
required value of the target point.
Figure 2.86 The Iteration / PID property page when calculation mode is chosen
as Iterative Solution.
142 MIKE 11
Tabular view: Structures
General
Most dambreak setups consist of a single or several channels, a reservoir,
the dam structure and perhaps auxiliary dam structures such as spillways,
bottom outlets etc. Further downstream the river may be crossed by
bridges, culverts etc. It is important to describe the river setup accurately
in order to obtain reasonable results. There is no limit to the number of
dam structures in a MIKE 11 model.
The description of the reservoir storage is carried out directly in the proc-
essed data. The only columns which contain 'real' data are those contain-
ing the water level and the additional flooded area.
In this way the surface storage area of the dam is described as a function
of the water level. The lowest water level should be somewhere below the
final breach elevation of the dam, and should be associated with some
finite flooded area. (This first value, hence, describes a type of 'slot' in the
reservoir).
The cross-sectional area is set to a large finite value. It is only used when
calculating the inflow headloss into the breach.
The dam
At the Q-point where the dambreak structure is located, the momentum
equation is replaced by an equation which describes the flow through the
structure. This may be either critical or sub-critical. A check on the energy
levels at the structure and at the next downstream h-point is first carried
out to determine which description is applicable. Refer to the MIKE 11
HD Reference Manual, Dambreak Section.
As the momentum equation is not used at the Q-point, the ∆X step used
between the adjoining h-points is of no consequence. The maximum ∆X
step should, however, be greater than the difference between given chain-
ages to prevent the insertion of interpolated cross- sections.
144 MIKE 11
Tabular view: Structures
At the node where the two branches meet, the surface flooded area is taken
as the sum of the individual flooded areas specified at each h-point.
Hence, if the reservoir storage has already been described in the reservoir
h-point, the spillway h-point should contain no additional surface areas. In
this case both the width and the additional flooded areas should be set to
zero. The cross-sectional area, hydraulic radii, etc. can be given as for the
reservoir branch.
Location
– River Name: Name of the river branch in which the dambreak is
located.
– Chainage: Chainage at which the dambreak is located.
– ID: String identification of the structure. It is used to identify the
structure if there are multiple structures at the same location. It is
recommended always to give the structure an ID.
– Type:
The lcation type may be Regular, Side Structure or Side Structure +
Reservoir. See 2.3.3 Structure Types definition (p. 77) for details.
Dam Geometry
– Crest Level: The crest level of the dam before failure.
– Crest Length: The crest length (perpendicular to the flow) of the
before failure.
Yes: The development of the breach is limited (e.g. solid rock below
the dam). The shape of the limitation should be specified in the
Cross Section Editor (p. 179).
146 MIKE 11
Tabular view: Structures
– Time after failure when changing the time step: Time for the
increase of time step relative to the failure time and specified in
hours.
Initial Conditions
In many cases dam failures occur on a dry river bed downstream. How-
ever, such initial conditions should be treated with caution in MIKE 11.
1 A small lateral inflow is added at the first h-point in the river down-
stream of the dam. This will ensure some depth of water in the river
from which a steady-state can be reached
2 The inflow into the reservoir can be non-zero, if desired.
3 The dambreak structure should be specified not to fail, i.e. to ensure
that the maximum calculated reservoir level is greater than the speci-
fied failure reservoir level (i.e. failure will not occur during the genera-
tion of the steady- state hot start file).
Initial conditions (h and Q) for this 'hot start' simulation must be specified
in the supplementary data, including the reservoir level.
With the hot start file ready, the dambreak simulation can now be carried
out. It is suggested that a DELTA value of slightly more than the default of
0.5 be used to damp out short waves which may lead to numerical instabil-
ities. A time step of the order 1-10 minutes is suggested.
148 MIKE 11
Tabular view: Structures
This dialog (Figure 2.89) is accessed from the Dambreak Str. (p. 143)
property page in the Tabular view: Structures (p. 76) by pressing the
. button (can only be accessed if the Failure Mode is set to
Erosion Based). The dialog can only be used to specify erosion based fail-
ure modes.
Purpose
The breach depth relationship is calculated using the Engelund-Hansen
sediment transport formula. Breach width is determined from the product
of breach depth and the side erosion index specified by user.
Dambreak Geometry
z Upstream slope: Slope (horizontal: vertical) of the upstream face of the
dam structure.
z Downstream slope: Slope (horizontal: vertical) of the downstream face
of the of the dam structure.
z Top Width: The top width of the dam crest.
Material Properties
z Grain Diameter: Representative grain diameter of the dam core mate-
rial.
z Specific Gravity 2.5 - 2.7: Relative density of the dam core material.
z Porosity 0.3 - 0.5: Porosity of the dam core material.
z Crit. Shear Stress 0.03 - 0.06: Critical shear stress of dam core material
used for sediment transport estimation (Shields criteria).
z Side Erosion Index: Multiplication factor used to calculate breach
width erosion rates from breach depth predictions.
z Final bottom level: The minimum level to which the breach is allowed
to develop.
z Final bottom width: The maximum width to which the breach is
allowed to develop.
z Breach slope: Slope (horizontal: vertical) on either side of the breach.
Initial Failure
The failure of the dam can initially take place in two ways:
Breach Failure
z Initial Level: The level of the breach develops in one time step as an
initial breach shape.
z Initial Width: The width of the breach develops in one time step as an
initial breach shape.
Piping Failure
z Starting Level: The level at which piping failure begins to occur.
150 MIKE 11
Tabular view: Structures
z Volume Loss Ratio 0 - 1: When the dam collapses some of the material
may be carried out without depositing on the bed of the breach. The
volume loss ratio is the fraction of the material to be washed out imme-
diately after collapse.
z Calibration Coef. > 0: Calibration multiplication factor used to adjust
the calculated change in pipe radius.
When activated, the user defined structure will access a DLL (Dynamic
Link Library) written by the user. The Network Editor interface contains a
number of variables that can be used in the DLL. In addition, the DLL can
access any variable in MIKE11 through several records.
The Tabulated Structure property page is used for defining a structure reg-
ulated by a user defined relation between the discharge through the struc-
ture and the up- and downstream water level. The relation is defined in a
table. The property page consists of a number of dialog boxes (see
Figure 2.91) whose functionality is described below:
Details
River Name: Name of the river branch in which the structure is
located.
Type: The lcation type may be Regular, Side Structure or Side Struc-
ture + Reservoir. See section 2.3 Tabular view: Structures for details
Calculation Mode:
152 MIKE 11
Tabular view: Structures
The upstream water level must increase in the right direction and
the downstream water level must increase in the downward direc-
tion. The discharge can not increase in the right direction and it can
not decrease in the downward direction.
H U/S= f(h D/S, Q): The upstream water level is given as a func-
tion of the discharge and downstream water level. The downstream
water level (h D/S) must be tabulated in first column and the dis-
charge must be tabulated in the first row in the table. Then the cor-
responding upstream water levels must be tabulated.
The discharge must increase in the right direction and the down-
stream water level must increase in the downward direction. The
upstream water level must increase in the right and the downward
direction.
Water level datum: The water level datum is added to the up- and
downstream water level in the table.
The Energy Loss property page is used to define energy losses associated
with local flow obstructions such as sudden flow contractions or expan-
sions and gradual or abrupt changes in the river alignment. Moreover, a
user defined energy loss coefficient can be defined.
Q 2
∆H = ζ -----------2- (2.5)
2gA
154 MIKE 11
Tabular view: Routing
Details:
River name: Name of the river in which the Energy Loss point is
located.
ID: String identification of the Energy Loss point. The specified ID has
no influence on the simulation.
Positive flow: Denotes the energy loss coefficient in the case of posi-
tive flow across the Energy Loss point in question. Applies to user
defined loss, contraction loss and expansion loss.
Negative flow: Denotes the energy loss coefficient in the case of nega-
tive flow across the Energy Loss point in question. Applies to user
defined loss, contraction loss and expansion loss.
To allow for the insertion of routing components into a branch the branch
type must be defined as either “Kinematic Routing” or “Routing”. See
section 2.2.2 Branches (p. 64).
156 MIKE 11
Tabular view: Routing
K1, P1, Q1, K2, P2, Tl, Tlz: Parameters for the calculation. See technical
reference for more details.
Alpha, P, ALR, SR and ANR: Parameters for kinematic wave. See tech-
nical reference for more details.
ID: Name of the routing component. Does not influence the simulation.
Type: The user should select the actual type of flood control.
158 MIKE 11
Tabular view: Routing
Initial water level: If checked the water level specified will be applied,
otherwise the initial water level will be equal to the water level giving an
outflow equal to the initial inflow.
ID: Name of the routing component. Does not influence the simulation.
2.4.5 Diversions
The dialog for specifying the parameters for a diversion is shown in
Figure 2.97. (Only visible and applicable is a special M11-Japanese mod-
ule is included in the actual MIKE 11 license).
160 MIKE 11
Tabular view: Routing
Normally when applying the routing facilities the network does not split
the flow as a proper calculation of the split requires a water level to be cal-
culated. However, using the diversion facility the user is allowed to spec-
ify how a branch splits into two branches. This is done by pre-defining the
split of flow, i.e. for a range of inflow discharges the amount continuing
along the main branch and along the tributary branch should be specified.
ID: Name of the routing component. Does not influence the simulation.
Whether the main river downstream carries the majority of the flow does
not matter.
This facility does not allow a routing branch to split be into more than two
branches. If this is required an artificial routing branch with no routing
elements has to be applied. Figure 2.98 shows how this is done when a
branch splits into three branches.
Artificial branch
162 MIKE 11
Tabular view: Routing
Details
Location
River name: Name of the river in which the Kinematic Routing point
is located.
Discharge Computation
Muskingum method: A routing method that requires the following
input parameters:
– K: Time scale describing the travel time of the water through the
Kinematic Routing element in question.
– x: A weighting factor greater than zero and smaller than 0.5.
164 MIKE 11
Tabular view: Runoff / Groundwater Links
Observe that the Include all branches feature will overwrite existing spec-
ifications.
Location
Leakage
Exchange Type
3 different River - Aquifer types of exchange between surface water and
aquifer are available (described in the technical documentation of the
MIKE SHE User Manual); Aquifer Only, Aquifer + Bed, River Bed only.
z Gaining Reach
In Gaining Reaches, the river receives baseflow from the baseflow res-
ervoirs in MIKE SHE
z Losing Reach
In Losing Reaches, the river will exchange water to the baseflow reser-
voirs in MIKE SHE as a function of the Leakage coefficient, Water
depth in river, Bank width and length of the coupling reach.
166 MIKE 11
Tabular view: Runoff / Groundwater Links
numerical problems, the flow is reduced linearly to zero when the water
depth is below the minimum upstream height threshold.
Inundation
z No Flooding
The no flooding option is equivalent to the old formulation in MIKE
SHE where rivers are considered a line between two adjacent model
grids. If this option is used one of the three river-aquifer exchange for-
mulations will be adopted. River-Overland exchange is always one-
way, namely overland to river. Over-bank spilling is not possible when
the No flooding formulation is adopted. The river water level may rise
above the topographic elevation of the adjacent grids without flooding
the grids.
Flood Code
Specification needed when the automatic or manual flood area option is
chosen.
168 MIKE 11
Tabular view: Runoff / Groundwater Links
Bed Topography
Specification needed when the automatic or manual flood area option is
chosen.
The MIKE SHE ground surface elevation can be re-defined in flood area
grid points, depending on the bed topography option. It should be empha-
sised that the flood mapping and dynamic flooding during the simulation
requires a good consistency between the MIKE 11 cross-sections and the
ground surface elevations of the corresponding MIKE SHE flood grid
points.
z Use Cross-section
When this option is specified the ground surface elevations of the
actual flood grid points are substituted with values directly interpolated
from the MIKE 11 cross-sections of the actual coupling reach. The set-
up program performs an inverse-distance-weighted interpolation, using
points (elevations) on the MIKE 11 cross-sections as discrete input
points. When the distance between individual MIKE 11 cross-sections
is higher than ½ Dx (grid size) extra discrete points are generated by
linear interpolation between the MIKE 11 cross-sections before the
grid interpolation is made. This is done to ensure that an approximate
river cross-sectional topography is incorporated in all MIKE SHE grids
along the river and not only where a MIKE 11 cross-section is located.
Please note that the interpolated grid values are only used inside the
area delineated by the MIKE 11 cross-sections used for interpolation.
When the manual flood area option is used, the user defined flood area
is not necessarily identical with the flood area covered by the MIKE 11
cross-sections. If the automatic flood area option is used the area cov-
ered by the MIKE 11 cross-sections and the flood area will always be
consistent, as the flood-area is generated (automatically) based on the
MIKE 11 cross-sections.
Bed Leakage
Specification needed when the automatic or manual flood area option is
chosen.
170 MIKE 11
Tabular view: Runoff / Groundwater Links
The substitution is made in all flood grid points of the actual coupling
reach.
Catchment Definitions
Name: Name of input catchment.
Connection to Branches
Branch Name: Name of the river branch for catchment inflow.
Overview
The dialog supplies a tabular overview of the catchment and connection
details, which have been entered manually one by one.
Purpose
The page has two specific purposes:
172 MIKE 11
Tabular View: Grid Points
Control of Output
When reduced output is selected only those grid points highlighted with a
check mark in the right hand side tree view, will be saved. The three levels
in the tree view are model setup, model branch and model grid points.
These are described below.
Table 2.4
Table 2.4
174 MIKE 11
Tabular View: Grid Points
Table 2.4
176 MIKE 11
CROSS SECTION EDITOR
177
178 MIKE 11
Raw data View
There are two types of cross section data; the raw survey data and the
derived processed data. The raw data describes the shape of the cross sec-
tion and typically comes from a section survey of the river. The processed
data is derived from the raw data and contains all information used by the
computer model (e.g. level, cross section area, flow width, hydrau-
lic/resistance radius). The processed data can be calculated by the cross
section editor or entered manually.
z River Name: The name given to the river branch. String of any length.
z Topo ID: Topographical identification name. String of any length.
z Chainage: River chainage of cross section.
Refer to one of the following sections for more Information:
The raw data editor is made up by three views plus a number of additional
dialog boxes:
z Tree view: Provides a list of all cross sections in the file. The list is dis-
played using a tree structure with three levels. The upper level contains
river names, the second contains the Topo-IDs, and the third contains
cross section chainage.
z Tabular view: Selecting a cross section with the left mouse button will
display the section information in the tabular view.
z Graphical view: An x-z-plot of the cross sectional data with markers
and vegetation zones indicated (the latter only for the quasi two dimen-
sional steady state solver with vegetation).
180 MIKE 11
Raw data View
Cross section ID
An individual identification tag for each cross section may be entered
here. This tag is subsequently displayed in MIKEView and does not influ-
ence the calculations.
Section Type
The type of section is set here. Four options are available:
– Open section:
The typical setting for river cross sections.
– Closed irregular:
Closed sections with arbitrary shape.
– Closed circular:
Closed circular section shape where only the diameter must to be
defined.
– Closed rectangular:
Closed rectangular section, where Width and height must be
defined.
Radius Type
The type of hydraulic radius formulation is set here. The choices are:
– Resistance Radius:
A resistance radius formulation is used.
– Effective Area, Hydraulic Radius:
A hydraulic radius formulation where the area is adjusted to the
effective area according to the relative resistance variation.
– Total area, Hydraulic Radius:
A hydraulic radius formulation where the total area is equal to the
physical cross sectional area.
Datum
A datum value may be entered here. The Datum is normally used for
adjusting the levels of the cross sections such that they conform with a
specific reference datum in the model area. The Datum value is added to
all vertical coordinates in the tabular view and the levels of the calculated
processed data hence, includes the Datum level as well.
Coordinates
Plan coordinates may be entered here for the left and right end points of
the cross section. If non-zero values are defined the values are used in the
graphical view of the network to display the cross section width and addi-
tionally - if the ‘Apply coordinates’ tick-box is enabled - determines also
the extent of cross sections in the HD Mapping routine (see section 6.11
Maps (p. 354) in HD Parameter Editor section).
Correction of X-coor
This feature is used for determining the correction angle for the X-coordi-
nates in the profile. The correction may be used for situations where the
cross section profile isn’t perpendicular to the centre line of the river.
The correction angle can be automatically calculated from the river center
line defined in the graphical network editor and the section plan coordi-
nates by activating the ‘Calculate angle’ button. Note that this button
requires that the simulation file is open, a network file and a cross section
file are defined in the simulation editor and that georeferenced coordinates
have been applied as described above.
Thalweg
θ
Figure 3.2 Definition sketch of the correction angle.
182 MIKE 11
Raw data View
Morphological Model
A level of divide can be entered. This level of divide has two purposes:
1 To identify a level where flood plains and the main channel are seperat-
ing. This feature will be used in a morphological sediment transport
simulation to distinguish the morphological active main channel part
from the flood plain part of the cross section where morphological
changes does not occur. In such a simulation, an internal division in a
main channel part and a flood plain part of the section will be made at
the x-coordinates corresponding to the Divide level z-value and even-
tual morphological changes will only be applied in the main channel
part of the section.
2 The second functionality of the Divide level is that it can be used for
activating the flood plain Resistance feature in a hydrodynamic simula-
tion, where a uniform resistance value can be applied in a section
above the level of divide as described and specified in section 6.15
Flood Plain Resistance (p. 371).
Resistance numbers
In this section of the raw data window the user chooses to how to deal
with bed resistance. Two choices have to be made:
Transversal distribution
This transversal distribution defines the description of resistance across
the cross section. There are three choices:
1 Uniform:
A single resistance number will be applied throughout the cross section
2 High/Low flow zones:
Three resistance numbers are to be specified: 1) Left high flow resist-
ance applying between marker 1 and 4. 2) Right high flow resistance
applying between marker 5 and 3, and 3) Low flow resistance applying
between marker 4 and 5. If marker 4 and 5 do not exist the low flow
resistance number will apply throughout
3 Distributed:
The resistance number is to be specified for each X,Z data set in the
raw data table.
Resistance type
Multiple options exists for defining the desired type of resistance method
in cross sectons.
z Relative resistance.
The resistance is given relative to the resistance number specified in
the .hd11 file (HD parameter file). Resistance numbers higher than one
always corresponds to higher physical resistance than specified in the
.hd11 file (i.e. independent of the choice of resistance number type in
the .hd11file).
z Manning’s n:
The resistance number is specified as Manning’s n in the unit s/m(1/3).
z Manning’s M:
The resistance number is specified as Manning’s M in the unit m(1/3)/s.
(Manning’s M = 1/Mannig’s n)
z Chezy number:
The resistance number is specified as Chezy number in the unit m(1/2)/s.
z Darcy-Weisbach (k):
The resistance is specified in the form of an equivalent grain diameter.
184 MIKE 11
Raw data View
Figure 3.3 The setting of the Darcy-Weisbach resistance. The ‘Edit’ button is
only active if the tick box next to it has been ticked.
The vegetation parameters for the Pasche resistance formula are set within
the ‘additional resistance due to vegetation’ dialog. The parameters which
are required are (fig.Figure 3.4)
Figure 3.4 The additional parameters to be entered for the Pasche resistance
formula.
The import and export functionality from ASCII files (see section 3.3)
also functions with the additional vegetation parameters.
To import the additional parameters the txt file should have the line with
the RESISTANCE NUMBERS augmented with the additional 8 parame-
ters specified in the order:
Selecting a river branch, Topo ID or cross section with the mouse and
pressing the right mouse button will open context sensitive pop-up
menus with different features for editing content of either individual cross
sections or the entire cross section file presents the pop-up menu.
Figure 3.5 presents the context sensitive pop-up menu’s appearing for
Branches/Topo-ID’s and single chainages respectively.
Figure 3.5 Tree-view context sensitive pop-up menus activated through the
right mouse button click. A menu similar to that on the left (but also
including the ‘Combine’ command) appears if opened at the Branch
or Topo-ID level, whereas the right menu appears when opened
from a chainage value (individual cross section).
The following editing facilities are available from the context sensitiive
pop-up menus:
Insert...
Once the ‘Insert...’ facility is selected from either the Branch, Topo-ID or
Section part of the tree view, the Insert branch dialog as presented in
Figure 3.20 is activated. Through this dialog it is possible to insert either a
new River Name, a new Topo-ID and/or a new sections defined by a
chainage of a cross section.
186 MIKE 11
Raw data View
Insert Interpolated
This feature is only available when activating the pop-up menu from the
cross section/chainage level of the tree view. The cross section editor
gives the user the possibility of inserting interpolated cross sections in a
given set-up. When selecting this feature a seperate dialog appears.
Figure 3.6 The Insert Interpolated Cross section facility from the Raw Data,
Tree view.
The user can either choose to interpolate a single cross section at a given
chainage or multiple cross sections. In the latter case a maximum distance
between the interpolated cross sections must be defined, and additionally,
the range of the interpolation need to be specified.
Delete...
From the pop-up menu it is possible to delete either an entire river branch,
a single Topo-ID within a river or a single cross section. Once activated,
the delete-feature will present a confirmation message-box where the user
must confirm the deleting of selected item by pressing Yes.
Figure 3.8 Delete sections option dialog from Raw data, tree view.
Rename...
It is possible to rename River names and Topo-ID’s or even
rename/change the chainage of a section through the ‘Rename’ facility.
When activating the Rename dialog from the pop-up menu in one of the
tree view levels it is possible to change either name strings for the river
name and/or Topo-ID or chainage-values for a river section, see
Figure 3.9.
188 MIKE 11
Raw data View
Figure 3.9 Rename feature dialog from Raw data, tree view.
Copy...
The Copy facility makes it possible to copy a single cross section, a Top-
ID or a an entire River branch to a different name/location. The copy dia-
log requests a Topo-ID, branch name and chainage before copying the
cross section.
If one or more sections have been selected it must be selected in the Copy
dialog whether the copy-function shall include all sections or only the
selected ones as illustrated in Figure 3.10.
Figure 3.10 Copy facility dialog from Raw data, tree view
Combine...
The combine dialog is used to combine two river branches of the same
name but with differing Topo-ID. The combination is saved as a new river
branch of the same name and a specified Topo-ID. The facility is designed
for combining cross sections at chainage locations where two sources of
cross section data exist. The dialog which appears once Combine is
selected is presented in Figure 3.11 .
Figure 3.11 Combine cross section profiles dialog from Raw Data, tree view.
The top non-editable edit-field presents which river branch has been
selected for combining sections. Following information shall be defined in
the dialog.
z Maximal difference:
The Maximal difference is the tolerance limit within which sections in
DEM and SUR Topo-ID are considered to correspond.
190 MIKE 11
Raw data View
z Synchronize to:
Specifies the method for combining sections. Currently only one
option is available: Centre (Mark 2) as illustrated in Figure 3.12.
Figure 3.12 Centre (mark 2) The plot shows a DEM section (indicated by the
blue/dark dots) and a SUR section (indicated by the grey dots). The
combined section will include the DEM section from x=0 to the left
marker (the diamond) of the SUR profile, then the SUR profile to the
right marker (third diamond from the left) of the SUR profile and
finally the DEM profile for the remaining part of the section.
Select / Unselect
The tree view in the raw data dialog provides a feature for selecting cross
sections. Most features such as deleting, renaming, copying, processing
and plotting can be applied on either all cross sections in the file or on
selected cross sections only. Cross sections which are selected are marked
in the tree view with chainage in bold. Individual cross section can be
selected in five different ways:
X
This column contains the transversal coordinates of the raw data.
Z
The vertical coordinates of the raw data.
X and Z pairs can be typed in manually, or can be imported from a text file
as described in section 3.3 below. A third alternative is to Copy and Paste
the two columns from a spreadsheet (either together or one by one); you
need to highlight the column header X or Z before pasting.
Resist.
This column is used for setting relative resistance or it may be used for
setting local values of Manning’s M or n, Chezy number or the equivalent
grain diameter. The selection depends on the setting in the resistance
combo box (see Figure 3.13).
Mark
The column is used for setting the markers 1 to 7 and eventual additional
user defined marks. Clicking an element in the ‘Mark’ column opens a
marker dialog as shown below from which a requested marker number can
be assigned for the selected raw data coordinate set (x,z).
192 MIKE 11
Raw data View
– User marker:
Any number above 7 may be used as a user marker. User markers
do not impact the simulation results. They are an option for indicat-
ing a specific point in a cross section e.g. the location of the meas-
urement gauge.
To remove a user marker set the numeric value to 0 and deselect the
apply flag.
Zone
This field and the following are only of concern in conjunction with the
quasi two dimensional steady flow with vegetation module.
Veg. h.
If a zone is set to either vegetation or bank vegetation this field becomes
active. The vegetation height is set here and the average vegetation height
for the corresponding panel is displayed in the graphical view.
Toolbar icons
Two toolbar-icons for graphical editing are available once the raw data
graphical view is activated, See Figure 3.15. Simply activate the graphical
view by clicking with the mouse in the view.
194 MIKE 11
Raw data View
Listed from left the icons presented in Figure 3.15 are as follows:
– Zoom In :
Enables zoom in of graphical view
– Zoom Out :
Zoom out to full extent of all sections visible in view
– Previous zoom :
Zoom to previous
– Next zoom :
Zoom to next (active if Previous zoom has been activated one or
more times)
– Draw Grid :
On/Off switch for the drawing of a grid in view
– Undo :
Undo facility for the last graphically edited point
– Select :
Selector of point in the graphical view. When clicking on a point in
the graphical view the line in the tabular view containing values for
the actual point will be highlighted
– Move points :
Moves a point in the graphical view. x- and z-values in the tabular
view is automatically updated from the moving of the point.
– Parallel move :
Moves a selected number of points graphically. To use this feature,
first select a group of points by clicking at the first point to move
and the last point to move. Thereafter the group of points within the
selection can be moved by click and drag in the graphical view.
– Insert points :
Insert additional points in a section. By point to a line segment of a
cross section and activating the insert points a new point will be
inserted in between two existing points on that section reach.
– Delete points :
Deletes points from the graphical view, and automatically points
will be removed from the tabular view as well.
Figure 3.16 Raw data, Graphical view: Right mouse pop-up menu
196 MIKE 11
Raw data View
Clear : Clears the view such that only the active section is presented
and additionaly zoomes out such that the actual cross section uses the
entire view.
Font : opens a seperate standard Font settings dialog where type, size
and colours of fonts can be modified.
3 The third group of features are different facilities for making graphical
editing of the cross section raw data.
These facilities include;
Select points,
Move points,
Move points parallel,
Insert points
Delete points.
This list of facilities is identical to the graphical editing toolbar icon
facilities as described above. Similar functionlity if a feature is acti-
vated through the pop-up menu as compared to activating a toolbar
icon.
4 The last group of features in the menu is a version manager facility as
described in more detail below.
Version Manager
The version manager is a facility which enables a storing of multiple ver-
sions of one or more cross sections. It is possible through the version man-
ager to save different versions of the same cross section which eventually
undergoes different changes in e.g. a design-optimisation project.
The version manager as presented in Figure 3.17 saves the different sec-
tions in a tree-view format which enables the user to keep track of differ-
ent changes the section has undergone. The Version manager data can be
saved to an Ascii file by use of the ‘Save’ buttons - and can be re-loaded
by the ‘Load’ buttons.
198 MIKE 11
Raw data View
Additional buttons
Additional buttons with different functionality is present in the lower part
of the editor under the raw data table, see Figure 3.19.
Update markers
This button updates markers 1, 2 and 3 in the actual section as the
extremes of the cross section (left bank limit, lowest point and right bank
limit respectively).
Note, this facility overwrites eventual user defined settings of these three
markers unless the appropriate boxes in the dialog under ‘Settings –>
Cross section.... –> Update Markers’ have been unticked. That is, markers
are only updated with this button if they are activated in this settings dia-
log.
Info
This simply gives an overview of data in the cross section data base:
z Number of Rivers
z Number of Topo IDs
z Number of cross sections in actual Topo ID
z Number of X, Z in actual profile.
200 MIKE 11
Raw data View
To change any of the crpss section parameters listed in this dialog the
respective ‘change’ tick box must be selected.
Chainages
Activate the ‘Calculate’ tick-box to enable the Chainages calculation tool.
Markers
An option for either deleting or changing/updating marker positions glo-
bally from one of the suggested options in the drop-down selection box.
The ‘apply’ tick-boxes determines which markers should be updated as
only selected marker boxes will be updated when processing the sections.
Action to be done
A number of options to be applied to all the cross sections in the set-up are
available:
The OK button
Press the OK button to execute the section settings selected in the dialog.
202 MIKE 11
Raw data View
Please note: Applying any of the above will overwrite any user edited set-
tings/data. There is no undo feature so make sure to save the cross section
data before activating the OK button.
204 MIKE 11
Raw data View
Draw GIS marks: Marks the locations on a cross section where data has
been extracted from GIS images.
Drawing style
The drawing style controls the Z axis display in the graphical view. There
are three options available:
Axis label
The axis label feature allows for user-defined axis labelling and selection
whether the graphical plot should include units or not.
Resistance scale
The resistance scale option defines whether a fixed or automatic scaling of
the Resistance axis is presented. In case of a fixed scale selected, the min-
imum and maximum value of the axis scale must be defined.
206 MIKE 11
Raw data View
Checks
A feature for checking raw data to ensure that sections are open or closed.
If data does not pass the chosen check-option then a message will occur.
Confirmations
The user can specify whether a confirmation dialog box should appear
when deleting points or clearing history in the graphical view.
Align
A snap to grid feature in the cross section editors graphical view.
Figure 3.28 Settings -> Cross Sections dialog; Update Markes page
This dialog is used for defining which of the markers 1, 2 and 3 should be
automatically updated. Only the selected markers in this dialog will be
Note on exception: when utilizing the quasi two dimensional steady state
with vegetation module the processed data does not reflect the values
used in the calculation. In the calculation kernel of this module the X, Z -
coordinates of the individual cross sections are used for determining the
hydrodynamic parameters of the individual panels.
The processed data view is similar to the raw data display. A tree view
exists on the left where the required cross sections can be selected. A tab-
ular view presents all processed data and a graphical view on the right
hand side displays the processed data graphically (see Figure 3.29).
208 MIKE 11
Processed data view
Level
Levels for which processed data are calculated in the cross section.
Default levels definition range from the lowest z-value and up to the high-
est z-value in the raw data table for the acutal section
The method of determining processed data levels can be selected using the
Levels Dialog activated by the Levels button ( section 3.2.4).
Radius
A resistance or hydraulic radius depending on the selected type in Raw
Data Editor (Resistance Radius or Hydraulic Radius).
Storage width
Top width of the cross section at each processed data level.
Resistance factor
This factor can be used to apply a level dependent, variable resistance in
the section cross section.The Resistance factor can contain two types of
values depending on the Resistance Type definition in the raw data;
210 MIKE 11
Processed data view
Conveyance
The Conveyance values are not used in the simulation but is primarily dis-
played as part of the processed data for the purposes of checking and veri-
fying that the conveyance relationship is monotonously increasing with
increasing water level, which is one of the key assumptions for the open
water hydraulics.
Toolbar Icons
A set of toolbar icons are acitve with the graphical view, see Figure 3.30
Listed from left the icons presented in Figure 3.30 are as follows:
– Zoom Out : Zoom out to full extent of all sections visible in view
– Next zoom : Zoom to next (active if Previous zoom has been acti-
vated one or more times)
212 MIKE 11
Processed data view
Figure 3.31 Right mouse pop-up menu for Processed data, Graphical view
Clear history: Clears the view such that only the active graph is pre-
sented and additionaly zoomes out such that the actual cross section
data graph uses the entire view.
Figure 3.32 Drop-down selection box for plot of processed data types
214 MIKE 11
Processed data view
The Levels dialog controls the method for calculating processed data in
sections as well as user-defined indications for number of processed data
levels as well as eventual definitions of minimum and maximum levels in
the processed data.
1 Automatic:
The levels are selected automatically.
In case Resistance Radius is applied in the Raw Data Editor, levels are
selected according to variations in section flow width.
In case Hydraulic Radius is applied, levels are selected according to
variation in the section conveyance.
2 Equidistant:
The levels are selected with equidistant level difference determined
from the number of levels specified. If required, the Minimum level
can be adjusted to control the bottom level of the processed data table.
However, it is generally recommended to include the full cross section
definition and hence, leave the Minimum and Maximum tick-boxes
un-checked.
3 User defined:
The levels can be fully or partially user-defined. The selected levels are
entered to the levels table on the dialog. If the number of defined levels
is less than required by the Number of Levels specification the remain-
ing levels will be selected automatically.
Minimum Level
The minimum calculation level to be presented as the first level in the
processed data table. Default is the lowest point in the section defined in
the Raw Data table.
Maximum Level
The maximum calculation level to be presented in the processed data.
However, it is strongly recommended not to apply this feature as it is
important to have the cross section properly processed up to the highest
point as defined in the Raw data.
Number of Levels
The desired number of calculation levels. The automatic level selection
method may not use the full number of level specified. This will occur
when a fewer number of levels is sufficient to describe the variation of
cross sectional parameters.
Table of Levels
This section of the dialog is only applicable if the level selection method is
user defined. The required levels are entered into the table manually. Lev-
els can be added by pressing the Tab key while positioned at the bottom of
the table. Levels can be deleted by selecting the row number and pressing
the Delete key.
216 MIKE 11
Importing cross sections using File Import
Figure 3.34 Cross section Main menu File dialog including options for Import,
Export and Print of cross sections.
From the text file MIKE 11 can load the data and change them to MIKE
11's internal data base format. The text file formats must correspond to
one of two types, depending on whether raw or processed data is to be
read.
Figure 3.35 File format of ASCII file used for importing data into MIKE11.
X x-coordinate
Z z-coordinate
rr relative resistance
218 MIKE 11
Importing cross sections using File Import
<#1> marker 1
<#2> marker 2
<#4> marker 3
<#8> marker 4
<#16> marker 5
<#32> marker 6
<#64> marker 7
Note that if a point has two or more markers the number after # is found as
a summation, for example: ‘<#6>’ indicates that the point represents
marker 2 and 3.
Each type of information in the ASCII file starts with an explanatory text
line followed by one or more lines containing numerical information. This
text line must start with three fixed characters, depending on the type of
data:
z Horizontal coordinates
Text line: COORDINATES
Numerical line: 1 27.43 13.293
where the first value is a coordinate data type indicator, which can be
defined as either ‘0’, ‘1’ or ‘2’.
where the first value is a flow direction type indicator, which can be
defined as either ‘0’ or‘1’.
z Datum adjustment
Text line: DATUM
z Closed section
Text line: CLOSED SECTION.
If this text line does not occur the section will be taken as open
z Radius formulation
Text line: RADIUS TYPE
Numerical line: 0
z x-z coordinates
Text line: PROFILE
z End of a cross-section
220 MIKE 11
Importing cross sections using File Import
If there is an error in the text file, the loading will be terminated and infor-
mation will be given regarding the erroneous line.
If data for a particular cross-section already exists in the data base, the
data in the text file will be ignored.
Selecting File –> Import –> Import Raw Data and Recompute it is possi-
ble to import raw data into MIKE 11’s cross section data base and recom-
pute the processed data automatically.
The first three lines are the same as for raw data: Topo-ID, River Name
and River Chainage. As for raw data format, it is hereinafter possible to
specify information about:
z processed data
The explanatory text line (see raw data) initiating the processed data
must start with PROCESSED DATA.
After this line, two text lines (headings), followed by M number of
lines with individual levels and the hydraulic parameters defined.
The processed data for each cross-section must finish up with a line
containing: *********.
Selecting File –> Import –> Import and overwrite Processed Data it is pos-
sible to import processed data into MIKE 11’s cross section data base and
overwrite the existing processed data. This facility is often used if for
example additional storage areas have been added to the processed data
and these data are copied into another data base.
The format of the ASCII text-file containing Levee marks coordinates is:
River Name, Topo-ID, Chainage, Left X, Right X, Left Y, Right Y (items
can be divided by 2 or more spaces or 1 or more tabs). One line for each
series of coordinates.
Example:
Donau<Tab>2005<Tab>0.00<Tab>LeftX1<Tab>RightX1<Tab>LeftY1<Tab>RightY1
Donau<Tab>2005<Tab>750.00<Tab>LeftX2<Tab>RightX2<Tab>LeftY2<Tab>RightY2
Donau<Tab>2005<Tab>2000.00<Tab>LeftX3<Tab>RightX3<Tab>LeftY3<Tab>RightY3
--
--
Export of cross sections can be done either for all sections through ‘Export
All Sections’, which will export the entire content of cross seciton file to
the text-file, or for Selected cross sections only (‘Export Selected Sec-
tions’).
222 MIKE 11
Plotting Multiple Cross Sections
1 Export All...:
Both raw and Processed data is exported to the same text file.
2 Export Raw...:
Only the raw data is exported to a text file.
3 Export Processed...:
Only processed data is exported to a text file.
Figure 3.37 Dialog with settings for multiple cross section plotting
Margins
The horizontal and vertical margins and the horizontal and vertical dis-
tances between the plots can be controlled in this section.
Fixed for all sections: User defined values for minimum and maximum of
the axis will be applied for all cross section plots.
224 MIKE 11
Plotting Multiple Cross Sections
Fixed scale: The scaling of the axis will be selected according to the user
defined ratio between the physical cross section size and the printed size.
I.e. all cross sections will be plotted on an axis with the same width (max-
imum minus minimum). The minimum of the scale can be controlled as
either a fixed offset below the data minimum for each cross section or as a
fixed value applied for all cross sections.
Design profile
Each plot will normally contain one cross section data set. However,
selecting this option it is possible to have another data set from the same
location (river name and chainage) drawn in each plot. The Topo ID for
the second cross section to be drawn is specified by the user. The second
cross section will be drawn using the graphical settings for “passive” cross
sections, and a legend for both cross section lines can optionally be drawn.
Output
The plots can either be routed to the printer or saved as meta graphics in a
number of metafiles. Each metafile can only contain one page. I.e. the
number of selected cross sections requires more than one page several
metafiles will be written. The file names are generated automatically by
adding _01, _02, _03 and so on to the file name specified by the user.
226 MIKE 11
BOUNDARY EDITOR
227
228 MIKE 11
Overview of the Boundary Editor
4 BOUNDARY EDITOR
The boundary editor is used to specify boundary conditions to a MIKE11
Model. It is used not only to specify common boundary conditions such as
water levels and inflow hydrographs but also for the specification of lat-
eral flows along river reaches, solute concentrations of the inflow
hydrographs, various meteorological data and certain boundary conditions
used in connection with structures applied in a MIKE 11 model.
The top split window contains the overall details of boundary conditions
defined in the model. Each boundary condition is defined as one row in
the Boundary Table and the table therefore, contains all boundaries
included in a model set up. There is no limit to the number of boundaries
(= table rows) that can be included in a model.
The view and contents of the second and the third split window depends
on the specifications of the selected boundary identified by the highlighted
row in the upper window (active row is highlighted) in the Boundary
Table. Additional information needed to specify the boundary conditions
are entered in the second and third split windows.
Figure 4.1 Layout of the boundary file when opened for the first time in an
application. Note that the third split window is empty.
– Boundary Description,
– Boundary Type,and
– Boundary Location (River Name and Chainage).
– In addition a Boundary ID can be entered, although this is optional.
Specifying an ID can be convenient for identifying the boundary,
but it has no effect on the calculation.
Select the actual Boundary Description and the Boundary Type by placing
the cursor at the right end of the edit field in question and left-clicking the
mouse. A drop-down list appears from which the appropriate boundary
type can be selected, see figure 4.2.
Figure 4.2 Drop down list for selection of the Boundary Description
To insert new boundaries (= rows in the Boundary Table in the first split
window), press the “Insert” button on the keyboard or use the <Tab> key.
A boundary (row) can be moved up or down in the table by selecting a
row (clicking in the left column) and dragging to the desired location. It is
also possible to sort the boundaries alphabetically by double clicking the
column headers. This operates for all the column headers in the top split
window.
z Open;
z Point Source;
z Distributed Source,
z Global;
z Structures;
230 MIKE 11
Overview of the Boundary Editor
z Closed.
The boundary description options are explained in detail below.
Boundary Type specifies the kind of data required for the boundary. For
each Boundary Description there are a number of valid Boundary Types.
Once a Boundary Description has been selected only the valid choices of
Boundary Type are displayed.
232 MIKE 11
Overview of the Boundary Editor
Figure 4.5 shows an example in which both global and local wind bounda-
ries are applied. The globally defined wind stress will be applied all over
the model except in the branch Main between chainage 0 - 10000 where a
different time series (wind speed and direction) has been applied
234 MIKE 11
Overview of the Boundary Editor
Figure 4.5 Example of the application of both global and local boundary condi-
tons.
equidistant time axis can both be applied). Further please ensure that the
extent of this file covers the full the simulation period.
When using the NWS methods the terminal breach level is defined by the
last value in the corresponding time series.
The content of the third split window will again depend on the specifica-
tions given in the second split window. The third split window deals pri-
marily with boundaries for the AD and ST modules.
236 MIKE 11
Overview of the Boundary Editor
z Include HD Calculation.
This box must be checked if the discharge time series is to be included
in the water balance in the HD calculation.
z Include AD calculation.
This box must be checked if the discharge is to be used with a concen-
tration to compute the mass inflow of a component in an AD simula-
tion. When checked, the associated concentrations are entered in the
third split window.
z Mike 12.
If this check box is checked, the boundary is applied to a two layered
branch.
– TS-Defined
This option is used when the specified concentration defined as the
boundary condition (constant or time-varying) shall be applied
directly as the concentration value at the boundary node.
– Open, Concentration.
This option is used at locations where outflow from the model area
takes place. When an outflow boundary becomes an inflow bound-
ary during a model simulation (eg due to tidal conditions) the
boundary condition is adjusted according to:
– t mix k mix
C = C bf + ( C out – C bf )e
– Open, Transport.
This type of boundary should be used where only inflow takes
place. The transport into the model area is computed using the spec-
ified boundary concentration and the discharge computed by the
HD model and hence, a well-defined amount (mass) of pollutant is
introduced into the model. Concequently, the computed concentra-
tion in the boundary node can differ from the concentration speci-
fied in the boundary file. It may be used where appreciable storage
and hence dilution of the inflow can take place close to the bound-
ary.
238 MIKE 11
Overview of the Boundary Editor
If the ‘Mike 12’ box is checked the layout of the boundary file changes as
shown in figure 4.8. It is now possible to define a discharge for both the
upper and the lower layer. Further there are now four possible AD bound-
ary types as each of the two layers can be a closed boundary. If the Bound-
ary Description were chosen as closed then both of the layers would be
regarded as closed boundaries. By specifying an Open boundary in the
Boundary Description it is still possible to set one of the layers as closed
e.g. the top layer can be of the Open, Transport type and the bottom layer
Closed. This combination is often used at upstream boundaries. In the
lower window, each component concentration needs to be defined for both
top and bottom layers.
240 MIKE 11
Overview of the Boundary Editor
The second split window is similar to that displayed for Open inflow
boundaries, with one additional facility:
z AD-RR.
If this check box is checked AD components can be included with the
inflow generated by the rainfall runoff models integrated in MIKE11.
In Figure 4.10 only the ‘Include HD calculation’ box is checked. It thus
represents a standard lateral inflow used in a HD simulation. Only the dis-
charge need be specified, as either a constant value or a time series.
If the ‘Include AD calculation’ box is also checked then the third split
window becomes editable and boundaries for the different AD compo-
nents can be entered (see Figure 4.11).The discharge specified in the sec-
ond split window is used both in the water balance and in the AD
calculation. In the AD calculation it is multiplied with the concentrations
in order to calculate the mass inflow for the different components.
Figure 4.11 Specification of a point source boundary for both HD and AD simu-
lations.
If the ‘Mike 12’ box is now also checked a new data section appears in the
second split window, see figure 4.12. (Note that the ‘AD-RR’ check box is
now hidden as this facility is not available in combinations with MIKE12
simulations.) In the second split window, a discharge time series must be
specified together with the level at which the inflow occurs.
242 MIKE 11
Overview of the Boundary Editor
Figure 4.13 shows the layout of the boundary file if we now uncheck the
‘Include HD calculation’ box. It is actually the same layout as if the
‘Include HD calculation’ box is left checked. This is because information
on the discharge is still required in order to calculate the mass flux of AD
components into the river branch. The only difference is that the discharge
is no longer used in the water balance.
Figure 4.13 Specification of a Point Source inflow boundary for a combined AD-
MIKE12 simulation. The discharge is not included in the water bal-
ance when the ‘Include HD calculation’ is left unchecked.
In figure 4.14 the ‘AD-RR’ box is checked. This facility can be used
where the concentration of rainfall or runoff from a NAM model are to be
used in an AD simulation. All other check boxes are now invisible.
Instead a section where information on the Catchment Name, Catchment
Area and Runoff Type must be specified are now present.
The runoff type must be selected between Total Runoff, Surface Runoff,
Root Zone Runoff, Groundwater Runoff or Rainfall. Note that Interflow is
only available for the NAM model. The AD-RR facility is not available
with MIKE12 branches. Also note that the Rainfall Runoff model must
run in parallel with the HD and AD models for the AD-RR facility to
operate.
244 MIKE 11
Overview of the Boundary Editor
Figure 4.17 Specification of a water level boundary for a combined HD-AD simu-
lation. Open concentration boundaries with a suitable kmix value are
used as this is an outflow boundary.
Figure 4.18 A water level boundary with ‘Include AD calculation’ and ‘Mike 12’
check boxes checked. The water level must be specified at the sec-
ond split window together with information on the AD boundary type.
Boundaries for the AD components are specified in the third split
window.
246 MIKE 11
Overview of the Boundary Editor
The Q-h relation is given in the table in the second split window. The Q-h
relation can either be entered from a known rating curve (eg. copied and
pasted from Excel) or automatically generated by selecting Tools -> ‘Auto
Calculation of Q/h Table...’ in the top menu bar. If this later option is
selected a new dialog appears, see figure 4.20.
Figure 4.20 Dialog used when making automatic calculation of Q-h relation.
Both methods will take H-values from the processed data of the specific
cross section (Levels column) and Q-values is calculated from either
Equation (4.1) or Equation (4.2).
A(h)
Q ( h ) = A ( h ) g ------------- (4.1)
W( h)
Q ( h ) = Conv ( h ) I (4.2)
where:
A(h) the level depending Area (from Cross section processed data)
W(h) the level depending Width (from Cross section processed data)
248 MIKE 11
Overview of the Boundary Editor
Conv ( h ) = C ( h ) ⋅ A ( h ) ⋅ R ( h )
where:
M(h), n(h) and C(h) are the respective Resistance number definitions
extracted from the resistance factor column in the in the cross section
processed data view.
As data from the specific cross section is required for calculating the Q/h-
relation, this facility is only available when the simulation editor is open
and the paths to the cross section and boundary files have been specified.
Please note that the location defined for the boundary and the chainage of
the cross section used to calculate the Q/h-relation must be within 1mm
distance of each other. If not, the Q/h-calculation tool can not find a rele-
vant cross section and no calculations takes place.
The second split window indicates that three time series must be specified:
Dam Breach Level; Dam Breach Width; and Dam Breach Slope.
Figure 4.22 Specification of a Dam Break boundary. Note the absence of the
third split window, which is not necessary for this, and the other
combinations of Boundary Description and Boundary Type shown.
Dam Boundary
Dam Boundaries are used in connection with stratified branches (MIKE
Reservoir model) when extraction from the dam needs to be specified. The
discharge value or time series is given in the second split window, while
the Level, Width and Height of the discharge point (extraction) should be
given in the third split window, see figure 4.23.
250 MIKE 11
Overview of the Boundary Editor
Figure 4.23 Specification of a Dam Boundary for a stratified branch (MIKE Res-
ervoir model). Discharge is specified in the second split window and
the geometrical data of the extraction point is specified in the third
split window.
Figure 4.24 The layout of the boundary file when the Boundary Type is chosen
to be rainfall. The second split window now contains a check box
used to specify if AD components should be included. AD compo-
nents must be specified in the third split window.
The second split window holds information on the time series needed and
also requires information about the type of sediment transport included in
the computations: single or multiple (graded) sediment fractions.
In figure 4.25, data for only one sediment fraction is needed and the rele-
vant time series is assigned in the second split window.
252 MIKE 11
Overview of the Boundary Editor
Figure 4.26 shows the layout of the boundary file for a graded sediment
model boundary. The third split window now prompts the user for fraction
numbers. The fraction numbers refer to the different fractions defined in
the Sediment Transport Editor.
dow the user specifies whether the boundary data are the absolute bottom
level or the change in bottom level. The time series for both types are
specified in the second split window. If the data is selected as change in
bottom level the absolute bottom level is calculated during the simulation
based its initial value.
254 MIKE 11
Boundary Editor Tools
top menu bar and is available when the lower split window is active. If the
Boundary Type is chosen as Sediment Transport a dialog will appear, see
figure 4.28. This dialog can be used to specify several boundaries simulta-
neously. The boundaries can be either constant values or time series,
depending on the selection made in the TS Type edit field.
Figure 4.28 Dialog for quick specification of graded sediment inflow boundaries
If the “Make List of Fractions” Tool is used for Boundary Types equal to
Bottom Level the dialog is slightly changed, see figure 4.29. Now the user
must specify if the boundaries should be ‘Fraction Value’ or ‘Change in
Fraction Value’.
Figure 4.29 Dialog for quick specification of graded sediment bottom level-
boundaries
select from the top menu bar ‘Make List of Components’. The dialog that
appears is shown in figure 4.30. This dialog works in a similar manner to
the dialog used to set up graded sediment boundaries. If constant boundary
values are requested the user must select the Data Type and enter the
number of components and the boundary value. If time varying bounda-
ries are requested the user must select the appropriate file. All legal (ie
concentration) time series items in the time series file will then be used as
boundaries. The first legal time series is used for component number one
the second legal time series for component number 2 etc.
If the boundary is open and used for a MIKE 12 simulation the tool oper-
ates slightly differently. If time varying boundaries are requested then the
two first legal time series are used for component number 1: The first is
used for the top layer and the second for the bottom layer. If constant
boundaries are requested the user will have to specify whether the values
are for the top or the bottom layer, see figure 4.31.
Figure 4.31 Dialog for quick specification of MIKE 12 boundaries. The user must
define if it is valid for the top layer or the bottom layer.
256 MIKE 11
Boundary Editor Tools
When this option is selected a dialog appears, see figure 4.32. This dialog
reflects a HD point source with a location (branch name and chainage), a
boundary ID and a discharge boundary. In this dialog each row represents
a new boundary. The user specifies the branch name and chainage for the
new boundary in the first two columns, and optionally the boundary ID.
The last two columns (with the common header Discharge) are used to
specify the discharge. If they are left empty the same discharge boundary
will be used for the new boundaries. If the new boundaries should use
other discharges the necessary information is entered here. If constant val-
ues are requested only the File /Value edit field should be filled out. The
dialog can also be filled by copying data from an Excel spreadsheet. When
the dialog is closed the user is asked if the new boundaries should be
pasted into the boundary editor.
If a time series is requested from a .dfs0 file with more than one item,
click on the button to specify which item you want.
If the point source being copied also includes AD boundaries the dialog
will also offer a possibility to change these boundaries, see figure 4.33, in
which a point source for three AD components are being copied. The col-
umns for the components work in the same way as the columns used to
copy the discharge.
The Scale Factor field is hidden by default but can be made visible by
right clicking the mouse in the File/Value edit field and then selecting
Scale factor from the pop up menu. The specified Scale Factor will be
multiplied with the boundary conditions value (constant or time series).
Figure 4.29 shows an AD point source with three components in which the
second component is reduced by 20 percent.
Figure 4.34 AD point source where the second component is reduce with 20
percent.
The user specifies the branch name and chainage interval to which the
change in scale factor(s) applies. Leaving the branch name blank corre-
sponds to selecting all branches, but the chainage interval must still be
specified. A new scale factor for the discharge boundary can be entered in
the appropriate field. If this is left blank no changes will take place.
258 MIKE 11
Users Upgrading from MIKE 11 Version 2002 or Previous Versions
260 MIKE 11
RAINFALL-RUNOFF EDITOR
261
262 MIKE 11
5 RAINFALL-RUNOFF EDITOR
The Rainfall Runoff Editor (RR-editor) provides the following facilities:
Simulation
The Rainfall Runoff Editor builds a file containing all the specified data
with extension .RR11. Once the catchments have been defined and the
rainfall-runoff, and the model parameters specified in the rainfall-runoff
editor, the Simulation is started from the MIKE 11 Run (or simulation)
Editor. It should be noticed that:
z Time step: It is recommended to use a time step not larger than the time
step in the rainfall series and not larger than the time constant for rout-
ing of overland flow. See example on Figure 5.2.
z Simulated catchment results can be linked with the River Network.
Catchment runoff/discharges and be inputted as lateral inflows and
summed to Normal and Routing river branch types, see sections 2.5.2
and 2.4 in the River Network Editor guide.
Results
MIKE 11 generates a variety of output types from a Rainfall Runoff simu-
lation ready to be used for model calibration and result presentation. These
are described in Section 5.11
Figure 5.1 Input page to the rainfall-runoff simulation in the Simulation Editor
264 MIKE 11
Specifying model Catchments
Figure 5.3 The Catchment page. Additional catchments are prepared via the
Insert Catchment dialog. The Example includes 2 sub-catchments
and a combined catchment which includes the 2 sub-catchments
Inserting Catchments
New catchments are defined via the Insert Catchment dialog (see
Figure 5.4). The insert catchment dialog is automatically activated for the
first catchment, when creating a new RR-parameter. A new RR parameter
File is created from the MIKEZero File dialog. Additional catchments are
defined when pressing the button: Insert catchment.
Catchments Definitions
A catchment is defined by:
Catchment Name
Simulations can be carried out for several catchments at the same time.
The catchment name could reflect e.g. the location of the outflow point.
266 MIKE 11
The NAM Rainfall-runoff model
Catchment Area
Defined as the upstream area at the outflow point from a catchment.
Calibration plot
A calibration plot will automatically be prepared for catchments, where
the time series for observed discharge have been specified on the Time
series Page and the selection of calibration plot has been ticked off. The
calibration can be loaded from the Plot composed and is saved in the sub-
directory RRCalibration with the file name: Catchment-name.plc. The
time series in these plots are also available in DFS0 format in the subdi-
rectory RRcalibration with the file name: Catchmentname.dfs0.
Figure 5.32 shows an example on a calibration plot.
Calculated Areas
The Calculated area shown in the Catchment Overview is based on the
digitised catchment boundaries in the Graphical display. The calculated
area is activated when the Basin View has been selected, see section 5.10.
The Catchment Area is shown in the edit fields for Area and Calculated
Area, when transferring a catchment from the Basin View to the catch-
ment page. The Area which is used in the model calculation can after-
wards be modified manually.
5.2.1 Surface-rootzone
Parameters used in the surface and the root zone are described below (see
Figure 5.5).
268 MIKE 11
The NAM Rainfall-runoff model
The Ground Water parameters are described below (see Figure 5.6).
Overall Parameters
270 MIKE 11
The NAM Rainfall-runoff model
Abstraction
Ground water abstraction or pumping may be specified in a time series
input file, in millimetres, or given as monthly values in mm by clicking
Edit Abstraction.
The simple degree-day approach uses only the two overall parameters: a
constant degree-day coefficient and a base temperature.
The Snow melt module uses a temperature input time series, usually mean
daily temperature, which is specified on the Timeseries page.
The Snow Melt parameters are described below (see Figure 5.7).
272 MIKE 11
The NAM Rainfall-runoff model
Overall Parameters
Radiation coefficient
May be introduced when time series data for incoming radiation is availa-
ble. The timeseries input file is specified separately on the time series
page. The total snow melt is calculated as a contribution from the tradi-
tional snow melt approach based on Csnow (representing the convective
term) plus a term based on the radiation.
Elevation Zones
Elevations zones are prepared in the elevation zone dialog (see Figure 5.8)
274 MIKE 11
The NAM Rainfall-runoff model
Correction of precipitation
Specifies the lapse rate for adjustment of precipitation. Precipitation in the
actual elevation zone is calculated based on a linear transformation of the
precipitation at the reference station to the actual zone defined as precipi-
tation lapse rate (C/100m) multiplied by the difference in elevation
between the reference station and the actual zone.
Area of each zone is specified in the table. The total area of the elevation
zones must equal the area of the catchment.
5.2.4 Irrigation
Minor irrigation schemes within a catchment will normally have negligi-
ble influence on the catchment hydrology, unless transfer of water over the
catchment boundary is involved. Large schemes, however, may signifi-
cantly affect the runoff and ground water recharge through local increases
in evaporation and infiltration. If the effect of an irrigation area within a
catchment is to be simulated, separate NAM catchments are defined for
the irrigated area and the remaining area and a combined catchment
defined to accumulate the runoff.
276 MIKE 11
The NAM Rainfall-runoff model
Include irrigation
Ticked for a sub-catchments with irrigation included.
Infiltration Parameters
Irrigation sources
Can be local ground water, a local river, an external river, or a combina-
tion of these. Local ground water will be taken from the NAM ground
water storage and irrigation water taken from a local river will be sub-
tracted from the simulated runoff. If all the water is abstracted from an
external source, outside the catchment, no subtractions are made.
Ground water
Initial values for baseflow must always be specified. When lower base-
flow are included a value for the initial lower baseflow must also be spec-
ified.
Snow melt
Initials values of the snow storage are specified when the snow melt rou-
tine is used. When the catchment are delineated into elevation zones, the
snow storage and the water content in each elevation zones are specified.
278 MIKE 11
The NAM Rainfall-runoff model
5.2.6 Autocalibration
Automatic calibration is possible for the most important parameters in the
NAM model. A detailed description of the automatic calibration is given
in the Rainfall-runoff reference manual.
Include Autocalibration
Ticked for a sub-catchment with autocalibration included.
Calibration parameters
The automatic calibration routine includes the 9 model parameters:
280 MIKE 11
The NAM Rainfall-runoff model
Objective Function
In automatic calibration, the calibration objectives have to be formulated
as numerical goodness-of fit measures that are optimised automatically.
For the four calibration objectives defined above the following numerical
performance measures are used:
Stopping Criteria
The automatic calibration will stop either when the optimisation algorithm
ceases to give an improvement in the calibration objective or when the
maximum number of model evaluation is reached.
5.3 UHM
Introduction
The UHM (Unit Hydrograph) module constitutes an alternative to the
NAM model for flood simulation in areas, where no streamflow records
are available or where unit hydrograph techniques are already well estab-
lished.
282 MIKE 11
UHM
Hydrograph
The distribution of the runoff in time can be described using different
methods:
Six other methods for describing the hydrograph are available. These are:
– Storage Function
– Quasi Linear Storage Function
– Nakayasu
– Rational method
– Kinematic Wave (rectangular basin)
– Kinematic Wave (Non-uniform slope length)
For each of these a number of parameters are to be given. These parame-
ters are described in more details in the reference manual.
Loss model
Constant loss
The infiltration is described as an initial loss at the beginning of the storm
followed by a constant infiltration:
Proportional loss
A runoff coefficient is specified as the ratio of runoff to the rainfall.
– Nakayasu
– f1-Rsa
– No loss
Lag time
Can be specified directly in hours or calculated by the standard SCS for-
mula:
SCS formular
Three parameters are specified: Hydraulic Length, Slope and Curve
Number. Use the ‘calculate’ button to calculate the actual lag time.
5.4 SMAP
Introduction
SMAP is simple rainfall runoff model of the lumped conceptual type.
It has been designed to work on the basis of monthly input data and there-
fore constitutes an economic alternative to the NAM model in scenarios
where a daily resolution of the results is not required. This is often the
case in overall water resources planning or for analyses of longterm reser-
voir operations. In such situations data preparation time may be saved if
simulations are carried out with monthly time steps only.
The SMAP model has been tested by DHI on various dry tropical and sub-
tropical catchments and has shown almost the same degree of accuracy on
the simulated monthly flow as the NAM model. The model does not
include a snow melt routine and is not recommended to be used in areas
where snow melt has significant influence on the hydrographs.
Model Parameters
The model accounts for the water storage in two linear reservoirs repre-
senting the root zone and the groundwater reservoirs respectively.
284 MIKE 11
SMAP
Similar to the NAM model many of the process descriptions in the SMAP
model depends on the current saturation fraction of the root zone storage.
I.e. the current storage of water (RSOL) divided by the max. possible stor-
age (MAX).
Evaporation Exponent
The actual evaporation (EA) is calculated as a fraction of the potential
Evapotranspiration (EP). It depends on the current saturation degree of the
root zone and the exponent E1. Small E1 will increase the Evaporation.
The parameter influences the total amount of base flow generated by the
model.
Autocalibration Option
Not yet implemented!
In addition to the above parameters the root zone content (in mm) at the
start of the simulation and the initial base flow (in m3/s) needs to be spec-
ified.
286 MIKE 11
Urban
Please note, however, that the discharge output in the main result file is in
m3/s and represent an instantaneous value at by the end of the last calcula-
tion time step. I.e. it is not the average discharge during the storing inter-
val.
Values of specific discharge (in mm) accumulated over the storing interval
are available in the file for additional results. This file also includes time
series of other relevant parameters such as groundwater recharge, base
flow and root zone moisture.
5.5 Urban
5.5.1 Introduction
Two different urban runoff computation concepts are available in the
Rainfall Runoff Module as two different runoff models:
The Model type (A/B) is selected in the first group box Model Parameters
-> Model (see Figure 5.16 and Figure 5.17)
The Parameters for Model A are described below (see Figure 5.16)
Impervious Area
The Impervious area represents the reduced catchment area, which con-
tributes to the surface runoff
Time of Concentration
Defines the time, required for the flow of water from the most distant part
of the catchment to the point of outflow
Initial Loss
Defines the precipitation depth, required to start the surface runoff. This is
a one-off loss, comprising the wetting and filling of catchment depres-
sions.
Reduction factor
Runoff reduction factor, accounts for water losses caused by e.g. evapo-
transpiration, imperfect imperviousness, etc. on the contributing area.
288 MIKE 11
Urban
Time/Area Curve
Accounts for the shape of the catchment lay-out, determines the choice of
the available T/A curve to be used in the computations.
1) rectangular catchment
The Parameters for Model B are described below (see Figure 5.17)
Length
Conceptually, definition of the catchment shape, as the flow channel. The
model assumes a prismatic flow channel with rectangular cross section.
The channel bottom width is computed from catchment area and length.
Slope
Average slope of the catchment surface used in the Manning equation for
the runoff.
Area (percentages)
The area distribution percentages allow the catchment area to be divided
in up to five sub-catchments with identical geometrical, but distinct hydro-
logical properties. The five sub catchment types are:
– impervious steep
– impervious flat
290 MIKE 11
Urban
Wetting loss
One-off loss, accounts for wetting of the catchment surface.
Storage loss
One-off loss, defines the precipitation depth required for filling the
depressions on the catchment surface prior to occurrence of runoff. Not
applicable to ‘Impervious Steep’ type of catchment.
Start infiltration
Defines the maximum rate of infiltration (Horton) for the specific surface
type.
End infiltration
Defines the minimum rate of infiltration (Horton) for the specific surface
type.
Horton's Exponent
Time factor "characteristic soil parameter" with unit of s-1. Determines
how quickly during a rainfall event the infiltration capacity rate reduces
over time from the ‘Start Infiltration’ to the ‘End Infiltration’. The actual
infiltration capacity is made dependent of time since the rainfall start only.
Manning's number
Describes roughness of the conceptual catchment surface, used in hydrau-
lic routing of the runoff (Manning's formula).
Note that the Mannings number here is specified as Manning M - the
reciprocal of Manning’s n.
Additional runoff
Additional runoff Evaporation check box - controls if the evapo-transpira-
tion process will be included in the runoff computations can be specified
as a constant flow or specified as load based on inhabitants (PE). An
additional time series for load (qload) is specified on the time series, when
the flow is based on load based on inhabitants (PE>0). The flow is calcu-
lated as:
Evaporation
Evaporation check box - controls if the evapo-transpiration shall be calcu-
lated based on a time series (when checked the time series is specified on
the Time series page) or based on a constant loss (equal to 0.05 mm/hour).
Snow melt
Snow melt check box - controls if snow melt is included in the calculation.
The content of the snow storage melts at a rate defined by the degree-day
coefficient CSnow multiplied with the temperature deficit above 0 Degree
Celsius. Typical values for Csnow is 2-4 mm/day. When snow melt is
checked a time series for temperature is specified on the Time series Page.
292 MIKE 11
Flood Estimation Handbook (FEH)
1. Overview
5. Catchment Descriptors
294 MIKE 11
Flood Estimation Handbook (FEH)
296 MIKE 11
Flood Estimation Handbook (FEH)
· Rainfall generation
· CWI estimation
· Contribution of Snowmelt
298 MIKE 11
Flood Estimation Handbook (FEH)
Very simply, the time to peak of the instantaneous unit hydrograph, Tp(0),
for a PMF computation is assume to be 0.67 times the standard value. This
affects both the peak (Up) and time base (TB). See Vol. 4, Chap 4.2.1.
Rainfall Generation
Contribution of Snowmelt
Snowmelt may contribute to both the storm depth and antecedent rainfall,
and therefore the CWI. The user should define a snow melt rate (mm/h)
from which both these effects can be computed. See Vol. 4, Chap 4.3.4
and example 4.1f. Output from this part is an adjusted CWI (denoted
CWI') as well as a modified storm profile in dfs0 format (as for step 10,
above).
If using a winter PMP, the SPR is set to a minimum of 53% to account for
frozen ground. See Vol. 4 Chap 4.2.2. In addition, a revised formulation
for PR is made, Vol. 4, Eq. 4.12.
Initially, the SPR is computed as for step 12, above, and subsequently
checked to ensure that the value is greater than or equal to 0.53 (if a winter
profile is used).
Catchment rainfall is provided by the user as a dfs0 file. Please note that
mean area rainfall computation is done on the time series page in the Rain-
fall-Riunoff editor. Following the specification of a catchment rainfall file,
the period (start and end) covered by the time series will be shown. Using
this information the user is required to set the design storm period, which
defines the storm duration and the rainfall depth.
Storm Profile
An option should be available to allow the user to use the measured rain-
fall time series distribution to generate the hydrograph, or else one of the
two standard profiles (summer or winter).
The CWI calculation is based on the observed rainfall record in the 5 days
prior to the start of the event and the observed soil moisture deficit (SMD).
The procedure is described in detail in Vol. 4, Appendix A, section A.4.2.
SMD is defined by the user. Subsequently, CWI is computed by MIKE
FEH.
5.6.6 Results
The user may set the origin of the time axis of all the result files computed
in MIKE FEH. This may be desirable if a hydraulic analysis (using MIKE
11 HD) is to be done afterwards.
300 MIKE 11
DRiFt
If, for a T-Year event, multiple return periods have been specified, the
result files contain multiple columns, one column for each event.
5.6.7 Validation
Once the user presses the Compute button on the Results page, MIKE
FEH starts the validation of the provided input. If the input is accepted, the
model proceeds, otherwise the validation error messages are shown in the
interface.
In the interface, the user may indicate for each catchment whether or not a
log should be created.
5.7 DRiFt
Introduction
The DRiFt module (DRiFt = Discharge River Forecast) is a semi-distrib-
uted rainfall-runoff model based on a morphological approach. The model
is able to consider the topography of each site analyzed and the spatial
variability of soil characteristics and rainfall patterns. Input data for the
DRiFt model is divided into three groups:
Geo-morphological Parameters
DEM
The DEM (Digital Elevation Model) of the basin - prepared and saved in a
two-dimensional grid-file (*.dfs2). For a better interpretation of the
derived files (like draining network files, etc.), it is opportune to mask the
useless parts of the territory (e.g. the sea) setting them to "no data value".
302 MIKE 11
DRiFt
Geo-morphological exponent, k
The Geo-morphological exponent, k is the exponent used in the expres-
sion for the threshold value ASk. Default value is 1.7.
Other files are automatically generated together with of the draining net-
work file (e.g. catchment_name.pnt.dfs2 representing the slope orienta-
tion, catchment_name.area.dfs2 representing the drained upstream area).
These files are used in the simulation even though they do not appear
within the interface, therefore they always must be located together in the
same directory. The user must pay attention to the operations which could
alter this conditions (renaming, moving, deleting these files).
Surface Parameters
Activating the ‘Average value’ tick-mark will activate the average value
field and a constant value must be defined. De-activating the tick-mark
will require a dfs2-file selection in the ‘Distributed’ field by use of the
browse button.
Default value for flow velocity on channels is 1 m/s. The normal range of
the variable is from 0,1 m/s to 10 m/s.
Default value for flow velocity on hillslopes is 0,1 m/s. The normal range
of the variable is from 0,001 m/s to 1 m/s.
304 MIKE 11
DRiFt
5.7.3 Rainfall
The input for the rainfall-runoff simulation is supplied as rainfall (or pre-
cipitation) data. Rainfall data can be specified either as a constant value,
as a timeseries or as a time-varying distributed rainfall pattern.
Precipitation Rate
DRiFt utilises a spatial distributed map of precipitation as input to the
computational part. It is therefore required to preprocess a rainfall map
from either spatial and/or temporal distributed sources.
Spatial distribution
Spatial distribution of precipitation can be made either Uniform or Dis-
tributed. Select the required option from the Spatial Distribution combo-
box.
Temporal distribution
Temporal distribution of precipitation can be made either Constant (= con-
stant value in space and time) or Time Varying.
306 MIKE 11
DRiFt
TS-File
With a selection of uniform spatial distribution and time varying temporal
distribution as presented in Figure 5.21, it is required to select a time
series file (dfs0-file) with rainfall data in the ‘TS-file’ filename field.
Rainfall file
With a selection of distributed precipitation maps as presented in
Figure 5.22, it is required to select a time varying grid-file (dfs2-file) with
rainfall data in the ‘Rainfall file’ filename field.
308 MIKE 11
DRiFt
Interpolation type
Interpolation options for generation of time varying distributed precipita-
tion maps are: ‘Thiessen’ and ‘Inverse squared distance’. Select the
desired interpolation type from the Interpolation type combo-box.
Precipitation timestep
The precipitation time step is the temporal resolution of the new distrib-
uted maps [seconds]. In order to exploit the whole available information,
it is recommended to set the precipitation time step on the same value of
raingauges measurements resolution.
310 MIKE 11
Time Series
Rainfall
A time series, representing the average catchment rainfall. The time inter-
val between values, may vary through the input series. The rainfall speci-
fied at a given time should be the rainfall depth accumulated since the
previous value.
Evaporation
The potential evaporation is typically given as monthly values. Like rain-
fall, the time for each potential evaporation value should be the accumu-
lated depth at the end of the period it represents. The monthly potential
evaporation in June should be dated 30 June or 1 July.
Temperature
A time series of temperature, usually mean daily values, is required only if
snow melt calculations are included in the simulations.
Irrigation
An input time series is required to provide information on the amount of
irrigation water applied, if the irrigation module is included in a NAM
simulation
Abstraction
Groundwater abstraction can be included in NAM simulations for areas,
where this is expected to influence e.g. the baseflow. The data should be
given in mm.
Radiation
A time series of incoming solar radiation can be used as input to the
extended snow melt routine.
Degree-day coefficient
A time series of seasonal variation of the degree-day coefficient can be
specified as input to the extended snow melt routine.
Observed Discharge
A time series of observed discharge values can be specified and used for
model calibration. The observed discharge must be specified when auto-
matic calibration is included.
312 MIKE 11
Time Series
1 Directly within the Rainfall Runoff Editor. From the top toolbar menu
select Basin Work Area and the Calculate mean precipitation. The cal-
culation is made without requiring a model run.
2 During a simulation. A new simulation is started in the Simulation Edi-
tor: If the weighted time series is ticked, the Mean Area weighting cal-
culation is carried out as part of the model run.
It is recommended to use option 1.This will ensure that the available peri-
ods of the input files known in the simulations editor.
After having calculated the weighted time series once the calculation can
be disconnected when removing the tick mark for weighted time series.
Alternatively, the user may specify the weight to be used for specific com-
bination of missing data. For each such catchment, a suitable weight
should be specified for the reporting stations and a weight of “-1.0” given
for the non-reporting station(s), including missing data.
Distribution in time
If data is available from stations reporting at different frequencies, e.g.
both daily and hourly stations, the Distribution in time of the average
catchment rainfall may be determined using a weighted average of the
high-frequency stations. You may, for example, use all daily and hourly
stations to determine the daily mean rainfall over the catchment and sub-
sequently use the hourly stations to the distribute (desegregate) this daily
rainfall in time. Different weight combinations for different cases of miss-
ing values may be applied also to this calculation of the distribution in
time.
Deleting stations
Stations which are not longer valid in the weight combinations are
removed from the editor by deleting the station number in the editor.
Delete values
The delete value used in the time series indicating periods with missing
data is usually specified with the default delete value ‘1e-30’. The default
delete value can be changed via the MIKEZero Data Utility tool.
314 MIKE 11
Basin View
When opening a new Basin View at least one catchment (usually the
default) must exist in the Rainfall Runoff Tabular View, which must be
open same time as the Basin View. This initializes the Rainfall Runoff
Editor. The default catchment can afterwards be deleted from the catch-
ment page in the Tabular View, such as the catchments in the Basin View
and on the Tabular View are the same.
316 MIKE 11
Basin View
1 n1
X1 Y1
X2 Y2
2 n2
X1 Y1
X2 Y2
LIne 5 to 8 are repeated for the following sections. The # Marks the of
all sections
1 478.2 98.0
2 488.5 110.1
3 462.5 113.2
4 425.0 151.9
Thiessen Options
Preparation of Thiessen Weights takes place from the “Thiessen Option”-
dialog. Select number 1 for the first combination and press OK (see
Figure 5.26). Thiessen weights have now been prepared on the Time
series page (see Figure 5.23, Time series page in the Rainfall Runoff Edi-
tor).
Apply the weight “-1.00” (for stations with missing data) on the timeseries
page before calculating of other combinations.
Isohyetal Options
The Isohyetal Option is used as a post processing tool to calculate average
catchment rainfall for a fixed period based on isohyetal lines. The tool has
no link to data on the Timeseries page in the Rainfall Runoff Editor. It
should therefore be noticed that the Isohyetal Option can not be used to
prepare weights and time series of mean area rainfall used as input to the
rainfall-runoff calculation. Select the Isohyetal Options to activate the Iso-
hyetal Option dialog (see Figure 5.27). The dialog has the following
pages:
1 Preparation of periods.
2 Grid Interpolation
3 Isoline Options
318 MIKE 11
Basin View
Combination Definitions
Options used to View different Thiessen Polygons on the Basin View.
Graphical Settings
Graphical Settings can be modified from the Graphical Settings Dialog
(see Figure 5.28).
The “Graphics” page is used to adjust display options for the following
graphical objects:
– Thiessen Objects
The page “Mouse” is used to adjust the digitizing distance and the Mouse
sensitivity for digitizing on the screen.
Create Polygons
After having digitized the catchment boundaries this option is used to cre-
ate catchment polygons (alternatively press the Create Polygon Catch-
ments icon: ). Each catchment will be created in the Rainfall Runoff
Editor, including an automatic calculation of the area.
320 MIKE 11
Basin View
Testing catchment
After having prepared the catchment boundaries the “Test fill catchment”-
icon ( ) can be used to test the validity of the digitized catchment poly-
gons.
Create Polygons
Catchments are created within the Rainfall Runoff Editor using the “Cre-
ate Polygon Catchments”-icon ( ) after having digitized the catchment
boundaries. Each catchment will be created in the Rainfall Runoff Editor,
including automatic calculation of the area. Catchment names can be mod-
ified in the Rainfall Runoff Editor.
Defining Stations
New rainfall stations are created with the “Create New Stations” icon
( ). Click in the Basin View on the Station Location and use the “Edit
Station”-dialog to select the time series and select the name for the Rain-
fall station (see Figure 5.29).
Deleting Stations
Rainfall stations are deleted from the Basin View as follows:
Editing Stations
Stations are modified in the “Edit Station”-dialog as follows:
322 MIKE 11
Result Presentation
Figure 5.30 shows an example of a Basin View for two catchment show-
ing the catchment boundaries, 7 rainfall stations and the Thiessen poly-
gons for all 7 stations.
Figure 5.30 Basin View with catchment boundaries, rainfall stations and Thies-
sen polygons.
Results
MIKE11 generates two Rainfall Runoff Result files. The first result file
contains simulated runoff and net precipitation. The second, additional
result file (RRAdd) contains time series of all calculated variables, such as
the moisture contents in all storages, the baseflow etc., and can be very
useful during model calibration. The results of the simulation can be gen-
erated in two formats, either as RES11 or DFS0 filetype. The format of the
result file should be selected before running the simulation. Three facili-
ties are available to plot and analyze the results of a rainfall-runoff simu-
lation:
2. MikeZero Time series Editor The time series editor can also be used
to view and compare simulated and measured results and to export results
to e.g. a spreadsheet for further processing. The result file should then be
given a DFS0 extension.
Summarised output
MIKE 11 generates as standard a table with yearly summarised values of
simulated discharge. The table is stored as the textfile “RRStat.txt” in the
current simulation directory. The table is extended with observed dis-
charge for catchments, where the time series for observed discharge have
been specified on the Timeseries Page. This includes a comparison
between observed and simulated discharge with calculation of the water
balance error and the coefficient of determination.
324 MIKE 11
Result Presentation
Calibration Plot
A calibration plot will automatically be prepared for catchments, where
the time series for observed discharge have been specified on the Time
series Page and the selection of calibration plot has been ticked off on the
catchment page. The calibration can be loaded from the Plot composed
and is saved in the subdirectory RRCalibration with the file name: Catch-
ment-name.plc. The time series in these plots are also available in DFS0
format in the subdirectory RRcalibration with the file name: Catchment-
name.dfs0. The plot shows following results (see Figure 5.32):
326 MIKE 11
A Step-by-step procedure for using the RR-Editor
328 MIKE 11
HYDRODYNAMIC EDITOR
329
330 MIKE 11
6 HYDRODYNAMIC PARAMETERS EDITOR
The Hydrodynamic parameters editor (HD-editor) is used for setting sup-
plementary data used for the simulation. Most of the parameters in this
editor have default values and in most cases these values are sufficient for
obtaining satisfactory simulation results. The HD parameter editor con-
tains a number of tabs as present in Figure 6.1 below.
A list of the tab-page available in the Editor are presented below and
described in the following: paragraphs.
6.1 Initial
Initial conditions for the hydrodynamic model are specified on this page.
The initial values may be specified as discharge and as either water level
or water depth. The radio button determines whether the specifications are
interpreted as water level or depth.The global values are applied over the
entire network at the start of the computation. Specific local values can be
332 MIKE 11
Initial
New lines in the Local Values table are created either by placing the cursor
in the last column and pressing <TAB> or through the ‘Grid’ option in the
main menu bar;.
This definition wil result in a local varying initial water depth from 5.70 to
5.00 with a linear relationship between chainage 0 and 3000 and the dis-
charge also varies between 1.000 and 1.400 with a linear relationship
between chainage 0 and 3000. At any other locations in the model, Global
values will be applied as the initial value.
every single calculation point in the model. Local Values should only bee
defined at locations where change in values should be applied.
6.2 Wind
Wind friction on the water surface can be accounted for by adding wind
shear stress in the simulation. Wind fields can be applied to the entire
model network using the wind property page.
The wind factors defined in the present page is applied to the simulations
in the following way:
– Topographical Factor:
The topographical factor is a multiplication factor which is multi-
plied to the Wind Velocity, and hence by eventual use of Local Fac-
tors it is possible to adjust the Wind Field applied with different
values at different locations.
– Friction Factor:
Friction Factor is a global parameter used in the calculation of
exchange between vertical layers and as such, the Friction factor is
only applicable to Stratified flows or MIKE 12 applications.
334 MIKE 11
Bed Resistance
level depending ‘Resistance factor’ which is specified for the cross sec-
tions in the cross section editor (.xns11 files).
336 MIKE 11
Bed Resistance
Zone separator lines must be defined in the User Defined Markers page
(see description in Activation of Bed resistance Triple Zone Approach
(p. 372)).
Global and local values of bed resistance for each zone can be specified as
described for the Uniform approach (see Figure 6.6).
Due to the special description in the friction term in the higher order fully
dynamic wave description the triple zone approach is only available for
fully dynamic and diffusive wave descriptions.
338 MIKE 11
Bed Resistance Toolbox
Five options are available in the Bed Resistance Equation combo box:
z Not Active
Bed resistance values used in the computation are those specified in the
Bed Resistance page (Uniform or Triple zone approach)
z n = 1/M = a*ln(VR)^b
The bed resistance is calculated as a function of ln(velocity*Hydraulic
Radius).
z n = 1/M = a*D^b
The bed resistance is calculated as a function of the depth.Please note,
that the Depth is calculated as an average depth (or ‘Hydraulic Depth’)
to avoid local deep parts of a section to control the resistance for the
entire cross section. D is therefore calculated as: D = Area / Width.
z n = 1/M = a*V^b
The bed resistance is calculated as a function of the velocity.
Note. To define the first line in the table, click the ‘Velocity’ bar in the
upper half of the page. Thereafter, press the <TAB> button and a new line
will be present in the grid in the upper part of the page.
All features (equations and table) can be defined both globally and locally.
If one of the equations has been applied, the user must define values for
the coefficient, a, and exponent, b. Additionally, a minimum and a maxi-
mum value must be specified to control, that bed resistance values calcu-
lated from the equations are inside the interval considered reasonable by
the user for the specific setup.
Note that when using a Chezy or Manning description the maximum bed
resistance requires the smallest Manning’s M or Chezy’s C. Similar for the
minimum bed resistance requiring the highest resistance number.
Further due to the special description in the friction term in the higher
order fully dynamic wave description. The bed resistance toolbox is only
available for fully dynamic and diffusive wave descriptions.
There are four possible flow description available in MIKE 11. The flow
descriptions can be selected globally for the system and/or locally for indi-
vidual branches. Locally specified flow descriptions must be specified for
the whole branch.
340 MIKE 11
Wave Approx
used when the simplifications/assumptions upon which they are based are
valid (see below).
The ‘high order fully dynamic’ flow description contains specific high
order and upstream centred friction terms in the momentum equation. This
modification typically allows simulations to be performed at longer time
steps than the ‘fully dynamic’ description.
The default value property page contains various parameters related to the
computational scheme. These parameters are essential for the simulation
and have been given default values. The parameters can be modified if
required.
Delta
The time-centring of the gravity term in the momentum equation. A velue
of 0.5 will produce the most accurate calculations provided there are no
numerical instabilities, but instabilites are less likely with higher Delta
values. For MIKE FLOOD models incorporating MIKE 11 channels, the
smaller timesteps are usual and Delta of about 0.85 generally provised
accurate calculations and stable computations. Large value of Delta
(towards 1.0) has a significant dissipative effect which can significantly
influence model dynamics and high values of Delta should therefore be
342 MIKE 11
Default values
Delhs
The minimum allowable water level difference across a weir. To obtain a
steady solution for differences below this limit a linear flow description is
used.
Delh
The Delh factor controls the dimensions of an artificial ‘slot’, which is
introduced to a cross section to prevent ‘drying out’ of the section. The
artificial slot is a small void introduced at the base of the section and
allows a small volume of water to remain in the section preventing com-
putational instabilities at low flows. The slot is inserted at height Delh
above the river bottom and extends to a depth of 5.Delh below this level.
Alpha
The velocity distribution coefficient used in the convective acceleration
term of the momentum equation.
Theta
A weighting factor used in the quadratic part of the convective accelera-
tion term of the momentum equation.
Eps
The water surface slope used in the diffusive wave approximation. If the
water surface slope becomes greater than EPS, the computational scheme
will become fully forwarded upstream. The parameter can be used to con-
trol the stability of the computation.
Dh node
Not used in the present version of the Engine.
Zeta min
The minimum head loss coefficient allowed in the computation of flow
over structures.
Struc Fac
Not used in the present version of the Engine
6.6.2 Switches
Node Compatibility
This switch should be set to water level since the energy compatibility is
not implemented and available in the present calculation engine.
344 MIKE 11
Quasi Steady
Relax
Weighting parameter used in the quasi-steady solution. For single
branches without bifurcation the value should be 1. In more complex sys-
tems the value should be less than 1.
Target_Branch
Computed water levels/discharges are shown on the screen at each itera-
tion for branch number equal to ‘Target Branch’. No computations are
shown if ‘Target Branch’ is negative.
Beta_Limit
Factor used to avoid underflow in ‘horizontal’ branches.
Fac_0
Factor used to control the stop criteria for the discharge convergence test.
Qconv_factor
Q convergence factor used in the stop criterion for the backwater compu-
tation iterations.
Hconv_factor
H convergence factor used in the stop criterion for the backwater compu-
tation iterations.
Min_Hconv_In_Branch
Minimum stop criterion to avoid underflow.
Q_struc_factor
Q structure factor: Used to determine the discharge at structures where a
slot description is introduced due to zero flow conditions.
H_stop
Stop criteria in the water level convergence test. Used also by the quasi
two dimensional steady state solver with vegetation as the convergence
criteria in the outer loop.
346 MIKE 11
Heat Balance
The information needed for the heat exchange calculation are (information
is also needed in the boundary-editor):
Latitude (N pos.)
Latitude of the considered area. Used in solar radiation calculation.
Longitude (W pos.)
Longitude of the considered area.
Displacement in time
Summertime correction: +1hour if the clock is 1hour ahead.
Light attenuation
Attenuation of solar radiation in the water column. Used to distribute the
incoming solar radiation over the different layers.
348 MIKE 11
Stratification
I sun
E fac = -------- ( 1 – β )exp(– a ( D – z )) (6.1)
I tot
where Isun is the solar radiation, β is constant in Beer's law, (D-z) is dis-
tance from surface and a is light absorption.
Radiation Parameter A
Daily radiation under cloudy skies is determined by:
H n
------ = A + B ------ (6.2)
H0 Nd
Radiation Parameter B
See above.
Vaporization Parameter A
Vaporative heat loss is determined by:
q e = LC e ( A + BW 2 ) ( Q W – Q a ) (6.3)
Vaporization Parameter B
See above.
6.9 Stratification
The property page used for setting up stratified flow simulations is illus-
trated in Figure 6.12.
Note that if stratified flow is to be simulated then the specific branch must
be defined as being stratified.
No. of layers
Number of layers in the stratified branches. The same number of layers is
assumed in all stratified branches. The thickness of a layer is equal to the
local depth divided by the number of layers.
Density calculated
Tick means yes and no tick means no. If densities are calculated it is done
on the basis of the simulated water temperatures, and if not density is
assumed to be 1000 kg/m3.
Turbulence model
Viscosity
In the case one chose a constant viscosity under turbulence model it is the
used viscosity in the calculations.
350 MIKE 11
Stratification
Corrections (reductions)
352 MIKE 11
Time Series Output
z Dfs Timeseries Bridge: This will generate a .dfs0 file in the standard
DHI file format for time series. Files can be loaded into for instance the
time series editor, plot composer, MIKE View or used as simulated
data for the auto calibration tool.
z ASCII format 1 Timeseries Bridge. This will generate an ASCII file
with a column based format. There will be one time column and one
column for grid point item selected for output. The number of lines
will equal the number of saved time steps.
z ASCII format 2 Timeseries Bridge. This will generate an ASCII file
with a table based format. There will be one table for each time step
saved. Each table will have rows corresponding to the number of
selected grid points, and columns corresponding to the number of
selected items.
z Minimum: A separate table at the bottom of the file will show the min-
imum value of the output items.
z Maximum: A separate table at the bottom of the file will show the
maximum value of the output items.
z Time of minimum: A separate table at the bottom of the file will show
the time of minimum value of the output items.
z Time of maximum: A separate table at the bottom of the file will show
the time of maximum value of the output items.
z Delimiter in output file: Between each column a special character like
semicolon, comma etc can be requested. This makes import to for
instance Excel easier.
z Width of columns: The desired minimum width of each column can be
specified.
quency for the regular .res11 output file as specified in the simulation edi-
tor.
A user defined name for the item can be specified. If not specified, the
default name will apply.
The number of decimals will only apply for output in ASCII files.
6.11 Maps
MIKE 11 may be used to produce two dimensional maps based on the
one-dimensional simulations. The maps are made as dfs2-files (rectangu-
lar grids) and constructed through interpolation in space of the grid point
results. Thus the maps constructed in this way should be viewed as a two
dimensional interpretation of results from a one -dimensional model. For
MIKE 11 simulations producing a map as output, the maps will therefore
be calculated only within the extent of cross sections defined in the map-
ping area as illustrated in Figure 6.14.
354 MIKE 11
Maps
Figure 6.14 MIKE 11 Mapping feature; Illustration of area included in map calcu-
lation. Thick black line represents river line, cross section lines
along the river and the dashed polygon represents the area for the
MIKE 11 simulation and hence, the active area where simulation
results will be transferred into a map representation.
z water level
z water depth
z velocity
z velocity times depth
z advection dispersion component
z (h, p, q)
z minimum values
the overall minimum value throughout the simulation period for each
gridcell included in the map.
z maximum values
the overall maximum value throughout the simulation period for each
gridcell included in the map.
z dynamic
Time-varying grid-file - for animations etc.
Only dynamic maps require a time span to be specified.
Through the mapping page, the user specifies the area that is to be
mapped, the type of map to be produced and the desired result item.
Depending on the user selections additional information may be required
such as the advection dispersion component, the time span for which the
map is to be produced and the storing frequency.
The option named ‘(h, p, q)’ produces a map-file containing the three
results-items; Depth, P-flux and Q-flux. This map-type is primarily tar-
geted the creation of ‘seamless dfs2 result-files’ with results from a MIKE
Flood simulation (coupled 2D- and 1D-simulation) such that inside river
cross sections coverage, the results in the map is calculated from MIKE 11
and the results outside the river sections coverage is calculation results
obtained from MIKE 21 and transferred to the mapping creation routine
during simulation of the coupled model..
356 MIKE 11
Maps
used. This option may further be used for introducing surveyed cross sec-
tions into a DEM covering the surrounding area. Figure 6.16 shows an
example of a DEM generated using the map capabilities of MIKE 11. The
DEM consists of the coarse background DEM (input topography) which
has been superimposed with the topography found through interpolation
of the cross sections in the MIKE 11 set-up.
Figure 6.16 An example of generating a DEM using MIKE 11. The DEM is con-
structed based on a background DEM in the form of a dfs2 file and
surveyed cross sections found in the MIKE 11 set-up. To visualize
the DEM MIKE Animator has been used.
To activate the Map Generation in a simulation you must select the ‘Gen-
erate maps’ tick-box. This will open the property page for editing of data.
358 MIKE 11
Maps
z ‘Grid Size’
The grid cell size of the maps produced. Maps generated can only be
squared (equal spacing of dx and dy of the grid)
z ‘J cells no.’ and ‘K Cells. no.’
Number of grid cells in the J- and K-direction respectively. J-direction
is equivalent to the X-axis and K-direction is equivalent to the Y-axis)
z ‘Filename’
The filename for the map you are creating.
Note: The map-file generated is a dfs2-file and hence, you should
ensure that the filename you specify has an extension named ‘.dfs2’.
use the ‘Browse’ button (‘...’) to search for existing filename or define
a new name for the map to be generated:
z ‘Type’
The map-type to be generated; Minimum, Maximum or Dynamic
z ‘Item’
The Resultt Item to be included in the Map (Water Level, Water Depth,
Velocity, Velocity * Depth, AD Component, DEM, (h,p,q)).
z ‘Component’
If AD Component has been selected as the Map-Item it is required to
define which component number (or State Variable) the maps shall be
produced.
z ‘Storing Frequency’
The frequency between timesteps saved in a dynamic map.
The timestep in the generated map-file will be defined as the ‘Storing
Frequency’ multiplied with the specified result saving frequency in the
simulation file. (example; a simulation storing frequency of 2 hours
and a Map Storing Frequency of 6 will result in a map grid-file which
will contain map-results for every 12 hours).
z ‘Default Period’
The Default Period check-box relates only to Dynamic type of maps.
Activate the Default Period box to create dynamic maps throughout the
entire simulation period.
Deactivate the check-box to create dynamic maps for only a limited
period within the simulation period.
Reducing the time in which the maps are created have two main pur-
poses: a) reducing the overall simulation time as map processes are
only activated in a limited period of time, and b) reducing the file-size
of a dynamic map.
Activate the ‘Apply DEM input data for ground elevations’ tick-box to
activate this feature.
360 MIKE 11
Maps
– During simulation and map generation, flood maps for the area cov-
ered by the helping grid are generated by using the calculated water
level in the actual cross section and eventual flooding in the addi-
tional area is included in the generated map.
Figure 6.18 Map area including additional Area outside of cross section defined
mapping area.
Figure 6.19 Definition of grid used for mapping. The rotation angle is positive
clockwise. The definition of origo is also shown. The two blue lines
indicate the location of the river system.
Figure 6.20 Overview of the location of the map in the network file. Note that the
map is rotated with respect to the working area coordinates.
362 MIKE 11
Maps
projection as the working area and contain an item describing the ele-
vation.
5 Do not select the ‘Apply DEM data between cross sections’ tick box.
6 Browse for the file(s) containing the additional topography and set the
item number of the elevation. Please refer to the documentation for the
bathymetry editor on how to generate a dfs2 bathymetry file.
7 If there are multiple background topography files then the ranking
order needs to be specified. The lowest order number is given the high-
est priority.
8 Run the simulation.
9 The output DEM may now be viewed using the Result Viewer or the
grid editor. If installed MIKE Animator may also be used to view the
generated DEM (as presented in Figure 6.16).
364 MIKE 11
Reach Lengths
On the Reach lengths page the user specifies the reach lengths of the left
and the right overflow banks in the downstream direction (increasing
chainage). Note that the reach lengths are to be specified at all cross sec-
tion locations except for the cross section with the highest chainage in a
branch (the area downstream of such sections is either a node/junction or
beyond the model area). The reach lengths are defined as the distances to
the next downstream cross section (next chainage). Please refer to fig. 1
below. As an aid to the user the "Load branch and chainages" button will
populate the table with all rivers and h-point locations. Please note that the
functionality of this button is only available if the simulation file is open
and the data from the cross section and network editors are accessible.
Since the use of reach lengths is based on raw data the method requires
that there exist cross sections at all h-points in the grid. If this is not the
case please insert cross sections or increase max-dx in the network editor
to avoid the code generating interpolated h-points at run time. The grid
generating function in the network editor can be used to check whether
cross sections are present at all h-points
366 MIKE 11
Add. Output
Structures
Structure flow, area and velocity. In case of control structures the gate
level is also stored.
Velocity
Velocities are calculated as the discharge divided by the cross sectional
area.
Discharge
The discharge calculated at h-points is a weighting of up- and downstream
discharges calculated at Q-points.
Slope
The free water surface slope.
Top width
The channel width at the free surface level.
Radius
The resistance radius.
Resistance
The cross-sectional resistance (resistance number multiplied by the resist-
ance factor).
Conveyance
The conveyance
Froude number
Defined as:
Q
F = ---------------- (6.4)
A
A g ----
bs
Where F is the Froude number, Q the discharge, A the cross sectional area,
g the acceleration due to gravity and bs the channel width at surface.
Volume
The volume calculated around the H-grid point.
368 MIKE 11
Add. Output
Flooded Area
H-points: The flooded area of the water surface between two neighbouring
Q-points.
Total: The total surface water area for the river system.
Mass Error
The mass error is defined as the difference between the volume calculated
in the model and the true volume. At nodal points with more than two con-
nections the mass error is distributed uniformly between each connection.
Lateral Inflows
Lateral inflows due to boundary conditions, catchment runoff, Flood fore-
casting updating or coupling to MIKE SHE.
Energy level
Energy level at water level points.
An alternative option for Bed Shear Stress saved to the additional result
file es available. The MIKE11.ini file contains a parameter;
‘AddOutput_ShearStress_UseEnergyLevelSlope’ which controls whether
the bed shear stress shall be calculated in the most common form using the
friction slope as indicated above - or if the bed shear stress should be cal-
culated fro using the Energy level slope.
dE
τ = ρgR -------
dx
where E is the energy level and x is the gridspacing between calculation
points along the river.
CREATE_QSSVEG_VELOCITY_FILE=ON
With this setting 8 .txt-files are generated and saved in the working direc-
tory i.e. where the simulation file is stored. The files are titled:
370 MIKE 11
Flood Plain Resistance
The global resistance number is applied on all flood plains unless local
values are specified. Local values are linearly interpolated at intermediate
chainage values. The resistance number value -99 indicates that the flood
plain resistance should be calculated from the raw data in the cross-section
data-base.
Note. To define the first Marker in an empty page, click the ‘Mark title’
bar in the upper half of the page. Thereafter, press the <TAB> button and a
new line will be present in the grid in the upper part of the page as well as
a new column is introduced in the ‘location grid’ in the lower half of the
page. Write the name of the marker in the empty line in the upper grid, and
this name will automatically be transferred as the name of the column.
Markers can be defined as single points only and as markers defined along
a river stretch. The ’Interpolate’ column must be checked in case a linear
interpolation is requested on stretches between chainages and marker lev-
els defined in this page. In case a user defined mark should be presented
on the longitudinal profile as a single point (e.g. a bridge location or flood
mark indicator) the Interpolate check-box must be un-checked.
Please Note: In case the Triple Zone Approach has been activated and
zone separator lines are not defined for the entire setup, MIKE 11 uses the
uniform bed resistance values in the points where separator lines are not
defined. The resistance value used at these points is the value (global or
local) defined for the lower zone.
372 MIKE 11
Encroachment
rator lines has been defined in RIVER1 in the reach from chainage 0.0 to
5000.
6.17 Encroachment
Encroachment simulations are setup through this page. An encroachment
simulation consist of two or more simulations. The first simulation acts as
a reference simulation to which all other results are compared. The refer-
ence simulation is set up as any ordinary steady state simulation.
Note that only MIKE 11’s default steady state solver may be used. Further
since the encroachment module utilizes the steady state solver, the instal-
6.17.1 Iteration
6.17.2 Location
The location of the encroachment station is entered here through a river
name and a chainage. If a location is entered for which no corresponding
cross section exists a warning is issued at run time and the station will be
ignored in the subsequent simulation.
374 MIKE 11
Encroachment
Method
A total of five different methods are available:
Sides
It is possible for the encroachment to take place on both sides of the main
channel or only on one of the sides. For this purpose the sides combo box
may be used.
Note: If the method chosen for encroachment is ’Fixed width’ then the
sides switch is automatically set to both sides. Since a fixed width
encroachment only makes sense if both sides are to be encroached.
Reduction type
The way that the conveyance reduction should be accomplished is speci-
fied here. Three possibilities are available:
Total reduction
If the reduction type is set to either ’Equal’ or ’Relative’ this field
becomes active. The total required conveyance reduction is entered here in
percentage.
376 MIKE 11
Encroachment
Water level target: The encroachment may be carried out so that only
a water level target is considered. This strategy is achieved by setting
the water level change to a non-zero value and the energy level change
to zero.
Water level target and energy level target: The encroachment is car-
ried out so that a water level target is met. Once this has been achieved
the energy level is checked. If the energy level is above the energy
level target the code will reconsider the encroachment and try to satisfy
the energy level request instead. This strategy is achieved by setting
both the water level change and the energy level change to non-zero
values.
Please note that the position of the encroachments are found through an
iterative procedure. This procedure considers each cross section individu-
ally starting downstream and working upstream. To ensure that this
method is successful do not use method 5 for river reaches which form
part of a loop in a network. Further method 5 is designed for encroaching
river reaches where the discharge distribution can be determined a priori,
thus the method will be less successful for networks having river bifurca-
tions in a downstream direction as opposed to bifurcations in upstream
directions. Finally it should be mentioned that not all user specified targets
can be reached. If this is the case the code will issue a warning and return
the encroachment which is closest to the requested target.
For encroachment simulations only the initial start time in the simulation
editor is used. This start time is used for determining the boundary values
in the river set-up. Note that constant boundary conditions in MIKE 11 are
specified by the use of non-varying boundary conditions in the boundary
editor.
The present page contains variables for specifically the model parameters
related to entrainment and exchange between the two layers in MIKE 12.
378 MIKE 11
MIKE 12 Initial
380 MIKE 11
Mixing Coefficients
Used only in conjunction with the quasi two dimensional steady state veg-
etation module. This menu is used for setting the mixing coefficients
between adjacent panels in the river cross sections. Both global and local
values may be set here.
Important! The table should start with b/B=0 and end with b/B=1 and all
intermediate values of b/B must be monotonically increasing. If the table
does not meet this criteria a warning is issued and the default settings are
used.
Note. To define the first line in the table, click the ‘b/B’ bar in the upper
half of the page. Thereafter, press the <TAB> button and a new line will
be present in the grid in the upper part of the page.
Expansion/Contraction f
The mixing coefficients at a water/water boundary at a dead water inter-
face.
6.20.2 Location
The river name and location (chainage) is displayed here.
The tab is illustrated in Figure 6.30 with all the different features all of
which are described below.
382 MIKE 11
W. L. Incr.- Curve
6.21.1 General
At the bottom of the editor a table is displayed with river name, chainage
and the four parameters appropriate for the determination of the water
level increment. The parameters which are not greyed may be edited.
1 Average Range.
2 Curvature Radius.
3 Water Surface Width.
4 Velocity.
Note! If either 2 or 3 is ticked the velocity is also automatically ticked.
In the column "Average Range" the user can control the calculation of the
curvature radius. If the average range is set to "None" no water level incre-
ment due to curvature applies. For other values of average range a curva-
ture radius is initially calculated or assigned (depending on the what's
selected in the group box "System Definition") individually in each h-
point. If the average range equals "Single" the curvature radius is kept
unchanged, otherwise this is averaged over a number of h-points. Consec-
utive h-points with the same average range setting is lumped together
when calculating an average curvature radius. "Multiple 1/2" is used for h-
points to be included in both the upstream and downstream averaging
reach.
The tab is illustrated in Figure 6.31 with all the different features all of
which are described below.
384 MIKE 11
W. L. Incr.- Sand Bars
6.22.1 General
At the bottom of the editor a table is displayed with river name, chainage
and parameters appropriate for the determination of the water level incre-
ment.
z Bed slope, low water channel width and water area of annual maxi-
mum discharge
or
386 MIKE 11
ADVECTION-DISPERSION EDITOR
387
388 MIKE 11
Module Description
7 ADVECTION-DISPERSION EDITOR
The AD Editor is used in conjunction with the following modules:
∂AC ∂QC ∂ ∂C
----------- + ------------ – ----- AD ------- = – AKC + C 2 q (7.1)
∂t ∂x ∂x ∂x
where
C : concentration
D : dispersion coefficient
A : cross-sectional area
C2 : source/sink concentration
q : lateral inflow
x : space coordinate
t : time coordinate
390 MIKE 11
Components
The module can also be used when resuspension of sediment affects water
quality. This is because the resuspension of cohesive sediment often gives
rise to oxygen depletion due to the high organic content and associated
oxygen demand (COD) in the cohesive sediment. Likewise resuspension
of cohesive sediment can give rise to heavy metal pollution since heavy
metals adhere to the sediment.
7.2 Components
Component names and numbers must be specified in this dialog.
Please note, that the combination presented in the figure is primarily for
demonstration purpose. In a real application it would not be recommended
for several reasons include two different types of Cohesive Sediment com-
ponent types in the same simulation.
Please Note: For Water Quality simulations using the ECO Lab Frame-
work it is not required to define any Components in this property page.
Component definitions (names) in an ECO Lab (Water Quality) simula-
tion are taken from the ‘State Variables’ definition in the applied
Ecolab11-file.
Component
Here one or more component names are defined by a user-defined string.
Units
The requested unit for each of the components is selected through the
drop-down selection box. Options available:
392 MIKE 11
Components
Type
– Normal: A component used for standard Advection-Dispersion
simulations.
– Single layer cohesive: A component used only in the single layer
cohesive sediment transport model (the CST module).
– Multi cohesive: A component used in the multi layer cohesive sed-
iment transport model (the ACST module).
– Non-cohesive: A component defined as a Non-Cohesive is actually
not simulated using the Advection-Dispersion equation - but rather
using traditional sediment transport equations.
Please be informed that this option for simulation a non-cohesive
component within the AD Parameter Editor Framework is actually
a leftover from earlier versions of MIKE 11 where combined Sedi-
ment Water Quality interaction was simulated through this inter-
face. However, now WQ-simulations are made using the ECO Lab
Framework and the non-cohesive component type is therefore
rather un-useful in the present release.
Simulating a Non-cohesive sediment component through the AD
Parameter Editor does Not include morphological updating and as
such it is not very userful anymore (and will be removed in a future
release).
For Non-cohesive sediment transport simulations, please apply the
ST-Module instead.
7.3 Dispersion
The dispersion coefficient, D, is described as a function of the mean flow
velocity, V, as shown below.
b
D = aV (7.2)
Both the ‘dispersion factor’ and the ‘dispersion exponent’ can be speci-
fied. If the dispersion exponent is zero then the dispersion coefficient D
becomes constant (equal to the dispersion factor). By default the disper-
sion is zero (i.e. there is only advective transport and no dispersion). The
‘Minimum dispersion coefficient’ and the ‘Maximum dispersion coeffi-
cient’ parameters are used to control the range of the calculated dispersion
coefficients.
394 MIKE 11
Dispersion
Global values
The Global dispersion values are defined in this group of fields.These val-
ues will be defined in all branches and calculation points unless Local Val-
ues are defined additionally.
Dispersion factor
The dispersion factor equivalent to parameter a in equation (7.2) is
defined here.
Exponent
The dispersion exponent equivalent to parameter b in equation (7.2) is
deifned here.
Local Values
Mike11 will use the values specified under global values except for those
places were local values have been specified.
River Name
Name of the river with local dispersion values.
Chainage
Chainage in river with local dispersion values
Dispersion factor
Local value of the dispersion factor
Exponent
Local value of the dispersion exponent
Example
In Figure 7.2 both global and local values are entered. In ‘RIVER 1’ the
dispersion coefficient is globally set to 10 m2/s (independent of the flow
velocity because b equals 0). Local values are defined in the reach
between chainages 1000 m and 2500 m, where the dispersion coefficient
is dependent on the velocity (D= 15V, 5 < D < 25)
396 MIKE 11
Init. cond.
Component
Here the component in question is selected. It is possible to choose
between the components defined in the Components property page by
selecting from the drop-down selection box present in the ‘Component’
column.
Concentration
The Initial concentration-value for the specific component is defined here
Global
This box must be checked if the value entered in the Concentration field
should be used as a global value. If it is left unchecked the value will be
used as a local value.
River name
The name of the river with the local initial value.
Chainage
The chainage in the river with the local value.
Example
In Figure 7.3 initial conditions is defined for the component ‘SALINITY’.
The initial, Global concentration of SALINITY is set to 35 PSU for the
entire river network. However, the initial concentration of SALINITY var-
ies linearly from 20 PSU to 25 PSU to 27 PSU up to 33 PSU along the
reach defined as river branch named ‘RIVER 1’ at the chainages 5000m to
10000m to 12000 m to 15000 m.
Component
Here the component in question is selected. Presently only temperature
can be selected.
Conc. S
Temperature at the surface.
Conc. 2
Temperature at layer k2 above the bottom.
Conc. 3
Temperature at layer k3 above the bottom.
Conc. B
Temperature at the bottom.
k2
Layer number above the bed
k3
Layer number above the bed.
Global
This box must be checked if the value entered in the Concentration field
should be used as a global value. If it is left unchecked the value will be
used as a local value.
398 MIKE 11
Decay
River name
The name of the river with the local initial value.
Chainage
The chainage in the river with the local value.
7.5 Decay
This page contains information for non-conservative components. These
components are assumed to decay according to a first-order expression:
dC
------- = KC (7.3)
dt
Component
Here the component in question is selected. Wither type in the Component
name manually or select from the Components drop-down box present in
the Components column
Decay const
Here the value of the decay constant are entered (default unit is [per hour])
Global
This box must be checked if the value entered in the Decay const. field
should be used as a global value. If it is left unchecked the value will be
used as a local value.
River name
The name of the river with the local initial value.
Chainage
The chainage in the river with the local value.
400 MIKE 11
Cohesive ST
Example
In Figure 7.4 the component COMP1 has been selected to be non-conserv-
ative. The decay constant is 1.00 globally in the river network and has a
value of 2.00 in RIVER 1 between the chainages 10000 m and 20000 m.
7.6 Cohesive ST
Data used for each of the two cohesive sediment transport models availa-
ble are entered on this page. When using the cohesive sediment transport
models (either the simple or the advanced) all components specified in the
AD editor must be defined as ‘Single layer cohesive’ or ‘Multi layer cohe-
sive’ in the Components dialog.
Figure 7.5 The Cohesive sediment property page when a single layer model is
selected.
Below the parameters that apply to the ‘Single layer cohesive’ sediment
transport model are described.
Fall Velocity
w0
The free settling velocity.
Deposition
Time centring
This centring factor used in the deposition formula. Typical range is: 0.5-
1.0.
Erosion
Erosion coefficient
The erosion coefficient is applied linearly in the erosion expression. Typi-
cal range: 0.20 - 0.50 g/m2/s.
Erosion exponent
The erosion exponent describes the degree of non-linearity in the rate of
erosion. Typical range: 1-4.
Overview
At the bottom of the property page a overview table is shown.
Global
If this box is checked the entered parameters are used globally.
River Name
The name of the river for which the data applies.
402 MIKE 11
Cohesive ST
Chainage
The chainage of the river for which the entered data applies.
Figure 7.6 The cohesive sediment property page when a multi layer model is
selected.
Below the parameters that apply to the ‘Multi layer cohesive’ sediment
transport model are described.
Fall velocity
C-offset
Concentration limit for flocculation affected settling velocity. For higher
concentrations the settling velocity is affected by hindered settling.
g
Exponent used in the settling velocity expression. Typical range: 3 - 5.
m
Exponent in the settling velocity expression for concentrations below C-
offset.
w0
Free settling velocity. Typical range: 0.0025 - 0.01 m/s.
swi
Sediment volume index used in the settling velocity expression.
Deposition
Time centring
This centring factor used in the deposition formula. Typical range is: 0.5-
1.0.
Erosion
Erosion coefficient
The erosion coefficient is applied linearly in the erosion expression. Typi-
cal range: 0.20 - 0.50 g/m2/s.
Erosion exponent
The erosion exponent describes the degree of non-linearity in the rate of
erosion. In case that ‘Instantaneous erosion’ of layer 1 is selected the ero-
sion exponent is not applicable for layer one. Typical range: 1-4.
Consolidation
Transition rates
The consolidation of the sediment layers is described by transition rates
between the layers. The transition rates include hindered settling and con-
solidation. Typical ranges:
404 MIKE 11
Cohesive ST
Overview
At the bottom of the property page a overview table is shown.
Global
If this box is checked the entered parameters are used globally.
River Name
The name of the river for which the data applies.
Chainage
The chainage of the river for which the entered data applies.
WC τ
- 1 – ------- , τ ≤ τ cd
S = -------- (7.4)
* τ
h cd
where,
The bed shear stress can be given by the Manning formula (as an exam-
ple):
2
V
τ = ρg -----------------
2 1⁄3
(7.5)
M h
where,
Substituting the bed shear stress into the deposition equation results in the
following equation:
WC- V 2
S = --------
*
1 – -------- , V ≤ V cd (7.6)
h V cd
where,
*
M τ n
S = ------- 1 – ------ , τ ≤ τ ce (7.7)
h τ ce
where
406 MIKE 11
Cohesive ST
n erosion exponent
*
M V 2 n
S = ------- 1 – ------- , V ≥ V ce (7.8)
h V ce
where
τb
S d = W s 1 – -------- c , τ < τ c, d (7.9)
τ c, d
where
n n+1
c = ( 1 – θ )c j + θc j (7.10)
where:
j spatial index
n time index
n
S E = E o ( τ b – τ c, e ) , τ b > τ c, e (7.11)
where,
408 MIKE 11
Sediment layers
Location
The Location specification includes the definition of locations of initial
values of sediment layers as well as the values for sediment layer defini-
tions.
Component
Select from the drop-down box one of the allowed component types: Sin-
gle Layer Cohesive, Multi Layer Cohesive and Non-Cohesive. If no com-
ponents have been selected as one of these types, the selection box will
appear empty.
Layers
Only available when Component is chosen as a Multi Layer Cohesive
component. The user can select between Upper, Middle and Lower repre-
senting the three layers in the Multi Layer Cohesive model. Parameters
must be specified for each of the layers.
Table
Only applicable for Multi Layer model components. Instead of giving the
initial conditions in Height (p. 410), a more detailed initial condition can
be specified using a width-Height table, see Parameters (p. 410).
Height
Although the header says 'Height' the initial data should be entered as vol-
ume of sediment per length of river. In order to convert this initial data
into an amount MIKE 11 uses the porosity and the relative density speci-
fied in the Non-cohesive ST (p. 414) property page.
Density
The density of the layer.
Global
If this box is checked the entered parameters are used globally.
River Name
The name of the river for which the data applies.
Chainage
The chainage of the river for which the entered data applies.
Parameters
For multi layer components a volume width relation can be entered. The
width in this relation is the width of the cross section, the volume is the
volume of sediment per length of the river. It is hereby possible to vary the
thickness of the sediment layer along the transverse of the cross section.
410 MIKE 11
Sediment layers
– Filename are ‘LayerInp.00?’ where the last ‘?’ in the file extension
refers to the Component number for which the layer thickness shall
be defined. That is, if an ACST component is defined e.g. as com-
ponent number 2 in the list of defined components, then the
filename must be; ‘LayerInp.002’
– The ASCII files must be located in the simulation folder. That is,
the folder where the .sim11 file for the actual simulation is located.
– If the simulation engine detects a file with the requested name, a
message box will appear as presented in Figure 7.8 below. From the
dialog it can now be chosed whether the initial layer depths should
be read from the external file (Press Yes) or using the definitions
from the AD Parameter file (Press No).
Figure 7.8 Confirmation dialog for reading initial sediment layer depths
from external ASCII file
Figure 7.9File format of ASCII input file for initial cohesive sediment layer
depths.
Note that the LayerOut.00? file(s) always contains the sediment layer val-
ues as defined in the AD Parameter, so the content of this file will not
412 MIKE 11
Sediment layers
change regardless of which layer input file you eventually use for the sim-
ulation.
Different to the LayerOut.00? file, this CST-file will contain the values for
the sediment layers at the time where the simulation ended, which could
eventually be useful initial conditions for consecutive simulations.
If Hotstart is used as initial conditions for the AD-simulation then the Lay-
erInp.00? file option is not used as initial conditions will be applied from
the hotstart result-file and concequently, the message-box (Figure 7.8) will
not appear for AD-simulations using hotstart.
7.8 Non-cohesive ST
414 MIKE 11
Non-cohesive ST
Model constants
Model Type
A pop down menu provides a choice from two types of sediment transport
formulations; the Engelund-Fredsøe and the van Rijn model.
Fac.1
Calibration factor for bed load transport. The calculated bed load is multi-
plied by the calibration factor.
Fac.2
Calibration factor for suspended load transport. The calculated suspended
load is multiplied by the calibration factor.
Beta
Dynamic friction factor used in the Engelund-Fredsøe model.
Kin.visc.
The kinematic viscosity of water.
Porosity
The porosity of the sediment.
Rel.dens.
The relative density of the sediment.
Thetac
Shield’s critical parameter. Typical range: 0.04 - 0.06.
Data
Component
Here a Non-cohesive component is selected.
grain size
The D50 value.
st dev.
Standard deviation in the grain size distribution.
Global
If this box is checked the entered parameters are used globally.
River Name
The name of the river for which the data applies.
Chainage
The chainage of the river for which the entered data applies.
Mass
The mass in the system. Given in the units specified on the ‘Components’
property page. Total and total accumulated as well as grid and grid accu-
mulated values can be selected.
Mass balance
The mass balance is given in o/oo (per thousands). Total and total accumu-
lated as well as grid and grid accumulated values can be selected.
1. order decay
The 1st order decay is given in the units specified on the ‘Components’
property page, per second. Total and total accumulated as well as grid and
grid accumulated values can be selected.
416 MIKE 11
Dispersion, MIKE 12
Mass in branches
The mass in river branches given in the unit specified on the ‘Compo-
nents’ property page. Total and total accumulated values can be selected.
Transport, total
The total transport is given in the unit specified on the ‘Components’
property page, per second. Grid and grid accumulated values can be
selected.
Dispersive transport
The dispersive transport is given in the unit specified on the ‘Components’
property page per second. Grid and grid accumulated values can be
selected.
Convective transport
The convective transport is given in the unit specified on the ‘Compo-
nents’ property page, per seconds. Grid and grid accumulated values can
be selected.
Please note, that this property page are only visible if a MIKE 12 module
is included in your license and most MIKE 11 users will therefore not see
this tab-page when working with the AD Parameter editor.
Please note, that this property page are only visible if a MIKE 12 module
is included in your license and most MIKE 11 users will therefore not see
this tab-page when working with the AD Parameter editor.
418 MIKE 11
AD Simulations, Additional features
The ULTIMATE method is very strong and efficient to resolve very sharp
concentration gradients at neighbouring calculation points, where the
standard Finite Difference method typically make a solution with signifi-
cant numerical ‘wiggles’. The Ultimate method is able to keep the sharp
gradients without generating these wiggles as illustrated in Figure 7.14.
Please note, that it is recommended to copy the mike11.ini file from the
bin-folder to your project simulation folder (must be located where you
have your sim11-file). The simulation engine searches for the MIKE11.Ini
file locally first, and secondly for the global file in the bin-folder. By cop-
ying the file to the local folder, you ensure, that the changes you eventu-
420 MIKE 11
AD Simulations, Additional features
ally make in the ini-file only concerns the present project, which is often
required as changes in the ini-file are possibly required and beneficial for
one project, whereas different values or parameters must be altered in
other projects.
AD-Variables no 1:
‘ADLOG=On’
AD-Variables no 3:
‘MORPH_UPD_IN_COHE_3LMODEL=OFF’
422 MIKE 11
WQ ECO LAB EDITOR
423
424 MIKE 11
Model Definition
ECO Lab utilises so-called Template files with the filename ‘.Ecolab’ (in
the following named ECO Lab File). The ECO Lab template files are cus-
tomized collections of equations and parameters required for a specific
type of Water Quality simulations. The user may use one of the predefined
ECO Lab Templates, which is installed together with the MIKE Zero
installation - or may choose to develop own model concepts. The module
can describe dissolved substances, particulate matter of dead or living
material, living biological organisms and other components (all referred to
as state variables in this context).
Having selected the .ecolab file, a brief summary of the contents of the
model is shown in the Dialog as presented in Figure 8.1.
Please note that every time a new .ecolab file is selected, the specifications
of all the remaining ECO Lab Dialogs are reset to default values.
426 MIKE 11
State Variables
The selection of the Time Step of the ECO Lab model, and hereby the
Update Frequency, has to be based on considerations of the time scales of
the processes involved. Please notice that this selection can be rather deci-
sive for the precision of the numerical solution as well as for the CPU time
of the simulation. A large Update Frequency will decrease the precision as
well as the CPU time. It is therefore advisable to perform a sensitivity
analysis on the Update Frequency before making the final selection.
For each State Variable its initial value within the model area should be
specified . It can be specified in one of two ways: As a constant value
applied to all points in the area through the ‘Value’ column or with local
exceptions. Local exceptions for parameter values are defined by opening
the local variables sub-dialog. Press the button in the ‘Local’ column
and a local values sub-dialog will open as presented in
To generate table-lines where local values can be defined, use the <TAB>
key at the last column of the table or press the <Insert> key. Local values
can be defined for one or multiple of the state variables in this page. Select
428 MIKE 11
Constants
the variable for which the local variables must be defined from the drop-
down selection box containing all the model state variables in the
‘Description’ column and specify river location (Rivername and Chain-
age) and the local value of the specific parameter.
Note, that if you want to specify local variations at the same locations for
different state variables, you must define this in individual lines for each
State Variable.
8.3 Constants
The Constants are defined as any input parameter (physical constant, coef-
ficient, rate, etc.) in the ECO Lab model, which is constant in time. The
Constants are essentially divided into two groups:
Constants can be declared as spatial varying in the ECO Lab template file
and for these, it is possible to specify local variation of the parameter val-
ues through the ‘Local’ column. Definition of local values are made in the
same was as described for defining local values of State Variables.
8.4 Forcings
The Forcings are defined as any input parameter (physical property, rate,
etc.) in the ECO Lab model, which is varying in time. Examples of a Forc-
ing are: Temperature, salinity, solar radiation and water depth. The Forc-
ings are essentially divided into two groups:
The built-in Forcings are automatically provided by the model system dur-
ing execution, whereas the user-specified Forcings have to be specified in
the present Dialog. Depending on the Spatial Variation of the Forcing, as
defined in the ECO Lab model, it can be specified as a “Constant value”
or a “Type 0 data file”.
Alto note that the forcings called "Water depth " and "Water layer height"
both hold the actual hydraulic radius as a measure for the water depth
430 MIKE 11
Auxiliary Variables
rather than for instance the maximum or the average depth across the cross
section.
8.6 Processes
Processes defined in the ECO Lab file, and at the same time marked as
OUTPUT, can be saved in the additional output file named '<AD
filename>WQAdd.res11'. Also for the processes the author of the ECO
Lab file has decided which of the processes that you can store as addi-
tional output. Simply tick the processes you want to save.
432 MIKE 11
SEDIMENT TRANSPORT EDITOR
433
434 MIKE 11
9 SEDIMENT TRANSPORT EDITOR
The MIKE 11 non-cohesive sediment transport module (NST) permits the
computation of non-cohesive sediment transport capacity, morphological
changes and alluvial resistance changes of a river system.
Some of the sediment transport formulas and other features of the Non
Cohesive Sediment Transport module have been developed in cooperation
with CTI Engineering CO., Ltd., Japan.
The explicit sediment transport mode is active if the check box; ‘Calcula-
tion of Bottom Level’ is un-checked (in the ‘Transport model’ page).
436 MIKE 11
Sediment grain diameter
No general guidelines can be given for the preference of one model over
another, as the applicability of each depends on a number of factors. Fur-
ther details can be found by consulting the NST Reference Manual.
Features and usage of the ST Parameter Editor pages are described below.
The Standard Deviation defined in this page is only applied for two of the
Transport Models available;
– the Van Riin model utilises the standard deviaiont for calculating
the mean diameter and particle size diameter used in the equation
set of the Van Rijn model, and
– the Engelund-Fredsoe model utilises the standard deviation in the
Suspended Load calculation to determination the grain size diame-
ter applied in calculating the effective fall velocity
– The other transport models available does not include the standard
deviation in the respective set of equations and hence, adjusting the
standard deviation when applying any of these other transport mod-
els will not affect the simulation results.
Figure 9.2 shows an example where the sediment grain diameter is glo-
bally set to 0.5 mm. This value will be used in the entire river network
except for the reaches ‘RIVER1’ between 1000 m and 2500 m, where the
local grain diameter varies linearly between 1.2 and 1.5 mm. and between
2500 m and 4400 m where the grain diameter varies linearly between 1.5
and 1.1 mm. At the same chainages, the standard deviation varies linearly
between 1.2, 1.2 and 1.0.
438 MIKE 11
Transport model
the model type (Total load or bed load/suspended load model) and the
appropriate transport model(s) as well as adjusting transport parameters if
required prior to making any ST-simulation.
Figure 9.3 shows an example of how to set the transport model type and
appropriate parameters in the dialog. In this example, the bed load trans-
port will be calculated using the ‘Engelund and Fredsoe’ model and the
suspended load transport calculated using Van Rijn formula. Morphologi-
cal computation is selected as the check box for ‘Bottom Level’ is acti-
vated, but there will be no computing of the bed shear stress.
Figure 9.3 Sediment transport model parameters. Selection of bed load and
suspended load models. Morphological updating is included but Bed
Shear stress updating excluded in this example.
the formulation of - and division between - suspended and bed load char-
acteristics.
Rel. density
Relative density (or Specific gravity) of the sediment. Relative dentisy
determined as; ρ sediment ⁄ ρ water .
Kin. Viscosity
Kinematic viscosity of water.
Note, that the Kinematic Viscosity must be specified as 'value . 10-6' m2/s
(for SI-Units). That is, if a value of 0.000001 m2/s should be used, in the
dialog, you must specify 1.0.
Beta
Dynamic friction coefficient used in the Engelund-Fredsoe model. Please
consult the MIKE 11 Reference Manual for further details.
Theta Critical
Critical Shields' parameter for incipient motion of sediment.
Gamma
Calibration parameter applied to suspended load with the Engelund-Fred-
soe model when calculating the height of sand dunes. For further details
please consult the MIKE 11 Reference Manual.
Ackers-White
Switch used in the Ackers-White model indicating whether the applied
grain size represents d35 or d65.
440 MIKE 11
Transport model
Storing
– Bed / Suspended load
Storing of suspended load and bed load as individual result items in
the ST result file from a simulation. This feature is only applicable
for those of the transport models which separates the sediment
transport into bed load and/or suspended load components.
dH/dZ
Calculation parameter for the morphological model.
PSI
Centring of the morphological computation scheme in space.
FI
Centring of the morphological computation scheme in time.
FAC
Calibration parameter for computation of derivatives in the morphological
model.
Note that this parameter implicitly defines the step length for a number of
numerical derivatives. For this reason the parameter must be greater than
unity. If this is not the case MIKE 11 sets its value equal to 1.01 internally.
Porosity
Porosity of the sediment (pore volume/total volume).
Omega
Calibration parameter for the resistance number. (ResistanceST =
OMEGA * ResistanceHD). Note that Omega is applied to the resistance
number, which is internally in the Mike11 engine calculated as Mannings
M. Hence, an increased value of Omega will result in a smaller overall
resistance (higher value of M) and a reduced value of Omega will result in
a higher overall resistance.
Minimum/Maximum
Minimum/maximum limits for the calculated resistance number in the
computations.
Also note, that only few of the authors of sediment transport formulas has
considered formulations for updated alluvial resistance. Hence, MIKE 11a
applies alternative formulations for these formulas as presented in
Table 9.1
442 MIKE 11
Transport model
Table 9.1 Methods applied for Updating bed shear stress in sediment trans-
port formulas available.
Please note, that for most of the methods there are no specific formulation
for a bed resistance update method and the Engelund-Hansen method has
been adopted in these formulations.
Engelund-Fredsoe model
When selecting the Engelund-Fredsoe transport model, dune height and
dune length are computed - if calculation of Bed Shear Stress is included.
Therefore, an additional property page; ‘Initial Dune Dimensions’ is made
visible in the ST Editor when either the bed load or suspended load trans-
port model is chosen as Engelund and Fredsoe (see Section 9.8).
Smart-Jaeggi model
When selecting the Smart-Jaeggi transport model, the model parameters
must be edited as for all other transport models. Additionally, user-defined
coefficients and exponents used in the Smart-Jaeggi formulation can be
specified. Therefore, when selecting the transport model for Total Load as
‘Smart and Jaeggi’ values for coefficients and exponents can be edited in a
separate dialog.
Figure 9.4 Additional dialog for defining Smart and Jaeggi model factors.
Coefficients and exponents are essential for the Smart and Jaeggi transport
model and a simulation should therefore not be performed until this dialog
has been edited.
z Method no 1.
Deposition in horizontal layers from the bottom. Erosion proportional
with depth below bank level
z Method no 2
Deposition and erosion uniformly distributed below the water surface.
No deposition and erosion above.
z Method no 3
444 MIKE 11
Transport model
z Method no 4
Deposition and erosion uniformly distributed over the whole cross sec-
tion (i.e. below the bank level).
z Method no 5
Deposition and erosion proportional with depth below bank level
Figure 9.5 Options available for updating of bed levels and cross sections
parameters.
Please note, that it is not unfortunately not possible to extract and present
the updated cross sections shapes due to morphological updating. from
any result-file. Analysis of impacts on cross section parameters (Area,
Width, Radius) can be made if these parameters are saved through the HD
additional output feature. If saved to the HD Additional Result-file then
the Area, Width and Radius will all contain the updated parameter values
as calculated during the sediment transport calculation.
446 MIKE 11
Data for graded ST
It is possible to specify a lower limit for the active layer depth ('Min. depth
active layer') and an initial depth for the passive layer.
Please note, that Graded Sediment Transport option is not available for
the Meyer, Peter & Mueller sediment transport formula.
(n + 1) n Qm
QSt m = ∑ QSt ⋅ -------------------- (9.1)
up ∑ Qi
down
Where:
(n + 1)
QSt m : Sediment transport at timestep (n+1) for outflow branch ‘m’
448 MIKE 11
Preset distribution of sediment in nodes
n
∑ QS t : Acc.sediment transport in upstream branch-ends
up (= outflow from node)
Bm
(n + 1) nAm ⋅ Qm
QSt m = ∑ QSt ⋅ -------------------------------
-
Bi
(9.2)
up ∑ i i A ⋅ Q
down
Where
Index ‘m’ refers to branch ‘m’ and index ‘i’ refers to inflow branches to
the specific node.
The coefficients and exponents are given for each branch, specified by its
upstream and downstream chainage as linked to the node as illustrated in
the figure below, where ‘River 1’ discharges into a node (defined by the
‘Downstream’ type) and ‘River 2’ and ‘River 3’ respectively receives
water from the node (= outflow branches in the noe. Defined by the
‘Upstream. type).
Sediment can be transported into a passive branch, but the sediment even-
tually transported into a passive branch is ‘trapped’, as no sediment can be
transported inside and out of a passive branch.
450 MIKE 11
Non-Scouring Bed Level
The Thickness of active layer is used in the Graded sediment transport cal-
culations. Default formulations in MIKE 11 define the thickness of active
layer as half the dune height, but this value can also be user-defined
through this property page. The value must be given as a depth. That is; a
height above bottom of river bed.
The Non scouring bed level item gives a possibility for the user to define
levels (global and/or locally) where a non-erodible surface is present. If,
during a morphological simulation, bed erosion occurs and the bottom of
the bed reaches the defined Non scouring bed level, no further bed erosion
will take place.
452 MIKE 11
Initial dune dimensions (Engelund Fredsoe model)
Figure 9.11 shows an example where the global dune height has been set
to 0.25 m, and the global dune length has been set to 12.50 m. These val-
ues will be used in the entire river network, except in the reach ‘RIVER1’,
between chainage 5.000 km and 10.000 km, where the dune height varies
linearly between 0.25 m and 0.40 m.
If no dune dimensions are given, or the dune height and length equals
zero, then the dune height will be calculated as the water depth divided by
6 with a dune length of 15 times the water depth.
The default MIKE 11 formulation for defining the transport layer depth
(active layer) is the water depth divided by (D/6). However, this formula-
tion only gives realistic results in case of slowly varying flow and similar
water depths in the entire modelled river reaches.
454 MIKE 11
ST Simulations; Additional features
Figure 9.13 Fileformat of ASCII file ‘Shield.Inp’ for defining alternative options
for shielding (hiding/exposure) of grains.
ln ( sf ) 2
Shield Type = 1 : τ c, shield = τ c ⋅ ----------------------------------
ln ( sf ⋅ d i ⁄ d )
m
d m sf
Shield Type = 2: τ c, shield = τ c ⋅ ------
di
where:
Please note, that it is recommended to copy the mike11.ini file from the
bin-folder to your project simulation folder (must be located where you
have your sim11-file). The simulation engine searches for the MIKE11.Ini
file locally first, and secondly for the global file in the bin-folder. By cop-
ying the file to the local folder, you ensure, that the changes you eventu-
ally make in the ini-file only concerns the present project, which is often
required as changes in the ini-file are possibly required and beneficial for
one project, whereas different values or parameters must be altered in
other projects.
ST-Variables no 1:
BOTTOM_LEVEL_UPDATE_METHODS = 4
ST-variable 13-17:
COUPLED_ST_AD=OFF
CORRECT_ST_AD=OFF
NO_AD_ST_FRACTIONS=1
AD_COMP_NO=1
ST_FRAC_NO=1
456 MIKE 11
FLOOD FORECASTING EDITOR
457
458 MIKE 11
Basic definitions
Real-time mode
Real time mode defines a condition where MIKE 11 FF is used to execute
simulations applying real-time hydrometeorological data as boundary
conditions. The common time span of the boundary data defines the hind-
.
Figure 10.2 Definition of Hindcast Period and ToF
Historical mode
While real-time telemetry data form the boundary conditions in an opera-
tional forecasting mode, historical hydrometeorological data are applied
as boundary conditions in the calibration and validation phase of forecast
modelling.
In the example shown below in Figure 10.3 the hindcast period starts on
the 4 January 1999 at 12:00 and last up to 7 January 1999 at 12:00.
460 MIKE 11
Forecast
10.2 Forecast
The main forecast parameters are specified in the Forecast Menu,
Figure 10.4.
10.2.3 Accuracy
The Boundary Conditions estimated after the Time of Forecast are obvi-
ously uncertain. The effect of a specified uncertainty level can be included
in the simulations.
Specify either global and/or local values for the deviation. Global values
are applied to all catchments or HD boundary conditions, except those
which are listed in the ‘Local Values’ fields.
Estimated boundary conditions with Upper and Lower levels are stored in
the ‘Boundary Estimates’ directory as described in Section 10.3.4.
462 MIKE 11
Forecast
Seasonal forecasting
Not yet implemented
Storage timestep
The storage frequency of forecast results can be more or less frequent than
the general MIKE 11 HD storage frequency specified in the Results menu
in the sim11 editor.
Boundary conditions applied during the forecast period are in this manual
described as Estimated boundary conditions.
464 MIKE 11
Boundary estimates
10.3.1 Setup
Specify catchment name (RR) or River name and Chainage (HD) to locate
the actual boundary
Type
Specify the appropriate data type:
Filename
Press the “...” button to select the appropriate time series file.
Filetype
The Axis type for the dfs0 files applied in the forecast period can be either
‘Calendar axis’ or ‘Relative axis’. If a dfs0 file is based on a ‘Relative
time axis’ the start time of that particular time series will be interpreted as
ToF.
10.3.2 Editing
All files included in the setup menu will be listed in the ‘Editing’ menu as
seen in Figure 10.9 above. Pressing the “Edit” button will start the MIKE
Zero time series editor with the actual time series loaded. In this manner it
is possible to view and edit the boundary estimate time series.
Table 10.1
466 MIKE 11
Boundary estimates
Table 10.1
468 MIKE 11
Update specifications
10.4.1 Comparison
Station
The location of the update point is defined via its River name and Chain-
age. If the specified chainage does not correspond to the computational
network it is shifted to the nearest h- or Q-point by the FF module and a
warning message is issued.
Data type
The Data type can be specified as water level or discharge. In general,
water level data should be specified at all sites where level forecasts are to
be issued, and discharge at reservoir inflow points.
470 MIKE 11
Update specifications
Method
Iterations
See No. of Iterations
Implicit solution
The specified time series are applied as internal boundary conditions in
the model. In the Continuity Equation hn+1 is substituted by the observed
water level and the lateral inflow qn+½ is calculated and applied as the
updating discharge.
No. of iterations
If a river branch includes a number of update points the specified No. of
iterations should be equal to or larger than this number. For large rivers
with few update points it may increase the update efficiency to use an
even larger number of iterations. Different numbers of iterations should be
tested before operational forecasting is initiated. A larger number will
increase the accuracy but also increase the required calculation time.
Frequency
Frequency of updating, i.e. the number of MIKE 11 HD time steps
between data observations in the time series used for updating.
10.4.2 Correction
The updating routine will calculate a correction discharge to be routed into
the river system along the correction branch. The correction branch is
specified by River name, First chainage and Last chainage.
10.4.3 Parameters
Table 10.2
472 MIKE 11
DATA ASSIMILATION EDITOR
473
474 MIKE 11
11 DATA ASSIMILATION EDITOR
The Data Assimilation editor is used for specifying the parameters needed
when carrying out
z Uncertainty estimation
z Model updating
z Forecasting
The uncertainty assessment and updating methods are applicable to hydro-
dynamic and/or advection dispersion simulations. Uncertainty assessment
and updating for NAM rainfall runoff will become available in a later
release. Thus any reference to catchment modelling in the present inter-
face should be neglected. One exception being that the uncertainty may be
applied to output from catchments.
Uncertainty assessment
The uncertainty assessment is a powerful tool for evaluating the effect of
uncertainties on the boundary conditions in a river network. To carry out
an uncertainty assessment a minimum of input parameters is required
from the user.
Model updating
If reliable measurements are available within the model domain, these
may be utilized to improve the model results of MIKE 11 using updating.
Two different updating methods are available:
Forecasting
With the appropriate data the updating routines of MIKE 11 can improve
model results prior to a time of forecast. After the time of forecast the
model is corrected using forecasts of errors identified prior to the time of
forecast. The type of error correction to be applied in the forecast period
depends on the option chosen for model updating. If the Kalman filter is
applied, estimated errors in the boundary conditions at the time of forecast
are phased out according to an exponential decay. If the combined weight-
ing function and error correction procedure is applied, errors at measure-
ment points are forecasted. Thus, the updating algorithm has an effect also
after the time of forecast. Furthermore, by using the Kalman filter proce-
dure a stochastic or ensemble forecast may be generated.
Depending on the method applied different menus within the editor are
required.
11.1 General
Figure 11.1 The general data assimilation model parameter tab. Note that
depending on the choice made in the Module selection section one
or more of the other boxes will be made inactive.
The first option requires the user to choose which updating technique that
should be applied. The choices are:
z Kalman filter: The model uses the ensemble Kalman filter based on
Monte Carlo simulation techniques to estimate the updates to be
applied.
476 MIKE 11
General
Ensemble size
If applying either the uncertainty assessment simulation method or the
Kalman filter updating method, the user is required to supply the ensemble
size. The ensemble size is defined as the number of simultaneous runs that
are to be carried out to evaluate the statistical properties needed for the
uncertainty assessment output and also for determining the updating
parameters applied in the Kalman filter.
The above considerations should be taken into account when choosing the
ensemble size. Recommended values are 50-200. If only reliable estimates
of standard deviations are of importance, a smaller value may be chosen.
When producing confidence intervals an estimate of the full uncertainty
distribution is needed, thus an increase in the ensemble size is recom-
mended for such cases.
11.1.3 Forecast
The forecasting parameters are only applicable to the updating methods.
The time of forecast should be within the simulation period specified in
the Simulation editor.
When using the Kalman filter, the user may select either a deterministic or
a stochastic forecast. The first uses the ensemble mean of the state varia-
bles including corrections at the time of forecast for a deterministic run,
whereas the latter continues the stochastic run omitting the updating step.
11.2 Measurements
Figure 11.2 The Measurements tab. Depending on the choice made in the gen-
eral tab some fields are made inactive.
The measurements page is used for specifying the location and parameters
of the measurements that are to be used for the updating. Note that the
page is only used when applying the updating option.
To add a new measurement place the cursor in the last row in the overview
and press the right tab until an additional row appears.
478 MIKE 11
Measurements
Branch Name
The branch name where the measurement is located is specified.
Chainage
The chainage of the measurement is specified.
Variable Type
Presently three types of measurements may be applied:
z Water level
z Discharge
z Concentration
If a concentration measurement is selected, the component number is
required. The component number is given by the order found in the Com-
ponents page of the AD editor.
File
The location of the measurement data is required. The dfs0 file may be
selected by browsing.
Note that for the percentage option the percentage is taken based on the
absolute value for discharge and concentration. For water level the per-
centage is interpreted as being with respect to the water depth.
For the weighting function method the standard deviation box is not appli-
cable.
Type
Three different types of weighting functions are available:
Amplitude
The Amplitude specifies the fraction of the observed error at the measure-
ment location that should be applied as error correction at that point. The
Amplitude should reflect the confidence of the observation as compared to
the model forecast. That is, if the amplitude equals one, the measurement
is assumed to be perfect, whereas for smaller amplitudes less emphasis is
put on the measurement as compared to the model forecast.
Lower chainage
The Lower chainage specifies the lower point (smallest chainage value) in
the branch where the measurement point is located for which the weight-
ing function is applied.
480 MIKE 11
Measurements
Note: Lower chainage must be smaller than the Upper chainage (regard-
less of branch flow direction definition).
Upper chainage
The Upper chainage specifies the upper point (highest chainage value) in
the branch where the measurement point is located for which the weight-
ing function is applied.
Note: Upper chainage must be larger than the Lower chainage (regardless
of branch flow direction definition).
Fade up
The user specifies the number of time steps used to fade up the corrections
initially, starting from zero corrections at the first time step up to full cor-
rections after the fade-up period. Employment of this function ensures that
model instabilities are avoided because the desired error correction value
is obtained gradually rather than abruptly.
The Error forecast model tick box is not applicable for the Kalman filter
method.
The Equation editor is used to define the error forecast models used in the
weighting function update procedure. Error forecast models can be
defined as general linear or non-linear functions with a one-step ahead
prediction that depends on the previous errors at the measurement loca-
tion, the model state at any grid points in the river network, and the bound-
ary conditions (or any other time series defined in a dfs0 file).
11.3.1 General
To add a new equation place the cursor in the last row in the Overview of
Error Forecast Model Equations and press the right tab until an additional
row appears.
Name
Unique name of the defined equation. This name is shown in the pull-
down menu of available equations when defining an error forecast model
for the Weighting function method on the Measurements page.
482 MIKE 11
Equation Editor
Equation
The equation editor is based on an equation parser that uses the general
arithmetic operators:
z Addition (+)
z Subtraction (-)
z Multiplication (*)
z Division (/)
In addition, a number of mathematical functions are supported, see
Table 11.1
Table 11.1 Mathematical functions used by the equation parser (X and Y are
variable names).
Syntax Function
SQR(X) Square function
SQRT(X) Square root function
SIN(X) Sine function. SIN returns the sine of the angle X in
radians.
COS(X) Cosine function. COS returns the cosine of the angle X
in radians.
TAN(X) Tangent function. TAN returns the cosine of the angle
X in radians.
COTAN(X) Cotangent function. COTAN returns the cosine of the
angle X in radians.
ATAN(X) ArcTangent function
EXP(X) Exponential function
LN(X) Natural logarithmic function
LOG(X) 10 based logarithmic function
SINH(X) Sinus Hyperbolic function
COSH(X) Cosine Hyperbolic function
INTPOW(X,Y) The INTPOW function raises X to an integer power Y,
e.g. INTPOW(2, 3) = 8. Note that the result of
INTPOW(2, 3.4) = 8 as well.
Table 11.1 Mathematical functions used by the equation parser (X and Y are
variable names).
Syntax Function
POW(X,Y) The POW function raises X to any power Y
ABS(X) Absolute value
SIGN(X) SIGN(X) returns -1 if X<0; +1 if X>0, 0 if X=0.
TRUNC(X) Discards the fractional part of a number, e.g.
TRUNC(3.2) is 3.
MIN(X,Y) Minimum of X and Y, e.g. MIN(2, 3) is 2.
MAX(X,Y) Maximum of X and Y, e.g. Max(2,3) is 3.
E(L) Error function with lag L (L<0), e.g. E(-1) is the error
at the previous time step.
5*X1-SQR(X2)+2*LN(X3)
If one wants to apply a second order auto regressive model, the Equation
field should be written:
A*E(-1)+B*E(-2)
Estimation period
For parameters defined as values in the equation, automatic parameter
estimation can be applied based on the record of observed errors. The
period of the record to be used for the parameter estimation can be speci-
fied relative to the time of forecast.
484 MIKE 11
Equation Editor
Name
Unique name of parameter specified in the equation.
Type
Three different parameter types are available:
Estimated
For a parameter that is assigned a numerical value, automatic parameter
estimation can be applied to estimate the value from the time series of
errors. Otherwise the value should be specified by the user in the Value
field.
Value
Numerical value assigned to the parameter. If the Estimated tick box is
active, this field is made inactive.
Minimum
Allowable lower bound on the parameter when automatic parameter esti-
mation is adopted.
Maximum
Allowable upper bound on the parameter when automatic parameter esti-
mation is adopted.
Variable type
If the parameter is defined as a state variable, the variable type should be
defined:
z Water level
z Discharge
z Concentration
Component number
If a concentration variable is selected, the component number is required.
The component number is given by the order found in the Components
page of the AD editor.
River name
The river name of the grid point of the state variable.
Chainage
The chainage of the grid point of the state variable.
File
If the parameter is defined as a time series, a dfs0 file item must be speci-
fied.
486 MIKE 11
Boundary Statistics
This menu is only applicable to the Kalman filter updating method and the
uncertainty assessment simulation modes. The structure of this page has
been chosen to facilitate that the same uncertainty statistic can be applied
to multiple boundaries. For instance, it may be that the user would like to
apply a standard deviation of say 10 percent to all discharge boundary
conditions in the set-up.
11.4.1 Details
The details for the boundaries where uncertainty is to be applied are spec-
ified here.
Boundary Type
The boundary types that uncertainty may be applied to are:
z Water levels
z Discharges
z Lateral source point inflows
z Q/h-relations
z Output from catchments
z Concentrations
z Wind fields
Note that at present uncertainty cannot be attributed to distributed lateral
discharge sources. Also note that the interface has been prepared for
including updates on rainfall runoff models, though the present release
does not include any rainfall runoff update capabilities.
Once the boundary type has been selected the locations where the uncer-
tainty is to be applied is entered in the second of the upper tables. Use the
tab key on the keyboard to supply additional locations. A location is spec-
ified through a branch name and a chainage. If a concentration is to be
given, the component number is also required. The component number is
given by the order found in the Components page of the AD editor. For
rainfall/runoff, only the catchment name needs to be specified. Wind
fields in MIKE 11 can be of two types either distributed or global. If a dis-
tributed wind field is to be augmented with uncertainty, then the chainage
where the distributed wind source starts is required. If uncertainty is to be
described for the global wind field, then the river name should be chosen
as ‘GLOBAL WIND’ without a chainage.If the user-supplied location
does not coincide with a boundary/forcing term, a warning is issued and
the boundary information is ignored in the subsequent simulation.
For wind fields and catchments two standard deviations must be supplied
since both of these types consist of two input variables. That is, wind
velocity and wind direction for a wind field, and runoff and net rainfall for
a catchment.
The standard deviation editor is used for defining the different standard
deviations that may be applied at the boundaries.
488 MIKE 11
Standard deviation editor
The errors applied at the boundaries may be described through a first order
auto regressive process:
ξ n + 1 = φξ n + ε (11.1)
where
ε white noise
∆tln 2
φ = exp – --------------- (11.2)
T½
where
The time constant TC should be interpreted as the time it takes for the cor-
rection to drop to half the initial magnitude (exponential decay). For wind
fields and Catchment runoff, two values must be supplied. The numerical
value of the regression coefficient must be less than unity to ensure that
the variance of the model error is limited. A negative time constant results
in a regression coefficient which is greater than unity. Therefore, if a neg-
ative time constant is entered, a warning is issued and the regression coef-
ficient is set to 0 in the simulation.
As an option, the time constant used when forecasting may differ from the
time constant used for updating. Thus, all in all up to four time constants
should be supplied by the user.
The model errors being described through the auto regressive process are
used in the forecast period to generate errors to be added to the bounda-
490 MIKE 11
Output
ries. If the forecast is deterministic, the model errors are described through
the auto regressive process without adding white noise. The model errors
fade out according to an exponential decay.
11.6 Output
The Output page is used for specifying additional output from an uncer-
tainty assessment analysis or a Kalman filter update. The output page is
split in two; one for hydrodynamic output and one for advection disper-
sion output. For both types the user may select the following output:
z Standard deviation
z Corrections (the corrections of the state variables in the river network
caused by updating with each of the measurements)
z Confidence intervals
The latter requires the user to specify a percentage for each confidence
interval. For example in Figure 11.6 as additional output two confidence
intervals (75% and 90%) along with the standard deviation have been
selected for the hydrodynamic variables consisting of water level and dis-
charge. Thus, in the additional HD result file these confidence intervals
and the standard deviation in the whole network will be stored. Further, in
Figure 11.6 two confidence intervals have been chosen for the advection
dispersion components. These confidence intervals are stored in the addi-
tional AD result file.
Figure 11.7 Selecting the Data assimilation mode in the simulation editor.
492 MIKE 11
A step by step guide to uncertainty assessment
7 Choose the ensemble size. Remember that the larger this number the
more accurate the uncertainty assessment, the trade-off being longer
run times. A value of 100 is recommended if feasible.
8 Consider boundary conditions applied to the model of the type
– Water level
– Discharge
– Q/H-relation
– Wind
– Source point lateral discharge
– Concentration
– Output from Catchment
For each of these boundaries consider which are to be assessed with
respect to uncertainty.
9 Based on the selection above the boundaries are added in the DA edi-
tor. Note that the same boundary statistics may be applied to multiple
boundaries of the same type. Thus, for instance a 10 percent uncer-
tainty can be added to all discharge boundaries by selecting all the
locations in the second box as shown below.
Figure 11.9 Multiple locations with the same boundary statistics. The same sta-
tistics will be applied to the upstream boundaries in main stem,
STtrib1 and STtrib2. To add a line place the cursor in the last line or
if no lines are present in the upper left hand corner and use the right
tab key on the keyboard.
Once the uncertainties have been decided upon the Standard Deviation
Editor is used to define the individual standard deviation items. The
standard deviation is given an ID. This ID is subsequently used for
10 Finally, the time constant for the temporal development of the errors
are set. Only the time constant prior to the time of forecast is used. The
time constants describe the “colouring” of the noise. If pure white
noise is to be applied at the boundaries, the time constant should be set
to zero. Typically, the noise at the boundaries will not be independent
from time step to time step (white noise), thus a non-zero time-constant
(coloured noise) is recommended. Please refer to section 11.5 for a
description of the time constants.
11 Decide what additional output is required. Tick the appropriate boxes
for HD and/or AD and supply the confidence intervals percentage.
Note that the confidence intervals are found by sorting each of the sim-
ulation values in the ensemble at every grid point. The confidence
internal of say 90% is found by locating the values in the lowest 5%
range of the ensemble and the values in the highest 5% range. Thus, as
an absolute minimum the ensemble size must be sufficiently large to
represent this fraction. For the 90% case an ensemble size of at least 20
(1/20 = 5%) is required. It is recommended to use a larger ensemble
size than the minimum dictated by the confidence interval.
494 MIKE 11
A step by step guide to updating using the Weighting function method
11.10 Examples
With the installation of MIKE 11 follows a number of examples illustrat-
ing the use of the Data assimilations methods.
The examples are all based on a simple set-up with three branches. The
examples are designed to illustrate the main features of the data assimila-
tion. The are 5 examples available all located in the folder:
496 MIKE 11
Examples
‘...\examples\m11\DataAssimilation\example1’
Figure 11.10 The discharge through the structure in Trib2. The upper and lower
limits of the 90% confidence interval are also illustrated.
Figure 11.11 Updating using the Kalman filtering technique. The discharge at the
upstream end of the main stem is illustrated. The blue curve is the
result obtained using MIKE 11 without updating the black curve with
crosses is the result obtained using the updating method. The
updating is based on a measurement located further downstream in
the main stem.
Figure 11.12 shows the temporal variation of the concentration at the junc-
tion of the main stem and Trib2 along with a 90% confidence interval
band.
498 MIKE 11
Examples
Figure 11.12 The concentration at the junction of the main stem and trib2. The
upper and lower limits of the 90% confidence interval are also illus-
trated.
Figure 11.13 The concentration at the junction of the main stem and trib2. The
red curve illustrates the original simulation without erroneous
boundary condition. The blue curve illustrates the updated results.
Finally for comparison the erroneous results (no update) are illus-
trated as a horizontal line at 0.5 g/m3.
500 MIKE 11
BATCH SIMULATION EDITOR
501
502 MIKE 11
12 BATCH SIMULATION EDITOR
The Batch Simulation Editor offers a possibility for setting up a batch sim-
ulation from the MIKEZero shell. That is, the Batch Simulation Editor is
used to pre-define a number of simulations, where all items included in a
simulation (input-files, simulation parameters, output files etc.) can be
changed from simulation to simulation and multiple simulations are then
performed automatically when starting the Batch simulation.
Figure 12.1 above shows an example on how to use the Batch Simulation
Editor. The figure presents a batch simulation setup, where two different
network files are combined with two different HD Parameter files. A setup
like this could be used to investigate the impact of variations in bed resist-
ance values (Manning numbers) at locations where a hydraulic structure
(weir) has been planned. The two different network files will then be iden-
tical except one of the files will contain description on the new proposed
weir, and the two HD Parameter files will only differ in the local variation
of the Manning numbers.
Output from the four different batch simulations has also been defined
such, that results from each simulation are saved in different result-files.
The Batch Simulation Editor has been developed in cooperation with CTI
Engineering, CO., Ltd., Japan.
Filename and path to the base Sim11 file must be defined in the ‘Base
Simulation File’ field (Use the ‘...’ button to browse for the Base Sim11
file on your computer).
Each line in the ‘Selected Parameters’ grid shall only contain specifica-
tions of the parameters or input files which should be different from the
definitions in the base simulation file. Parameters which should differ
from the base simulation file is selected in the tree-view on the left-part of
the Batch Simulation Editor, (see Figure 12.2). Open any tree-view items
by clicking the ‘+’ and select the item/parameter, which should be modi-
fied in the batch simulation by double-clicking the empty square in front
of the specific item. After double-clicking the item, a new column will be
introduced in the ‘Selected Parameters’ grid and a red check-mark is pre-
sented in the square in front of the selected item in the tree-view. Within
the new column in the Selected Parameters , it is now possible for the user
to select different input files or define variations in input parameters for
the different simulation performed in the batch simulation.
504 MIKE 11
Setting up a Batch Simulation
Figure 12.2 Tree view from the Batch Simulation Editor dialog for selecting batch
simulation parameters
If e.g. the Network file should be different in some simulations, open the
‘Input files’ item in the tree-view and double-click the Network square.
After this a Network column is presented in the ‘Selected Parameters’ grid
and network-files can now be specified in this column - either manually or
by pressing the ‘...’ button to browse for the required file. If e.g. the net-
work file in one simulation should be the same as in the base simulation
file - but other parameters are changed - the ‘base network file’ must be
defined in the network field, as it is not allowed to have any blank cells in
the ‘Selected Parameters’ grid.
After all files and parameters for the batch simulation have been specified,
it is required to save the data to a Batch Simulation file (*.BS11).
506 MIKE 11
CORRELATION ANALYSIS AND
GAP FILLING TOOL
507
508 MIKE 11
Introduction
13.1 Introduction
This tool can, as it’s name indicates, do the following:
The tool itself consists of one window being divided into a Tree view, a
Data view and an Error Message view. See Figure 13.1. When applying
the tool you would have to step through all pages by selecting the items in
the tree view from top to bottom. The content of the tree view may be less
than shown in Figure 13.1.
This tool has been developed in cooperation with CTI Engineering, CO.,
Ltd., Japan.
Correlation Analysis
The choice of including correlation analysis or not has to be made. If
selected an output file name has to be specified. The output file is an
ASCII/PFS file.
510 MIKE 11
CorrAnalysis page
Gap Filling
The choice of including gap filling or not has to be made. If selected a
regressive formula input and an gap filling output file name has to be
specified. Most often the user will specify the regressive formula input file
as being equal to the output file from the correlation analysis. The output
file is a .dfs0 time series file which will contain the result of the gap fill-
ing.
Figure 13.2 shows the page when all data selections and choices has been
made.
Regression Formula
On this page the user selects what regression formula to use. Four types of
regression formulas all being linear are available:
1 Y=AX+B (Y regression)
2 Y=AX+B (XY regression)
3 Y=AX (Y regression)
Lower Limit
For some data types, for instance rainfall, is may be desirable to exclude
values less than a lower limit in the analysis.
Run
Once the run button is pressed the correlation analysis in carried out, and
the result is written to the output file specified on the first page.
512 MIKE 11
RegEquationOutput page
This page contains a table with number of rows and columns equal to the
number of items in the time series input file. For each item in the file a
row is shown holding the result of the correlation analysis against all items
in the time series file. The result of each analysis is shown with two num-
bers: the correlation coefficient and number of values applied in the analy-
sis. The number of values applied may differ from one analysis to the next
as the number of missing values may differ from one item to the next. In
the diagonal of the table the correlation coefficient is always 1.000 as this
is the result of correlating the item to itself.
This page contains a table with number of rows and columns equal to the
number of items in the time series input file. For each item in the file a
row is shown holding the coefficients in the regression equation corre-
sponding to the other items in the time series file. The two numbers shown
are the coefficients A (upper) and B (lower) in the equation Y=AX+B. In
the diagonal of the table the A and B is not shown as these always are 1
and 0 respectively.
514 MIKE 11
Gap Filling page
Copy
Pressing the button copies the plot to the clipboard as Windows meta file
information. This allows for subsequent paste into for instance a word
processing tool.
Print
Pressing this button will send the plot to a printer.
The gap filling works such that for each missing value in each item an
estimated value will be calculated by using the value at the same time step
from another item transformed using the corresponding regressing equa-
tion. For each item the other items are ranked in the order of decreasing
correlation coefficient. When a missing value occurs a non-missing value
is searched for in the other items at the same time step, and the non-miss-
ing value in the item with the highest rank (lowest order) is selected.
516 MIKE 11
Detailed Result page
Detailed condition
In the table the item specific conditions for the gap filling is given. For
each item to be gap filled a row with specifications is available. The first
column in the row indicate if the specifications are to be achieved auto-
matically (equal to the output of the correlation analysis) or to be user
defined. In the remaining columns the items are shown in order with the
following information for each item:
Run
Pressing this button starts the gap filling process, and the output .dfs0 file
is generated. The result is shown also in two pages Details Results and
Outlined Results. See Figure 13.8.
The detailed results shows all the time series data after gap filling. Artifi-
cial values that has been calculated by the gap filling has an integer
number to the left of the number. This number is the item number used for
the gap filling. There may still be gaps in the data as no gap filling is done
if all items on a particular time step has missing values.
518 MIKE 11
Outlined Output page
520 MIKE 11
FLOW RESISTANCE AND
VEGETATION A
521
522 MIKE 11
Flow Channels in Halkaer AA
Results and findings from the Danish gauging programmes and investiga-
tions on the weeds influence on flow resistance are described in the fol-
lowing.
Fig A.1.1 Q-h curves determined for varying flow channel width
524 MIKE 11
Laboratory measurements using Bur Reed
n = aln(VR) + b (A.1.1)
β
n = aD η (A.1.2)
A
D η = --- (A.1.3)
B
where A is the flow area and B is the width of the section at water surface.
526 MIKE 11
Experiments in ‘Arn Å’
A.1.5 References
/1/ Bakry, M.F.; T.K.Gates; A.F.Khattab:
528 MIKE 11
ADDITIONAL TOOLS B
529
530 MIKE 11
Merging .pfs files
where
...\ denotes the full path to the application located in the bin director of
the installation.
Note that the above syntax is based on a call from the data directory
(the directory where the pfsfiles are located).
Note: When converting v.3.2 network-files (.RDF) all relevant cross sec-
tion files (.pst, .ix0, .ix1) must be located in the same directory as the
.RDF file.
To execute the program, open a Dos Prompt (on Windows XP the Dos
Prompt can normally be found under; Start-> Programs -> Accessories ->
Command Prompt).
532 MIKE 11
Converting simulation results to text files (Res11Read)
where :
‘...\’ denotes the full path to the application located in the bin directory
of the MIKE Zero installation.
The different Options available for usage with the application is listed and
described briefly below:
z ‘-xy’
Exports X- and Y-coordinates as well as levels for all calculation points
(h- and Q-points) in the model.
(example: ‘Res11Read -xy xx.res11 out.txt’)
z ‘-xyh’
Exports X- and Y-coordinates and levels for all the h-points in the
model.
(example: ‘Res11Read -xyh xx.res11 out.txt’)
z ‘-xyq’
Exports X- and Y-coordinates and levels for all the Q-points in the
model.
(example: ‘Res11Read -xyq xx.res11 out.txt’)
z ‘-xyxsec’
Exports X- and Y-coordinates and levels for all h-points defined by
cross sections in the model.
(example: ‘Res11Read -xyxsec xx.res11 out.txt’)
z ‘-raw’
Exports Raw data for all cross sections used in the simulation.
(example: ‘Res11Read -raw xx.res11 out.txt’)
z ‘-sim’
Content of the .sim11 file used for the simulation is saved as static data
in the Res11-file. This can be extracted using this option (useful for
extracting details on simulation settings).
(example: ‘Res11Read -sim xx.res11 out.txt’)
z ‘-minX’
Exports the absolute Minimum values obtained in all grid points
throughout the simulation for item no X.
The item-number can be found by using the option ‘-Items’ initially.
(For a standard HD-Result file Water Level is Item number 1 and Dis-
charge is Item-number 2).
(example: ‘Res11Read -min2 xx.res11 out.txt’)
z ‘-maxX’
Exports the absolute Maximum values obtained in all grid points
throughout the simulation for item no X
(example: ‘Res11Read -max1 xx.res11 out.txt’)
z ‘-xsecids’
Exsports the Cross section ID’s from cross sections active in the simu-
lation.
(example: ‘Res11Read -xsecids xx.res11 out.txt’)
z ‘-usermarks’
Exports User defined marks from the result file.
(example: ‘Res11Read -usermarks xx.res11 out.txt’)
z ‘-items’
Exports a List of all dynamic items present in the result file.
(example: ‘Res11Read -items xx.res11 out.txt’)
z ‘-allres‘
Exports all simulation results (all timesteps) for all calculation points
in the model.
Note: This option can occasionally result in a large ascii output file!
(example: ‘Res11Read -allres xx.res11 out.txt’)
z ‘-someresFILE’
Simulation results from selected calculation points only are written to
the output file. The selected points to be exported is defined in a con-
figuration file with a user-defined filename. In this example the file is
named; ‘FILE’ (see format of FILE below).
(example: ‘Res11Read -SomeresInp.txt xx.res11 out.txt’
where’inp.txt’ is the input configuration filename)
534 MIKE 11
Converting simulation results to text files (Res11Read)
z ‘-compareFILE’
This option allows a comparion of two Res11-files from the same
model, and the option will issue an error if differences in result-values
between simulation points in the compared results are larger than a
user-defined threshold value.
The FILE is the configuration file for this option and contains informa-
tion on which threshold value shall be active for which result items.
(see format of FILE below).
(example: ‘Res11Read -CompareInpFile.txt xx.res11 yy.res11 out.txt’
where ‘InpFile.txt is the configuration filename)
z ‘-silent’
Writing to prompt is cancelled. Used in conjunction with one or more
of the other options.
(example: ‘Res11Read -allres -silent xx.res11 out.txt’)
z ‘-MessageCompare’
Writes additional message on Return code in output file on result of
comparison ( 0 or 2 for the Compare results option).
(example: ‘Res11Read -CompareInp.txt -MessageCompare xx.res11
yy.res11 out.txt’
where ‘Inp.txt’ is the configuration file for file comparison)
z ‘-DHIASCII’
Option for suppressing header information in Outputfile in order to
convert e.g. results/Timeseries into a DHI standard ascii format.
Should be used in conjunction with one or more of the above.
(example: ‘Res11Read -max1 -dhiascii xx.res11 out.txt’)
z ‘-FloodWatch‘
Special option for converting results from a MIKE 11 simulation into a
specific format used by the FLOOD Watch real-time data handling sys-
tem from DHI.
z ‘-MakeDfs0’
Option to automatically create a timeseries containing results from spe-
cific calculation points.
Note: This option can only be used in conjunction with the -Someres-
FILE option
(example: ‘Res11Read -SomeResFILE -MakeDfs0 xx.res11 out.txt’)
z ‘-SeperatorSTRING’
Option for defining the seperator string between output-results.
Note: This option is only applicable in conjunction with the flood-
Watch option
(example: ‘Res11Read -Floodwatch -Separator, xx.res11 out.txt’)
z ‘-RRResults’
Extraction of Results from Rainfall Runoff simulation specifically ded-
icated toward a FLOOD Watch system. Hence, this option is only
applicable with the FloodWatch option.
(example: ‘Res11Read -Floodwatch -RRResults rr.res11 out.txt’)
z ‘-ArcView’
Exports results into an ArcView table import format. Used in conjunc-
tion with other options. Typically RRResults as a sub-option to Flood-
watch.
(example: ‘Res11Read -Floodwatch -RRResults -ArcView rr.res11
out.txt’)
Some of the above listed options includes a ‘FILE’ in the calling conven-
sion. FILE is here a name of a configuration file used by the Res11Read
application to detect details on what results shall be converted.
When FILE is used in connection with the option ‘someresFILE’, the for-
mat of the configuration file is as presented below:
Figure 13.10 Format for use in the file used for the someresFILE option.
Note, that if results are extracted from a Rainfall Runoff simulation results
file, the format of the FILE is slightly different.
Instead of Chainage and Rivername you must instead specify the Catch-
ment-Area and Catchment-Name as parameters as presented below:
536 MIKE 11
Converting simulation results to text files (Res11Read)
Figure 13.11 Format of FILE when using option SomeRes and extacting results
from a Rainfall-Runoff result-file
This can be done using the ">" character at the end of a line in FILE:
Figure 13.12 Format of file used for the someresFILE option, with alternate output
file option.
Figure 13.13
The ‘Limit’ parameter defined in the FILE is the threshold for the compar-
ison between files, such that comparison errors will not be issued if differ-
ences are smaller than the define Limit-value.
...\M11QHCalcLauncher.exe simulationfilename.sim11
where
...\ denotes the full path to the application located in the bin directory
of the MIKE Zero installation.
538 MIKE 11
INDEX
539
Index
540 MIKE 11
Index
541
Index
542 MIKE 11