[Section] Physical processes swave [-] Include short waves (1), exclude short waves (0) lwave [-] Include short wave forcing on NLSW equations and boundary conditions (1), or exclude (0) flow [-] Include flow calculation (1), or exclude (0) sedtrans [-] Include sediment transport (1) or exclude (0) morphology [-] Include morphology (1) or exclude (0) avalanching [-] Include avalanching (1) or exclude (0) nonh [-] (advanced) Non-hydrostatic pressure option: 0 = NSWE, 1 = NSW + non-hydrostatic pressure compensation Stelling & Zijlema, 2003 gwflow [-] (advanced) Turn on (1) or off (0) groundwater flow module q3d [-] (advanced) Turn on (1) or off (0) quasi-3D sediment transport module swrunup [-] (advanced) Turn on (1) or off (0) short wave runup ships [-] (advanced) Turn on (1) or off (0) ship waves vegetation [-] (advanced) Turn on (1) or off (0) interaction of waves and flow with vegetation snells [-] (advanced) Turn on (1) or off (0) Snell's law for wave refraction bchwiz [-] Use beachwizard, 0 = beachwizard off, 1 beachwizard on, bed update off, 2 beachwiz on, bed update on. (also requires morphology == 1) setbathy [-] Provide timeseries of prescribed bathy input, 0 = off (default), 1 = on [Section] Grid parameters depfile [-] Name of the input bathymetry file posdwn [-] Bathymetry is specified positive down (1) or positive up (-1) nx [-] Number of computiation cell corners in x-direction ny [-] Number of computiation cell corners in y-direction alfa [deg] Angle of x-axis from East vardx [-] Switch for variable grid spacing: 1 = irregular spacing, 0 = regular grid spacing dx [m] Regular grid spacing in x-direction dy [m] Regular grid spacing in y-direction xfile [name] Name of the file containing x-coordinates of the calculation grid yfile [name] Name of the file containing y-coordinates of the calculation grid xori [m] X-coordinate of origin of axis yori [m] Y-coordinate of origin of axis thetamin [deg] Lower directional limit (angle w.r.t computational x-axis) thetamax [deg] Higher directional limit (angle w.r.t computational x-axis) dtheta [deg] Directional resolution thetanaut [-] Thetamin,thetamax in cartesian (0) or nautical convention (1) gridform [name] Swicth to read in grid bathy files with 'XBeach' or 'Delft3D' format respectively xyfile [name] Name of the file containing (Delft3D) xy-coordinates of the calculation grid [Section] Model time tstop [s] Stop time of simulation, in morphological time CFL [-] Maximum Courant-Friedrichs-Lewy number tunits [-] (advanced) Units can be defined in udunits format (seconds since 1970-01-01 00:00:00.00 +1:00) [Section] Physical constants g [ms^-2] Gravitational acceleration rho [kgm^-3] Density of water depthscale [-] (advanced) depthscale of (lab)test simulated. 1 = default, which corresponds to teh real world (nature) the following (numerical) parameters are scaled with the depth scale (see Brandenburg, 2010): eps = eps_default/depthscale hmin = hmin_default/depthscale hswitch = hswitch/depthscale dzmax = dzmax/depthscale**1.5d0 Brandenburg concluded that also the following parameters potentially need to be scaled: wetslp, turb (suggested to turn off at depthscales<20) & ucr (distinguish ucr_bed load and ucr_sus at depthscales<20) [Section] Initial conditions zs0 [m] Inital water level zsinitfile [name] Name of inital condition file zs hotstartflow [-] (advanced) Switch to hotstart flow conditions with pressure gradient balanced by wind and bed stress [Section] Wave boundary condition parameters instat [-] Wave boundary condtion type taper [s] Spin-up time of wave boundary conditions, in morphological time Hrms [m] Hrms wave height for instat = 0,1,2,3 Tm01 [s] (deprecated) Old name for Trep Trep [s] Representative wave period for instat = 0,1,2,3 Tlong [s] Wave group period for case instat = 1 dir0 [deg] Mean wave direction (Nautical convention) for instat = 0,1,2,3 nmax [-] (advanced) maximum ratio of cg/c fro computing long wave boundary conditions m [-] Power in cos^m directional distribution for instat = 0,1,2,3 lateralwave [-] (deprecated) Switch for lateral boundary at left, 'neumann' = E Neumann, 'wavefront' = along wave front leftwave [-] old name for lateralwave rightwave [-] old name for lateralwave [Section] Wave-spectrum boundary condition parameters bcfile [-] Name of spectrum file random [-] (advanced) Random seed on (1) or off (0) for instat = 4,5,6 boundary conditions fcutoff [Hz] (advanced) Low-freq cutoff frequency for instat = 4,5,6 boundary conditions nspr [-] (advanced) nspr = 1 long wave direction forced into centres of short wave bins, nspr = 0 regular long wave spreadin trepfac [-] (advanced) Compute mean wave period over energy band: par%trepfac*maxval(Sf) for instat 4,5,6; converges to Tm01 for trepfac = 0.0 and sprdthr [-] (advanced) Threshold ratio to maxval of S above which spec dens are read in (default 0.08*maxval) oldwbc [-] (deprecated) (1) Use old version wave boundary conditions for instat 4,5,6 correctHm0 [-] (advanced) Turn off or on Hm0 correction oldnyq [-] (advanced) Turn off or on old nyquist switch Tm01switch [-] (advanced) Turn off or on Tm01 or Tm-10 switch rt [s] Duration of wave spectrum at offshore boundary, in morphological time dtbc [s] (advanced) Timestep used to describe time series of wave energy and long wave flux at offshore boundary (not affected by morfac) dthetaS_XB [deg] (advanced) The (counter-clockwise) angle in the degrees needed to rotate from the x-axis in SWAN to the x-axis pointing East nspectrumloc [-] (advanced) Number of input spectrum locations wbcversion [-] (advanced) Version of wave boundary conditions nonhspectrum [-] (advanced) Switch between spectrum format for wave action balance of nonhydrostatic waves [Section] Flow boundary condition parameters front [-] Switch for seaward flow boundary: 0 = radiating boundary(Ad), 1 = Van Dongeren, 1997 left [-] Switch for lateral boundary at ny+1, 'neumann' = vv computed from NSWE, 'wall' = reflective wall; vv=0 right [-] Switch for lateral boundary at right, 0 = vv computed from NSWE, 1 = reflective wall; vv=0 back [-] Switch for boundary at bay side, 0 = radiating boundary (Ad), 1 = reflective boundary; uu=0 ARC [-] (advanced) Switch for active reflection compensation at seaward boundary: 0 = reflective, 1 = weakly (non) reflective order [-] (advanced) Switch for order of wave steering, 1 = first order wave steering (short wave energy only), 2 = second oder wave steering (bound long wave corresponding to short wave forcing is added) freewave [-] (advanced) Switch for free wave propagation 0 = use cg (default); 1 = use sqrt(gh) in instat = 3 carspan [-] (deprecated) Switch for Carrier-Greenspan test 0 = use cg (default); 1 = use sqrt(gh) in instat = 3 for c&g tests epsi [-] (advanced) Ratio of mean current to time varying current through offshore boundary nc [-] (advanced) Smoothing distance (defined as nr of cells) for estimating umean tidetype [-] (advanced) Switch for offfshore boundary, velocity boundary or instant water level boundary (default) [Section] Tide boundary conditions zs0file [-] Name of tide boundary condition series tideloc [-] Number of corner points on which a tide time series is specified paulrevere [-] Specifies tide on sea and land ('land') or two sea points ('sea') if tideloc = 2 if tideloc =>2, then this indicates where the time series are to be applied. Input for tidal information to xbeach options (3): 1. one tidal record --> specify tidal record everywhere 2. two tidal records --> Need to specify keyword 'paulrevere' paulrevere=='land' implies to apply one tidal record to both sea corners and one tidal record to both land corners paulrevere=='sea' implies to apply the first tidal record (column 2 in zs0input.dat) to the (x=1,y=1) sea corner and the second tidal record (third column) to the (x=1,y=N) sea corner 3. four tidal records --> Need to list tidal records in zs0input.dat in order of: (x=1,y=1) (x=1,y=N) (x=N,y=N) (x=N,y=1) NOTE: clockwise from (1,1) corner [Section] Flow discharge boundary conditions ndischarge [-] (advanced) Number of discharge locations ntdischarge [-] (advanced) Length of discharge time series disch_loc_file [-] (advanced) Name of discharge locations file disch_timeseries_file [-] (advanced) Name of discharge timeseries file [Section] Wave breaking parameters break [-] Type of breaker formulation (1=roelvink, 2=baldock, 3=roelvink adapted, 4=roelvink on/off breaking) gamma [-] Breaker parameter in Baldock or Roelvink formulation gamma2 [-] End of breaking parameter in break = 4 formulation alpha [-] (advanced) Wave dissipation coefficient in Roelvink formulation n [-] (advanced) Power in Roelvink dissipation model gammax [-] (advanced) Maximum ratio wave height to water depth delta [-] (advanced) Fraction of wave height to add to water depth fw [-] (advanced) Bed friction factor fwcutoff [-] Depth greater than which the bed friction factor is NOT applied breakerdelay [-] (advanced) Turn on (1) or off (0) breaker delay model shoaldelay [-] (advanced) Turn on (1) or off (0) shoaling delay facsd [-] (advanced) fraction of the local wave length to use for shoaling delay depth facrun [-] (advanced) calibration coefficient for short wave runup [Section] Roller parameters roller [-] (advanced) Turn on (1) or off(0) roller model beta [-] (advanced) Breaker slope coefficient in roller model rfb [-] (advanced) If rfb = 1 then maximum wave surface slope is feeded back in roller energy balance; else rfb = par%Beta [Section] Wave-current interaction parameters wci [-] Turns on (1) or off (0) wave-current interaction hwci [m] (advanced) Minimum depth until which wave-current interaction is used cats [Trep] (advanced) Current averaging time scale for wci, in terms of mean wave periods [Section] Flow parameters bedfriction [-] Bed friction formulation: 'chezy','white-colebrook' bedfricfile [-] Bed friction file (only valid with values of C) C [m^0.5s^-1] Chezy coefficient cf [-] (advanced) Friction coefficient flow nuh [m^2s^-1] Horizontal background viscosity nuhfac [-] (advanced) Viscosity switch for roller induced turbulent horizontal viscosity nuhv [-] (advanced) Longshore viscosity enhancement factor, following Svendsen (?) smag [-] (advanced) Switch for smagorinsky subgrid model for viscocity [Section] Coriolis force parameters wearth [hour^-1] (advanced) Angular velocity of earth calculated as: 1/rotation_time (in hours), later changed in calculation code to rad/s lat [deg] (advanced) Latitude at model location for computing Coriolis [Section] Wind parameters rhoa [kgm^-3] (advanced) Air density Cd [-] (advanced) Wind drag coefficient windv [ms^-1] Wind velocity, in case of stationary wind windth [deg] Nautical wind direction, in case of stationary wind windfile [-] Name of file with non-stationary wind data [Section] Groundwater parameters kx [ms^-1] (advanced) Darcy-flow permeability coefficient in x-direction [m/s] ky [ms^-1] (advanced) Darcy-flow permeability coefficient in y-direction [m/s] kz [ms^-1] (advanced) Darcy-flow permeability coefficient in z-direction [m/s] dwetlayer [m] (advanced) Thickness of the top soil layer interacting more freely with the surface water aquiferbot [m] (advanced) Level of uniform aquifer bottom aquiferbotfile [-] (advanced) Name of the aquifer bottom file gw0 [m] (advanced) Level initial groundwater level gw0file [-] (advanced) Name of initial groundwater level file gwnonh [-] (advanced) Switch to turn on or off non-hydrostatic pressure for groundwater gwfastsolve [-] (advanced) Reduce full 2D non-hydrostatic solution to quasi-explicit in longshore direction gwscheme [-] (advanced) Scheme for momentum equation (laminar, turbulent) gwReturb [-] (advanced) Reynolds number for start of turbulent flow in case of gwscheme = turbulent gwheadmodel [-] (advanced) Model to use for vertical groundwater head: 'parabolic' (default), or 'exponential' gwhorinfil [-] (advanced) switch to include horizontal infiltration from surface water to groundwater (default = 0) [Section] Q3D sediment transport parameters vonkar (advanced) von Karman constant vicmol (advanced) molecular viscosity kmax [-] (advanced) Number of sigma layers in Quasi-3D model; kmax = 1 (default) is without vertical structure of flow and suspensions sigfac [-] (advanced) dsig scales with log(sigfac). Default = 1.3 [Section] Non-hydrostatic correction parameters solver [-] (advanced) Solver used to solve the linear system, sip, or tridiag (only for 1d) solver_maxit [-] (advanced) Maximum number of iterations in the linear sip solver solver_acc [-] (advanced) accuracy with respect to the right-hand side used in the following termination criterion: ||b-Ax || < acc*||b|| solver_urelax [-] (advanced) Underrelaxation parameter kdmin [-] (advanced) Minimum value of kd ( pi/dx > minkd ) dispc [?] (advanced) Coefficient in front of the vertical pressure gradient, Default = 1. Topt [s] (advanced) Absolute period to optimize coefficient nhbreaker [-] (advanced) Turn on or off nonhydrostatic breaker model breakviscfac [-] (advanced) Factor to increase viscosity during breaking maxbrsteep [-] (advanced) Maximum wave steepness criterium secbrsteep [-] (advanced) Secondary maximum wave steepness criterium reformsteep [-] (advanced) Wave steepness criterium to reform after breaking breakvisclen [-] (advanced) Ratio between local depth and length scale in extra breaking viscosity [Section] Bed composition parameters rhos [kgm^-3] Solid sediment density (no pores) ngd [-] Number of sediment classes nd [-] (advanced) Number of computational layers in the bed dzg1 [m] (advanced) Thickness of top sediment class layers dzg2 [m] (advanced) Nominal thickness of variable sediment class layer dzg3 [m] (advanced) Thickness of bottom sediment class layers por [-] Porosity D50 [m] D50 grain size per grain type D90 [m] D90 grain size per grain type sedcal [-] (advanced) Sediment transport calibration coefficient per grain type ucrcal [-] (advanced) Critical velocity calibration coefficient per grain type [Section] Sediment transport parameters waveform [-] Option for waveshape model: 1 = Ruessink & Van Rijn, 2 = Van Thiel de Vries, 2009 form [-] Equilibrium sed. conc. formulation: 1 = Soulsby van Rijn, 1997, 2 = Van Rijn 2008 with modifications by Van Thiel sws [-] (advanced) 1 = short wave & roller stirring and undertow, 0 = no short wave & roller stirring and undertow lws [-] (advanced) 1 = long wave stirring, 0 = no long wave stirring BRfac [-] (advanced) Calibration factor surface slope facsl [-] (advanced) Factor bedslope effect z0 [m] (advanced) Zero flow velocity level in Soulsby van Rijn (1997) sed.conc. expression smax [-] (advanced) Being tested: maximum Shields parameter for ceq Diane Foster tsfac [-] (advanced) Coefficient determining Ts = tsfac * h/ws in sediment source term facua [-] (advanced) Calibration factor time averaged flows due to wave skewness and asymmetry facSk [-] (advanced) Calibration factor time averaged flows due to wave skewness facAs [-] (advanced) Calibration factor time averaged flows due to wave asymmetry turbadv [-] (advanced) Switch 0/1 to activate turbulence advection model for short and or long wave turbulence turb [-] (advanced) Switch to include short wave turbulence: none = no turbulence, wave_averaged = wave averaged turbulence bore_averaged = maximum turbulence Tbfac [-] (advanced) Calibration factor for bore interval Tbore: Tbore = Tbfac*Tbore Tsmin [s] (advanced) Minimum adaptation time scale in advection diffusion equation sediment lwt [-] (advanced) Switch 0/1 long wave turbulence betad [-] (advanced) Dissipation parameter long wave breaking turbulence swtable [-] (deprecated)Name of intra short wave assymetry and skewness table sus [-] (advanced) Calibration factor for suspensions transports [0..1] bed [-] (advanced) Calibration factor for bed transports [0..1] bulk [-] (advanced) Option to compute bedload and suspended load seperately; 0 = seperately, 1 = bulk (as in previous versions) facDc [-] (advanced) Option to control sediment diffusion coefficient [0..1] jetfac [-] (advanced) Option to mimic turbulence production near revetments [0..1] [Section] Morphology parameters morfac [-] Morphological acceleration factor morfacopt [-] (advanced) Option indicating whether times should be adjusted (1) or not(0) for morfac morstart [s] Start time morphology, in morphological time morstop [s] Stop time morphology, in morphological time wetslp [-] Critical avalanching slope under water (dz/dx and dz/dy) dryslp [-] Critical avalanching slope above water (dz/dx and dz/dy) hswitch [m] (advanced) Water depth at which is switched from wetslp to dryslp dzmax [m/s/m] (advanced) Maximum bedlevel change due to avalanching struct [-] Switch for hard structures ne_layer [name] Name of file containing depth of hard structure [Section] Output variables timings [-] (advanced) Switch to turn on (1) or off (0) progress output to screen tstart [s] Start time of output, in morphological time tint [s] (deprecated) Interval time of global output (replaced by tintg) tintg [s] Interval time of global output tintp [s] Interval time of point and runup gauge output tintc [s] (advanced) Interval time of cross section output tintm [s] Interval time of mean,var,max,min output tsglobal [-] (advanced) Name of file containing timings of global output tspoints [-] (advanced) Name of file containing timings of point output tscross [-] (advanced) Name of file containing timings of cross section output tsmean [-] (advanced) Name of file containing timings of mean, max, min and var output nglobalvar [-] Number of global output variables (as specified by user) globalvars [-] (advanced) Mnems of global output variables, not per se the same sice as nglobalvar (invalid variables, defaults) nmeanvar [-] Number of mean,min,max,var output variables meanvars [-] (advanced) Mnems of mean output variables (by variables) npointvar [-] Number of point output variables pointvars [-] (advanced) Mnems of point output variables (by variables) npoints [-] Number of output point locations nrugauge [-] Number of output runup gauge locations nrugdepth [-] (advanced) Number of depths to compute runup in runup gauge rugdepth [m] (advanced) Minimum depth for determination of last wet point in runup gauge ncross [-] (advanced) Number of output cross sections outputformat [-] (advanced) Choice of output file format: 'netcdf', 'fortran', or 'debug' ncfilename [-] (advanced) xbeach netcdf output file name Projection units (not to be used, only pass to output, this limit is too short for WKT....) This could be the proj4 string +init=epsg:28992 [Section] Output projection projection [-] (advanced) projection string rotate [-] Rotate (1) postprocess output with the rotate function. [Section] Drifters parameters ndrifter [-] Number of drifers drifterfile [-] Name of drifter data file [Section] Ship parameters shipfile [-] Name of ship data file nship [-] (advanced) Number of ships [Section] Vegetation parameters vegiefile [-] Name of vegie species list file vegiemapfile [-] Name of vegie species map file nveg [-] Number of vegetation species [Section] Wave numerics parameters scheme [-] (advanced) Use first-order upwind (upwind_1), second order upwind (upwind_2) or Lax-Wendroff (lax_wendroff) for wave propagation wavint [s] Interval between wave module calls (only in stationary wave mode) maxerror [m] (advanced) Maximum wave height error in wave stationary iteration maxiter [-] (advanced) Maximum number of iterations in wave stationary [Section] Flow numerics parameters eps [m] Threshold water depth above which cells are considered wet eps_sd [m/s] Threshold velocity difference to determine conservation of energy head vs momentum umin [m/s] Threshold velocity for upwind velocity detection and for vmag2 in eq. sediment concentration hmin [m] Threshold water depth above which Stokes drift is included secorder [-] (advanced) Use second order corrections to advection/non-linear terms based on MacCormack scheme oldhu [-] (advanced) Turn on / off old hu calculation [Section] Sediment transport numerics parameters thetanum [-] (advanced) Coefficient determining whether upwind (1) or central scheme (0.5) is used. sourcesink [-] (advanced) In suspended transport use source-sink terms to calculate bed level change (1) or sus transport gradients (0) cmax [-] (advanced) Maximum allowed sediment concentration [Section] Bed update numerics parameters frac_dz [-] (advanced) Relative thickness to split time step for bed updating nd_var [-] (advanced) Index of layer with variable thickness split [-] (advanced) Split threshold for variable sediment layer (ratio to nominal thickness) merge [-] (advanced) Merge threshold for variable sediment layer (ratio to nominal thickness) Prescribed bed evolution parameters nsetbathy [-] (advanced) Number of prescribed bed updates setbathyfile [-] (advanced) Name of prescribed bed update file [Section] MPI parameters mpiboundary [-] (advanced) Fix mpi boundaries along y-lines ('y'), x-lines ('x'), use manual defined domains ('man') or find shortest boundary ('auto') mmpi [-] (advanced) Number of domains in m direction (cross-shore) when manually specifying mpi domains nmpi [-] (advanced) Number of domains in n direction (alongshore) when manually specifying mpi domains