# Generated on 08/23/2016 17:24:00
# Deltares, FM-Suite DFlowFM Model Version 1.1.5.35222, DFlow FM Version 1.1.171.46655
[model]
Program = D-Flow FM
Version = 1.1.171.46655
MDUFormatVersion = 1.02
GuiVersion = 1.1.5.35222
AutoStart = 0
[geometry]
NetFile = river1D_net.nc # Unstructured grid file *_net.nc
BathymetryFile = # Bathymetry points file *.xyb
DryPointsFile = # Dry points file *.xyz (third column dummy z values), or dry areas polygon file *.pol (third column 1/-1: inside/outside)
WaterLevIniFile = # Initial water levels sample file *.xyz
LandBoundaryFile = # Land boundaries file *.ldb, used for visualization
ThinDamFile = # Polyline file *_thd.pli, containing thin dams
FixedWeirFile = # Polyline file *_fxw.pliz, containing fixed weirs with rows x, y, crest level, left ground level, right ground level
StructureFile = # File *.ini containing list of structures (pumps, weirs, gates and general structures)
VertplizFile = # Vertical layering file *_vlay.pliz with rows x, y, Z, first Z, nr of layers, second Z, layer type
ProflocFile = # Channel profile location file *_proflocation.xyz with rows x, y, z, profile number ref
ProfdefFile = # Channel profile definition file *_profdefinition.def with definition for all profile numbers
ProfdefxyzFile = # Channel profile definition file _profdefinition.def with definition for all profile numbers
Uniformwidth1D = 2 # Uniform width for channel profiles not specified by profloc
ManholeFile = # File *.ini containing manholes
WaterLevIni = 8 # Initial water level at missing s0 values
Bedlevuni = -5 # Uniform bed level used at missing z values if BedlevType > 2
Bedslope = 0 # Bed slope inclination if BedlevType > 2
BedlevType = 3 # Bathymetry specification (1: at cell centers (from BathymetryFile), 2: at faces (from BathymetryFile), 3: at nodes, face levels mean of node values, 4: at nodes, face levels min. of node values, 5: at nodes, face levels max. of node values, 6: at nodes, face levels max. of cell center values)
Blmeanbelow = -999 # If not -999d0, below this level the cell center bed level is the mean of surrouding net nodes
Blminabove = -999 # If not -999d0, above this level the cell center bed level is the min. of surrouding net nodes
PartitionFile = # Domain partition polygon file *_part.pol for parallel run
AngLat = 0 # Angle of latitude S-N (deg), 0: no Coriolis
AngLon = 0 # Angle of longitude E-W (deg), 0: Greenwich
Conveyance2D = -1 # -1: R,HU,0: R,H, 1: R,A/P, 2: K,analytic-1D conv, 3: K,analytic-2D conv
Nonlin2D = 0 # Non-linear 2D volumes, only used if ibedlevtype=3 and Conveyance2D>=1
Sillheightmin = 0.5 # Weir treatment only if both sills larger than this value (m)
Makeorthocenters = 0 # Switch from circumcentres to orthocentres in geominit (1: yes, 0: no)
Dcenterinside = 1 # Limit cell center (1.0: in cell, 0.0: on c/g)
Bamin = 1E-06 # Minimum grid cell area, in combination with cut cells
OpenBoundaryTolerance= 3 # Search tolerance factor between boundary polyline and grid cells, in cell size units
RenumberFlowNodes = 0 # Renumber the flow nodes (1: yes, 0: no)
Kmx = 0 # Maximum number of vertical layers
Layertype = 1 # Vertical layer type (1: all sigma, 2: all z, 3: use VertplizFile)
Numtopsig = 0 # Number of sigma layers in top of z-layer model
SigmaGrowthFactor = 1 # Layer thickness growth factor from bed up
[numerics]
CFLMax = 0.7 # Maximum Courant number
AdvecType = 33 # Advection type (0: none, 1: Wenneker, 2: Wenneker q(uio-u), 3: Perot q(uio-u), 4: Perot q(ui-u), 5: Perot q(ui-u) without itself
TimeStepType = 2 # Time step handling (0: only transport, 1: transport + velocity update, 2: full implicit step-reduce, 3: step-Jacobi, 4: explicit)
Limtyphu = 0 # Limiter type for waterdepth in continuity eqn. (0: none, 1: minmod, 2: van Leer, 3: Kooren, 4: monotone central)
Limtypmom = 4 # Limiter type for cell center advection velocity (0: none, 1: minmod, 2: van Leer, 3: Kooren, 4: monotone central)
Limtypsa = 4 # Limiter type for salinity transport (0: none, 1: minmod, 2: van Leer, 3: Kooren, 4: monotone central)
TransportMethod = 1 # Transport method (0: Herman's method, 1: Sander's method)
Vertadvtypsal = 5 # Vertical advection type for salinity (0: none, 1: upwind explicit, 2: central explicit, 3: upwind implicit, 4: central implicit, 5: central implicit but upwind for neg. stratif., 6: higher order explicit, no Forester)
Icgsolver = 4 # Solver type (1: sobekGS_OMP, 2: sobekGS_OMPthreadsafe, 3: sobekGS, 4: sobekGS + Saadilud, 5: parallel/global Saad, 6: parallel/Petsc, 7: parallel/GS)
Maxdegree = 6 # Maximum degree in Gauss elimination
FixedWeirScheme = 0 # Fixed weir scheme (0: none, 1: compact stencil, 2: whole tile lifted, full subgrid weir + factor)
FixedWeirContraction= 1 # Fixed weir flow width contraction factor
FixedWeirfrictscheme= 1 # Fixed weir friction scheme (0: friction based on hu, 1: friction based on subgrid weir friction scheme)
Fixedweirtopwidth = 3 # Uniform width of the groyne part of fixed weirs
Fixedweirtopfrictcoef= -999 # Uniform friction coefficient of the groyne part of fixed weirs
Fixedweirtalud = 0.25 # Uniform talud slope of fixed weirs
Izbndpos = 0 # Position of z boundary (0: D3Dflow, 1: on net boundary, 2: on specifiend polyline)
Tlfsmo = 3600 # Fourier smoothing time (s) on water level boundaries
Slopedrop2D = 0 # Apply drop losses only if local bed slope > Slopedrop2D, (<=0: no drop losses)
Chkadvd = 0.1 # Check advection terms if depth < chkadvdp, => less setbacks
Teta0 = 0.55 # Theta of time integration (0.5 < theta < 1)
Qhrelax = 0.01 # Relaxation on Q-h open boundaries
Jbasqbnddownwindhs= 0 # Water depth scheme at discharge boundaries (0: original hu, 1: downwind hs)
cstbnd = 0 # Delft-3D type velocity treatment near boundaries for small coastal models (1: yes, 0: no)
Maxitverticalforestersal= 100 # Forester iterations for salinity (0: no vertical filter for salinity, > 0: max nr of iterations)
Maxitverticalforestertem= 0 # Forester iterations for temperature (0: no vertical filter for temperature, > 0: max nr of iterations)
Jaorgsethu = 1 # Velocity reconstruction scheme (0 : setumod, sethu, setau sequence, 1 : sethu, setau, setumod sequence (standard))
Turbulencemodel = 1 # Turbulence model (0: none, 1: constant, 2: algebraic, 3: k-epsilon, 4: k-tau)
Turbulenceadvection= 3 #
AntiCreep = 0 # Include anti-creep calculation (0: no, 1: yes)
Maxwaterleveldiff = 0 # Upper bound (in m) on water level changes (<=0: no bounds)
Maxvelocitydiff = 0 # Upper bound (in m/s) on velocity changes (<=0: no bounds)
Epshu = 0.0001 # Threshold water depth for wet and dry cells
SobekDFM_umin = 0
[physics]
UnifFrictCoef = 0.255 # Uniform friction coefficient (0: no friction)
UnifFrictType = 2 # Uniform friction type (0: Chezy, 1: Manning, 2: White-Colebrook, 3: idem, WAQUA style)
UnifFrictCoef1D = 0.023 # Uniform friction coefficient in 1D links (0: no friction)
UnifFrictCoefLin = 0 # Uniform linear friction coefficient for ocean models (m/s) (0: no friction)
Umodlin = 0 # Linear friction umod, for ifrctyp=4,5,6
Vicouv = 1 # Uniform horizontal eddy viscosity (m2/s)
Dicouv = 1 # Uniform horizontal eddy diffusivity (m2/s)
Vicoww = 5E-05 # Uniform vertical eddy viscosity (m2/s)
Dicoww = 5E-05 # Uniform vertical eddy diffusivity (m2/s)
Vicwminb = 0 # Minimum viscosity in prod and buoyancy term (m2/s)
Smagorinsky = 0 # Smagorinsky factor in horizontal turbulence
Elder = 0 # Elder factor in horizontal turbulence
Irov = 0
wall_ks = 0 # Nikuradse roughness for side walls
Rhomean = 1000 # Average water density (kg/m3)
Idensform = 1 # Density calulation (0: uniform, 1: Eckard, 2: Unesco, 3: baroclinic case)
Ag = 9.81 # Gravitational acceleration
TidalForcing = 0 # Tidal forcing, if jsferic = 1 (0: no, 1: yes)
Doodsonstart = 55.565 # TRIWAQ: 55.565, D3D: 57.555
Doodsonstop = 375.575 # TRIWAQ: 375.575, D3D: 275.555
Doodsoneps = 0 # TRIWAQ = 0.0 400 cmps , D3D = 0.03 60 cmps
Salinity = 0 # Include salinity, (0: no, 1: yes)
InitialSalinity = 0 # Uniform initial salinity concentration (ppt)
Sal0abovezlev = -999 # Vertical level (m) above which salinity is set 0
DeltaSalinity = -999 # for testcases
Backgroundsalinity= 30 # Background salinity for eqn. of state (ppt)
Temperature = 0 # Include temperature (0: no, 1: only transport, 3: excess model of D3D, 5: composite (ocean) model)
InitialTemperature= 6 # Uniform initial water temperature (degC)
Secchidepth = 1 # Water clarity parameter (m)
Stanton = -1 # Coefficient for convective heat flux
Dalton = -1 # Coefficient for evaporative heat flux
Backgroundwatertemperature= 6 # Background water temperature for eqn. of state (deg C)
SecondaryFlow = 0 # Secondary flow (0: no, 1: yes)
EffectSpiral = 0 # Weight factor of the spiral flow intensity on transport angle
BetaSpiral = 0 # Weight factor of the spiral flow intensity on flow dispersion stresses
[wind]
ICdtyp = 2 # Wind drag coefficient type (1: constant, 2: S&B 2 breakpoints, 3: S&B 3 breakpoints, 4: Charnock constant, 5: Whang)
Cdbreakpoints = 0.00063 0.00723 # Wind drag coefficient break points
Windspeedbreakpoints= 0 100 # Wind speed break points (m/s)
Rhoair = 1.2 # Air density (kg/m3)
PavBnd = 0 # Average air pressure on open boundaries (N/m2) (only applied if > 0)
PavIni = 0 # Average air pressure for initial water level correction (N/m2) (only applied if > 0)
[waves]
Wavemodelnr = 0 # Wave model nr. (0: none, 1: fetch/depth limited hurdlestive, 2: Young-Verhagen, 3: SWAN, 4: wave group forcing)
WaveNikuradse = 0.01 # Wave friction Nikuradse ks coefficient (m), used in Krone-Swart
Rouwav = FR84 # Friction model for wave induced shear stress
Gammax = 1 # Maximum wave height/water depth ratio
[time]
RefDate = 20120101 # Reference date (yyyymmdd)
Tzone = 0 # Time zone assigned to input time series
DtUser = 600 # Time interval (s) for external forcing update
DtNodal = # Time interval (s) for updating nodal factors in astronomical boundary conditions
DtMax = 600 # Maximal computation timestep (s)
DtInit = 60 # Initial computation timestep (s)
Tunit = S # Time unit for start/stop times (H, M or S)
TStart = 0 # Start time w.r.t. RefDate (in TUnit)
TStop = 864000 # Stop time w.r.t. RefDate (in TUnit)
[restart]
RestartFile = # Restart netcdf-file, either *_rst.nc or *_map.nc
RestartDateTime = # Restart date and time (YYYYMMDDHHMMSS) when restarting from *_map.nc
[external forcing]
ExtForceFile = # Old format for external forcings file *.ext, link with tim/cmp-format boundary conditions specification
ExtForceFileNew = river1D.ext # New format for external forcings file *.ext, link with bc-format boundary conditions specification
[trachytopes]
TrtRou = Y # Include alluvial and vegetation roughness (trachytopes) (Y: yes, N: no)
TrtDef = ruw.ttd # File (*.ttd) containing trachytope definitions
TrtL = river1D.arl # File (*.arl) containing distribution of trachytope definitions
TrtDt = 60 # Trachytope roughness update time interval (s)
[output]
OutputDir = output # Output directory of map-, his-, rst-, dat- and timings-files, default: DFM_OUTPUT_. Set to . for current dir.
FlowGeomFile = # Flow geometry NetCDF *_flowgeom.nc
ObsFile = river1D.xyn # Points file *.xyn with observation stations with rows x, y, station name
CrsFile = river1D.crs # Polyline file *_crs.pli defining observation cross sections
HisFile = # History output file *_his.nc in NetCDF format
HisInterval = 600 259200 # History output times, given as 'interval' 'start period' 'end period' (s)
XLSInterval = # Interval (s) between XLS history
MapFile = # Map file output file *_map.nc in NetCDF format
MapInterval = 3600 172800 864000 # Map file output times, given as 'interval' 'start period' 'end period' (s)
RstInterval = 0 # Restart file output times, given as 'interval' 'start period' 'end period' (s)
S1incinterval = # Interval (m) in incremental file for water levels s1
WaqInterval = # Interval (in s) between DELWAQ file outputs
StatsInterval = 3600 # Interval (in s) between simulation statistics output
Writebalancefile = 0 # Write balance file (1: yes, 0: no)
TimingsInterval = # Timings statistics output interval
MapFormat = 1 # Map file format (1: NetCDF, 2: Tecplot, 3: NetCFD and Tecplot)
Wrihis_balance = 1
Wrihis_structure_gen= 0
Wrihis_structure_dam= 0
Wrihis_structure_pump= 0
Wrihis_structure_gate= 0
Wrimap_waterlevel_s0= 0
Wrimap_waterlevel_s1= 1
Wrimap_velocity_component_u0= 1
Wrimap_velocity_component_u1= 1
Wrimap_velocity_vector= 1
Wrimap_upward_velocity_component= 0
Wrimap_density_rho= 0
Wrimap_horizontal_viscosity_viu= 0
Wrimap_horizontal_diffusivity_diu= 0
Wrimap_flow_flux_q1= 1
Wrimap_spiral_flow= 0
Wrimap_numlimdt = 1
Wrimap_taucurrent = 0
Wrimap_chezy = 1
Wrimap_turbulence = 0
Wrimap_wind = 0
Wrimap_heat_fluxes= 0
Richardsononoutput= 0 # Write Richardson numbers (1: yes, 0: no)
MapOutputTimeVector= # File (*.mpt) containing fixed map output times (s) w.r.t. RefDate
FullGridOutput = 0 # Full grid output mode (0: compact, 1: full time-varying grid data)
EulerVelocities = 0 # Euler velocities output (0: GLM, 1: Euler velocities)