!--------------------------------------------------------------------------
   !---------------------------- pvort_only ----------------------------------
   DO l = ll_begin, ll_end
      !DIR$ SIMD
      DO ij=ij_begin_ext, ij_end_ext
         etav = 0.d0
         etav = etav + ne_rup*u(ij+u_rup,l)
         etav = etav + ne_left*u(ij+t_rup+u_left,l)
         etav = etav + (-ne_lup)*u(ij+u_lup,l)
         hv=0.
         hv = hv + Riv2(ij,vup)*rhodz(ij,l)
         hv = hv + Riv2(ij+t_rup,vldown)*rhodz(ij+t_rup,l)
         hv = hv + Riv2(ij+t_lup,vrdown)*rhodz(ij+t_lup,l)
         qv(ij+z_up,l) = (etav + fv(ij+z_up)*Av(ij+z_up) )/(hv*Av(ij+z_up))
         etav = 0.d0
         etav = etav + (-ne_rdown)*u(ij+u_rdown,l)
         etav = etav + ne_right*u(ij+t_ldown+u_right,l)
         etav = etav + ne_ldown*u(ij+u_ldown,l)
         hv=0.
         hv = hv + Riv2(ij,vdown)*rhodz(ij,l)
         hv = hv + Riv2(ij+t_ldown,vrup)*rhodz(ij+t_ldown,l)
         hv = hv + Riv2(ij+t_rdown,vlup)*rhodz(ij+t_rdown,l)
         qv(ij+z_down,l) = (etav + fv(ij+z_down)*Av(ij+z_down) )/(hv*Av(ij+z_down))
      END DO
   END DO
   DO l = ll_begin, ll_end
      !DIR$ SIMD
      DO ij=ij_begin, ij_end
         qu(ij+u_right,l)=0.5d0*(qv(ij+z_rdown,l)+qv(ij+z_rup,l))
         qu(ij+u_lup,l)=0.5d0*(qv(ij+z_up,l)+qv(ij+z_lup,l))
         qu(ij+u_ldown,l)=0.5d0*(qv(ij+z_ldown,l)+qv(ij+z_down,l))
      END DO
   END DO
   !---------------------------- pvort_only ----------------------------------
   !--------------------------------------------------------------------------