!  DO NOT EDIT THIS FILE
!  But edit spaceparams.tmpl and/or makeincludes.F90
!  compile makeincludes.F90 and run the compilate
!  Compiling and running is taken care of by the Makefile
        case(  1)
          pointr2  => 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 comp. grid'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(  2)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr0  => 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)
          pointr0  => 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)
          pointr0  => 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 coordinates'
        case( 28)
          pointr0  => 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)
          pointr0  => 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)
          pointr0  => 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)
          pointr0  => 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)
          pointr0  => 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)
          pointi0  => 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)
          pointi0  => 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)
          pointr0  => 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)
          pointr0  => 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)
          pointr0  => 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)
          pointr0  => 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)
          pointr1  => 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)
          pointr1  => 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)
          pointr1  => 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)
          pointr1  => 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)
          pointi0  => 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)
          pointi0  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr1  => 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)
          pointr1  => 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)
          pointi0  => 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)
          pointi0  => 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 for wave_directions only'
        case( 51)
          pointr0  => 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)
          pointr0  => 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 for wave_directions only'
        case( 53)
          pointr0  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr3  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr1  => 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)
          pointr2  => 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)
          pointr1  => 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)
          pointr1  => 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)
          pointr1  => s%winddirts
          t%rank =   1
          t%type = 'r'
          t%name = 'winddirts'
          t%btype = 'b'
          t%units= 'deg'
          t%standardname= 'nautical'
          t%description= '[]       input wind direction'
          t%dimensions(1:  1) = (/ 's%windlen           '/)
        case( 93)
          pointr1  => 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)
          pointr1  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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)
          pointr2  => 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(134)
          pointi2  => 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(135)
          pointi2  => 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(136)
          pointi2  => 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(137)
          pointi2  => 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(138)
          pointr2  => 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(139)
          pointr2  => 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(140)
          pointr2  => 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(141)
          pointr2  => 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(142)
          pointr2  => 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(143)
          pointr3  => 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(144)
          pointr2  => 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(145)
          pointr2  => 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(146)
          pointr2  => 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(147)
          pointr2  => 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(148)
          pointr1  => 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(149)
          pointr2  => 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(150)
          pointr2  => 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(151)
          pointr2  => 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(152)
          pointr2  => 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(153)
          pointi0  => 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(154)
          pointr1  => 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, Lodewijk'
          t%dimensions(1:  1) = (/ 'par%ngd             '/)
        case(155)
          pointr1  => 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(156)
          pointr1  => 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(157)
          pointr1  => 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(158)
          pointr1  => 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(159)
          pointr3  => 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(160)
          pointr3  => 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(161)
          pointr3  => 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(162)
          pointr3  => 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(163)
          pointr3  => 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(164)
          pointr3  => 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(165)
          pointr3  => 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(166)
          pointr2  => 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(167)
          pointr2  => 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(168)
          pointr2  => s%kb
          t%rank =   2
          t%type = 'r'
          t%name = 'kb'
          t%btype = 'd'
          t%units= 'm^2/s^2'
          t%standardname= ''
          t%description= 'near bed turbulence intensity due to depth induces breaking'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(169)
          pointr2  => 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(170)
          pointr2  => 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(171)
          pointr2  => 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(172)
          pointr2  => 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(173)
          pointr2  => 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(174)
          pointr2  => 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(175)
          pointr2  => 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(176)
          pointr2  => 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(177)
          pointr2  => 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(178)
          pointr2  => 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(179)
          pointr2  => 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(180)
          pointr2  => 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(181)
          pointr2  => 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(182)
          pointr2  => 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(183)
          pointr2  => 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(184)
          pointr2  => 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(185)
          pointr2  => 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(186)
          pointr2  => 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(187)
          pointr2  => 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(188)
          pointr2  => 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(189)
          pointr2  => 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(190)
          pointr2  => s%kturb
          t%rank =   2
          t%type = 'r'
          t%name = 'kturb'
          t%btype = 'd'
          t%units= 'm^2/s^2'
          t%standardname= ''
          t%description= 'depth averaged turbulence intensity due to long wave breaking'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(191)
          pointr3  => 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(192)
          pointr3  => 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(193)
          pointr3  => 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(194)
          pointi2  => 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(195)
          pointr4  => s%pbbed
          t%rank =   4
          t%type = 'r'
          t%name = 'pbbed'
          t%btype = 'd'
          t%units= '-'
          t%standardname= ''
          t%description= 'Concentration of every sediment type, per bed sediment layer, per grid cell'
          t%dimensions(1:  4) = (/ 's%nx+1              ','s%ny+1              ','par%nd              ','par%ngd             '/)
        case(196)
          pointr3  => s%dzbed
          t%rank =   3
          t%type = 'r'
          t%name = 'dzbed'
          t%btype = 'd'
          t%units= '-'
          t%standardname= ''
          t%description= 'Vertical dimension of each bed sediment layer (nd), per grid cell (nx,ny)'
          t%dimensions(1:  3) = (/ 's%nx+1              ','s%ny+1              ','par%nd              '/)
        case(197)
          pointr2  => s%z0bed
          t%rank =   2
          t%type = 'r'
          t%name = 'z0bed'
          t%btype = 'd'
          t%units= '-'
          t%standardname= ''
          t%description= 'Reference level for the bed, located at the initial bottom of the bed sediment layers'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(198)
          pointr2  => 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(199)
          pointr2  => 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(200)
          pointr2  => 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(201)
          pointr2  => 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(202)
          pointr2  => 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(203)
          pointr2  => 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(204)
          pointr2  => 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(205)
          pointr2  => 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(206)
          pointr2  => s%pres
          t%rank =   2
          t%type = 'r'
          t%name = 'pres'
          t%btype = 'd'
          t%units= 'm^2/s^2'
          t%standardname= ''
          t%description= 'normalized dynamic pressure'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(207)
          pointr2  => s%dU
          t%rank =   2
          t%type = 'r'
          t%name = 'dU'
          t%btype = 'd'
          t%units= 'm^2/s^2'
          t%standardname= ''
          t%description= 'The u-velocity difference between two vertical layars (reduced 2-Layer non-hydrostatic model)'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(208)
          pointr2  => s%dV
          t%rank =   2
          t%type = 'r'
          t%name = 'dV'
          t%btype = 'd'
          t%units= 'm^2/s^2'
          t%standardname= ''
          t%description= 'The v-velocity difference between two vertical layars (reduced 2-Layer non-hydrostatic model)'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(209)
          pointr2  => 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(210)
          pointr2  => 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(211)
          pointr2  => 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(212)
          pointr2  => 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(213)
          pointr2  => s%taubx
          t%rank =   2
          t%type = 'r'
          t%name = 'taubx'
          t%btype = 'd'
          t%units= 'N/m^2'
          t%standardname= ''
          t%description= 'bed shear stress, x-component'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(214)
          pointr2  => s%tauby
          t%rank =   2
          t%type = 'r'
          t%name = 'tauby'
          t%btype = 'd'
          t%units= 'N/m^2'
          t%standardname= ''
          t%description= 'bed shear stress, y-component'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(215)
          pointr2  => s%Df
          t%rank =   2
          t%type = 'r'
          t%name = 'Df'
          t%btype = 'd'
          t%units= 'W/m^2'
          t%standardname= ''
          t%description= 'dissipation rate due to bed friction'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(216)
          pointr2  => s%Dp
          t%rank =   2
          t%type = 'r'
          t%name = 'Dp'
          t%btype = 'd'
          t%units= 'W/m^2'
          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(217)
          pointr2  => 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(218)
          pointr2  => 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(219)
          pointr2  => 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(220)
          pointr2  => 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(221)
          pointr2  => 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(222)
          pointr2  => 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(223)
          pointr2  => 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(224)
          pointr2  => s%cf
          t%rank =   2
          t%type = 'r'
          t%name = 'cf'
          t%btype = 'd'
          t%units= '-'
          t%standardname= ''
          t%description= 'Dimensionless friction coefficient flow'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(225)
          pointr2  => s%cfu
          t%rank =   2
          t%type = 'r'
          t%name = 'cfu'
          t%btype = 'd'
          t%units= '-'
          t%standardname= ''
          t%description= 'Dimensionless friction coefficient flow in u-points'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(226)
          pointr2  => s%cfv
          t%rank =   2
          t%type = 'r'
          t%name = 'cfv'
          t%btype = 'd'
          t%units= '-'
          t%standardname= ''
          t%description= 'Dimensionless friction coefficient flow in v-points'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(227)
          pointr2  => 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(228)
          pointr2  => 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(229)
          pointr2  => 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(230)
          pointr3  => 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(231)
          pointr1  => 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(232)
          pointi2  => s%pdisch
          t%rank =   2
          t%type = 'i'
          t%name = 'pdisch'
          t%btype = 'b'
          t%units= '-'
          t%standardname= ''
          t%description= 'Discharge locations'
          t%dimensions(1:  2) = (/ 'par%ndischarge      ','4                   '/)
        case(233)
          pointi1  => 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(234)
          pointr2  => s%qdisch
          t%rank =   2
          t%type = 'r'
          t%name = 'qdisch'
          t%btype = 'b'
          t%units= 'm^2/s'
          t%standardname= ''
          t%description= 'Discharges'
          t%dimensions(1:  2) = (/ 'par%ntdischarge     ','par%ndischarge      '/)
        case(235)
          pointr1  => 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(236)
          pointr1  => 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(237)
          pointr1  => 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(238)
          pointr1  => 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(239)
          pointr1  => 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(240)
          pointr1  => 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(241)
          pointr1  => 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(242)
          pointr1  => 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(243)
          pointr1  => 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(244)
          pointr1  => 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(245)
          pointr2  => 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(246)
          pointr2  => 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(247)
          pointi0  => 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(248)
          pointr2  => s%dobs
          t%rank =   2
          t%type = 'r'
          t%name = 'dobs'
          t%btype = 'd'
          t%units= 'W/m2'
          t%standardname= 'beachwizard'
          t%description= 'observed dissipation'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(249)
          pointr2  => s%sig2prior
          t%rank =   2
          t%type = 'r'
          t%name = 'sig2prior'
          t%btype = 'd'
          t%units= 'm2'
          t%standardname= 'beachwizard'
          t%description= 'prior std squared'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(250)
          pointr2  => s%zbobs
          t%rank =   2
          t%type = 'r'
          t%name = 'zbobs'
          t%btype = 'd'
          t%units= 'm'
          t%standardname= 'beachwizard'
          t%description= 'observed depth'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(251)
          pointr2  => s%shobs
          t%rank =   2
          t%type = 'r'
          t%name = 'shobs'
          t%btype = 'd'
          t%units= 'm'
          t%standardname= 'beachwizard'
          t%description= 'observed shoreline'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(252)
          pointr2  => s%bwalpha
          t%rank =   2
          t%type = 'r'
          t%name = 'bwalpha'
          t%btype = 'd'
          t%units= '-'
          t%standardname= 'beachwizard'
          t%description= 'weighting factor'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(253)
          pointr2  => s%dcmdo
          t%rank =   2
          t%type = 'r'
          t%name = 'dcmdo'
          t%btype = 'd'
          t%units= 'W/m2'
          t%standardname= 'beachwizard'
          t%description= 'computed minus observed dissipation'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(254)
          pointr2  => s%dassim
          t%rank =   2
          t%type = 'r'
          t%name = 'dassim'
          t%btype = 'd'
          t%units= 'm'
          t%standardname= 'beachwizard'
          t%description= 'depth change'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(255)
          pointr2  => s%cobs
          t%rank =   2
          t%type = 'r'
          t%name = 'cobs'
          t%btype = 'd'
          t%units= 'm/s'
          t%standardname= 'beachwizard'
          t%description= 'observed wave celerity'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(256)
          pointr1  => 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(257)
          pointr1  => 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(258)
          pointr1  => 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(259)
          pointr1  => s%shipFx
          t%rank =   1
          t%type = 'r'
          t%name = 'shipFx'
          t%btype = 'b'
          t%units= 'm'
          t%standardname= ''
          t%description= 'force on ship in x-direction'
          t%dimensions(1:  1) = (/ 'par%nship           '/)
        case(260)
          pointr1  => s%shipFy
          t%rank =   1
          t%type = 'r'
          t%name = 'shipFy'
          t%btype = 'b'
          t%units= 'm'
          t%standardname= ''
          t%description= 'force on ship in y-direction'
          t%dimensions(1:  1) = (/ 'par%nship           '/)
        case(261)
          pointr1  => s%shipFz
          t%rank =   1
          t%type = 'r'
          t%name = 'shipFz'
          t%btype = 'b'
          t%units= 'm'
          t%standardname= ''
          t%description= 'force on ship in z-direction'
          t%dimensions(1:  1) = (/ 'par%nship           '/)
        case(262)
          pointr1  => s%shipMx
          t%rank =   1
          t%type = 'r'
          t%name = 'shipMx'
          t%btype = 'b'
          t%units= 'm'
          t%standardname= ''
          t%description= 'moment on ship around x-axis'
          t%dimensions(1:  1) = (/ 'par%nship           '/)
        case(263)
          pointr1  => s%shipMy
          t%rank =   1
          t%type = 'r'
          t%name = 'shipMy'
          t%btype = 'b'
          t%units= 'm'
          t%standardname= ''
          t%description= 'moment on ship around y-axis'
          t%dimensions(1:  1) = (/ 'par%nship           '/)
        case(264)
          pointr1  => s%shipMz
          t%rank =   1
          t%type = 'r'
          t%name = 'shipMz'
          t%btype = 'b'
          t%units= 'm'
          t%standardname= ''
          t%description= 'moment on ship around z-axis'
          t%dimensions(1:  1) = (/ 'par%nship           '/)
        case(265)
          pointr1  => s%shipphi
          t%rank =   1
          t%type = 'r'
          t%name = 'shipphi'
          t%btype = 'b'
          t%units= 'm'
          t%standardname= ''
          t%description= 'turning angle arround x-axis'
          t%dimensions(1:  1) = (/ 'par%nship           '/)
        case(266)
          pointr1  => s%shipchi
          t%rank =   1
          t%type = 'r'
          t%name = 'shipchi'
          t%btype = 'b'
          t%units= 'm'
          t%standardname= ''
          t%description= 'turning angle arround y-axis'
          t%dimensions(1:  1) = (/ 'par%nship           '/)
        case(267)
          pointr1  => s%shippsi
          t%rank =   1
          t%type = 'r'
          t%name = 'shippsi'
          t%btype = 'b'
          t%units= 'm'
          t%standardname= ''
          t%description= 'turning angle arround z-axis'
          t%dimensions(1:  1) = (/ 'par%nship           '/)
        case(268)
          pointr2  => 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(269)
          pointr2  => 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(270)
          pointr2  => 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(271)
          pointr2  => s%ududx
          t%rank =   2
          t%type = 'r'
          t%name = 'ududx'
          t%btype = 'd'
          t%units= 'm^2/s^2'
          t%standardname= ''
          t%description= 'advection'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(272)
          pointr2  => s%vdvdy
          t%rank =   2
          t%type = 'r'
          t%name = 'vdvdy'
          t%btype = 'd'
          t%units= 'm^2/s^2'
          t%standardname= ''
          t%description= 'advection'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(273)
          pointr2  => s%udvdx
          t%rank =   2
          t%type = 'r'
          t%name = 'udvdx'
          t%btype = 'd'
          t%units= 'm^2/s^2'
          t%standardname= ''
          t%description= 'advection'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(274)
          pointr2  => s%vdudy
          t%rank =   2
          t%type = 'r'
          t%name = 'vdudy'
          t%btype = 'd'
          t%units= 'm^2/s^2'
          t%standardname= ''
          t%description= 'advection'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(275)
          pointr2  => s%viscu
          t%rank =   2
          t%type = 'r'
          t%name = 'viscu'
          t%btype = 'd'
          t%units= 'm^2/s^2'
          t%standardname= ''
          t%description= 'viscosity'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(276)
          pointr2  => s%viscv
          t%rank =   2
          t%type = 'r'
          t%name = 'viscv'
          t%btype = 'd'
          t%units= 'm^2/s^2'
          t%standardname= ''
          t%description= 'viscosity'
          t%dimensions(1:  2) = (/ 's%nx+1              ','s%ny+1              '/)
        case(277)
          pointr3  => 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(278)
          pointr1  => 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(279)
          pointi2  => 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              '/)
!directions for vi vim: filetype=fortran : syntax=fortran
