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



case(  1)
   if (associated(s%x)) deallocate(s%x)
   allocate(s%x(s%nx+1,s%ny+1))
case(  2)
   if (associated(s%y)) deallocate(s%y)
   allocate(s%y(s%nx+1,s%ny+1))
case(  3)
   if (associated(s%xz)) deallocate(s%xz)
   allocate(s%xz(s%nx+1,s%ny+1))
case(  4)
   if (associated(s%yz)) deallocate(s%yz)
   allocate(s%yz(s%nx+1,s%ny+1))
case(  5)
   if (associated(s%xu)) deallocate(s%xu)
   allocate(s%xu(s%nx+1,s%ny+1))
case(  6)
   if (associated(s%yu)) deallocate(s%yu)
   allocate(s%yu(s%nx+1,s%ny+1))
case(  7)
   if (associated(s%xv)) deallocate(s%xv)
   allocate(s%xv(s%nx+1,s%ny+1))
case(  8)
   if (associated(s%yv)) deallocate(s%yv)
   allocate(s%yv(s%nx+1,s%ny+1))
case(  9)
   if (associated(s%dsu)) deallocate(s%dsu)
   allocate(s%dsu(s%nx+1,s%ny+1))
case(  10)
   if (associated(s%dsv)) deallocate(s%dsv)
   allocate(s%dsv(s%nx+1,s%ny+1))
case(  11)
   if (associated(s%dsz)) deallocate(s%dsz)
   allocate(s%dsz(s%nx+1,s%ny+1))
case(  12)
   if (associated(s%dsc)) deallocate(s%dsc)
   allocate(s%dsc(s%nx+1,s%ny+1))
case(  13)
   if (associated(s%dnu)) deallocate(s%dnu)
   allocate(s%dnu(s%nx+1,s%ny+1))
case(  14)
   if (associated(s%dnv)) deallocate(s%dnv)
   allocate(s%dnv(s%nx+1,s%ny+1))
case(  15)
   if (associated(s%dnz)) deallocate(s%dnz)
   allocate(s%dnz(s%nx+1,s%ny+1))
case(  16)
   if (associated(s%dnc)) deallocate(s%dnc)
   allocate(s%dnc(s%nx+1,s%ny+1))
case(  17)
   if (associated(s%dsdnui)) deallocate(s%dsdnui)
   allocate(s%dsdnui(s%nx+1,s%ny+1))
case(  18)
   if (associated(s%dsdnvi)) deallocate(s%dsdnvi)
   allocate(s%dsdnvi(s%nx+1,s%ny+1))
case(  19)
   if (associated(s%dsdnzi)) deallocate(s%dsdnzi)
   allocate(s%dsdnzi(s%nx+1,s%ny+1))
case(  20)
   if (associated(s%alfaz)) deallocate(s%alfaz)
   allocate(s%alfaz(s%nx+1,s%ny+1))
case(  21)
   if (associated(s%alfau)) deallocate(s%alfau)
   allocate(s%alfau(s%nx+1,s%ny+1))
case(  22)
   if (associated(s%alfav)) deallocate(s%alfav)
   allocate(s%alfav(s%nx+1,s%ny+1))
case(  23)
   if (associated(s%sdist)) deallocate(s%sdist)
   allocate(s%sdist(s%nx+1,s%ny+1))
case(  24)
   if (associated(s%ndist)) deallocate(s%ndist)
   allocate(s%ndist(s%nx+1,s%ny+1))
case(  25)
case(  26)
case(  27)
case(  28)
case(  29)
case(  30)
case(  31)
case(  32)
case(  33)
case(  34)
case(  35)
case(  36)
case(  37)
case(  38)
case(  39)
   if (associated(s%xyzs01)) deallocate(s%xyzs01)
   allocate(s%xyzs01(2))
case(  40)
   if (associated(s%xyzs02)) deallocate(s%xyzs02)
   allocate(s%xyzs02(2))
case(  41)
   if (associated(s%xyzs03)) deallocate(s%xyzs03)
   allocate(s%xyzs03(2))
case(  42)
   if (associated(s%xyzs04)) deallocate(s%xyzs04)
   allocate(s%xyzs04(2))
case(  43)
case(  44)
case(  45)
   if (associated(s%zb)) deallocate(s%zb)
   allocate(s%zb(s%nx+1,s%ny+1))
case(  46)
   if (associated(s%zb0)) deallocate(s%zb0)
   allocate(s%zb0(s%nx+1,s%ny+1))
case(  47)
   if (associated(s%theta)) deallocate(s%theta)
   allocate(s%theta(s%ntheta))
case(  48)
   if (associated(s%theta_s)) deallocate(s%theta_s)
   allocate(s%theta_s(s%ntheta_s))
case(  49)
case(  50)
case(  51)
case(  52)
case(  53)
case(  54)
   if (associated(s%thetamean)) deallocate(s%thetamean)
   allocate(s%thetamean(s%nx+1,s%ny+1))
