44 real(dp),
intent(in) :: dxi, deta, dzeta_c
46 integer(i4b) :: i, j, kc
47 real(dp) :: dxi_inv, deta_inv
48 real(dp) :: dvx_dxi, dvx_deta, dvy_dxi, dvy_deta
49 real(dp) :: aqxy1(0:kcmax)
50 real(dp) :: cvis1(0:kcmax)
55 deta_inv = 1.0_dp/deta
58 aqxy1(kc) = deform/(ea-1.0_dp)*eaz_c(kc)*dzeta_c
66 if (maske(j,i)==0_i2b.and..not.flag_sf(j,i))
then
70 vis_int_g(j,i) = (h_c(j,i)+h_t(j,i)) / flui_ave_sia(j,i)
72 else if ((maske(j,i)==1_i2b).or.(maske(j,i)==2_i2b))
then
76 vis_int_g(j,i) = 0.0_dp
83 dvx_dxi = (vx_m(j,i)-vx_m(j,i-1))*dxi_inv
84 dvy_deta = (vy_m(j,i)-vy_m(j-1,i))*deta_inv
86 dvx_deta = 0.25_dp*deta_inv &
87 *(vx_m(j+1,i)+vx_m(j+1,i-1)-vx_m(j-1,i)-vx_m(j-1,i-1))
88 dvy_dxi = 0.25_dp*dxi_inv &
89 *(vy_m(j,i+1)+vy_m(j-1,i+1)-vy_m(j,i-1)-vy_m(j-1,i-1))
91 de_ssa(j,i) = sqrt( dvx_dxi*dvx_dxi &
94 + 0.25_dp*(dvx_deta+dvy_dxi)*(dvx_deta+dvy_dxi) )
100 cvis1(kc) = aqxy1(kc)*h_c(j,i) &
102 temp_c(kc,j,i), temp_c_m(kc,j,i), 0.0_dp, &
103 enh_c(kc,j,i), .true.)
109 vis_int_g(j,i) = 0.0_dp
112 vis_int_g(j,i) = vis_int_g(j,i)+0.5_dp*(cvis1(kc+1)+cvis1(kc))