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 t(klon,klev), 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 t(klon,klev), q(klon,klev), u(klon,klev), v(klon,klev)
22	REAL utran_e(klon), utran_q(klon), vtran_e(klon), vtran_q(klon)
23	c
24	INTEGER i, l
25	c ------------------------------------------------------------------
26	REAL geom(klon,klev), e
27	c ------------------------------------------------------------------
28	DO i = 1, klon
29	utran_e(i) = 0.0
30	utran_q(i) = 0.0
31	vtran_e(i) = 0.0
32	vtran_q(i) = 0.0
33	ENDDO
34	c
35	DO l = 1, klev
36	DO i = 1, klon
37	e = RCPDt(i,l) + RLVTTq(i,l) + geom(i,l)
38	utran_e(i)=utran_e(i)+ u(i,l)e(paprs(i,l)-paprs(i,l+1))/RG
39	utran_q(i)=utran_q(i)+ u(i,l)*q(i,l)
40	. *(paprs(i,l)-paprs(i,l+1))/RG
41	vtran_e(i)=vtran_e(i)+ v(i,l)e(paprs(i,l)-paprs(i,l+1))/RG
42	vtran_q(i)=vtran_q(i)+ v(i,l)*q(i,l)
43	. *(paprs(i,l)-paprs(i,l+1))/RG
44	ENDDO
45	ENDDO
46	c
47	RETURN
48	END