case(  55)
   if (associated(s%Fx)) deallocate(s%Fx)
   allocate(s%Fx(s%nx+1,s%ny+1))
case(  56)
   if (associated(s%Fy)) deallocate(s%Fy)
   allocate(s%Fy(s%nx+1,s%ny+1))
case(  57)
   if (associated(s%Sxy)) deallocate(s%Sxy)
   allocate(s%Sxy(s%nx+1,s%ny+1))
case(  58)
   if (associated(s%Syy)) deallocate(s%Syy)
   allocate(s%Syy(s%nx+1,s%ny+1))
case(  59)
   if (associated(s%Sxx)) deallocate(s%Sxx)
   allocate(s%Sxx(s%nx+1,s%ny+1))
case(  60)
   if (associated(s%n)) deallocate(s%n)
   allocate(s%n(s%nx+1,s%ny+1))
case(  61)
   if (associated(s%H)) deallocate(s%H)
   allocate(s%H(s%nx+1,s%ny+1))
case(  62)
   if (associated(s%cgx)) deallocate(s%cgx)
   allocate(s%cgx(s%nx+1,s%ny+1,s%ntheta))
case(  63)
   if (associated(s%cgy)) deallocate(s%cgy)
   allocate(s%cgy(s%nx+1,s%ny+1,s%ntheta))
case(  64)
   if (associated(s%cx)) deallocate(s%cx)
   allocate(s%cx(s%nx+1,s%ny+1,s%ntheta))
case(  65)
   if (associated(s%cy)) deallocate(s%cy)
   allocate(s%cy(s%nx+1,s%ny+1,s%ntheta))
case(  66)
   if (associated(s%ctheta)) deallocate(s%ctheta)
   allocate(s%ctheta(s%nx+1,s%ny+1,s%ntheta))
case(  67)
   if (associated(s%ee)) deallocate(s%ee)
   allocate(s%ee(s%nx+1,s%ny+1,s%ntheta))
case(  68)
   if (associated(s%thet)) deallocate(s%thet)
   allocate(s%thet(s%nx+1,s%ny+1,s%ntheta))
case(  69)
   if (associated(s%costh)) deallocate(s%costh)
   allocate(s%costh(s%nx+1,s%ny+1,s%ntheta))
case(  70)
   if (associated(s%sinth)) deallocate(s%sinth)
   allocate(s%sinth(s%nx+1,s%ny+1,s%ntheta))
case(  71)
   if (associated(s%sigt)) deallocate(s%sigt)
   allocate(s%sigt(s%nx+1,s%ny+1,s%ntheta))
case(  72)
   if (associated(s%rr)) deallocate(s%rr)
   allocate(s%rr(s%nx+1,s%ny+1,s%ntheta))
case(  73)
   if (associated(s%cgx_s)) deallocate(s%cgx_s)
   allocate(s%cgx_s(s%nx+1,s%ny+1,s%ntheta_s))
case(  74)
   if (associated(s%cgy_s)) deallocate(s%cgy_s)
   allocate(s%cgy_s(s%nx+1,s%ny+1,s%ntheta_s))
case(  75)
   if (associated(s%ctheta_s)) deallocate(s%ctheta_s)
   allocate(s%ctheta_s(s%nx+1,s%ny+1,s%ntheta_s))
case(  76)
   if (associated(s%ee_s)) deallocate(s%ee_s)
   allocate(s%ee_s(s%nx+1,s%ny+1,s%ntheta_s))
case(  77)
   if (associated(s%thet_s)) deallocate(s%thet_s)
   allocate(s%thet_s(s%nx+1,s%ny+1,s%ntheta_s))
case(  78)
   if (associated(s%costh_s)) deallocate(s%costh_s)
   allocate(s%costh_s(s%nx+1,s%ny+1,s%ntheta_s))
case(  79)
   if (associated(s%sinth_s)) deallocate(s%sinth_s)
   allocate(s%sinth_s(s%nx+1,s%ny+1,s%ntheta_s))
case(  80)
   if (associated(s%k)) deallocate(s%k)
   allocate(s%k(s%nx+1,s%ny+1))
case(  81)
   if (associated(s%c)) deallocate(s%c)
   allocate(s%c(s%nx+1,s%ny+1))
case(  82)
   if (associated(s%cg)) deallocate(s%cg)
   allocate(s%cg(s%nx+1,s%ny+1))
case(  83)
   if (associated(s%sigm)) deallocate(s%sigm)
   allocate(s%sigm(s%nx+1,s%ny+1))
case(  84)
   if (associated(s%wm)) deallocate(s%wm)
   allocate(s%wm(s%nx+1,s%ny+1))
case(  85)
   if (associated(s%hh)) deallocate(s%hh)
   allocate(s%hh(s%nx+1,s%ny+1))
