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