!--------------------------------------------------------------------------
   !---------------------------- remap_theta ----------------------------------
   ! IN : thetarhodz, eta
   ! TMP : thetarhodz_cum, new_thetarhodz_cum
   ! OUT : thetarhodz
   !$OMP DO SCHEDULE(STATIC)
   DO ij=1,primal_num
      thetarhodz_cum(1,ij)=0.
      DO l = 1,llm
         thetarhodz_cum(l+1,ij) = thetarhodz_cum(l,ij) + thetarhodz(l,ij)
      END DO
      DO l = 1,llm+1
         X = eta(l,ij)
         level = MIN(llm,FLOOR(X)) ! eta=llm+1 => level=llm, X=1
         X = X-level
         new_thetarhodz_cum(l,ij) = thetarhodz_cum(level,ij)+X*thetarhodz(level,ij)
      END DO
      DO l = 1,llm
         thetarhodz(l,ij) = new_thetarhodz_cum(l+1,ij) - new_thetarhodz_cum(l,ij)
      END DO
   END DO
   !$OMP END DO
   !---------------------------- remap_theta ----------------------------------
   !--------------------------------------------------------------------------