!* 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] [projection_x_coordinate] x-coord. original comp. grid real*8 2 y s%nx+1 s%ny+1 d [m] [projection_y_coordinate] y-coord. original comp. grid real*8 2 xz s%nx+1 s%ny+1 d [m] [projection_x_coordinate] x-coord. comp. grid (positive shoreward, perp. to coastline) real*8 2 yz s%nx+1 s%ny+1 d [m] [projection_y_coordinate] y-coord. comp. grid real*8 2 xu s%nx+1 s%ny+1 d [m] [projection_x_coordinate] x-coord. comp. grid u-points real*8 2 yu s%nx+1 s%ny+1 d [m] [projection_y_coordinate] y-coord. comp. grid u-points real*8 2 xv s%nx+1 s%ny+1 d [m] [projection_x_coordinate] x-coord. comp. grid v-points real*8 2 yv s%nx+1 s%ny+1 d [m] [projection_y_coordinate] y-coord. comp. grid v-points real*8 2 dsu s%nx+1 s%ny+1 d [m] [] grid distance in s-direction, centered around u-point real*8 2 dsv s%nx+1 s%ny+1 d [m] [] grid distance in s-direction, centered around v-point real*8 2 dsz s%nx+1 s%ny+1 d [m] [] grid distance in s-direction, centered around z-point (=eta-point) real*8 2 dsc s%nx+1 s%ny+1 d [m] [] grid distance in s-direction, centered around c-point real*8 2 dnu s%nx+1 s%ny+1 d [m] [] grid distance in n-direction, centered around u-point real*8 2 dnv s%nx+1 s%ny+1 d [m] [] grid distance in n-direction, centered around v-point real*8 2 dnz s%nx+1 s%ny+1 d [m] [] grid distance in n-direction, centered around z-point (=eta-point) real*8 2 dnc s%nx+1 s%ny+1 d [m] [] grid distance in n-direction, centered around c-point real*8 2 dsdnui s%nx+1 s%ny+1 d [1/m2] [] inverse of grid cell surface, centered around u-point real*8 2 dsdnvi s%nx+1 s%ny+1 d [1/m2] [] inverse of grid cell surface, centered around v-point real*8 2 dsdnzi s%nx+1 s%ny+1 d [1/m2] [] inverse of grid cell surface, centered around z-point real*8 2 alfaz s%nx+1 s%ny+1 d [rad] [angle_of_rotation_from_east_to_x] grid orientation at z-point real*8 2 alfau s%nx+1 s%ny+1 d [rad] [angle_of_rotation_from_east_to_x] grid orientation at u-point real*8 2 alfav s%nx+1 s%ny+1 d [rad] [angle_of_rotation_from_east_to_x] grid orientation at v-point real*8 2 sdist s%nx+1 s%ny+1 d [m] [] cum. distance from offshore boundary along s-direction real*8 2 ndist s%nx+1 s%ny+1 d [m] [] cum. distance from right boundary along n-direction 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] [projection_x_coordinate] x-origin of grid in world coordinates real*8 0 yori b [m] [projection_y_coordinate] y-origin of grid in world coordinates real*8 0 alfa b [rad] [angle_of_rotation_from_east_to_x] (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] [angle_of_rotation_from_east_to_x] minimum angle of computational wave grid (cart. in rad) real*8 0 thetamax b [rad] [angle_of_rotation_from_east_to_x] minimum angle of computational wave grid (cart. 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 0 zs01 b [m] [sea_surface_height_above_sea_level] Initial water level first sea boundary real*8 0 zs02 b [m] [sea_surface_height_above_sea_level] Initial water level second sea boundary real*8 0 zs03 b [m] [sea_surface_height_above_sea_level] Initial water level first land boundary real*8 0 zs04 b [m] [sea_surface_height_above_sea_level] Initial water level second land boundary real*8 1 xyzs01 2 b [-] [] global xy coordinates of corner (x=1,y=1) real*8 1 xyzs02 2 b [-] [] global xy coordinates of corner (x=1,y=N) real*8 1 xyzs03 2 b [-] [] global xy coordinates of corner (x=N,y=N) real*8 1 xyzs04 2 b [-] [] global xy coordinates of corner (x=N,y=1) integer 0 tidelen b [-] [] length of tide time series integer 0 windlen b [-] [] length of tide time series real*8 2 zb s%nx+1 s%ny+1 d [m] [altitude] bed level real*8 2 zb0 s%nx+1 s%ny+1 d [m] [altitude] initial bed level real*8 1 theta s%ntheta b [rad] [sea_surface_wind_wave_to_direction] 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] [sea_surface_wind_wave_to_direction] mean incident wave angle real*8 2 thetamean s%nx+1 s%ny+1 d [rad] [sea_surface_wind_wave_to_direction] mean wave angle real*8 2 Fx s%nx+1 s%ny+1 d [N/m2] [] wave force, x-component real*8 2 Fy s%nx+1 s%ny+1 d [N/m2] [] wave force, y-component real*8 2 Sxy s%nx+1 s%ny+1 d [N/m] [] radiation stress, y-component real*8 2 Syy s%nx+1 s%ny+1 d [N/m] [] radiation stress, y-component real*8 2 Sxx s%nx+1 s%ny+1 d [N/m] [] radiation stress, x-component 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] [] Hrms wave height based on instantaneous wave energy real*8 3 cgx s%nx+1 s%ny+1 s%ntheta d [m/s] [] group velocity, x-component real*8 3 cgy s%nx+1 s%ny+1 s%ntheta d [m/s] [] group velocity, y-component real*8 3 cx s%nx+1 s%ny+1 s%ntheta d [m/s] [] wave celerity, x-component real*8 3 cy s%nx+1 s%ny+1 s%ntheta d [m/s] [] wave celerity, y-component 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] [sea_surface_wind_wave_to_direction] wave angles real*8 3 costh s%nx+1 s%ny+1 s%ntheta d [-] [] cos of wave angles relative to grid direction real*8 3 sinth s%nx+1 s%ny+1 s%ntheta d [-] [] sin of wave angles relative to grid direction 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] [sea_surface_height_above_sea_level] water level real*8 2 zs0 s%nx+1 s%ny+1 d [m] [sea_surface_height_above_sea_level] water level due to tide alone real*8 1 tideinpt s%tidelen b [s] [time] input time of input tidal signal real*8 2 tideinpz s%tidelen par%tideloc b [m] [sea_surface_height_above_sea_level] input tidal signal real*8 1 windinpt s%windlen b [s] [time] input time of input wind signal real*8 1 windvelts s%windlen b [m/s] [] input wind velocity real*8 1 winddirts s%windlen b [deg nautical] [] input wind direction real*8 1 windxts s%windlen b [m/s] [] time series of input wind velocity (not S direction), x-component real*8 1 windyts s%windlen b [m/s] [] time series of input wind velocity (not N direction), y-component real*8 2 windsu s%nx+1 s%ny+1 d [m/s] [] wind velocity in S direction in u point at current time step real*8 2 windnv s%nx+1 s%ny+1 d [m/s] [] wind velocity in N direction in v point at current time step 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] [sea_water_x_velocity] GLM velocity in u-points, x-component real*8 2 vu s%nx+1 s%ny+1 d [m/s] [sea_water_y_velocity] GLM velocity in u-points, y-component real*8 2 uv s%nx+1 s%ny+1 d [m/s] [sea_water_x_velocity] GLM velocity in v-points, x-component real*8 2 vv s%nx+1 s%ny+1 d [m/s] [sea_water_y_velocity] GLM velocity in v-points, y-component real*8 2 qx s%nx+1 s%ny+1 d [m2/s] [] discharge in u-points, x-component real*8 2 qy s%nx+1 s%ny+1 d [m2/s] [] discharge in u-points, y-component 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 in, x-component real*8 2 vi s%nx+1 s%ny+1 d [m/s] [] incident bound wave velocity in, y-component 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] [sea_water_x_velocity] Eulerian velocity in u-points, x-component real*8 2 vev s%nx+1 s%ny+1 d [m/s] [sea_water_y_velocity] Eulerian velocity in u-points, y-component 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] [] Eulerian 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] [] Eulerian velocity magnitude v-points real*8 2 u s%nx+1 s%ny+1 d [m/s] [sea_water_x_velocity] GLM velocity in cell centre, x-component real*8 2 v s%nx+1 s%ny+1 d [m/s] [sea_water_y_velocity] GLM velocity in cell centre, y-component real*8 2 ue s%nx+1 s%ny+1 d [m/s] [sea_water_x_velocity] Eulerian velocity in cell centre, x-component real*8 2 ve s%nx+1 s%ny+1 d [m/s] [sea_water_y_velocity] Eulerian velocity in cell centre, y-component 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] [] Stokes drift, x-component real*8 2 vwf s%nx+1 s%ny+1 d [m/s] [] Stokes drift, y-component 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 s%nx+1 s%ny+1 d [m/s] [sea_water_x_velocity] longterm mean velocity at bnds in u-points, x-component real*8 2 vmean s%nx+1 s%ny+1 d [m/s] [sea_water_y_velocity] longterm mean velocity at bnds in u-points, y-component real*8 2 ur s%nx+1 s%ny+1 d [m/s] [] reflected velocity at bnds in u-points integer 0 vardx b [-] [] 0 = uniform grid size, 1 = variable grid size 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 for each sediment class (excluding pores), x-component real*8 3 Svsg s%nx+1 s%ny+1 par%ngd d [m2/s] [] suspended sediment transport for each sediment class (excluding pores), y-component real*8 3 Subg s%nx+1 s%ny+1 par%ngd d [m2/s] [] bed sediment transport for each sediment class (excluding pores), x-component real*8 3 Svbg s%nx+1 s%ny+1 par%ngd d [m2/s] [] bed sediment transport for each sediment class (excluding pores), y-component 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 [m] [] wave length (used in dispersion relation) 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 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] [sea_water_x_velocity] representative flow velocity for sediment advection and diffusion, x-component real*8 2 vreps s%nx+1 s%ny+1 d [m/s] [sea_water_y_velocity] representative flow velocity for sediment advection and diffusion, y-component real*8 2 urepb s%nx+1 s%ny+1 d [m/s] [sea_water_x_velocity] representative flow velocity for sediment advection and diffusion, x-component real*8 2 vrepb s%nx+1 s%ny+1 d [m/s] [sea_water_y_velocity] representative flow velocity for sediment advection and diffusion, y-component real*8 2 ustar s%nx+1 s%ny+1 d [m/s] [sea_water_x_velocity] bed shear stress velocity (based on Nielsen 1992,2002,2006), x-component real*8 2 vstar s%nx+1 s%ny+1 d [m/s] [sea_water_y_velocity] bed shear stress velocity (based on Nielsen 1992,2002,2006), y-component real*8 2 umwci s%nx+1 s%ny+1 d [m/s] [sea_water_x_velocity] velocity (time-averaged) for wci, x-component real*8 2 vmwci s%nx+1 s%ny+1 d [m/s] [sea_water_y_velocity] velocity (time-averaged) for wci, y-component real*8 2 rolthick s%nx+1 s%ny+1 d [m] [] long wave roller thickness 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] [] bed shear stress, x-component real*8 2 tauby s%nx+1 s%ny+1 d [N/m^2] [] bed shear stress, y-component 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, x-component 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, y-component 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) waves 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 real*8 3 zs0fac s%nx+1 s%ny+1 2 d [-] [] relative weight of offshore boundary and bay boundary for each grid point is stored in zs0fac real*8 1 tdisch par%ntdischarge b [-] [] Discharge time series real*8 2 pdisch par%ndischarge 4 b [-] [] Discharge locations integer 1 pntdisch par%ndischarge b [-] [] Point discharge locations (no momentum) real*8 2 qdisch par%ntdischarge par%ndischarge b [m^2/s] [] Discharges real*8 1 idrift par%ndrifter b [-] [] Drifter x-coordinate in grid space real*8 1 jdrift par%ndrifter b [-] [] Drifter y-coordinate in grid space real*8 1 tdriftb par%ndrifter b [s] [] Drifter release time real*8 1 tdrifte par%ndrifter b [s] [] Drifter retrieval time real*8 1 runup s%ny+1 d [m] [] Short wave runup height real*8 1 Hrunup s%ny+1 d [m] [] Short wave height used in runup formulation real*8 1 xHrunup s%ny+1 d [m] [] Location at which short wave height for runup is taken real*8 2 Dc s%nx+1 s%ny+1 d [m2/s] [] diffusion coefficient real*8 2 ph s%nx+1 s%ny+1 d [m] [] pressure head due to ship real*8 2 zsmean s%nx+1 s%ny+1 d [m] [] time-average water level