1 | !-------------------------------------------------------------------------- |
---|
2 | !---------------------------- pvort_only ---------------------------------- |
---|
3 | DO l = ll_begin, ll_end |
---|
4 | !DIR$ SIMD |
---|
5 | DO ij=ij_begin_ext, ij_end_ext |
---|
6 | etav = 0.d0 |
---|
7 | etav = etav + ne_rup*u(ij+u_rup,l) |
---|
8 | etav = etav + ne_left*u(ij+t_rup+u_left,l) |
---|
9 | etav = etav + (-ne_lup)*u(ij+u_lup,l) |
---|
10 | hv=0. |
---|
11 | hv = hv + Riv2(ij,vup)*rhodz(ij,l) |
---|
12 | hv = hv + Riv2(ij+t_rup,vldown)*rhodz(ij+t_rup,l) |
---|
13 | hv = hv + Riv2(ij+t_lup,vrdown)*rhodz(ij+t_lup,l) |
---|
14 | qv(ij+z_up,l) = (etav + fv(ij+z_up)*Av(ij+z_up) )/(hv*Av(ij+z_up)) |
---|
15 | etav = 0.d0 |
---|
16 | etav = etav + (-ne_rdown)*u(ij+u_rdown,l) |
---|
17 | etav = etav + ne_right*u(ij+t_ldown+u_right,l) |
---|
18 | etav = etav + ne_ldown*u(ij+u_ldown,l) |
---|
19 | hv=0. |
---|
20 | hv = hv + Riv2(ij,vdown)*rhodz(ij,l) |
---|
21 | hv = hv + Riv2(ij+t_ldown,vrup)*rhodz(ij+t_ldown,l) |
---|
22 | hv = hv + Riv2(ij+t_rdown,vlup)*rhodz(ij+t_rdown,l) |
---|
23 | qv(ij+z_down,l) = (etav + fv(ij+z_down)*Av(ij+z_down) )/(hv*Av(ij+z_down)) |
---|
24 | END DO |
---|
25 | END DO |
---|
26 | DO l = ll_begin, ll_end |
---|
27 | !DIR$ SIMD |
---|
28 | DO ij=ij_begin, ij_end |
---|
29 | qu(ij+u_right,l)=0.5d0*(qv(ij+z_rdown,l)+qv(ij+z_rup,l)) |
---|
30 | qu(ij+u_lup,l)=0.5d0*(qv(ij+z_up,l)+qv(ij+z_lup,l)) |
---|
31 | qu(ij+u_ldown,l)=0.5d0*(qv(ij+z_ldown,l)+qv(ij+z_down,l)) |
---|
32 | END DO |
---|
33 | END DO |
---|
34 | !---------------------------- pvort_only ---------------------------------- |
---|
35 | !-------------------------------------------------------------------------- |
---|