case(  86)
   if (associated(s%zs)) deallocate(s%zs)
   allocate(s%zs(s%nx+1,s%ny+1))
case(  87)
   if (associated(s%zs0)) deallocate(s%zs0)
   allocate(s%zs0(s%nx+1,s%ny+1))
case(  88)
   if (associated(s%tideinpt)) deallocate(s%tideinpt)
   allocate(s%tideinpt(s%tidelen))
case(  89)
   if (associated(s%tideinpz)) deallocate(s%tideinpz)
   allocate(s%tideinpz(s%tidelen,par%tideloc))
case(  90)
   if (associated(s%windinpt)) deallocate(s%windinpt)
   allocate(s%windinpt(s%windlen))
case(  91)
   if (associated(s%windvelts)) deallocate(s%windvelts)
   allocate(s%windvelts(s%windlen))
case(  92)
   if (associated(s%winddirts)) deallocate(s%winddirts)
   allocate(s%winddirts(s%windlen))
case(  93)
   if (associated(s%windxts)) deallocate(s%windxts)
   allocate(s%windxts(s%windlen))
case(  94)
   if (associated(s%windyts)) deallocate(s%windyts)
   allocate(s%windyts(s%windlen))
case(  95)
   if (associated(s%windsu)) deallocate(s%windsu)
   allocate(s%windsu(s%nx+1,s%ny+1))
case(  96)
   if (associated(s%windnv)) deallocate(s%windnv)
   allocate(s%windnv(s%nx+1,s%ny+1))
case(  97)
   if (associated(s%dzsdt)) deallocate(s%dzsdt)
   allocate(s%dzsdt(s%nx+1,s%ny+1))
case(  98)
   if (associated(s%dzsdx)) deallocate(s%dzsdx)
   allocate(s%dzsdx(s%nx+1,s%ny+1))
case(  99)
   if (associated(s%dzsdy)) deallocate(s%dzsdy)
   allocate(s%dzsdy(s%nx+1,s%ny+1))
case(  100)
   if (associated(s%dzbdx)) deallocate(s%dzbdx)
   allocate(s%dzbdx(s%nx+1,s%ny+1))
case(  101)
   if (associated(s%dzbdy)) deallocate(s%dzbdy)
   allocate(s%dzbdy(s%nx+1,s%ny+1))
case(  102)
   if (associated(s%dzbdt)) deallocate(s%dzbdt)
   allocate(s%dzbdt(s%nx+1,s%ny+1))
case(  103)
   if (associated(s%dzbnow)) deallocate(s%dzbnow)
   allocate(s%dzbnow(s%nx+1,s%ny+1))
case(  104)
   if (associated(s%uu)) deallocate(s%uu)
   allocate(s%uu(s%nx+1,s%ny+1))
case(  105)
   if (associated(s%vu)) deallocate(s%vu)
   allocate(s%vu(s%nx+1,s%ny+1))
case(  106)
   if (associated(s%uv)) deallocate(s%uv)
   allocate(s%uv(s%nx+1,s%ny+1))
case(  107)
   if (associated(s%vv)) deallocate(s%vv)
   allocate(s%vv(s%nx+1,s%ny+1))
case(  108)
   if (associated(s%qx)) deallocate(s%qx)
   allocate(s%qx(s%nx+1,s%ny+1))
case(  109)
   if (associated(s%qy)) deallocate(s%qy)
   allocate(s%qy(s%nx+1,s%ny+1))
case(  110)
   if (associated(s%sedero)) deallocate(s%sedero)
   allocate(s%sedero(s%nx+1,s%ny+1))
case(  111)
   if (associated(s%dcbdx)) deallocate(s%dcbdx)
   allocate(s%dcbdx(s%nx+1,s%ny+1))
case(  112)
   if (associated(s%dcbdy)) deallocate(s%dcbdy)
   allocate(s%dcbdy(s%nx+1,s%ny+1))
case(  113)
   if (associated(s%dcsdx)) deallocate(s%dcsdx)
   allocate(s%dcsdx(s%nx+1,s%ny+1))
case(  114)
   if (associated(s%dcsdy)) deallocate(s%dcsdy)
   allocate(s%dcsdy(s%nx+1,s%ny+1))
case(  115)
   if (associated(s%ui)) deallocate(s%ui)
   allocate(s%ui(s%nx+1,s%ny+1))
case(  116)
   if (associated(s%vi)) deallocate(s%vi)
   allocate(s%vi(s%nx+1,s%ny+1))
case(  117)
   if (associated(s%E)) deallocate(s%E)
   allocate(s%E(s%nx+1,s%ny+1))
case(  118)
   if (associated(s%R)) deallocate(s%R)
   allocate(s%R(s%nx+1,s%ny+1))
case(  119)
   if (associated(s%urms)) deallocate(s%urms)
   allocate(s%urms(s%nx+1,s%ny+1))
