!--------------------------------------------------------------------------
   !---------------------------- caldyn_mil ----------------------------------
   !$OMP DO SCHEDULE(STATIC)
   DO ij = 1, primal_num
      l=1
      m_il(l,ij) = .5*rhodz(l,ij)
      !DIR$ SIMD
      DO l = 2, llm
         m_il(l,ij) = .5*(rhodz(l,ij)+rhodz(l-1,ij))
      END DO
      l=llm+1
      m_il(l,ij) = .5*rhodz(l-1,ij)
   END DO
   !$OMP END DO
   !---------------------------- caldyn_mil ----------------------------------
   !--------------------------------------------------------------------------