!--------------------------------------------------------------------------
   !---------------------------- compute_divergence ----------------------------------
   DO l = 1, llm
      !DIR$ SIMD
      DO ij=ij_begin_ext, ij_end_ext
         div_ij = 0.d0
         div_ij = div_ij + ne_rup*ue(ij+u_rup,l)*le(ij+u_rup)
         div_ij = div_ij + ne_lup*ue(ij+u_lup,l)*le(ij+u_lup)
         div_ij = div_ij + ne_left*ue(ij+u_left,l)*le(ij+u_left)
         div_ij = div_ij + ne_ldown*ue(ij+u_ldown,l)*le(ij+u_ldown)
         div_ij = div_ij + ne_rdown*ue(ij+u_rdown,l)*le(ij+u_rdown)
         div_ij = div_ij + ne_right*ue(ij+u_right,l)*le(ij+u_right)
         div(ij,l) = div_ij / Ai(ij)
      END DO
   END DO
   !---------------------------- compute_divergence ----------------------------------
   !--------------------------------------------------------------------------