1 | ! |
---|
2 | ! $Header$ |
---|
3 | ! |
---|
4 | SUBROUTINE transp (paprs,tsol, & |
---|
5 | & t, q, u, v, geom, & |
---|
6 | & vtran_e, vtran_q, utran_e, utran_q) |
---|
7 | !c |
---|
8 | USE dimphy |
---|
9 | IMPLICIT none |
---|
10 | !c====================================================================== |
---|
11 | !c Auteur(s): Z.X.Li (LMD/CNRS) |
---|
12 | !c Date: le 25 avril 1994 |
---|
13 | !c Objet: Calculer le transport de l'energie et de la vapeur d'eau |
---|
14 | !c====================================================================== |
---|
15 | !c |
---|
16 | !cym#include "dimensions.h" |
---|
17 | !cym#include "dimphy.h" |
---|
18 | #include "YOMCST.h" |
---|
19 | !c |
---|
20 | REAL paprs(klon,klev+1), 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 = RCPD*t(i,l) + RLVTT*q(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 SUBROUTINE transp |
---|