! $Header$

SUBROUTINE transp_lay(paprs, tsol, t, q, u, v, geom, vtran_e, vtran_q, &
        utran_e, utran_q)

  USE dimphy
  USE lmdz_yomcst

  IMPLICIT NONE
  ! ======================================================================
  ! Auteur(s): Z.X.Li (LMD/CNRS)
  ! Date: le 25 avril 1994
  ! Objet: Calculer le transport de l'energie et de la vapeur d'eau
  ! ======================================================================

  REAL paprs(klon, klev + 1), tsol(klon)
  REAL t(klon, klev), q(klon, klev), u(klon, klev), v(klon, klev)
  REAL utran_e(klon, klev), utran_q(klon, klev)
  REAL vtran_e(klon, klev), vtran_q(klon, klev)

  INTEGER i, l
  ! ------------------------------------------------------------------
  REAL geom(klon, klev), esh
  ! ------------------------------------------------------------------
  DO l = 1, klev
    DO i = 1, klon
      utran_e(i, l) = 0.0
      utran_q(i, l) = 0.0
      vtran_e(i, l) = 0.0
      vtran_q(i, l) = 0.0
    END DO
  END DO

  DO l = 1, klev
    DO i = 1, klon
      esh = rcpd * t(i, l) + rlvtt * q(i, l) + geom(i, l)
      utran_e(i, l) = utran_e(i, l) + u(i, l) * esh * (paprs(i, l) - paprs(i, l + 1)) / &
              rg
      utran_q(i, l) = utran_q(i, l) + u(i, l) * q(i, l) * (paprs(i, l) - paprs(i, l + 1 &
              )) / rg
      vtran_e(i, l) = vtran_e(i, l) + v(i, l) * esh * (paprs(i, l) - paprs(i, l + 1)) / &
              rg
      vtran_q(i, l) = vtran_q(i, l) + v(i, l) * q(i, l) * (paprs(i, l) - paprs(i, l + 1 &
              )) / rg
    END DO
  END DO

END SUBROUTINE transp_lay