# 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)