Directory: | ./ |
---|---|
File: | phys/transp.f90 |
Date: | 2022-01-11 19:19:34 |
Exec | Total | Coverage | |
---|---|---|---|
Lines: | 18 | 18 | 100.0% |
Branches: | 6 | 6 | 100.0% |
Line | Branch | Exec | Source |
---|---|---|---|
1 | |||
2 | ! $Id: transp.F90 3257 2018-03-13 13:41:20Z fairhead $ | ||
3 | |||
4 | 480 | SUBROUTINE transp(paprs, tsol, t, q, ql, qs, u, v, geom, vtran_e, vtran_q, utran_e, & | |
5 | utran_q, vtran_w, utran_w) | ||
6 | |||
7 | USE dimphy | ||
8 | IMPLICIT NONE | ||
9 | ! ====================================================================== | ||
10 | ! Auteur(s): Z.X.Li (LMD/CNRS) | ||
11 | ! Date: le 25 avril 1994 | ||
12 | ! Objet: Calculer le transport de l'energie et de la vapeur d'eau | ||
13 | ! ====================================================================== | ||
14 | |||
15 | include "YOMCST.h" | ||
16 | |||
17 | REAL paprs(klon, klev+1), tsol(klon) | ||
18 | REAL t(klon, klev), q(klon, klev), ql(klon, klev), qs(klon, klev) | ||
19 | REAL u(klon, klev), v(klon, klev) | ||
20 | REAL utran_e(klon), utran_q(klon), vtran_e(klon), vtran_q(klon) | ||
21 | REAL utran_w(klon), vtran_w(klon) | ||
22 | |||
23 | INTEGER i, l | ||
24 | ! ------------------------------------------------------------------ | ||
25 | REAL geom(klon, klev), e | ||
26 | ! ------------------------------------------------------------------ | ||
27 |
2/2✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
|
477600 | DO i = 1, klon |
28 | 477120 | utran_e(i) = 0.0 | |
29 | 477120 | utran_q(i) = 0.0 | |
30 | 477120 | vtran_e(i) = 0.0 | |
31 | 477120 | vtran_q(i) = 0.0 | |
32 | 477120 | utran_w(i) = 0.0 | |
33 | 477600 | vtran_w(i) = 0.0 | |
34 | END DO | ||
35 | |||
36 |
2/2✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
|
19200 | DO l = 1, klev |
37 |
2/2✓ Branch 0 taken 18607680 times.
✓ Branch 1 taken 18720 times.
|
18626880 | DO i = 1, klon |
38 | ! e = rcpd*t(i, l) + rlvtt*q(i, l) + geom(i, l) | ||
39 | 18607680 | e = rcpd*t(i, l) + geom(i, l) | |
40 | 18607680 | utran_e(i) = utran_e(i) + u(i, l)*e*(paprs(i,l)-paprs(i,l+1))/rg | |
41 | 18607680 | utran_q(i) = utran_q(i) + u(i, l)*q(i, l)*(paprs(i,l)-paprs(i,l+1))/rg | |
42 | utran_w(i) = utran_w(i) + u(i, l)*(q(i, l)+ql(i, l)+qs(i, l)) & | ||
43 | 18607680 | *(paprs(i,l)-paprs(i,l+1))/rg | |
44 | 18607680 | vtran_e(i) = vtran_e(i) + v(i, l)*e*(paprs(i,l)-paprs(i,l+1))/rg | |
45 | 18607680 | vtran_q(i) = vtran_q(i) + v(i, l)*q(i, l)*(paprs(i,l)-paprs(i,l+1))/rg | |
46 | vtran_w(i) = vtran_w(i) + v(i, l)*(q(i, l)+ql(i, l)+qs(i, l)) & | ||
47 | 18626400 | *(paprs(i,l)-paprs(i,l+1))/rg | |
48 | END DO | ||
49 | END DO | ||
50 | |||
51 | 480 | RETURN | |
52 | END SUBROUTINE transp | ||
53 |