case(  120)
   if (associated(s%D)) deallocate(s%D)
   allocate(s%D(s%nx+1,s%ny+1))
case(  121)
   if (associated(s%Qb)) deallocate(s%Qb)
   allocate(s%Qb(s%nx+1,s%ny+1))
case(  122)
   if (associated(s%ust)) deallocate(s%ust)
   allocate(s%ust(s%nx+1,s%ny+1))
case(  123)
   if (associated(s%ueu)) deallocate(s%ueu)
   allocate(s%ueu(s%nx+1,s%ny+1))
case(  124)
   if (associated(s%vev)) deallocate(s%vev)
   allocate(s%vev(s%nx+1,s%ny+1))
case(  125)
   if (associated(s%vmagu)) deallocate(s%vmagu)
   allocate(s%vmagu(s%nx+1,s%ny+1))
case(  126)
   if (associated(s%vmageu)) deallocate(s%vmageu)
   allocate(s%vmageu(s%nx+1,s%ny+1))
case(  127)
   if (associated(s%vmagv)) deallocate(s%vmagv)
   allocate(s%vmagv(s%nx+1,s%ny+1))
case(  128)
   if (associated(s%vmagev)) deallocate(s%vmagev)
   allocate(s%vmagev(s%nx+1,s%ny+1))
case(  129)
   if (associated(s%u)) deallocate(s%u)
   allocate(s%u(s%nx+1,s%ny+1))
case(  130)
   if (associated(s%v)) deallocate(s%v)
   allocate(s%v(s%nx+1,s%ny+1))
case(  131)
   if (associated(s%ue)) deallocate(s%ue)
   allocate(s%ue(s%nx+1,s%ny+1))
case(  132)
   if (associated(s%ve)) deallocate(s%ve)
   allocate(s%ve(s%nx+1,s%ny+1))
case(  133)
   if (associated(s%ue_sed)) deallocate(s%ue_sed)
   allocate(s%ue_sed(s%nx+1,s%ny+1))
case(  134)
   if (associated(s%ve_sed)) deallocate(s%ve_sed)
   allocate(s%ve_sed(s%nx+1,s%ny+1))
case(  135)
   if (associated(s%hold)) deallocate(s%hold)
   allocate(s%hold(s%nx+1,s%ny+1))
case(  136)
   if (associated(s%wetu)) deallocate(s%wetu)
   allocate(s%wetu(s%nx+1,s%ny+1))
case(  137)
   if (associated(s%wetv)) deallocate(s%wetv)
   allocate(s%wetv(s%nx+1,s%ny+1))
case(  138)
   if (associated(s%wetz)) deallocate(s%wetz)
   allocate(s%wetz(s%nx+1,s%ny+1))
case(  139)
   if (associated(s%wete)) deallocate(s%wete)
   allocate(s%wete(s%nx+1,s%ny+1))
case(  140)
   if (associated(s%hu)) deallocate(s%hu)
   allocate(s%hu(s%nx+1,s%ny+1))
case(  141)
   if (associated(s%hv)) deallocate(s%hv)
   allocate(s%hv(s%nx+1,s%ny+1))
case(  142)
   if (associated(s%hum)) deallocate(s%hum)
   allocate(s%hum(s%nx+1,s%ny+1))
case(  143)
   if (associated(s%hvm)) deallocate(s%hvm)
   allocate(s%hvm(s%nx+1,s%ny+1))
case(  144)
   if (associated(s%vmag)) deallocate(s%vmag)
   allocate(s%vmag(s%nx+1,s%ny+1))
case(  145)
   if (associated(s%ccg)) deallocate(s%ccg)
   allocate(s%ccg(s%nx+1,s%ny+1,par%ngd))
case(  146)
   if (associated(s%uwf)) deallocate(s%uwf)
   allocate(s%uwf(s%nx+1,s%ny+1))
case(  147)
   if (associated(s%vwf)) deallocate(s%vwf)
   allocate(s%vwf(s%nx+1,s%ny+1))
case(  148)
   if (associated(s%ustr)) deallocate(s%ustr)
   allocate(s%ustr(s%nx+1,s%ny+1))
case(  149)
   if (associated(s%usd)) deallocate(s%usd)
   allocate(s%usd(s%nx+1,s%ny+1))
case(  150)
   if (associated(s%bi)) deallocate(s%bi)
   allocate(s%bi(s%ny+1))
case(  151)
   if (associated(s%DR)) deallocate(s%DR)
   allocate(s%DR(s%nx+1,s%ny+1))
case(  152)
   if (associated(s%umean)) deallocate(s%umean)
   allocate(s%umean(s%nx+1,s%ny+1))
case(  153)
   if (associated(s%vmean)) deallocate(s%vmean)
   allocate(s%vmean(s%nx+1,s%ny+1))
case(  154)
   if (associated(s%ur)) deallocate(s%ur)
   allocate(s%ur(s%nx+1,s%ny+1))
