!--------------------------------------------------------------------------
   !---------------------------- div ----------------------------------
   DO l = ll_begin, ll_end
      !DIR$ SIMD
      DO ij=ij_begin_ext, ij_end_ext
         div_ij=0.
         div_ij = div_ij + ne_rup*le_de(ij+u_rup)*u(ij+u_rup,l)
         div_ij = div_ij + ne_lup*le_de(ij+u_lup)*u(ij+u_lup,l)
         div_ij = div_ij + ne_left*le_de(ij+u_left)*u(ij+u_left,l)
         div_ij = div_ij + ne_ldown*le_de(ij+u_ldown)*u(ij+u_ldown,l)
         div_ij = div_ij + ne_rdown*le_de(ij+u_rdown)*u(ij+u_rdown,l)
         div_ij = div_ij + ne_right*le_de(ij+u_right)*u(ij+u_right,l)
         divu(ij,l) = div_ij / Ai(ij)
      END DO
   END DO
   !---------------------------- div ----------------------------------
   !--------------------------------------------------------------------------