libf/phylmd/transp.f

!
! $Header: /home/cvsroot/LMDZ4/libf/phylmd/transp.F,v 1.1.1.1 2004/05/19 12:53:09 lmdzadmin Exp $
!
      SUBROUTINE transp (paprs,tsol,
     e                   t, q, u, v, geom,
     s                   vtran_e, vtran_q, utran_e, utran_q)
c
      use dimens_m
      use dimphy
      use SUPHEC_M
      IMPLICIT none
c======================================================================
c Auteur(s): Z.X.Li (LMD/CNRS)
c Date: le 25 avril 1994
c Objet: Calculer le transport de l'energie et de la vapeur d'eau
c======================================================================
c
c
      REAL, intent(in):: paprs(klon,klev+1)
      real tsol(klon)
      REAL, INTENT (IN):: t(klon,klev)
      real q(klon,klev), u(klon,klev), v(klon,klev)
      REAL utran_e(klon), utran_q(klon), vtran_e(klon), vtran_q(klon)
c
      INTEGER i, l
c     ------------------------------------------------------------------
      REAL geom(klon,klev), e
c     ------------------------------------------------------------------
      DO i = 1, klon
         utran_e(i) = 0.0
         utran_q(i) = 0.0
         vtran_e(i) = 0.0
         vtran_q(i) = 0.0
      ENDDO
c
      DO l = 1, klev
      DO i = 1, klon
         e = RCPD*t(i,l) + RLVTT*q(i,l) + geom(i,l)
         utran_e(i)=utran_e(i)+ u(i,l)*e*(paprs(i,l)-paprs(i,l+1))/RG
         utran_q(i)=utran_q(i)+ u(i,l)*q(i,l)
     .                         *(paprs(i,l)-paprs(i,l+1))/RG
         vtran_e(i)=vtran_e(i)+ v(i,l)*e*(paprs(i,l)-paprs(i,l+1))/RG
         vtran_q(i)=vtran_q(i)+ v(i,l)*q(i,l)
     .                         *(paprs(i,l)-paprs(i,l+1))/RG
      ENDDO
      ENDDO
c
      RETURN
      END
1	!
2	! $Header: /home/cvsroot/LMDZ4/libf/phylmd/transp.F,v 1.1.1.1 2004/05/19 12:53:09 lmdzadmin Exp $
3	!
4	SUBROUTINE transp (paprs,tsol,
5	e t, q, u, v, geom,
6	s vtran_e, vtran_q, utran_e, utran_q)
7	c
8	use dimens_m
9	use dimphy
10	use SUPHEC_M
11	IMPLICIT none
12	c======================================================================
13	c Auteur(s): Z.X.Li (LMD/CNRS)
14	c Date: le 25 avril 1994
15	c Objet: Calculer le transport de l'energie et de la vapeur d'eau
16	c======================================================================
17	c
18	c
19	REAL, intent(in):: paprs(klon,klev+1)
20	real tsol(klon)
21	REAL, INTENT (IN):: t(klon,klev)
22	real q(klon,klev), u(klon,klev), v(klon,klev)
23	REAL utran_e(klon), utran_q(klon), vtran_e(klon), vtran_q(klon)
24	c
25	INTEGER i, l
26	c ------------------------------------------------------------------
27	REAL geom(klon,klev), e
28	c ------------------------------------------------------------------
29	DO i = 1, klon
30	utran_e(i) = 0.0
31	utran_q(i) = 0.0
32	vtran_e(i) = 0.0
33	vtran_q(i) = 0.0
34	ENDDO
35	c
36	DO l = 1, klev
37	DO i = 1, klon
38	e = RCPDt(i,l) + RLVTTq(i,l) + geom(i,l)
39	utran_e(i)=utran_e(i)+ u(i,l)e(paprs(i,l)-paprs(i,l+1))/RG
40	utran_q(i)=utran_q(i)+ u(i,l)*q(i,l)
41	. *(paprs(i,l)-paprs(i,l+1))/RG
42	vtran_e(i)=vtran_e(i)+ v(i,l)e(paprs(i,l)-paprs(i,l+1))/RG
43	vtran_q(i)=vtran_q(i)+ v(i,l)*q(i,l)
44	. *(paprs(i,l)-paprs(i,l+1))/RG
45	ENDDO
46	ENDDO
47	c
48	RETURN
49	END