case(  155)
case(  156)
   if (associated(s%D15)) deallocate(s%D15)
   allocate(s%D15(par%ngd))
case(  157)
   if (associated(s%D50)) deallocate(s%D50)
   allocate(s%D50(par%ngd))
case(  158)
   if (associated(s%D90)) deallocate(s%D90)
   allocate(s%D90(par%ngd))
case(  159)
   if (associated(s%sedcal)) deallocate(s%sedcal)
   allocate(s%sedcal(par%ngd))
case(  160)
   if (associated(s%ucrcal)) deallocate(s%ucrcal)
   allocate(s%ucrcal(par%ngd))
case(  161)
   if (associated(s%Tsg)) deallocate(s%Tsg)
   allocate(s%Tsg(s%nx+1,s%ny+1,par%ngd))
case(  162)
   if (associated(s%Susg)) deallocate(s%Susg)
   allocate(s%Susg(s%nx+1,s%ny+1,par%ngd))
case(  163)
   if (associated(s%Svsg)) deallocate(s%Svsg)
   allocate(s%Svsg(s%nx+1,s%ny+1,par%ngd))
case(  164)
   if (associated(s%Subg)) deallocate(s%Subg)
   allocate(s%Subg(s%nx+1,s%ny+1,par%ngd))
case(  165)
   if (associated(s%Svbg)) deallocate(s%Svbg)
   allocate(s%Svbg(s%nx+1,s%ny+1,par%ngd))
case(  166)
   if (associated(s%ceqbg)) deallocate(s%ceqbg)
   allocate(s%ceqbg(s%nx+1,s%ny+1,par%ngd))
case(  167)
   if (associated(s%ceqsg)) deallocate(s%ceqsg)
   allocate(s%ceqsg(s%nx+1,s%ny+1,par%ngd))
case(  168)
   if (associated(s%ua)) deallocate(s%ua)
   allocate(s%ua(s%nx+1,s%ny+1))
case(  169)
   if (associated(s%BR)) deallocate(s%BR)
   allocate(s%BR(s%nx+1,s%ny+1))
case(  170)
   if (associated(s%kb)) deallocate(s%kb)
   allocate(s%kb(s%nx+1,s%ny+1))
case(  171)
   if (associated(s%Tbore)) deallocate(s%Tbore)
   allocate(s%Tbore(s%nx+1,s%ny+1))
case(  172)
   if (associated(s%dzav)) deallocate(s%dzav)
   allocate(s%dzav(s%nx+1,s%ny+1))
case(  173)
   if (associated(s%maxzs)) deallocate(s%maxzs)
   allocate(s%maxzs(s%nx+1,s%ny+1))
case(  174)
   if (associated(s%minzs)) deallocate(s%minzs)
   allocate(s%minzs(s%nx+1,s%ny+1))
case(  175)
   if (associated(s%L1)) deallocate(s%L1)
   allocate(s%L1(s%nx+1,s%ny+1))
case(  176)
   if (associated(s%Sk)) deallocate(s%Sk)
   allocate(s%Sk(s%nx+1,s%ny+1))
case(  177)
   if (associated(s%As)) deallocate(s%As)
   allocate(s%As(s%nx+1,s%ny+1))
case(  178)
   if (associated(s%gwhead)) deallocate(s%gwhead)
   allocate(s%gwhead(s%nx+1,s%ny+1))
case(  179)
   if (associated(s%gwheadb)) deallocate(s%gwheadb)
   allocate(s%gwheadb(s%nx+1,s%ny+1))
case(  180)
   if (associated(s%gwlevel)) deallocate(s%gwlevel)
   allocate(s%gwlevel(s%nx+1,s%ny+1))
case(  181)
   if (associated(s%gwheight)) deallocate(s%gwheight)
   allocate(s%gwheight(s%nx+1,s%ny+1))
case(  182)
   if (associated(s%gwbottom)) deallocate(s%gwbottom)
   allocate(s%gwbottom(s%nx+1,s%ny+1))
case(  183)
   if (associated(s%gwu)) deallocate(s%gwu)
   allocate(s%gwu(s%nx+1,s%ny+1))
case(  184)
   if (associated(s%gwv)) deallocate(s%gwv)
   allocate(s%gwv(s%nx+1,s%ny+1))
case(  185)
   if (associated(s%gwqx)) deallocate(s%gwqx)
   allocate(s%gwqx(s%nx+1,s%ny+1))
case(  186)
   if (associated(s%gwqy)) deallocate(s%gwqy)
   allocate(s%gwqy(s%nx+1,s%ny+1))
case(  187)
   if (associated(s%gww)) deallocate(s%gww)
   allocate(s%gww(s%nx+1,s%ny+1))
case(  188)
   if (associated(s%gwcurv)) deallocate(s%gwcurv)
   allocate(s%gwcurv(s%nx+1,s%ny+1))
