! This file is generated by 'scripts/generate.py' using
!    /src/xbeachlibrary/variables.def
!    /src/xbeachlibrary/params.def
!    src/xbeachlibray/templates/indextos.mako
! Advice: do not edit this file, but above mentioned files.

!  DO NOT EDIT THIS FILE
!  But edit variable.f90 and scripts/generate.py
!  Compiling and running is taken care of by the Makefile




 case(  1)
   t%r2  => s%x
   t%rank =  2
   t%type = 'r'
   t%name = 'x'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'projection_x_coordinate'
   t%description= 'x-coord. original cmp. grid'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  2)
   t%r2  => s%y
   t%rank =  2
   t%type = 'r'
   t%name = 'y'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'projection_y_coordinate'
   t%description= 'y-coord. original comp. grid'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  3)
   t%r2  => s%xz
   t%rank =  2
   t%type = 'r'
   t%name = 'xz'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'projection_x_coordinate'
   t%description= 'x-coord. comp. grid (positive shoreward, perp. to coastline)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  4)
   t%r2  => s%yz
   t%rank =  2
   t%type = 'r'
   t%name = 'yz'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'projection_y_coordinate'
   t%description= 'y-coord. comp. grid'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  5)
   t%r2  => s%xu
   t%rank =  2
   t%type = 'r'
   t%name = 'xu'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'projection_x_coordinate'
   t%description= 'x-coord. comp. grid u-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  6)
   t%r2  => s%yu
   t%rank =  2
   t%type = 'r'
   t%name = 'yu'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'projection_y_coordinate'
   t%description= 'y-coord. comp. grid u-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  7)
   t%r2  => s%xv
   t%rank =  2
   t%type = 'r'
   t%name = 'xv'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'projection_x_coordinate'
   t%description= 'x-coord. comp. grid v-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  8)
   t%r2  => s%yv
   t%rank =  2
   t%type = 'r'
   t%name = 'yv'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'projection_y_coordinate'
   t%description= 'y-coord. comp. grid v-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  9)
   t%r2  => s%dsu
   t%rank =  2
   t%type = 'r'
   t%name = 'dsu'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'grid distance in s-direction, centered around u-point'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  10)
   t%r2  => s%dsv
   t%rank =  2
   t%type = 'r'
   t%name = 'dsv'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'grid distance in s-direction, centered around v-point'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  11)
   t%r2  => s%dsz
   t%rank =  2
   t%type = 'r'
   t%name = 'dsz'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'grid distance in s-direction, centered around z-point (=eta-point)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  12)
   t%r2  => s%dsc
   t%rank =  2
   t%type = 'r'
   t%name = 'dsc'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'grid distance in s-direction, centered around c-point'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  13)
   t%r2  => s%dnu
   t%rank =  2
   t%type = 'r'
   t%name = 'dnu'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'grid distance in n-direction, centered around u-point'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  14)
   t%r2  => s%dnv
   t%rank =  2
   t%type = 'r'
   t%name = 'dnv'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'grid distance in n-direction, centered around v-point'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  15)
   t%r2  => s%dnz
   t%rank =  2
   t%type = 'r'
   t%name = 'dnz'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'grid distance in n-direction, centered around z-point (=eta-point)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  16)
   t%r2  => s%dnc
   t%rank =  2
   t%type = 'r'
   t%name = 'dnc'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'grid distance in n-direction, centered around c-point'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  17)
   t%r2  => s%dsdnui
   t%rank =  2
   t%type = 'r'
   t%name = 'dsdnui'
   t%btype = 'd'
   t%units= '1/m2'
   t%standardname= ''
   t%description= 'inverse of grid cell surface, centered around u-point'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  18)
   t%r2  => s%dsdnvi
   t%rank =  2
   t%type = 'r'
   t%name = 'dsdnvi'
   t%btype = 'd'
   t%units= '1/m2'
   t%standardname= ''
   t%description= 'inverse of grid cell surface, centered around v-point'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  19)
   t%r2  => s%dsdnzi
   t%rank =  2
   t%type = 'r'
   t%name = 'dsdnzi'
   t%btype = 'd'
   t%units= '1/m2'
   t%standardname= ''
   t%description= 'inverse of grid cell surface, centered around z-point'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  20)
   t%r2  => s%alfaz
   t%rank =  2
   t%type = 'r'
   t%name = 'alfaz'
   t%btype = 'd'
   t%units= 'rad'
   t%standardname= 'angle_of_rotation_from_east_to_x'
   t%description= 'grid orientation at z-point'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  21)
   t%r2  => s%alfau
   t%rank =  2
   t%type = 'r'
   t%name = 'alfau'
   t%btype = 'd'
   t%units= 'rad'
   t%standardname= 'angle_of_rotation_from_east_to_x'
   t%description= 'grid orientation at u-point'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  22)
   t%r2  => s%alfav
   t%rank =  2
   t%type = 'r'
   t%name = 'alfav'
   t%btype = 'd'
   t%units= 'rad'
   t%standardname= 'angle_of_rotation_from_east_to_x'
   t%description= 'grid orientation at v-point'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  23)
   t%r2  => s%sdist
   t%rank =  2
   t%type = 'r'
   t%name = 'sdist'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'cum. distance from offshore boundary along s-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  24)
   t%r2  => s%ndist
   t%rank =  2
   t%type = 'r'
   t%name = 'ndist'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'cum. distance from right boundary along n-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  25)
   t%r0  => s%dx
   t%rank =  0
   t%type = 'r'
   t%name = 'dx'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= ''
   t%description= 'grid size x-direction'
 case(  26)
   t%r0  => s%dy
   t%rank =  0
   t%type = 'r'
   t%name = 'dy'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= ''
   t%description= 'grid size y-direction'
 case(  27)
   t%r0  => s%xori
   t%rank =  0
   t%type = 'r'
   t%name = 'xori'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= 'projection_x_coordinate'
   t%description= 'x-origin of grid in world coordinate'
 case(  28)
   t%r0  => s%yori
   t%rank =  0
   t%type = 'r'
   t%name = 'yori'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= 'projection_y_coordinate'
   t%description= 'y-origin of grid in world coordinates'
 case(  29)
   t%r0  => s%alfa
   t%rank =  0
   t%type = 'r'
   t%name = 'alfa'
   t%btype = 'b'
   t%units= 'rad'
   t%standardname= 'angle_of_rotation_from_east_to_x'
   t%description= '(deg on input) angle of grid w.r.t. East'
 case(  30)
   t%r0  => s%posdwn
   t%rank =  0
   t%type = 'r'
   t%name = 'posdwn'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'depths defined positive downwards (1) or upwards(-1)'
 case(  31)
   t%r0  => s%thetamin
   t%rank =  0
   t%type = 'r'
   t%name = 'thetamin'
   t%btype = 'b'
   t%units= 'rad'
   t%standardname= 'angle_of_rotation_from_east_to_x'
   t%description= 'minimum angle of computational wave grid (cart. in rad)'
 case(  32)
   t%r0  => s%thetamax
   t%rank =  0
   t%type = 'r'
   t%name = 'thetamax'
   t%btype = 'b'
   t%units= 'rad'
   t%standardname= 'angle_of_rotation_from_east_to_x'
   t%description= 'minimum angle of computational wave grid (cart. in rad)'
 case(  33)
   t%i0  => s%nx
   t%rank =  0
   t%type = 'i'
   t%name = 'nx'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'local number of grid cells x-direction'
 case(  34)
   t%i0  => s%ny
   t%rank =  0
   t%type = 'i'
   t%name = 'ny'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'local number of grid cells y-direction'
 case(  35)
   t%r0  => s%zs01
   t%rank =  0
   t%type = 'r'
   t%name = 'zs01'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= 'sea_surface_height_above_sea_level'
   t%description= 'Initial water level first sea boundary'
 case(  36)
   t%r0  => s%zs02
   t%rank =  0
   t%type = 'r'
   t%name = 'zs02'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= 'sea_surface_height_above_sea_level'
   t%description= 'Initial water level second sea boundary'
 case(  37)
   t%r0  => s%zs03
   t%rank =  0
   t%type = 'r'
   t%name = 'zs03'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= 'sea_surface_height_above_sea_level'
   t%description= 'Initial water level first land boundary'
 case(  38)
   t%r0  => s%zs04
   t%rank =  0
   t%type = 'r'
   t%name = 'zs04'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= 'sea_surface_height_above_sea_level'
   t%description= 'Initial water level second land boundary'
 case(  39)
   t%r1  => s%xyzs01
   t%rank =  1
   t%type = 'r'
   t%name = 'xyzs01'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'global xy coordinates of corner (x=1,y=1)'
   t%dimensions(1:1) = (/ '2                   ' /)
 case(  40)
   t%r1  => s%xyzs02
   t%rank =  1
   t%type = 'r'
   t%name = 'xyzs02'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'global xy coordinates of corner (x=1,y=N)'
   t%dimensions(1:1) = (/ '2                   ' /)
 case(  41)
   t%r1  => s%xyzs03
   t%rank =  1
   t%type = 'r'
   t%name = 'xyzs03'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'global xy coordinates of corner (x=N,y=N)'
   t%dimensions(1:1) = (/ '2                   ' /)
 case(  42)
   t%r1  => s%xyzs04
   t%rank =  1
   t%type = 'r'
   t%name = 'xyzs04'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'global xy coordinates of corner (x=N,y=1)'
   t%dimensions(1:1) = (/ '2                   ' /)
 case(  43)
   t%i0  => s%tidelen
   t%rank =  0
   t%type = 'i'
   t%name = 'tidelen'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'length of tide time series'
 case(  44)
   t%i0  => s%windlen
   t%rank =  0
   t%type = 'i'
   t%name = 'windlen'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'length of tide time series'
 case(  45)
   t%r2  => s%zb
   t%rank =  2
   t%type = 'r'
   t%name = 'zb'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'altitude'
   t%description= 'bed level'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  46)
   t%r2  => s%zb0
   t%rank =  2
   t%type = 'r'
   t%name = 'zb0'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'altitude'
   t%description= 'initial bed level'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  47)
   t%r1  => s%theta
   t%rank =  1
   t%type = 'r'
   t%name = 'theta'
   t%btype = 'b'
   t%units= 'rad'
   t%standardname= 'sea_surface_wind_wave_to_direction'
   t%description= 'wave angles directional distribution w.r.t. comp. x-axis'
   t%dimensions(1:1) = (/ 's%ntheta            ' /)
 case(  48)
   t%r1  => s%theta_s
   t%rank =  1
   t%type = 'r'
   t%name = 'theta_s'
   t%btype = 'b'
   t%units= 'rad'
   t%standardname= 'sea_surface_wind_wave_to_direction'
   t%description= 'wave angles directional distribution w.r.t. comp. x-axis'
   t%dimensions(1:1) = (/ 's%ntheta_s          ' /)
 case(  49)
   t%i0  => s%ntheta
   t%rank =  0
   t%type = 'i'
   t%name = 'ntheta'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'number of wave direction bins'
 case(  50)
   t%i0  => s%ntheta_s
   t%rank =  0
   t%type = 'i'
   t%name = 'ntheta_s'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'number of wave direction bins'
 case(  51)
   t%r0  => s%dtheta
   t%rank =  0
   t%type = 'r'
   t%name = 'dtheta'
   t%btype = 'b'
   t%units= 'rad'
   t%standardname= ''
   t%description= 'wave direction bin size'
 case(  52)
   t%r0  => s%dtheta_s
   t%rank =  0
   t%type = 'r'
   t%name = 'dtheta_s'
   t%btype = 'b'
   t%units= 'rad'
   t%standardname= ''
   t%description= 'wave direction bin size'
 case(  53)
   t%r0  => s%theta0
   t%rank =  0
   t%type = 'r'
   t%name = 'theta0'
   t%btype = 'b'
   t%units= 'rad'
   t%standardname= 'sea_surface_wind_wave_to_direction'
   t%description= 'mean incident wave angle'
 case(  54)
   t%r2  => s%thetamean
   t%rank =  2
   t%type = 'r'
   t%name = 'thetamean'
   t%btype = 'd'
   t%units= 'rad'
   t%standardname= 'sea_surface_wind_wave_to_direction'
   t%description= 'mean wave angle'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  55)
   t%r2  => s%Fx
   t%rank =  2
   t%type = 'r'
   t%name = 'Fx'
   t%btype = 'd'
   t%units= 'N/m2'
   t%standardname= ''
   t%description= 'wave force, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  56)
   t%r2  => s%Fy
   t%rank =  2
   t%type = 'r'
   t%name = 'Fy'
   t%btype = 'd'
   t%units= 'N/m2'
   t%standardname= ''
   t%description= 'wave force, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  57)
   t%r2  => s%Sxy
   t%rank =  2
   t%type = 'r'
   t%name = 'Sxy'
   t%btype = 'd'
   t%units= 'N/m'
   t%standardname= ''
   t%description= 'radiation stress, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  58)
   t%r2  => s%Syy
   t%rank =  2
   t%type = 'r'
   t%name = 'Syy'
   t%btype = 'd'
   t%units= 'N/m'
   t%standardname= ''
   t%description= 'radiation stress, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  59)
   t%r2  => s%Sxx
   t%rank =  2
   t%type = 'r'
   t%name = 'Sxx'
   t%btype = 'd'
   t%units= 'N/m'
   t%standardname= ''
   t%description= 'radiation stress, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  60)
   t%r2  => s%n
   t%rank =  2
   t%type = 'r'
   t%name = 'n'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'ratio group velocity/wave celerity'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  61)
   t%r2  => s%H
   t%rank =  2
   t%type = 'r'
   t%name = 'H'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'Hrms wave height based on instantaneous wave energy'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  62)
   t%r3  => s%cgx
   t%rank =  3
   t%type = 'r'
   t%name = 'cgx'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'group velocity, x-component'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta            ' /)
 case(  63)
   t%r3  => s%cgy
   t%rank =  3
   t%type = 'r'
   t%name = 'cgy'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'group velocity, y-component'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta            ' /)
 case(  64)
   t%r3  => s%cx
   t%rank =  3
   t%type = 'r'
   t%name = 'cx'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'wave celerity, x-component'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta            ' /)
 case(  65)
   t%r3  => s%cy
   t%rank =  3
   t%type = 'r'
   t%name = 'cy'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'wave celerity, y-component'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta            ' /)
 case(  66)
   t%r3  => s%ctheta
   t%rank =  3
   t%type = 'r'
   t%name = 'ctheta'
   t%btype = 'd'
   t%units= 'rad/s'
   t%standardname= ''
   t%description= 'wave celerity theta-direction (refraction)'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta            ' /)
 case(  67)
   t%r3  => s%ee
   t%rank =  3
   t%type = 'r'
   t%name = 'ee'
   t%btype = 'd'
   t%units= 'J/m2/rad'
   t%standardname= ''
   t%description= 'directionally distributed wave energy'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta            ' /)
 case(  68)
   t%r3  => s%thet
   t%rank =  3
   t%type = 'r'
   t%name = 'thet'
   t%btype = 'd'
   t%units= 'rad'
   t%standardname= 'sea_surface_wind_wave_to_direction'
   t%description= 'wave angles'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta            ' /)
 case(  69)
   t%r3  => s%costh
   t%rank =  3
   t%type = 'r'
   t%name = 'costh'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'cos of wave angles relative to grid direction'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta            ' /)
 case(  70)
   t%r3  => s%sinth
   t%rank =  3
   t%type = 'r'
   t%name = 'sinth'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'sin of wave angles relative to grid direction'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta            ' /)
 case(  71)
   t%r3  => s%sigt
   t%rank =  3
   t%type = 'r'
   t%name = 'sigt'
   t%btype = 'd'
   t%units= 'rad/s'
   t%standardname= ''
   t%description= 'relative frequency'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta            ' /)
 case(  72)
   t%r3  => s%rr
   t%rank =  3
   t%type = 'r'
   t%name = 'rr'
   t%btype = 'd'
   t%units= 'J/m2/rad'
   t%standardname= ''
   t%description= 'directionally distributed roller energy'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta            ' /)
 case(  73)
   t%r3  => s%cgx_s
   t%rank =  3
   t%type = 'r'
   t%name = 'cgx_s'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'group velocity, x-component'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta_s          ' /)
 case(  74)
   t%r3  => s%cgy_s
   t%rank =  3
   t%type = 'r'
   t%name = 'cgy_s'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'group velocity, y-component'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta_s          ' /)
 case(  75)
   t%r3  => s%ctheta_s
   t%rank =  3
   t%type = 'r'
   t%name = 'ctheta_s'
   t%btype = 'd'
   t%units= 'rad/s'
   t%standardname= ''
   t%description= 'wave celerity theta-direction (refraction)'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta_s          ' /)
 case(  76)
   t%r3  => s%ee_s
   t%rank =  3
   t%type = 'r'
   t%name = 'ee_s'
   t%btype = 'd'
   t%units= 'J/m2/rad'
   t%standardname= ''
   t%description= 'directionally distributed wave energy'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta_s          ' /)
 case(  77)
   t%r3  => s%thet_s
   t%rank =  3
   t%type = 'r'
   t%name = 'thet_s'
   t%btype = 'd'
   t%units= 'rad'
   t%standardname= 'sea_surface_wind_wave_to_direction'
   t%description= 'wave angles'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta_s          ' /)
 case(  78)
   t%r3  => s%costh_s
   t%rank =  3
   t%type = 'r'
   t%name = 'costh_s'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'cos of wave angles relative to grid direction'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta_s          ' /)
 case(  79)
   t%r3  => s%sinth_s
   t%rank =  3
   t%type = 'r'
   t%name = 'sinth_s'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'sin of wave angles relative to grid direction'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%ntheta_s          ' /)
 case(  80)
   t%r2  => s%k
   t%rank =  2
   t%type = 'r'
   t%name = 'k'
   t%btype = 'd'
   t%units= 'rad/m'
   t%standardname= ''
   t%description= 'wave number'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  81)
   t%r2  => s%c
   t%rank =  2
   t%type = 'r'
   t%name = 'c'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'wave celerity'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  82)
   t%r2  => s%cg
   t%rank =  2
   t%type = 'r'
   t%name = 'cg'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'group velocity'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  83)
   t%r2  => s%sigm
   t%rank =  2
   t%type = 'r'
   t%name = 'sigm'
   t%btype = 'd'
   t%units= 'rad/s'
   t%standardname= ''
   t%description= 'mean frequency'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  84)
   t%r2  => s%wm
   t%rank =  2
   t%type = 'r'
   t%name = 'wm'
   t%btype = 'd'
   t%units= 'rad/s'
   t%standardname= ''
   t%description= 'mean abs frequency'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  85)
   t%r2  => s%hh
   t%rank =  2
   t%type = 'r'
   t%name = 'hh'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'water depth'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  86)
   t%r2  => s%zs
   t%rank =  2
   t%type = 'r'
   t%name = 'zs'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'sea_surface_height_above_sea_level'
   t%description= 'water level'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  87)
   t%r2  => s%zs0
   t%rank =  2
   t%type = 'r'
   t%name = 'zs0'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'sea_surface_height_above_sea_level'
   t%description= 'water level due to tide alone'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  88)
   t%r1  => s%tideinpt
   t%rank =  1
   t%type = 'r'
   t%name = 'tideinpt'
   t%btype = 'b'
   t%units= 's'
   t%standardname= 'time'
   t%description= 'input time of input tidal signal'
   t%dimensions(1:1) = (/ 's%tidelen           ' /)
 case(  89)
   t%r2  => s%tideinpz
   t%rank =  2
   t%type = 'r'
   t%name = 'tideinpz'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= 'sea_surface_height_above_sea_level'
   t%description= 'input tidal signal'
   t%dimensions(1:2) = (/ 's%tidelen           ','par%tideloc         ' /)
 case(  90)
   t%r1  => s%windinpt
   t%rank =  1
   t%type = 'r'
   t%name = 'windinpt'
   t%btype = 'b'
   t%units= 's'
   t%standardname= 'time'
   t%description= 'input time of input wind signal'
   t%dimensions(1:1) = (/ 's%windlen           ' /)
 case(  91)
   t%r1  => s%windvelts
   t%rank =  1
   t%type = 'r'
   t%name = 'windvelts'
   t%btype = 'b'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'input wind velocity'
   t%dimensions(1:1) = (/ 's%windlen           ' /)
 case(  92)
   t%r1  => s%winddirts
   t%rank =  1
   t%type = 'r'
   t%name = 'winddirts'
   t%btype = 'b'
   t%units= 'deg_nautical'
   t%standardname= ''
   t%description= 'input wind direction'
   t%dimensions(1:1) = (/ 's%windlen           ' /)
 case(  93)
   t%r1  => s%windxts
   t%rank =  1
   t%type = 'r'
   t%name = 'windxts'
   t%btype = 'b'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'time series of input wind velocity (not S direction), x-component'
   t%dimensions(1:1) = (/ 's%windlen           ' /)
 case(  94)
   t%r1  => s%windyts
   t%rank =  1
   t%type = 'r'
   t%name = 'windyts'
   t%btype = 'b'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'time series of input wind velocity (not N direction), y-component'
   t%dimensions(1:1) = (/ 's%windlen           ' /)
 case(  95)
   t%r2  => s%windsu
   t%rank =  2
   t%type = 'r'
   t%name = 'windsu'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'wind velocity in S direction in u point at current time step'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  96)
   t%r2  => s%windnv
   t%rank =  2
   t%type = 'r'
   t%name = 'windnv'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'wind velocity in N direction in v point at current time step'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  97)
   t%r2  => s%dzsdt
   t%rank =  2
   t%type = 'r'
   t%name = 'dzsdt'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'rate of change water level'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  98)
   t%r2  => s%dzsdx
   t%rank =  2
   t%type = 'r'
   t%name = 'dzsdx'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'water surface gradient in x-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  99)
   t%r2  => s%dzsdy
   t%rank =  2
   t%type = 'r'
   t%name = 'dzsdy'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'water surface gradient in y-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  100)
   t%r2  => s%dzbdx
   t%rank =  2
   t%type = 'r'
   t%name = 'dzbdx'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'bed level gradient in x-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  101)
   t%r2  => s%dzbdy
   t%rank =  2
   t%type = 'r'
   t%name = 'dzbdy'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'bed level gradient in y-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  102)
   t%r2  => s%dzbdt
   t%rank =  2
   t%type = 'r'
   t%name = 'dzbdt'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'rate of change bed level'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  103)
   t%r2  => s%dzbnow
   t%rank =  2
   t%type = 'r'
   t%name = 'dzbnow'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'bed level change in current time step'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  104)
   t%r2  => s%uu
   t%rank =  2
   t%type = 'r'
   t%name = 'uu'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_x_velocity'
   t%description= 'GLM velocity in u-points, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  105)
   t%r2  => s%vu
   t%rank =  2
   t%type = 'r'
   t%name = 'vu'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_y_velocity'
   t%description= 'GLM velocity in u-points, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  106)
   t%r2  => s%uv
   t%rank =  2
   t%type = 'r'
   t%name = 'uv'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_x_velocity'
   t%description= 'GLM velocity in v-points, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  107)
   t%r2  => s%vv
   t%rank =  2
   t%type = 'r'
   t%name = 'vv'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_y_velocity'
   t%description= 'GLM velocity in v-points, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  108)
   t%r2  => s%qx
   t%rank =  2
   t%type = 'r'
   t%name = 'qx'
   t%btype = 'd'
   t%units= 'm2/s'
   t%standardname= ''
   t%description= 'discharge in u-points, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  109)
   t%r2  => s%qy
   t%rank =  2
   t%type = 'r'
   t%name = 'qy'
   t%btype = 'd'
   t%units= 'm2/s'
   t%standardname= ''
   t%description= 'discharge in u-points, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  110)
   t%r2  => s%sedero
   t%rank =  2
   t%type = 'r'
   t%name = 'sedero'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'cum. sedimentation/erosion'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  111)
   t%r2  => s%dcbdx
   t%rank =  2
   t%type = 'r'
   t%name = 'dcbdx'
   t%btype = 'd'
   t%units= 'kg/m3/m'
   t%standardname= ''
   t%description= 'bed concentration gradient x-dir.'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  112)
   t%r2  => s%dcbdy
   t%rank =  2
   t%type = 'r'
   t%name = 'dcbdy'
   t%btype = 'd'
   t%units= 'kg/m3/m'
   t%standardname= ''
   t%description= 'bed concentration gradient y-dir.'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  113)
   t%r2  => s%dcsdx
   t%rank =  2
   t%type = 'r'
   t%name = 'dcsdx'
   t%btype = 'd'
   t%units= 'kg/m3/m'
   t%standardname= ''
   t%description= 'suspended concentration gradient x-dir.'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  114)
   t%r2  => s%dcsdy
   t%rank =  2
   t%type = 'r'
   t%name = 'dcsdy'
   t%btype = 'd'
   t%units= 'kg/m3/m'
   t%standardname= ''
   t%description= 'suspended concentration gradient y-dir.'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  115)
   t%r2  => s%ui
   t%rank =  2
   t%type = 'r'
   t%name = 'ui'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'incident bound wave velocity in, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  116)
   t%r2  => s%vi
   t%rank =  2
   t%type = 'r'
   t%name = 'vi'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'incident bound wave velocity in, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  117)
   t%r2  => s%E
   t%rank =  2
   t%type = 'r'
   t%name = 'E'
   t%btype = 'd'
   t%units= 'Nm/m2'
   t%standardname= ''
   t%description= 'wave energy'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  118)
   t%r2  => s%R
   t%rank =  2
   t%type = 'r'
   t%name = 'R'
   t%btype = 'd'
   t%units= 'Nm/m2'
   t%standardname= ''
   t%description= 'roller energy'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  119)
   t%r2  => s%urms
   t%rank =  2
   t%type = 'r'
   t%name = 'urms'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'orbital velocity'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  120)
   t%r2  => s%D
   t%rank =  2
   t%type = 'r'
   t%name = 'D'
   t%btype = 'd'
   t%units= 'W/m2'
   t%standardname= ''
   t%description= 'dissipation'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  121)
   t%r2  => s%Qb
   t%rank =  2
   t%type = 'r'
   t%name = 'Qb'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'fraction breaking waves'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  122)
   t%r2  => s%ust
   t%rank =  2
   t%type = 'r'
   t%name = 'ust'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'Stokes drift'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  123)
   t%r2  => s%ueu
   t%rank =  2
   t%type = 'r'
   t%name = 'ueu'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_x_velocity'
   t%description= 'Eulerian velocity in u-points, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  124)
   t%r2  => s%vev
   t%rank =  2
   t%type = 'r'
   t%name = 'vev'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_y_velocity'
   t%description= 'Eulerian velocity in u-points, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  125)
   t%r2  => s%vmagu
   t%rank =  2
   t%type = 'r'
   t%name = 'vmagu'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'GLM velocity magnitude u-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  126)
   t%r2  => s%vmageu
   t%rank =  2
   t%type = 'r'
   t%name = 'vmageu'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'Eulerian velocity magnitude u-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  127)
   t%r2  => s%vmagv
   t%rank =  2
   t%type = 'r'
   t%name = 'vmagv'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'GLM velocity magnitude v-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  128)
   t%r2  => s%vmagev
   t%rank =  2
   t%type = 'r'
   t%name = 'vmagev'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'Eulerian velocity magnitude v-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  129)
   t%r2  => s%u
   t%rank =  2
   t%type = 'r'
   t%name = 'u'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_x_velocity'
   t%description= 'GLM velocity in cell centre, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  130)
   t%r2  => s%v
   t%rank =  2
   t%type = 'r'
   t%name = 'v'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_y_velocity'
   t%description= 'GLM velocity in cell centre, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  131)
   t%r2  => s%ue
   t%rank =  2
   t%type = 'r'
   t%name = 'ue'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_x_velocity'
   t%description= 'Eulerian velocity in cell centre, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  132)
   t%r2  => s%ve
   t%rank =  2
   t%type = 'r'
   t%name = 've'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_y_velocity'
   t%description= 'Eulerian velocity in cell centre, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  133)
   t%r2  => s%ue_sed
   t%rank =  2
   t%type = 'r'
   t%name = 'ue_sed'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_x_velocity'
   t%description= 'Advection velocity sediment in cell centre, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  134)
   t%r2  => s%ve_sed
   t%rank =  2
   t%type = 'r'
   t%name = 've_sed'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_x_velocity'
   t%description= 'Advection velocity sediment in cell centre, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  135)
   t%r2  => s%hold
   t%rank =  2
   t%type = 'r'
   t%name = 'hold'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'water depth previous time step'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  136)
   t%i2  => s%wetu
   t%rank =  2
   t%type = 'i'
   t%name = 'wetu'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'mask wet/dry u-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  137)
   t%i2  => s%wetv
   t%rank =  2
   t%type = 'i'
   t%name = 'wetv'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'mask wet/dry v-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  138)
   t%i2  => s%wetz
   t%rank =  2
   t%type = 'i'
   t%name = 'wetz'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'mask wet/dry eta-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  139)
   t%i2  => s%wete
   t%rank =  2
   t%type = 'i'
   t%name = 'wete'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'mask wet/dry wave-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  140)
   t%r2  => s%hu
   t%rank =  2
   t%type = 'r'
   t%name = 'hu'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'water depth in u-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  141)
   t%r2  => s%hv
   t%rank =  2
   t%type = 'r'
   t%name = 'hv'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'water depth in v-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  142)
   t%r2  => s%hum
   t%rank =  2
   t%type = 'r'
   t%name = 'hum'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'water depth in u-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  143)
   t%r2  => s%hvm
   t%rank =  2
   t%type = 'r'
   t%name = 'hvm'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'water depth in v-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  144)
   t%r2  => s%vmag
   t%rank =  2
   t%type = 'r'
   t%name = 'vmag'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'velocity magnitude in cell centre'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  145)
   t%r2  => s%qmag
   t%rank =  2
   t%type = 'r'
   t%name = 'qmag'
   t%btype = 'd'
   t%units= 'm2/s'
   t%standardname= ''
   t%description= 'depth-velocity product'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  146)
   t%r3  => s%ccg
   t%rank =  3
   t%type = 'r'
   t%name = 'ccg'
   t%btype = 'd'
   t%units= 'm3/m3'
   t%standardname= ''
   t%description= 'depth-averaged suspended concentration for each sediment fraction'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%ngd             ' /)
 case(  147)
   t%r2  => s%uwf
   t%rank =  2
   t%type = 'r'
   t%name = 'uwf'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'Stokes drift, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  148)
   t%r2  => s%vwf
   t%rank =  2
   t%type = 'r'
   t%name = 'vwf'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'Stokes drift, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  149)
   t%r2  => s%ustr
   t%rank =  2
   t%type = 'r'
   t%name = 'ustr'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'return flow due to roller'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  150)
   t%r2  => s%usd
   t%rank =  2
   t%type = 'r'
   t%name = 'usd'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'return flow due to roller after breaker delay'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  151)
   t%r1  => s%bi
   t%rank =  1
   t%type = 'r'
   t%name = 'bi'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'incoming bound long wave'
   t%dimensions(1:1) = (/ 's%ny+1              ' /)
 case(  152)
   t%r2  => s%DR
   t%rank =  2
   t%type = 'r'
   t%name = 'DR'
   t%btype = 'd'
   t%units= 'W/m2'
   t%standardname= ''
   t%description= 'roller energy dissipation'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  153)
   t%r2  => s%umean
   t%rank =  2
   t%type = 'r'
   t%name = 'umean'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_x_velocity'
   t%description= 'longterm mean velocity at bnds in u-points, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  154)
   t%r2  => s%vmean
   t%rank =  2
   t%type = 'r'
   t%name = 'vmean'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_y_velocity'
   t%description= 'longterm mean velocity at bnds in u-points, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  155)
   t%r2  => s%ur
   t%rank =  2
   t%type = 'r'
   t%name = 'ur'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'reflected velocity at bnds in u-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  156)
   t%i0  => s%vardx
   t%rank =  0
   t%type = 'i'
   t%name = 'vardx'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= '0 = uniform grid size, 1 = variable grid size'
 case(  157)
   t%r1  => s%D15
   t%rank =  1
   t%type = 'r'
   t%name = 'D15'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= ''
   t%description= 'D15 grain diameters for all sediment classses'
   t%dimensions(1:1) = (/ 'par%ngd             ' /)
 case(  158)
   t%r1  => s%D50
   t%rank =  1
   t%type = 'r'
   t%name = 'D50'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= ''
   t%description= 'D50 grain diameters for all sediment classses'
   t%dimensions(1:1) = (/ 'par%ngd             ' /)
 case(  159)
   t%r1  => s%D90
   t%rank =  1
   t%type = 'r'
   t%name = 'D90'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= ''
   t%description= 'D90 grain diameters for all sediment classses'
   t%dimensions(1:1) = (/ 'par%ngd             ' /)
 case(  160)
   t%r1  => s%sedcal
   t%rank =  1
   t%type = 'r'
   t%name = 'sedcal'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'equilibrium sediment concentartion factor for each sediment class'
   t%dimensions(1:1) = (/ 'par%ngd             ' /)
 case(  161)
   t%r1  => s%ucrcal
   t%rank =  1
   t%type = 'r'
   t%name = 'ucrcal'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'calibration factor for u critical for each sediment class'
   t%dimensions(1:1) = (/ 'par%ngd             ' /)
 case(  162)
   t%r3  => s%Tsg
   t%rank =  3
   t%type = 'r'
   t%name = 'Tsg'
   t%btype = 'd'
   t%units= 's'
   t%standardname= ''
   t%description= 'sediment response time for each sediment class'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%ngd             ' /)
 case(  163)
   t%r3  => s%Susg
   t%rank =  3
   t%type = 'r'
   t%name = 'Susg'
   t%btype = 'd'
   t%units= 'm2/s'
   t%standardname= ''
   t%description= 'suspended sediment transport for each sediment class (excluding pores), x-component'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%ngd             ' /)
 case(  164)
   t%r3  => s%Svsg
   t%rank =  3
   t%type = 'r'
   t%name = 'Svsg'
   t%btype = 'd'
   t%units= 'm2/s'
   t%standardname= ''
   t%description= 'suspended sediment transport for each sediment class (excluding pores), y-component'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%ngd             ' /)
 case(  165)
   t%r3  => s%Subg
   t%rank =  3
   t%type = 'r'
   t%name = 'Subg'
   t%btype = 'd'
   t%units= 'm2/s'
   t%standardname= ''
   t%description= 'bed sediment transport for each sediment class (excluding pores), x-component'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%ngd             ' /)
 case(  166)
   t%r3  => s%Svbg
   t%rank =  3
   t%type = 'r'
   t%name = 'Svbg'
   t%btype = 'd'
   t%units= 'm2/s'
   t%standardname= ''
   t%description= 'bed sediment transport for each sediment class (excluding pores), y-component'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%ngd             ' /)
 case(  167)
   t%r3  => s%ceqbg
   t%rank =  3
   t%type = 'r'
   t%name = 'ceqbg'
   t%btype = 'd'
   t%units= 'm3/m3'
   t%standardname= ''
   t%description= 'depth-averaged bed equilibrium concentration for each sediment class'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%ngd             ' /)
 case(  168)
   t%r3  => s%ceqsg
   t%rank =  3
   t%type = 'r'
   t%name = 'ceqsg'
   t%btype = 'd'
   t%units= 'm3/m3'
   t%standardname= ''
   t%description= 'depth-averaged suspended equilibrium concentration for each sediment class'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%ngd             ' /)
 case(  169)
   t%r2  => s%ua
   t%rank =  2
   t%type = 'r'
   t%name = 'ua'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'time averaged flow velocity due to wave assymetry'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  170)
   t%r2  => s%BR
   t%rank =  2
   t%type = 'r'
   t%name = 'BR'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'maximum wave surface slope used in roller dissipation formulation'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  171)
   t%r2  => s%kb
   t%rank =  2
   t%type = 'r'
   t%name = 'kb'
   t%btype = 'd'
   t%units= 'm2/s2'
   t%standardname= ''
   t%description= 'near bed turbulence intensity due to depth induces breaking'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  172)
   t%r2  => s%Tbore
   t%rank =  2
   t%type = 'r'
   t%name = 'Tbore'
   t%btype = 'd'
   t%units= 's'
   t%standardname= ''
   t%description= 'wave period interval associated with breaking induced turbulence'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  173)
   t%r2  => s%dzav
   t%rank =  2
   t%type = 'r'
   t%name = 'dzav'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'total bed level change due to avalanching'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  174)
   t%r2  => s%maxzs
   t%rank =  2
   t%type = 'r'
   t%name = 'maxzs'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'maximum elevation in simulation'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  175)
   t%r2  => s%minzs
   t%rank =  2
   t%type = 'r'
   t%name = 'minzs'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'minimum elevation in simulation'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  176)
   t%r2  => s%L1
   t%rank =  2
   t%type = 'r'
   t%name = 'L1'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'wave length (used in dispersion relation)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  177)
   t%r2  => s%Sk
   t%rank =  2
   t%type = 'r'
   t%name = 'Sk'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'skewness of short waves'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  178)
   t%r2  => s%As
   t%rank =  2
   t%type = 'r'
   t%name = 'As'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'asymmetry of short waves'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  179)
   t%r2  => s%gwhead
   t%rank =  2
   t%type = 'r'
   t%name = 'gwhead'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'groundwater head (differs from gwlevel)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  180)
   t%r2  => s%gwheadb
   t%rank =  2
   t%type = 'r'
   t%name = 'gwheadb'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'groundwater head at bottom (differs from gwlevel)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  181)
   t%r2  => s%gwlevel
   t%rank =  2
   t%type = 'r'
   t%name = 'gwlevel'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'groundwater table (min(zb,gwhead))'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  182)
   t%r2  => s%gwheight
   t%rank =  2
   t%type = 'r'
   t%name = 'gwheight'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'vertical size of aquifer through which groundwater can flow'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  183)
   t%r2  => s%gwbottom
   t%rank =  2
   t%type = 'r'
   t%name = 'gwbottom'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'level of the bottom of the aquifer'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  184)
   t%r2  => s%gwu
   t%rank =  2
   t%type = 'r'
   t%name = 'gwu'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'groundwater flow in x-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  185)
   t%r2  => s%gwv
   t%rank =  2
   t%type = 'r'
   t%name = 'gwv'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'groundwater flow in y-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  186)
   t%r2  => s%gwqx
   t%rank =  2
   t%type = 'r'
   t%name = 'gwqx'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'groundwater discharge in x-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  187)
   t%r2  => s%gwqy
   t%rank =  2
   t%type = 'r'
   t%name = 'gwqy'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'groundwater discharge in y-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  188)
   t%r2  => s%gww
   t%rank =  2
   t%type = 'r'
   t%name = 'gww'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'groundwater flow in z-direction (interaction between surface and ground water)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  189)
   t%r2  => s%gwcurv
   t%rank =  2
   t%type = 'r'
   t%name = 'gwcurv'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'curvature coefficient of groundwater head function'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  190)
   t%r2  => s%dinfil
   t%rank =  2
   t%type = 'r'
   t%name = 'dinfil'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'Infiltration layer depth used in quasi-vertical flow model for groundwater'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  191)
   t%r2  => s%infil
   t%rank =  2
   t%type = 'r'
   t%name = 'infil'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'Rate of exchange of water between surface and groundwater (positive from sea to groundwater)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  192)
   t%r2  => s%gw0back
   t%rank =  2
   t%type = 'r'
   t%name = 'gw0back'
   t%btype = '2'
   t%units= 'm'
   t%standardname= ''
   t%description= 'boundary condition back boundary for groundwater head'
   t%dimensions(1:2) = (/ '2                   ','s%ny+1              ' /)
 case(  193)
   t%r2  => s%Kx
   t%rank =  2
   t%type = 'r'
   t%name = 'Kx'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= '(Turbulent) Hydraulic conductivity in x-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  194)
   t%r2  => s%Ky
   t%rank =  2
   t%type = 'r'
   t%name = 'Ky'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= '(Turbulent) Hydraulic conductivity in y-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  195)
   t%r2  => s%Kz
   t%rank =  2
   t%type = 'r'
   t%name = 'Kz'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= '(Turbulent) Hydraulic conductivity in z-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  196)
   t%r2  => s%Kzinf
   t%rank =  2
   t%type = 'r'
   t%name = 'Kzinf'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= '(Turbulent) Hydraulic conductivity in z-direction for infiltration'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  197)
   t%r2  => s%kturb
   t%rank =  2
   t%type = 'r'
   t%name = 'kturb'
   t%btype = 'd'
   t%units= 'm2/s2'
   t%standardname= ''
   t%description= 'depth averaged turbulence intensity due to long wave breaking'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  198)
   t%r3  => s%ero
   t%rank =  3
   t%type = 'r'
   t%name = 'ero'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'bed erosion rate per fraction'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%ngd             ' /)
 case(  199)
   t%r3  => s%depo_im
   t%rank =  3
   t%type = 'r'
   t%name = 'depo_im'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'implicit bed deposition rate per fraction'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%ngd             ' /)
 case(  200)
   t%r3  => s%depo_ex
   t%rank =  3
   t%type = 'r'
   t%name = 'depo_ex'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'explicit bed deposition rate per fraction'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%ngd             ' /)
 case(  201)
   t%i2  => s%nd
   t%rank =  2
   t%type = 'i'
   t%name = 'nd'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'number of bed layers (can be different for each computational cell)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  202)
   t%r4  => s%pbbed
   t%rank =  4
   t%type = 'r'
   t%name = 'pbbed'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'NO DESCRIPTION'
   t%dimensions(1:4) = (/ 's%nx+1              ','s%ny+1              ','par%nd              ','par%ngd             ' /)
 case(  203)
   t%r3  => s%dzbed
   t%rank =  3
   t%type = 'r'
   t%name = 'dzbed'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'NO DESCRIPTION'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%nd              ' /)
 case(  204)
   t%r2  => s%z0bed
   t%rank =  2
   t%type = 'r'
   t%name = 'z0bed'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'NO DESCRIPTION'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  205)
   t%r2  => s%ureps
   t%rank =  2
   t%type = 'r'
   t%name = 'ureps'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_x_velocity'
   t%description= 'representative flow velocity for sediment advection and diffusion, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  206)
   t%r2  => s%vreps
   t%rank =  2
   t%type = 'r'
   t%name = 'vreps'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_y_velocity'
   t%description= 'representative flow velocity for sediment advection and diffusion, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  207)
   t%r2  => s%urepb
   t%rank =  2
   t%type = 'r'
   t%name = 'urepb'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_x_velocity'
   t%description= 'representative flow velocity for sediment advection and diffusion, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  208)
   t%r2  => s%vrepb
   t%rank =  2
   t%type = 'r'
   t%name = 'vrepb'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_y_velocity'
   t%description= 'representative flow velocity for sediment advection and diffusion, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  209)
   t%r2  => s%umwci
   t%rank =  2
   t%type = 'r'
   t%name = 'umwci'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_x_velocity'
   t%description= 'velocity (time-averaged) for wci, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  210)
   t%r2  => s%vmwci
   t%rank =  2
   t%type = 'r'
   t%name = 'vmwci'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= 'sea_water_y_velocity'
   t%description= 'velocity (time-averaged) for wci, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  211)
   t%r2  => s%rolthick
   t%rank =  2
   t%type = 'r'
   t%name = 'rolthick'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'long wave roller thickness'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  212)
   t%r2  => s%zswci
   t%rank =  2
   t%type = 'r'
   t%name = 'zswci'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'waterlevel (time-averaged) for wci'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  213)
   t%r2  => s%pres
   t%rank =  2
   t%type = 'r'
   t%name = 'pres'
   t%btype = 'd'
   t%units= 'm2/s2'
   t%standardname= ''
   t%description= 'normalized dynamic pressure'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  214)
   t%r2  => s%dU
   t%rank =  2
   t%type = 'r'
   t%name = 'dU'
   t%btype = 'd'
   t%units= 'm2/s2'
   t%standardname= ''
   t%description= 'u-velocity difference between two vertical layers (reduced 2-layer non-hydrostatic model)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  215)
   t%r2  => s%dV
   t%rank =  2
   t%type = 'r'
   t%name = 'dV'
   t%btype = 'd'
   t%units= 'm2/s2'
   t%standardname= ''
   t%description= 'v-velocity difference between two vertical layers (reduced 2-layer non-hydrostatic model)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  216)
   t%r2  => s%wb
   t%rank =  2
   t%type = 'r'
   t%name = 'wb'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'vertical velocity at the bottom'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  217)
   t%r2  => s%ws
   t%rank =  2
   t%type = 'r'
   t%name = 'ws'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'vertical velocity at the free surface'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  218)
   t%r2  => s%wscrit
   t%rank =  2
   t%type = 'r'
   t%name = 'wscrit'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'critial vertical velocity at the free surface for breaking'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  219)
   t%r2  => s%bedfriccoef
   t%rank =  2
   t%type = 'r'
   t%name = 'bedfriccoef'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'dimensional/dimensionless input bed friction coefficient; depends on value of parameter bedfriction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  220)
   t%r2  => s%taubx
   t%rank =  2
   t%type = 'r'
   t%name = 'taubx'
   t%btype = 'd'
   t%units= 'N/m2'
   t%standardname= ''
   t%description= 'bed shear stress, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  221)
   t%r2  => s%tauby
   t%rank =  2
   t%type = 'r'
   t%name = 'tauby'
   t%btype = 'd'
   t%units= 'N/m2'
   t%standardname= ''
   t%description= 'bed shear stress, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  222)
   t%r2  => s%Df
   t%rank =  2
   t%type = 'r'
   t%name = 'Df'
   t%btype = 'd'
   t%units= 'W/m2'
   t%standardname= ''
   t%description= 'dissipation rate due to bed friction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  223)
   t%r2  => s%Dp
   t%rank =  2
   t%type = 'r'
   t%name = 'Dp'
   t%btype = 'd'
   t%units= 'W/m2'
   t%standardname= ''
   t%description= 'dissipation rate in the swash due to transformation of kinetic wave energy to potential wave energy'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  224)
   t%r2  => s%Sutot
   t%rank =  2
   t%type = 'r'
   t%name = 'Sutot'
   t%btype = 'd'
   t%units= 'm2/s'
   t%standardname= ''
   t%description= 'Sediment transport integrated over bed load and suspended and for all sediment grains, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  225)
   t%r2  => s%Svtot
   t%rank =  2
   t%type = 'r'
   t%name = 'Svtot'
   t%btype = 'd'
   t%units= 'm2/s'
   t%standardname= ''
   t%description= 'Sediment transport integrated over bed load and suspended and for all sediment grains, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  226)
   t%r2  => s%cctot
   t%rank =  2
   t%type = 'r'
   t%name = 'cctot'
   t%btype = 'd'
   t%units= 'm3/m3'
   t%standardname= ''
   t%description= 'Sediment concentration integrated over bed load and suspended and for all sediment grains'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  227)
   t%r2  => s%wi
   t%rank =  2
   t%type = 'r'
   t%name = 'wi'
   t%btype = '2'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'Vertical velocity at boundary due to (short) waves'
   t%dimensions(1:2) = (/ '2                   ','s%ny+1              ' /)
 case(  228)
   t%r2  => s%dUi
   t%rank =  2
   t%type = 'r'
   t%name = 'dUi'
   t%btype = '2'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'Velocity difference at boundary due to (short) waves'
   t%dimensions(1:2) = (/ '2                   ','s%ny+1              ' /)
 case(  229)
   t%r2  => s%dVi
   t%rank =  2
   t%type = 'r'
   t%name = 'dVi'
   t%btype = '2'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'Velocity difference at boundary due to (short) waves'
   t%dimensions(1:2) = (/ '2                   ','s%ny+1              ' /)
 case(  230)
   t%r2  => s%zi
   t%rank =  2
   t%type = 'r'
   t%name = 'zi'
   t%btype = '2'
   t%units= 'm'
   t%standardname= ''
   t%description= 'Surface elevation at boundary due to (short) waves'
   t%dimensions(1:2) = (/ '2                   ','s%ny+1              ' /)
 case(  231)
   t%r2  => s%nuh
   t%rank =  2
   t%type = 'r'
   t%name = 'nuh'
   t%btype = 'd'
   t%units= 'm2/s'
   t%standardname= ''
   t%description= 'horizontal viscosity coefficient'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  232)
   t%r2  => s%cf
   t%rank =  2
   t%type = 'r'
   t%name = 'cf'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'Friction coefficient flow'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  233)
   t%r2  => s%cfu
   t%rank =  2
   t%type = 'r'
   t%name = 'cfu'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'Friction coefficient flow in u-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  234)
   t%r2  => s%cfv
   t%rank =  2
   t%type = 'r'
   t%name = 'cfv'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'Friction coefficient flow in v-points'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  235)
   t%r2  => s%D50top
   t%rank =  2
   t%type = 'r'
   t%name = 'D50top'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'Friction coefficient flow'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  236)
   t%r2  => s%D90top
   t%rank =  2
   t%type = 'r'
   t%name = 'D90top'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'Friction coefficient flow'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  237)
   t%r2  => s%structdepth
   t%rank =  2
   t%type = 'r'
   t%name = 'structdepth'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'Depth of structure in relation to instantaneous bed level'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  238)
   t%r3  => s%zs0fac
   t%rank =  3
   t%type = 'r'
   t%name = 'zs0fac'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'relative weight of offshore boundary and bay boundary for each grid point is stored in zs0fac'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','2                   ' /)
 case(  239)
   t%r1  => s%tdisch
   t%rank =  1
   t%type = 'r'
   t%name = 'tdisch'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'Discharge time series'
   t%dimensions(1:1) = (/ 'par%ntdischarge     ' /)
 case(  240)
   t%r2  => s%pdisch
   t%rank =  2
   t%type = 'r'
   t%name = 'pdisch'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'Discharge locations'
   t%dimensions(1:2) = (/ 'par%ndischarge      ','4                   ' /)
 case(  241)
   t%i1  => s%pntdisch
   t%rank =  1
   t%type = 'i'
   t%name = 'pntdisch'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'Point discharge locations (no momentum)'
   t%dimensions(1:1) = (/ 'par%ndischarge      ' /)
 case(  242)
   t%r2  => s%qdisch
   t%rank =  2
   t%type = 'r'
   t%name = 'qdisch'
   t%btype = 'b'
   t%units= 'm2/s'
   t%standardname= ''
   t%description= 'Discharges'
   t%dimensions(1:2) = (/ 'par%ntdischarge     ','par%ndischarge      ' /)
 case(  243)
   t%r1  => s%idrift
   t%rank =  1
   t%type = 'r'
   t%name = 'idrift'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'Drifter x-coordinate in grid space'
   t%dimensions(1:1) = (/ 'par%ndrifter        ' /)
 case(  244)
   t%r1  => s%jdrift
   t%rank =  1
   t%type = 'r'
   t%name = 'jdrift'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'Drifter y-coordinate in grid space'
   t%dimensions(1:1) = (/ 'par%ndrifter        ' /)
 case(  245)
   t%r1  => s%tdriftb
   t%rank =  1
   t%type = 'r'
   t%name = 'tdriftb'
   t%btype = 'b'
   t%units= 's'
   t%standardname= ''
   t%description= 'Drifter release time'
   t%dimensions(1:1) = (/ 'par%ndrifter        ' /)
 case(  246)
   t%r1  => s%tdrifte
   t%rank =  1
   t%type = 'r'
   t%name = 'tdrifte'
   t%btype = 'b'
   t%units= 's'
   t%standardname= ''
   t%description= 'Drifter retrieval time'
   t%dimensions(1:1) = (/ 'par%ndrifter        ' /)
 case(  247)
   t%r1  => s%runup
   t%rank =  1
   t%type = 'r'
   t%name = 'runup'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'Short wave runup height'
   t%dimensions(1:1) = (/ 's%ny+1              ' /)
 case(  248)
   t%r1  => s%Hrunup
   t%rank =  1
   t%type = 'r'
   t%name = 'Hrunup'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'Short wave height used in runup formulation'
   t%dimensions(1:1) = (/ 's%ny+1              ' /)
 case(  249)
   t%r1  => s%xHrunup
   t%rank =  1
   t%type = 'r'
   t%name = 'xHrunup'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'Location at which short wave height for runup is taken'
   t%dimensions(1:1) = (/ 's%ny+1              ' /)
 case(  250)
   t%r1  => s%istruct
   t%rank =  1
   t%type = 'r'
   t%name = 'istruct'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'location of revetments toe'
   t%dimensions(1:1) = (/ 's%ny+1              ' /)
 case(  251)
   t%r1  => s%iwl
   t%rank =  1
   t%type = 'r'
   t%name = 'iwl'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'location of water line (including long wave runup)'
   t%dimensions(1:1) = (/ 's%ny+1              ' /)
 case(  252)
   t%r1  => s%strucslope
   t%rank =  1
   t%type = 'r'
   t%name = 'strucslope'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'slope of structure'
   t%dimensions(1:1) = (/ 's%ny+1              ' /)
 case(  253)
   t%r2  => s%Dc
   t%rank =  2
   t%type = 'r'
   t%name = 'Dc'
   t%btype = 'd'
   t%units= 'm2/s'
   t%standardname= ''
   t%description= 'diffusion coefficient'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  254)
   t%r2  => s%ph
   t%rank =  2
   t%type = 'r'
   t%name = 'ph'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'pressure head due to ship'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  255)
   t%i0  => s%newstatbc
   t%rank =  0
   t%type = 'i'
   t%name = 'newstatbc'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= '(1) Use new stationary boundary conditions for instat is stat or stat_table'
 case(  256)
   t%i0  => s%isSet_Vbc
   t%rank =  0
   t%type = 'i'
   t%name = 'isSet_Vbc'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= '(1) V wave boundary condition vlocities in V direction set (1), or unset (0)'
 case(  257)
   t%r2  => s%dobs
   t%rank =  2
   t%type = 'r'
   t%name = 'dobs'
   t%btype = 'd'
   t%units= 'W/m2'
   t%standardname= ''
   t%description= 'beachwizard observed dissipation'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  258)
   t%r2  => s%sig2prior
   t%rank =  2
   t%type = 'r'
   t%name = 'sig2prior'
   t%btype = 'd'
   t%units= 'm2'
   t%standardname= ''
   t%description= 'beachwizard prior std squared'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  259)
   t%r2  => s%zbobs
   t%rank =  2
   t%type = 'r'
   t%name = 'zbobs'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'beachwizard observed depth'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  260)
   t%r2  => s%shobs
   t%rank =  2
   t%type = 'r'
   t%name = 'shobs'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'beachwizard observed shoreline'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  261)
   t%r2  => s%bwalpha
   t%rank =  2
   t%type = 'r'
   t%name = 'bwalpha'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'beachwizard weighting factor'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  262)
   t%r2  => s%dcmdo
   t%rank =  2
   t%type = 'r'
   t%name = 'dcmdo'
   t%btype = 'd'
   t%units= 'W/m2'
   t%standardname= ''
   t%description= 'beachwizard computed minus observed dissipation'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  263)
   t%r2  => s%dassim
   t%rank =  2
   t%type = 'r'
   t%name = 'dassim'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'beachwizard depth change'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  264)
   t%r2  => s%cobs
   t%rank =  2
   t%type = 'r'
   t%name = 'cobs'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'beachwizard observed wave celerity'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  265)
   t%r1  => s%shipxCG
   t%rank =  1
   t%type = 'r'
   t%name = 'shipxCG'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= ''
   t%description= 'x-coordinate of ship center of gravity'
   t%dimensions(1:1) = (/ 'par%nship           ' /)
 case(  266)
   t%r1  => s%shipyCG
   t%rank =  1
   t%type = 'r'
   t%name = 'shipyCG'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= ''
   t%description= 'y-coordinate of ship center of gravity'
   t%dimensions(1:1) = (/ 'par%nship           ' /)
 case(  267)
   t%r1  => s%shipzCG
   t%rank =  1
   t%type = 'r'
   t%name = 'shipzCG'
   t%btype = 'b'
   t%units= 'm'
   t%standardname= ''
   t%description= 'z-coordinate of ship center of gravity'
   t%dimensions(1:1) = (/ 'par%nship           ' /)
 case(  268)
   t%r1  => s%shipFx
   t%rank =  1
   t%type = 'r'
   t%name = 'shipFx'
   t%btype = 'b'
   t%units= 'N'
   t%standardname= ''
   t%description= 'force on ship in x-direction'
   t%dimensions(1:1) = (/ 'par%nship           ' /)
 case(  269)
   t%r1  => s%shipFy
   t%rank =  1
   t%type = 'r'
   t%name = 'shipFy'
   t%btype = 'b'
   t%units= 'N'
   t%standardname= ''
   t%description= 'force on ship in y-direction'
   t%dimensions(1:1) = (/ 'par%nship           ' /)
 case(  270)
   t%r1  => s%shipFz
   t%rank =  1
   t%type = 'r'
   t%name = 'shipFz'
   t%btype = 'b'
   t%units= 'N'
   t%standardname= ''
   t%description= 'force on ship in z-direction'
   t%dimensions(1:1) = (/ 'par%nship           ' /)
 case(  271)
   t%r1  => s%shipMx
   t%rank =  1
   t%type = 'r'
   t%name = 'shipMx'
   t%btype = 'b'
   t%units= 'Nm'
   t%standardname= ''
   t%description= 'moment on ship around x-axis'
   t%dimensions(1:1) = (/ 'par%nship           ' /)
 case(  272)
   t%r1  => s%shipMy
   t%rank =  1
   t%type = 'r'
   t%name = 'shipMy'
   t%btype = 'b'
   t%units= 'Nm'
   t%standardname= ''
   t%description= 'moment on ship around y-axis'
   t%dimensions(1:1) = (/ 'par%nship           ' /)
 case(  273)
   t%r1  => s%shipMz
   t%rank =  1
   t%type = 'r'
   t%name = 'shipMz'
   t%btype = 'b'
   t%units= 'Nm'
   t%standardname= ''
   t%description= 'moment on ship around z-axis'
   t%dimensions(1:1) = (/ 'par%nship           ' /)
 case(  274)
   t%r1  => s%shipphi
   t%rank =  1
   t%type = 'r'
   t%name = 'shipphi'
   t%btype = 'b'
   t%units= 'deg'
   t%standardname= ''
   t%description= 'turning angle arround x-axis'
   t%dimensions(1:1) = (/ 'par%nship           ' /)
 case(  275)
   t%r1  => s%shipchi
   t%rank =  1
   t%type = 'r'
   t%name = 'shipchi'
   t%btype = 'b'
   t%units= 'deg'
   t%standardname= ''
   t%description= 'turning angle arround y-axis'
   t%dimensions(1:1) = (/ 'par%nship           ' /)
 case(  276)
   t%r1  => s%shippsi
   t%rank =  1
   t%type = 'r'
   t%name = 'shippsi'
   t%btype = 'b'
   t%units= 'deg'
   t%standardname= ''
   t%description= 'turning angle arround z-axis'
   t%dimensions(1:1) = (/ 'par%nship           ' /)
 case(  277)
   t%i2  => s%vegtype
   t%rank =  2
   t%type = 'i'
   t%name = 'vegtype'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'vegetation type index'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  278)
   t%i2  => s%nsecveg
   t%rank =  2
   t%type = 'i'
   t%name = 'nsecveg'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'vegetation number of sections'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  279)
   t%i0  => s%nsecvegmax
   t%rank =  0
   t%type = 'i'
   t%name = 'nsecvegmax'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'maximum number of vegetation sections within domain'
 case(  280)
   t%r3  => s%Cdveg
   t%rank =  3
   t%type = 'r'
   t%name = 'Cdveg'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'vegetation drag coefficient'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%nsecvegmax        ' /)
 case(  281)
   t%r3  => s%ahveg
   t%rank =  3
   t%type = 'r'
   t%name = 'ahveg'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'vegetation height'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%nsecvegmax        ' /)
 case(  282)
   t%r3  => s%bveg
   t%rank =  3
   t%type = 'r'
   t%name = 'bveg'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'vegetation stem diameter'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%nsecvegmax        ' /)
 case(  283)
   t%r3  => s%Nveg
   t%rank =  3
   t%type = 'r'
   t%name = 'Nveg'
   t%btype = 'd'
   t%units= 'm-2'
   t%standardname= ''
   t%description= 'vegetation density'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','s%nsecvegmax        ' /)
 case(  284)
   t%r2  => s%Dveg
   t%rank =  2
   t%type = 'r'
   t%name = 'Dveg'
   t%btype = 'd'
   t%units= 'W/m2'
   t%standardname= ''
   t%description= 'dissipation due to short wave attenuation by vegetation'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  285)
   t%r2  => s%Fvegu
   t%rank =  2
   t%type = 'r'
   t%name = 'Fvegu'
   t%btype = 'd'
   t%units= 'N/m2'
   t%standardname= ''
   t%description= 'x-forcing due to long wave attenuation by vegetation'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  286)
   t%r2  => s%Fvegv
   t%rank =  2
   t%type = 'r'
   t%name = 'Fvegv'
   t%btype = 'd'
   t%units= 'N/m2'
   t%standardname= ''
   t%description= 'y-forcing due to long wave attenuation by vegetation'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  287)
   t%r2  => s%ududx
   t%rank =  2
   t%type = 'r'
   t%name = 'ududx'
   t%btype = 'd'
   t%units= 'm2/s2'
   t%standardname= ''
   t%description= 'advection'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  288)
   t%r2  => s%vdvdy
   t%rank =  2
   t%type = 'r'
   t%name = 'vdvdy'
   t%btype = 'd'
   t%units= 'm2/s2'
   t%standardname= ''
   t%description= 'advection'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  289)
   t%r2  => s%udvdx
   t%rank =  2
   t%type = 'r'
   t%name = 'udvdx'
   t%btype = 'd'
   t%units= 'm2/s2'
   t%standardname= ''
   t%description= 'advection'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  290)
   t%r2  => s%vdudy
   t%rank =  2
   t%type = 'r'
   t%name = 'vdudy'
   t%btype = 'd'
   t%units= 'm2/s2'
   t%standardname= ''
   t%description= 'advection'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  291)
   t%r2  => s%viscu
   t%rank =  2
   t%type = 'r'
   t%name = 'viscu'
   t%btype = 'd'
   t%units= 'm2/s2'
   t%standardname= ''
   t%description= 'viscosity'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  292)
   t%r2  => s%viscv
   t%rank =  2
   t%type = 'r'
   t%name = 'viscv'
   t%btype = 'd'
   t%units= 'm2/s2'
   t%standardname= ''
   t%description= 'viscosity'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  293)
   t%r3  => s%setbathy
   t%rank =  3
   t%type = 'r'
   t%name = 'setbathy'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'prescribed bed levels'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%nsetbathy       ' /)
 case(  294)
   t%r1  => s%tsetbathy
   t%rank =  1
   t%type = 'r'
   t%name = 'tsetbathy'
   t%btype = 'b'
   t%units= 's'
   t%standardname= ''
   t%description= 'points in time of prescibed bed levels'
   t%dimensions(1:1) = (/ 'par%nsetbathy       ' /)
 case(  295)
   t%i2  => s%breaking
   t%rank =  2
   t%type = 'i'
   t%name = 'breaking'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'indicator whether cell has breaking nonh waves'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  296)
   t%r2  => s%fw
   t%rank =  2
   t%type = 'r'
   t%name = 'fw'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'wave friction coefficient'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  297)
   t%r1  => s%sigz
   t%rank =  1
   t%type = 'r'
   t%name = 'sigz'
   t%btype = 'b'
   t%units= '-'
   t%standardname= ''
   t%description= 'vertical distribution of sigma layers Q3D'
   t%dimensions(1:1) = (/ 'par%nz              ' /)
 case(  298)
   t%r3  => s%uz
   t%rank =  3
   t%type = 'r'
   t%name = 'uz'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'velocity (Q3D) ksi-comp'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%nz              ' /)
 case(  299)
   t%r3  => s%vz
   t%rank =  3
   t%type = 'r'
   t%name = 'vz'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'velocity (Q3D) eta-comp'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%nz              ' /)
 case(  300)
   t%r3  => s%ustz
   t%rank =  3
   t%type = 'r'
   t%name = 'ustz'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'stokes velocity (Q3D)'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%nz              ' /)
 case(  301)
   t%r3  => s%nutz
   t%rank =  3
   t%type = 'r'
   t%name = 'nutz'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'turbulence viscosity'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%nz              ' /)
 case(  302)
   t%r2  => s%dzs0dn
   t%rank =  2
   t%type = 'r'
   t%name = 'dzs0dn'
   t%btype = 'd'
   t%units= '-'
   t%standardname= ''
   t%description= 'alongshore water level gradient due to tide alone'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  303)
   t%r3  => s%ccz
   t%rank =  3
   t%type = 'r'
   t%name = 'ccz'
   t%btype = 'd'
   t%units= 'm3/m3'
   t%standardname= ''
   t%description= 'concentration profile'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%nz              ' /)
 case(  304)
   t%r2  => s%refA
   t%rank =  2
   t%type = 'r'
   t%name = 'refA'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'reference level'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  305)
   t%r2  => s%ca
   t%rank =  2
   t%type = 'r'
   t%name = 'ca'
   t%btype = 'd'
   t%units= 'm3/m3'
   t%standardname= ''
   t%description= 'reference concentration'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  306)
   t%r2  => s%zs1
   t%rank =  2
   t%type = 'r'
   t%name = 'zs1'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= 'sea_surface_height_above_sea_level_without_tide'
   t%description= 'water level minus tide'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  307)
   t%r2  => s%taubx_add
   t%rank =  2
   t%type = 'r'
   t%name = 'taubx_add'
   t%btype = 'd'
   t%units= 'N/m2'
   t%standardname= ''
   t%description= 'additional bed shear stress due to boundary layer effects, x-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  308)
   t%r2  => s%tauby_add
   t%rank =  2
   t%type = 'r'
   t%name = 'tauby_add'
   t%btype = 'd'
   t%units= 'N/m2'
   t%standardname= ''
   t%description= 'additional bed shear stress due to boundary layer effects, y-component'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  309)
   t%r2  => s%hhw
   t%rank =  2
   t%type = 'r'
   t%name = 'hhw'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'water depth used in all wave computations, includes H*par%delta'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  310)
   t%r2  => s%hhws
   t%rank =  2
   t%type = 'r'
   t%name = 'hhws'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'water depth used in wave stationary computation (and single_dir wave directions)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  311)
   t%r2  => s%uws
   t%rank =  2
   t%type = 'r'
   t%name = 'uws'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'u-velocity used in wave stationary computation (and single_dir wave directions)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  312)
   t%r2  => s%vws
   t%rank =  2
   t%type = 'r'
   t%name = 'vws'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'v-velocity used in wave stationary computation (and single_dir wave directions)'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  313)
   t%r2  => s%hhwcins
   t%rank =  2
   t%type = 'r'
   t%name = 'hhwcins'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'water depth used in wave instationary computation in case of wci'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  314)
   t%r2  => s%uwcins
   t%rank =  2
   t%type = 'r'
   t%name = 'uwcins'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'u-velocity used in wave stationary computation in case of wci'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  315)
   t%r2  => s%vwcins
   t%rank =  2
   t%type = 'r'
   t%name = 'vwcins'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'v-velocity used in wave stationary computation in case of wci'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  316)
   t%r2  => s%ucan
   t%rank =  2
   t%type = 'r'
   t%name = 'ucan'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'in-canopy velocity in x-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  317)
   t%r2  => s%vcan
   t%rank =  2
   t%type = 'r'
   t%name = 'vcan'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'in-canopy velocity in y-direction'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  318)
   t%r3  => s%rainfallinput
   t%rank =  3
   t%type = 'r'
   t%name = 'rainfallinput'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'prescribed rainfall rate'
   t%dimensions(1:3) = (/ 's%nx+1              ','s%ny+1              ','par%nrainfallrate   ' /)
 case(  319)
   t%r1  => s%trainfallinput
   t%rank =  1
   t%type = 'r'
   t%name = 'trainfallinput'
   t%btype = 'b'
   t%units= 's'
   t%standardname= ''
   t%description= 'points in time of prescibed rainfall'
   t%dimensions(1:1) = (/ 'par%nrainfallrate   ' /)
 case(  320)
   t%r2  => s%rainfallrate
   t%rank =  2
   t%type = 'r'
   t%name = 'rainfallrate'
   t%btype = 'd'
   t%units= 'm/s'
   t%standardname= ''
   t%description= 'rainfall rate at current time'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)
 case(  321)
   t%r2  => s%hstokes
   t%rank =  2
   t%type = 'r'
   t%name = 'hstokes'
   t%btype = 'd'
   t%units= 'm'
   t%standardname= ''
   t%description= 'water depth used for Stokes drift, turbulence and wave forces'
   t%dimensions(1:2) = (/ 's%nx+1              ','s%ny+1              ' /)

!directions for vi vim: filetype=fortran : syntax=fortran


! vim: filetype=fortran