4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
FUNCTION qcheck(klon, klev, paprs, q, ql, aire) |
pure FUNCTION qcheck(paprs, q, ql) |
8 |
|
|
9 |
! From phylmd/physiq.F, v 1.22 2006/02/20 09:38:28 |
! From phylmd/physiq.F, v 1.22 2006/02/20 09:38:28 |
10 |
|
|
11 |
use YOMCST |
! Calculer et imprimer l'eau totale. A utiliser pour vérifier |
12 |
|
! la conservation de l'eau. |
13 |
|
|
14 |
! Calculer et imprimer l'eau totale. A utiliser pour verifier |
use comgeomphy, only: airephy |
15 |
! la conservation de l'eau |
use dimphy, only: klon, klev |
16 |
|
use SUPHEC_M, ONLY: rg |
17 |
|
|
18 |
INTEGER klon, klev |
REAL, intent(in):: paprs(:, :) ! (klon, klev + 1) |
19 |
REAL, intent(in):: paprs(klon, klev+1) |
real, intent(in):: q(:, :), ql(:, :) ! (klon, klev) |
20 |
real q(klon, klev), ql(klon, klev) |
|
21 |
REAL aire(klon) |
! Local: |
22 |
REAL qtotal, zx, qcheck |
REAL qtotal, zx, qcheck |
23 |
INTEGER i, k |
INTEGER i, k |
24 |
|
|
25 |
|
!--------------------------------------------------------- |
26 |
|
|
27 |
zx = 0.0 |
zx = 0.0 |
28 |
DO i = 1, klon |
DO i = 1, klon |
29 |
zx = zx + aire(i) |
zx = zx + airephy(i) |
30 |
ENDDO |
ENDDO |
31 |
qtotal = 0.0 |
qtotal = 0.0 |
32 |
DO k = 1, klev |
DO k = 1, klev |
33 |
DO i = 1, klon |
DO i = 1, klon |
34 |
qtotal = qtotal + (q(i, k)+ql(i, k)) * aire(i) & |
qtotal = qtotal + (q(i, k)+ql(i, k)) * airephy(i) & |
35 |
*(paprs(i, k)-paprs(i, k+1))/RG |
*(paprs(i, k)-paprs(i, k+1))/RG |
36 |
ENDDO |
ENDDO |
37 |
ENDDO |
ENDDO |
38 |
|
|
39 |
qcheck = qtotal/zx |
qcheck = qtotal / zx |
40 |
|
|
41 |
END FUNCTION qcheck |
END FUNCTION qcheck |
42 |
|
|