59 SUBROUTINE ana_vmix_tile (ng, tile, model, &
60 & LBi, UBi, LBj, UBj, &
61 & IminS, ImaxS, JminS, JmaxS, &
63 & h, z_r, z_w, zeta, Akv, Akt)
76 integer,
intent(in) :: ng, tile, model
77 integer,
intent(in) :: LBi, UBi, LBj, UBj
78 integer,
intent(in) :: IminS, ImaxS, JminS, JmaxS
79 integer,
intent(in) :: knew
82 real(r8),
intent(in) :: h(LBi:,LBj:)
83 real(r8),
intent(in) :: z_r(LBi:,LBj:,:)
84 real(r8),
intent(in) :: z_w(LBi:,LBj:,0:)
85 real(r8),
intent(in) :: zeta(LBi:,LBj:,:)
86 real(r8),
intent(out) :: Akv(LBi:,LBj:,0:)
87 real(r8),
intent(out) :: Akt(LBi:,LBj:,0:,:)
89 real(r8),
intent(in) :: h(LBi:UBi,LBj:UBj)
90 real(r8),
intent(in) :: z_r(LBi:UBi,LBj:UBj,N(ng))
91 real(r8),
intent(in) :: z_w(LBi:UBi,LBj:UBj,0:N(ng))
92 real(r8),
intent(in) :: zeta(LBi:UBi,LBj:UBj,:)
93 real(r8),
intent(out) :: Akv(LBi:UBi,LBj:UBj,0:N(ng))
94 real(r8),
intent(out) :: Akt(LBi:UBi,LBj:UBj,0:N(ng),NAT)
99 integer :: i, itrc, j, k
101#include "set_bounds.h"
111 akv(i,j,k)=1.0e-03_r8+95.0e-04_r8*exp(z_w(i,j,k)/50.0_r8)+ &
112 & 95.0e-04_r8*exp(-(z_w(i,j,k)+h(i,j))/50.0_r8)
116#elif defined CHANNEL_NECK
120 akv(i,j,k)=2.0e-04_r8+8.0e-04_r8*exp(z_w(i,j,k)/5.0_r8)
124#elif defined COUPLING_TEST
128 akv(i,j,k)=2.0e-03_r8+8.0e-03_r8*exp(z_w(i,j,k)/1500.0_r8)
132#elif defined ESTUARY_TEST
140#elif defined LAKE_SIGNELL
148#elif defined NJ_BIGHT
152 akv(i,j,k)=1.0e-03_r8+2.0e-04_r8*exp(z_r(i,j,k)/10.0_r8)
156#elif defined SED_TEST1
160 akv(i,j,k)=0.025_r8*(h(i,j)+z_w(i,j,k))* &
161 & (1.0_r8-(h(i,j)+z_w(i,j,k))/ &
162 & (h(i,j)+zeta(i,j,knew)))
163 akt(i,j,k,
itemp)=akv(i,j,k)*0.49_r8/0.39_r8
174 akv(i,j,k)=0.41_r8*0.01_r8*(h(i,j)+z_w(i,j,k))* &
175 & (1.0_r8-(h(i,j)+z_w(i,j,k))/ &
176 & (h(i,j)+zeta(i,j,knew)))
180#elif defined SHOREFACE
184 akv(i,j,k)=0.025_r8*(h(i,j)+z_w(i,j,k))* &
185 & (1.0_r8-(h(i,j)+z_w(i,j,k))/ &
186 & (h(i,j)+zeta(i,j,knew)))
190#elif defined TEST_CHAN
194 akv(i,j,k)=0.41_r8*0.0625_r8*(h(i,j)+z_w(i,j,k))* &
195 & (1.0_r8-(h(i,j)+z_w(i,j,k))/ &
196 & (h(i,j)+zeta(i,j,knew)))
200#elif defined UPWELLING
204 akv(i,j,k)=2.0e-03_r8+8.0e-03_r8*exp(z_w(i,j,k)/150.0_r8)
209 ana_vmix.h: no values provided for akv.
216 & lbi, ubi, lbj, ubj, 0, n(ng), &
222 & lbi, ubi, lbj, ubj, 0, n(ng), &
240#elif defined CHANNEL_NECK
244 akt(i,j,k,
itemp)=2.0e-06_r8+ &
245 & 8.0e-06_r8*exp(z_w(i,j,k)/5.0_r8)
249#elif defined COUPLING_TEST
260#elif defined ESTUARY_TEST
264 akt(i,j,k,
itemp)=akv(i,j,k)
266 akt(i,j,k,
isalt)=akv(i,j,k)
271#elif defined LAKE_SIGNELL
282#elif defined NJ_BIGHT
286 akt(i,j,k,
itemp)=1.0e-05_r8+ &
287 & 2.0e-06_r8*exp(z_r(i,j,k)/10.0_r8)
298 akt(i,j,k,
itemp)=akv(i,j,k)
300 akt(i,j,k,
isalt)=akv(i,j,k)
305#elif defined SHOREFACE
309 akt(i,j,k,
itemp)=akv(i,j,k)
311 akt(i,j,k,
isalt)=akv(i,j,k)
316#elif defined TEST_CHAN
320 akt(i,j,k,
itemp)=akv(i,j,k)*0.49_r8/0.39_r8
327#elif defined UPWELLING
339 ana_vmix.h: no values provided for akt.
347 & lbi, ubi, lbj, ubj, 0, n(ng), &
354 & lbi, ubi, lbj, ubj, 0, n(ng), 1, nat, &