42 real(dp),
intent(in) :: z_sl
43 real(dp),
intent(in) :: dxi, deta, dzeta_c
45 integer(i4b) :: i, j, kt, kc
46 real(dp),
dimension(0:JMAX,0:IMAX) :: vz_sl
47 real(dp),
dimension(0:KCMAX) :: zeta_c_sgz, eaz_c_sgz, eaz_c_quotient_sgz
48 real(dp) :: dxi_inv, deta_inv
49 real(dp) :: dvx_dxi, dvy_deta
52 deta_inv = 1.0_dp/deta
57 zeta_c_sgz(kc) = (kc+0.5_dp)*dzeta_c
58 eaz_c_sgz(kc) = exp(deform*zeta_c_sgz(kc))
59 eaz_c_quotient_sgz(kc) = (eaz_c_sgz(kc)-1.0_dp)/(ea-1.0_dp)
67 if (maske(j,i)==3_i2b)
then
71 dvx_dxi = insq_g11_g(j,i)*insq_g22_g(j,i) &
72 *(vx_m(j,i)*sq_g22_sgx(j,i)-vx_m(j,i-1)*sq_g22_sgx(j,i-1)) &
75 dvy_deta = insq_g11_g(j,i)*insq_g22_g(j,i) &
76 *(vy_m(j,i)*sq_g11_sgy(j,i)-vy_m(j-1,i)*sq_g11_sgy(j-1,i)) &
81 vz_b(j,i) = 0.5_dp*(vx_m(j,i)+vx_m(j,i-1))*dzb_dxi_g(j,i) &
82 +0.5_dp*(vy_m(j,i)+vy_m(j-1,i))*dzb_deta_g(j,i) &
83 +dzb_dtau(j,i)-q_b_tot(j,i)
88 vz_sl(j,i) = vz_b(j,i) - (z_sl-zb(j,i))*(dvx_dxi+dvy_deta)
92 vz_s(j,i) = vz_sl(j,i) - (zs(j,i)-z_sl)*(dvx_dxi+dvy_deta)
97 vz_c(kc,j,i) = vz_sl(j,i) &
98 -(zm(j,i)+eaz_c_quotient_sgz(kc)*h_c(j,i)-z_sl) &
102 if ((n_cts(j,i) == -1).or.(n_cts(j,i) == 0))
then
106 vz_t(kt,j,i) = vz_b(j,i)
112 vz_t(kt,j,i) = vz_sl(j,i) &
114 +0.5_dp*(zeta_t(kt)+zeta_t(kt+1))*h_t(j,i) &