case(  189)
   if (associated(s%dinfil)) deallocate(s%dinfil)
   allocate(s%dinfil(s%nx+1,s%ny+1))
case(  190)
   if (associated(s%infil)) deallocate(s%infil)
   allocate(s%infil(s%nx+1,s%ny+1))
case(  191)
   if (associated(s%gw0back)) deallocate(s%gw0back)
   allocate(s%gw0back(2,s%ny+1))
case(  192)
   if (associated(s%kturb)) deallocate(s%kturb)
   allocate(s%kturb(s%nx+1,s%ny+1))
case(  193)
   if (associated(s%ero)) deallocate(s%ero)
   allocate(s%ero(s%nx+1,s%ny+1,par%ngd))
case(  194)
   if (associated(s%depo_im)) deallocate(s%depo_im)
   allocate(s%depo_im(s%nx+1,s%ny+1,par%ngd))
case(  195)
   if (associated(s%depo_ex)) deallocate(s%depo_ex)
   allocate(s%depo_ex(s%nx+1,s%ny+1,par%ngd))
case(  196)
   if (associated(s%nd)) deallocate(s%nd)
   allocate(s%nd(s%nx+1,s%ny+1))
case(  197)
   if (associated(s%pbbed)) deallocate(s%pbbed)
   allocate(s%pbbed(s%nx+1,s%ny+1,max(par%nd,2),par%ngd))
case(  198)
   if (associated(s%dzbed)) deallocate(s%dzbed)
   allocate(s%dzbed(s%nx+1,s%ny+1,max(par%nd,2)))
case(  199)
   if (associated(s%z0bed)) deallocate(s%z0bed)
   allocate(s%z0bed(s%nx+1,s%ny+1))
case(  200)
   if (associated(s%ureps)) deallocate(s%ureps)
   allocate(s%ureps(s%nx+1,s%ny+1))
case(  201)
   if (associated(s%vreps)) deallocate(s%vreps)
   allocate(s%vreps(s%nx+1,s%ny+1))
case(  202)
   if (associated(s%urepb)) deallocate(s%urepb)
   allocate(s%urepb(s%nx+1,s%ny+1))
case(  203)
   if (associated(s%vrepb)) deallocate(s%vrepb)
   allocate(s%vrepb(s%nx+1,s%ny+1))
case(  204)
   if (associated(s%umwci)) deallocate(s%umwci)
   allocate(s%umwci(s%nx+1,s%ny+1))
case(  205)
   if (associated(s%vmwci)) deallocate(s%vmwci)
   allocate(s%vmwci(s%nx+1,s%ny+1))
case(  206)
   if (associated(s%rolthick)) deallocate(s%rolthick)
   allocate(s%rolthick(s%nx+1,s%ny+1))
case(  207)
   if (associated(s%zswci)) deallocate(s%zswci)
   allocate(s%zswci(s%nx+1,s%ny+1))
case(  208)
   if (associated(s%pres)) deallocate(s%pres)
   allocate(s%pres(s%nx+1,s%ny+1))
case(  209)
   if (associated(s%dU)) deallocate(s%dU)
   allocate(s%dU(s%nx+1,s%ny+1))
case(  210)
   if (associated(s%dV)) deallocate(s%dV)
   allocate(s%dV(s%nx+1,s%ny+1))
case(  211)
   if (associated(s%wb)) deallocate(s%wb)
   allocate(s%wb(s%nx+1,s%ny+1))
case(  212)
   if (associated(s%ws)) deallocate(s%ws)
   allocate(s%ws(s%nx+1,s%ny+1))
case(  213)
   if (associated(s%wscrit)) deallocate(s%wscrit)
   allocate(s%wscrit(s%nx+1,s%ny+1))
case(  214)
   if (associated(s%bedfriccoef)) deallocate(s%bedfriccoef)
   allocate(s%bedfriccoef(s%nx+1,s%ny+1))
case(  215)
   if (associated(s%taubx)) deallocate(s%taubx)
   allocate(s%taubx(s%nx+1,s%ny+1))
case(  216)
   if (associated(s%tauby)) deallocate(s%tauby)
   allocate(s%tauby(s%nx+1,s%ny+1))
case(  217)
   if (associated(s%Df)) deallocate(s%Df)
   allocate(s%Df(s%nx+1,s%ny+1))
case(  218)
   if (associated(s%Dp)) deallocate(s%Dp)
   allocate(s%Dp(s%nx+1,s%ny+1))
case(  219)
   if (associated(s%Sutot)) deallocate(s%Sutot)
   allocate(s%Sutot(s%nx+1,s%ny+1))
case(  220)
   if (associated(s%Svtot)) deallocate(s%Svtot)
   allocate(s%Svtot(s%nx+1,s%ny+1))
case(  221)
   if (associated(s%cctot)) deallocate(s%cctot)
   allocate(s%cctot(s%nx+1,s%ny+1))
