!* This serves as a template for the variables in spacepars !* from this, the declarations, string definitions and other stuff will !* be generated !* format: !* type ndim name dimensions(1..ndim) b/d [unit] description !* type is type of variable, !* ndim is number of dimensions, !* name is name of variable !* dimensions are the dimensions, use constants or names that are !* available in the program, for example s%nx+1 !* b/d = b, if variable is to be broadcasted at start of program !* b/d = d, if variable is to be scattered at start of program !* b/d = 2d, if if variable has dimensions (2,:) and the second dimension is to !* be scatterd !* b/d = 0 if nothing has to be done !* makes no difference for scalars, these are always broadcasted !* The unit should be enclosed between square brackets. If no unit is known or if the variable has no unit it should be marked as [-]. !* The description is everything from the unit onwards until an optional ! !* NOTES: only one variable definition per line, the variable should not continue on the next line !* each item a space-less word, for example !* !* real*8 2 x s%nx+1 s%ny+1 d [m] coordinate of the grid ! good !* real*8 2 x s%nx +1 s%ny+1 d [m] coordinate of the grid ! bad because of space after s%nx !* !* comment starts with ! !* lines starting with !* will be ignored real*8 2 x s%nx+1 s%ny+1 d [m] x-coord. comp. grid (positive shoreward, perp. to coastline) real*8 2 y s%nx+1 s%ny+1 d [m] y-coord. comp. grid real*8 1 xz s%nx+1 d [m] x-coord. comp. grid (positive shoreward, perp. to coastline) real*8 1 yz s%ny+1 d [m] y-coord. comp. grid real*8 1 xu s%nx+1 d [m] x-coord. in u points real*8 1 yv s%ny+1 d [m] y-coord. in v points real*8 2 xw s%nx+1 s%ny+1 d [m] world x-coordinates real*8 2 yw s%nx+1 s%ny+1 d [m] world y-coordinates real*8 0 dx b [m] grid size x-direction real*8 0 dy b [m] grid size y-direction real*8 0 xori b [m] x-origin of grid in world coordinates real*8 0 yori b [m] y-origin of grid in world coordinates real*8 0 alfa b [rad] (deg on input) angle of grid w.r.t. East real*8 0 posdwn b [-] depths defined positive downwards (1) or upwards(-1) real*8 0 thetamin b [rad] minimum angle of computational wave grid (carth. in rad) real*8 0 thetamax b [rad] minimum angle of computational wave grid (carth. in rad) integer 0 nx b [-] local number of grid cells x-direction integer 0 ny b [-] local number of grid cells y-direction real*8 2 zb s%nx+1 s%ny+1 d [m] bed level real*8 2 zb0 s%nx+1 s%ny+1 d [m] initial bed level real*8 1 theta s%ntheta b [rad] wave angles directional distribution w.r.t. comp. x-axis integer 0 ntheta b [-] number of wave direction bins real*8 0 dtheta b [rad] wave direction bin size real*8 0 theta0 b [rad] mean incident wave angle real*8 1 cxsth s%ntheta b [-] cos(theta) real*8 1 sxnth s%ntheta b [-] sin(theta) real*8 2 thetamean s%nx+1 s%ny+1 d [rad] mean wave angle real*8 2 Fx s%nx+1 s%ny+1 d [N/m2] wave force x-direction real*8 2 Fy s%nx+1 s%ny+1 d [N/m2] wave force y-direction real*8 2 Sxy s%nx+1 s%ny+1 d [N/m] radiation stress real*8 2 Syy s%nx+1 s%ny+1 d [N/m] radiation stress real*8 2 Sxx s%nx+1 s%ny+1 d [N/m] radiation stress real*8 2 n s%nx+1 s%ny+1 d [-] ratio group velocity/wave celerity real*8 2 H s%nx+1 s%ny+1 d [m] wave height real*8 3 cgx s%nx+1 s%ny+1 s%ntheta d [m/s] group velocity x-direction real*8 3 cgy s%nx+1 s%ny+1 s%ntheta d [m/s] group velocity y-direction real*8 3 cx s%nx+1 s%ny+1 s%ntheta d [m/s] wave celerity x-direction real*8 3 cy s%nx+1 s%ny+1 s%ntheta d [m/s] wave celerity y-direction real*8 3 ctheta s%nx+1 s%ny+1 s%ntheta d [rad/s] wave celerity theta-direction (refraction) real*8 3 ee s%nx+1 s%ny+1 s%ntheta d [J/m2/rad] directionally distributed wave energy real*8 3 thet s%nx+1 s%ny+1 s%ntheta d [rad] wave angles real*8 3 costhet s%nx+1 s%ny+1 s%ntheta d [-] cos of wave angles real*8 3 sinthet s%nx+1 s%ny+1 s%ntheta d [-] sin of wave angles real*8 3 sigt s%nx+1 s%ny+1 s%ntheta d [rad/s] relative frequency real*8 3 rr s%nx+1 s%ny+1 s%ntheta d [J/m2/rad] directionally distributed roller energy real*8 2 k s%nx+1 s%ny+1 d [rad/m] wave number real*8 2 c s%nx+1 s%ny+1 d [m/s] wave celerity real*8 2 cg s%nx+1 s%ny+1 d [m/s] group velocity real*8 2 sigm s%nx+1 s%ny+1 d [rad/s] mean frequency real*8 2 wm s%nx+1 s%ny+1 d [rad/s] mean abs frequency real*8 2 hh s%nx+1 s%ny+1 d [m] water depth real*8 2 zs s%nx+1 s%ny+1 d [m] water level real*8 2 zs0 s%nx+1 s%ny+1 d [m] water level due to tide alone real*8 1 tideinpt par%tidelen b [s] input time of input tidal signal real*8 2 tideinpz par%tidelen par%tideloc b [m] input tidal signal real*8 1 windinpt par%windlen b [s] input time of input wind signal real*8 1 windvel par%windlen b [m/s] input wind velocity real*8 1 winddir par%windlen b [deg nautical] input wind direction real*8 0 windvnow b [m/s] uniform wind velocity current time real*8 0 winddirnow b [m/s] uniform wind direction current time real*8 2 dzsdt s%nx+1 s%ny+1 d [m/s] rate of change water level real*8 2 dzsdx s%nx+1 s%ny+1 d [m/s] water surface gradient in x-direction real*8 2 dzsdy s%nx+1 s%ny+1 d [m/s] water surface gradient in y-direction real*8 2 dzbdx s%nx+1 s%ny+1 d [-] bed level gradient in x-direction real*8 2 dzbdy s%nx+1 s%ny+1 d [-] bed level gradient in y-direction real*8 2 dzbdt s%nx+1 s%ny+1 d [m/s] rate of change bed level real*8 2 uu s%nx+1 s%ny+1 d [m/s] (GLM) x-velocity in u-points real*8 2 vv s%nx+1 s%ny+1 d [m/s] (GLM) y-velocity in v-points real*8 2 qx s%nx+1 s%ny+1 d [m2/s] x-discharge in u-points real*8 2 qy s%nx+1 s%ny+1 d [m2/s] y-discharge in u-points real*8 2 sedero s%nx+1 s%ny+1 d [m] cum. sedimentation/erosion real*8 2 dcbdx s%nx+1 s%ny+1 d [kg/m3/m] bed concentration gradient x-dir. real*8 2 dcbdy s%nx+1 s%ny+1 d [kg/m3/m] bed concentration gradient y-dir. real*8 2 dcsdx s%nx+1 s%ny+1 d [kg/m3/m] suspended concentration gradient x-dir. real*8 2 dcsdy s%nx+1 s%ny+1 d [kg/m3/m] suspended concentration gradient y-dir. real*8 2 ui s%nx+1 s%ny+1 d [m/s] incident bound wave velocity real*8 2 E s%nx+1 s%ny+1 d [Nm/m2] wave energy real*8 2 R s%nx+1 s%ny+1 d [Nm/m2] roller energy real*8 2 urms s%nx+1 s%ny+1 d [m/s] orbital velocity real*8 2 D s%nx+1 s%ny+1 d [W/m2] dissipation real*8 2 Qb s%nx+1 s%ny+1 d [-] fraction breaking waves real*8 2 ust s%nx+1 s%ny+1 d [m/s] Stokes drift real*8 2 tm s%nx+1 s%ny+1 d [rad] mean wave direction real*8 2 ueu s%nx+1 s%ny+1 d [m/s] Eulerian mean velocity x-dir. real*8 2 vev s%nx+1 s%ny+1 d [m/s] Eulerian mean velocity y-dir. real*8 2 vmagu s%nx+1 s%ny+1 d [m/s] (GLM) velocity magnitude u-points real*8 2 vmageu s%nx+1 s%ny+1 d [m/s] (GLM) velocity magnitude u-points real*8 2 vmagv s%nx+1 s%ny+1 d [m/s] (GLM) velocity magnitude v-points real*8 2 vmagev s%nx+1 s%ny+1 d [m/s] (GLM) velocity magnitude v-points real*8 2 u s%nx+1 s%ny+1 d [m/s] (GLM) x-velocity cell centre (for output) real*8 2 v s%nx+1 s%ny+1 d [m/s] (GLM) y-velocity cell centre (for output) real*8 2 ue s%nx+1 s%ny+1 d [m/s] Eulerian mean x-velocity cell centre (for output) real*8 2 ve s%nx+1 s%ny+1 d [m/s] Eulerian mean y-velocity cell centre (for output) real*8 2 hold s%nx+1 s%ny+1 d [m] water depth previous time step integer 2 wetu s%nx+1 s%ny+1 d [-] mask wet/dry u-points integer 2 wetv s%nx+1 s%ny+1 d [-] mask wet/dry v-points integer 2 wetz s%nx+1 s%ny+1 d [-] mask wet/dry eta-points real*8 2 hu s%nx+1 s%ny+1 d [m] water depth in u-points real*8 2 hv s%nx+1 s%ny+1 d [m] water depth in v-points real*8 2 hum s%nx+1 s%ny+1 d [m] water depth in u-points real*8 2 hvm s%nx+1 s%ny+1 d [m] water depth in v-points real*8 2 vmag s%nx+1 s%ny+1 d [m/s] velocity magnitude in cell centre real*8 3 ccg s%nx+1 s%ny+1 par%ngd d [m3/m3] depth-averaged suspended concentration for each sediment fraction real*8 3 ccbg s%nx+1 s%ny+1 par%ngd d [m3/m3] depth-averaged bed concentration for each sediment fraction real*8 2 uwf s%nx+1 s%ny+1 d [m/s] x-comp. Stokes drift real*8 2 vwf s%nx+1 s%ny+1 d [m/s] y-comp. Stokes drift real*8 2 ustr s%nx+1 s%ny+1 d [m/s] return flow due to roller real*8 2 usd s%nx+1 s%ny+1 d [m/s] return flow due to roller after breaker delay real*8 1 bi s%ny+1 d [m] incoming bound long wave real*8 2 DR s%nx+1 s%ny+1 d [W/m2] roller energy dissipation real*8 2 umean 2 s%ny+1 2 [m/s] longterm mean velocity at bnds in x-direction real*8 2 vmean 2 s%ny+1 2 [m/s] longterm mean velocity at bnds in y-direction integer 0 vardx b [-] 0 = uniform grid size, 1 = variable grid size real*8 2 vu s%nx+1 s%ny+1 d [m/s] y velocity in u points real*8 2 uv s%nx+1 s%ny+1 d [m/s] x velocity in v points real*8 1 D50 par%ngd b [m] D50 grain diameters for all sediment classses real*8 1 D90 par%ngd b [m] D90 grain diameters for all sediment classses real*8 1 sedcal par%ngd b [-] equilibrium sediment concentartion factor for each sediment class real*8 1 ucrcal par%ngd b [-] calibration factor for u critical for each sediment class real*8 3 Tsg s%nx+1 s%ny+1 par%ngd d [s] sediment response time for each sediment class real*8 3 Susg s%nx+1 s%ny+1 par%ngd d [m2/s] suspended sediment transport x-dir. for each sediment class (excluding pores) real*8 3 Svsg s%nx+1 s%ny+1 par%ngd d [m2/s] suspended sediment transport y-dir. for each sediment class (excluding pores) real*8 3 Subg s%nx+1 s%ny+1 par%ngd d [m2/s] bed sediment transport x-dir. for each sediment class (excluding pores) real*8 3 Svbg s%nx+1 s%ny+1 par%ngd d [m2/s] bed sediment transport y-dir. for each sediment class (excluding pores) real*8 3 ceqbg s%nx+1 s%ny+1 par%ngd d [m3/m3] depth-averaged bed equilibrium concentration for each sediment class real*8 3 ceqsg s%nx+1 s%ny+1 par%ngd d [m3/m3] depth-averaged suspended equilibrium concentration for each sediment class real*8 2 ua s%nx+1 s%ny+1 d [m/s] time averaged flow velocity due to wave assymetry real*8 2 BR s%nx+1 s%ny+1 d [-] maximum wave surface slope used in roller dissipation formulation real*8 2 kb s%nx+1 s%ny+1 d [m^2/s^2] near bed turbulence intensity due to depth induces breaking real*8 2 Tbore s%nx+1 s%ny+1 d [s] wave period interval associated with breaking induced turbulence real*8 2 dzav s%nx+1 s%ny+1 d [m] total bed level change due to avalanching real*8 2 maxzs s%nx+1 s%ny+1 d [m] maximum elevation in simulation real*8 2 minzs s%nx+1 s%ny+1 d [m] minimum elevation in simulation real*8 2 L1 s%nx+1 s%ny+1 d for subroutine dispersion real*8 2 Sk s%nx+1 s%ny+1 d [-] skewness of short waves real*8 2 As s%nx+1 s%ny+1 d [-] asymmetry of short waves ! real*8 2 Fimpact s%nx+1 s%ny+1 d ! [N/m^2] Force on dune face real*8 2 gwhead s%nx+1 s%ny+1 d [m] groundwater head (differs from gwlevel) real*8 2 gwlevel s%nx+1 s%ny+1 d [m] groundwater table (min(zb,gwhead)) real*8 2 gwheight s%nx+1 s%ny+1 d [m] vertical size of aquifer through which groundwater can flow real*8 2 gwbottom s%nx+1 s%ny+1 d [m] level of the bottom of the aquifer real*8 2 gwu s%nx+1 s%ny+1 d [m/s] groundwater flow in x-direction real*8 2 gwv s%nx+1 s%ny+1 d [m/s] groundwater flow in y-direction real*8 2 gww s%nx+1 s%ny+1 d [m/s] groundwater flow in z-direction (interaction between surface and ground water) real*8 2 dinfil s%nx+1 s%ny+1 d [m] Infiltration layer depth used in quasi-vertical flow model for groundwater real*8 2 gw0back 2 s%ny+1 2 [m] boundary condition back boundary for groundwater head real*8 2 kturb s%nx+1 s%ny+1 d [m^2/s^2] depth averaged turbulence intensity due to long wave breaking real*8 3 ero s%nx+1 s%ny+1 par%ngd d [m/s] bed erosion rate per fraction real*8 3 depo_im s%nx+1 s%ny+1 par%ngd d [m/s] implicit bed deposition rate per fraction real*8 3 depo_ex s%nx+1 s%ny+1 par%ngd d [m/s] explicit bed deposition rate per fraction integer 2 nd s%nx+1 s%ny+1 d [-] number of bed layers (can be different for each computational cell) real*8 4 pbbed s%nx+1 s%ny+1 max(par%nd,2) par%ngd d [-] NO DESCRIPTION real*8 3 dzbed s%nx+1 s%ny+1 max(par%nd,2) d [-] NO DESCRIPTION real*8 2 z0bed s%nx+1 s%ny+1 d [-] NO DESCRIPTION real*8 2 ureps s%nx+1 s%ny+1 d [m/s] representative flow velocitie for sediment advection and diffusion in x-direction real*8 2 vreps s%nx+1 s%ny+1 d [m/s] representative flow velocitie for sediment advection and diffusion in y-direction real*8 2 urepb s%nx+1 s%ny+1 d [m/s] representative flow velocitie for sediment advection and diffusion in x-direction real*8 2 vrepb s%nx+1 s%ny+1 d [m/s] representative flow velocitie for sediment advection and diffusion in y-direction real*8 2 umwci s%nx+1 s%ny+1 d [m/s] u-velocity (time-averaged) for wci real*8 2 rolthick s%nx+1 s%ny+1 d [m] long wave roller thickness real*8 2 vmwci s%nx+1 s%ny+1 d [m/s] v-velocity (time-averaged) for wci real*8 2 zswci s%nx+1 s%ny+1 d [m] waterlevel (time-averaged) for wci real*8 2 pres s%nx+1 s%ny+1 d [m^2/s^2] normalized dynamic pressure real*8 2 wb s%nx+1 s%ny+1 d [m/s] vertical velocity at the bottom real*8 2 ws s%nx+1 s%ny+1 d [m/s] vertical velocity at the free surface real*8 2 taubx s%nx+1 s%ny+1 d [N/m^2] x-component of bed shear stress real*8 2 tauby s%nx+1 s%ny+1 d [N/m^2] y-component of bed shear stress real*8 2 Df s%nx+1 s%ny+1 d [W/m^2] dissipation rate due to bed friction real*8 2 Dp s%nx+1 s%ny+1 d [W/m^2] dissipation rate in the swash due to transformation of kinetic wave energy to potential wave energy real*8 2 Sutot s%nx+1 s%ny+1 d [m2/s] Sediment transport integrated over bed load and suspended and for all sediment grains real*8 2 Svtot s%nx+1 s%ny+1 d [m2/s] Sediment transport integrated over bed load and suspended and for all sediment grains real*8 2 cctot s%nx+1 s%ny+1 d [m3/m3] Sediment concentration integrated over bed load and suspended and for all sediment grains real*8 2 wi 2 s%ny+1 2 [m/s] Vertical velocity at boundary due to (short) waves2 real*8 2 zi 2 s%ny+1 2 [m] Surface elevation at boundary due to (short) waves real*8 2 nuh s%nx+1 s%ny+1 d [m2/s] horizontal viscosity coefficient real*8 2 cf s%nx+1 s%ny+1 d [-] Friction coefficient flow real*8 2 D50top s%nx+1 s%ny+1 d [-] Friction coefficient flow real*8 2 D90top s%nx+1 s%ny+1 d [-] Friction coefficient flow real*8 2 structdepth s%nx+1 s%ny+1 d [m] Depth of structure in relation to instantaneous bed level !* hints for vi vim: filetype=fortran : syntax=fortran