[model] Program = D-Flow FM # Program name Version = # FM kernel Version MDUFormatVersion = 1.02 # File format version (do not edit this) GuiVersion = # DeltaShell FM suite version AutoStart = 0 # Autostart simulation after loading MDU (0: no, 1: autostart, 2: autostartstop) [geometry] NetFile = ameland_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 = 0 # 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 = 52 # 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 = 600 # 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 = 1 # 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 = # 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 [physics] UnifFrictCoef = 65 # Uniform friction coefficient (0: no friction) UnifFrictType = 0 # 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 roughness type (0: free slip, 1: partial slip using wall_ks) 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 = 3 # 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 0.00723 # Wind drag coefficient break points Windspeedbreakpoints = 0 100 100 # Wind speed break points (m/s) Rhoair = 1.2 # Air density (kg/m3) PavBnd = -999 # 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) Gapres = 101325.0000000 # Global Atmospheric Pressure, (N/m2) [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 = 0.5 # Maximum wave height/water depth ratio [time] RefDate = 20120401 # Reference date (yyyymmdd) Tzone = 0 # Time zone assigned to input time series DtUser = 60 # Time interval (s) for external forcing update DtNodal = 0 # Time interval (s) for updating nodal factors in astronomical boundary conditions DtMax = 60 # Maximal computation timestep (s) DtInit = 60 # Initial computation timestep (s) Tunit = S # Time unit for start/stop times (H, M or S) TStart = 950400 # Start time w.r.t. RefDate (in TUnit) TStop = 1382400 # 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 = wad.ext # Old format for external forcings file *.ext, link with tim/cmp-format boundary conditions specification ExtForceFileNew = # New format for external forcings file *.ext, link with bc-format boundary conditions specification [trachytopes] TrtRou = N # Include alluvial and vegetation roughness (trachytopes) (Y: yes, N: no) TrtDef = # File (*.ttd) containing trachytope definitions TrtL = # File (*.arl) containing distribution of trachytope definitions TrtDt = 60 # Trachytope roughness update time interval (s) [output] OutputDir = # 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 = wad.xyn # Points file *.xyn with observation stations with rows x, y, station name CrsFile = # Polyline file *_crs.pli defining observation cross sections HisFile = # History output file *_his.nc in NetCDF format HisInterval = 600 # History output times, given as 'interval' 'start period' 'end period' (s) XLSInterval = 0 # Interval (s) between XLS history MapFile = # Map file output file *_map.nc in NetCDF format MapInterval = 1800 # Map file output times, given as 'interval' 'start period' 'end period' (s) RstInterval = 86400 S1incinterval = 0 # Interval (m) in incremental file for water levels s1 WaqInterval = 0 # Interval (in s) between DELWAQ file outputs StatsInterval = 0 # Interval (in s) between simulation statistics output Writebalancefile = 0 # Write balance file (1: yes, 0: no) TimingsInterval = 0 # Timings statistics output interval MapFormat = 1 # Map file format (1: NetCDF, 2: Tecplot, 3: NetCFD and Tecplot) Wrihis_balance = 1 Wrihis_structure_gen = 1 Wrihis_structure_dam = 1 Wrihis_structure_pump = 1 Wrihis_structure_gate = 1 Wrimap_waterlevel_s0 = 1 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 = 1 Wrimap_horizontal_viscosity_viu = 1 Wrimap_horizontal_diffusivity_diu= 1 Wrimap_flow_flux_q1 = 1 Wrimap_spiral_flow = 1 Wrimap_numlimdt = 1 Wrimap_taucurrent = 1 Wrimap_chezy = 1 Wrimap_turbulence = 1 Wrimap_wind = 1 Wrimap_heat_fluxes = 0 Richardsononoutput = 1 MapOutputTimeVector = # File (.mpt) containing fixed map output times (s) w.r.t. RefDate FullGridOutput = 0 # 0:compact, 1:full time-varying grid data EulerVelocities = 0 # Euler velocities output (0: GLM, 1: Euler velocities)