case(  222)
   if (associated(s%wi)) deallocate(s%wi)
   allocate(s%wi(2,s%ny+1))
case(  223)
   if (associated(s%dUi)) deallocate(s%dUi)
   allocate(s%dUi(2,s%ny+1))
case(  224)
   if (associated(s%zi)) deallocate(s%zi)
   allocate(s%zi(2,s%ny+1))
case(  225)
   if (associated(s%nuh)) deallocate(s%nuh)
   allocate(s%nuh(s%nx+1,s%ny+1))
case(  226)
   if (associated(s%cf)) deallocate(s%cf)
   allocate(s%cf(s%nx+1,s%ny+1))
case(  227)
   if (associated(s%cfu)) deallocate(s%cfu)
   allocate(s%cfu(s%nx+1,s%ny+1))
case(  228)
   if (associated(s%cfv)) deallocate(s%cfv)
   allocate(s%cfv(s%nx+1,s%ny+1))
case(  229)
   if (associated(s%D50top)) deallocate(s%D50top)
   allocate(s%D50top(s%nx+1,s%ny+1))
case(  230)
   if (associated(s%D90top)) deallocate(s%D90top)
   allocate(s%D90top(s%nx+1,s%ny+1))
case(  231)
   if (associated(s%structdepth)) deallocate(s%structdepth)
   allocate(s%structdepth(s%nx+1,s%ny+1))
case(  232)
   if (associated(s%zs0fac)) deallocate(s%zs0fac)
   allocate(s%zs0fac(s%nx+1,s%ny+1,2))
case(  233)
   if (associated(s%tdisch)) deallocate(s%tdisch)
   allocate(s%tdisch(par%ntdischarge))
case(  234)
   if (associated(s%pdisch)) deallocate(s%pdisch)
   allocate(s%pdisch(par%ndischarge,4))
case(  235)
   if (associated(s%pntdisch)) deallocate(s%pntdisch)
   allocate(s%pntdisch(par%ndischarge))
case(  236)
   if (associated(s%qdisch)) deallocate(s%qdisch)
   allocate(s%qdisch(par%ntdischarge,par%ndischarge))
case(  237)
   if (associated(s%idrift)) deallocate(s%idrift)
   allocate(s%idrift(par%ndrifter))
case(  238)
   if (associated(s%jdrift)) deallocate(s%jdrift)
   allocate(s%jdrift(par%ndrifter))
case(  239)
   if (associated(s%tdriftb)) deallocate(s%tdriftb)
   allocate(s%tdriftb(par%ndrifter))
case(  240)
   if (associated(s%tdrifte)) deallocate(s%tdrifte)
   allocate(s%tdrifte(par%ndrifter))
case(  241)
   if (associated(s%runup)) deallocate(s%runup)
   allocate(s%runup(s%ny+1))
case(  242)
   if (associated(s%Hrunup)) deallocate(s%Hrunup)
   allocate(s%Hrunup(s%ny+1))
case(  243)
   if (associated(s%xHrunup)) deallocate(s%xHrunup)
   allocate(s%xHrunup(s%ny+1))
case(  244)
   if (associated(s%istruct)) deallocate(s%istruct)
   allocate(s%istruct(s%ny+1))
case(  245)
   if (associated(s%iwl)) deallocate(s%iwl)
   allocate(s%iwl(s%ny+1))
case(  246)
   if (associated(s%strucslope)) deallocate(s%strucslope)
   allocate(s%strucslope(s%ny+1))
case(  247)
   if (associated(s%Dc)) deallocate(s%Dc)
   allocate(s%Dc(s%nx+1,s%ny+1))
case(  248)
   if (associated(s%ph)) deallocate(s%ph)
   allocate(s%ph(s%nx+1,s%ny+1))
case(  249)
case(  250)
   if (associated(s%dobs)) deallocate(s%dobs)
   allocate(s%dobs(s%nx+1,s%ny+1))
case(  251)
   if (associated(s%sig2prior)) deallocate(s%sig2prior)
   allocate(s%sig2prior(s%nx+1,s%ny+1))
case(  252)
   if (associated(s%zbobs)) deallocate(s%zbobs)
   allocate(s%zbobs(s%nx+1,s%ny+1))
case(  253)
   if (associated(s%shobs)) deallocate(s%shobs)
   allocate(s%shobs(s%nx+1,s%ny+1))
case(  254)
   if (associated(s%bwalpha)) deallocate(s%bwalpha)
   allocate(s%bwalpha(s%nx+1,s%ny+1))
case(  255)
   if (associated(s%dcmdo)) deallocate(s%dcmdo)
   allocate(s%dcmdo(s%nx+1,s%ny+1))
case(  256)
   if (associated(s%dassim)) deallocate(s%dassim)
   allocate(s%dassim(s%nx+1,s%ny+1))
