!--------------------------------------------------------------------------
   !---------------------------- compute_theta ----------------------------------
   IF(caldyn_eta==eta_mass) THEN ! Compute mass
      DO l = ll_begin, ll_end
         !DIR$ SIMD
         DO ij=ij_begin_ext, ij_end_ext
            m = MASS_DAK(ij,l)+(mass_col(ij)*g+ptop)*MASS_DBK(ij,l)
            rhodz(ij,l) = m/g
         END DO
      END DO
   END IF
   DO iq=1,nqdyn
      DO l = ll_begin, ll_end
         !DIR$ SIMD
         DO ij=ij_begin_ext, ij_end_ext
            theta(ij,l,iq) = theta_rhodz(ij,l,iq)/rhodz(ij,l)
         END DO
      END DO
   END DO
   !---------------------------- compute_theta ----------------------------------
   !--------------------------------------------------------------------------