case(  257)
   if (associated(s%cobs)) deallocate(s%cobs)
   allocate(s%cobs(s%nx+1,s%ny+1))
case(  258)
   if (associated(s%shipxCG)) deallocate(s%shipxCG)
   allocate(s%shipxCG(par%nship))
case(  259)
   if (associated(s%shipyCG)) deallocate(s%shipyCG)
   allocate(s%shipyCG(par%nship))
case(  260)
   if (associated(s%shipzCG)) deallocate(s%shipzCG)
   allocate(s%shipzCG(par%nship))
case(  261)
   if (associated(s%shipFx)) deallocate(s%shipFx)
   allocate(s%shipFx(par%nship))
case(  262)
   if (associated(s%shipFy)) deallocate(s%shipFy)
   allocate(s%shipFy(par%nship))
case(  263)
   if (associated(s%shipFz)) deallocate(s%shipFz)
   allocate(s%shipFz(par%nship))
case(  264)
   if (associated(s%shipMx)) deallocate(s%shipMx)
   allocate(s%shipMx(par%nship))
case(  265)
   if (associated(s%shipMy)) deallocate(s%shipMy)
   allocate(s%shipMy(par%nship))
case(  266)
   if (associated(s%shipMz)) deallocate(s%shipMz)
   allocate(s%shipMz(par%nship))
case(  267)
   if (associated(s%shipphi)) deallocate(s%shipphi)
   allocate(s%shipphi(par%nship))
case(  268)
   if (associated(s%shipchi)) deallocate(s%shipchi)
   allocate(s%shipchi(par%nship))
case(  269)
   if (associated(s%shippsi)) deallocate(s%shippsi)
   allocate(s%shippsi(par%nship))
case(  270)
   if (associated(s%vegtype)) deallocate(s%vegtype)
   allocate(s%vegtype(s%nx+1,s%ny+1))
case(  271)
   if (associated(s%Cdrag)) deallocate(s%Cdrag)
   allocate(s%Cdrag(s%nx+1,s%ny+1))
case(  272)
   if (associated(s%Dveg)) deallocate(s%Dveg)
   allocate(s%Dveg(s%nx+1,s%ny+1))
case(  273)
   if (associated(s%Fvegu)) deallocate(s%Fvegu)
   allocate(s%Fvegu(s%nx+1,s%ny+1))
case(  274)
   if (associated(s%Fvegv)) deallocate(s%Fvegv)
   allocate(s%Fvegv(s%nx+1,s%ny+1))
case(  275)
   if (associated(s%ududx)) deallocate(s%ududx)
   allocate(s%ududx(s%nx+1,s%ny+1))
case(  276)
   if (associated(s%vdvdy)) deallocate(s%vdvdy)
   allocate(s%vdvdy(s%nx+1,s%ny+1))
case(  277)
   if (associated(s%udvdx)) deallocate(s%udvdx)
   allocate(s%udvdx(s%nx+1,s%ny+1))
case(  278)
   if (associated(s%vdudy)) deallocate(s%vdudy)
   allocate(s%vdudy(s%nx+1,s%ny+1))
case(  279)
   if (associated(s%viscu)) deallocate(s%viscu)
   allocate(s%viscu(s%nx+1,s%ny+1))
case(  280)
   if (associated(s%viscv)) deallocate(s%viscv)
   allocate(s%viscv(s%nx+1,s%ny+1))
case(  281)
   if (associated(s%setbathy)) deallocate(s%setbathy)
   allocate(s%setbathy(s%nx+1,s%ny+1,par%nsetbathy))
case(  282)
   if (associated(s%tsetbathy)) deallocate(s%tsetbathy)
   allocate(s%tsetbathy(par%nsetbathy))
case(  283)
   if (associated(s%breaking)) deallocate(s%breaking)
   allocate(s%breaking(s%nx+1,s%ny+1))
case(  284)
   if (associated(s%fw)) deallocate(s%fw)
   allocate(s%fw(s%nx+1,s%ny+1))
case(  285)
   if (associated(s%sigz)) deallocate(s%sigz)
   allocate(s%sigz(par%nz))
case(  286)
   if (associated(s%uz)) deallocate(s%uz)
   allocate(s%uz(s%nx+1,s%ny+1,par%nz))
case(  287)
   if (associated(s%vz)) deallocate(s%vz)
   allocate(s%vz(s%nx+1,s%ny+1,par%nz))
case(  288)
   if (associated(s%ustz)) deallocate(s%ustz)
   allocate(s%ustz(s%nx+1,s%ny+1,par%nz))
case(  289)
   if (associated(s%nutz)) deallocate(s%nutz)
   allocate(s%nutz(s%nx+1,s%ny+1,par%nz))
case(  290)
   if (associated(s%dzs0dn)) deallocate(s%dzs0dn)
   allocate(s%dzs0dn(s%nx+1,s%ny+1))

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