source: LMDZ6/branches/Amaury_dev/libf/phylmd/qcheck.F90

Last change on this file was 5144, checked in by abarral, 3 months ago

Put YOMCST.h into modules

File size: 647 bytes
Line 
1FUNCTION qcheck(klon, klev, paprs, q, ql, aire)
2  USE lmdz_yomcst
3
4  IMPLICIT NONE
5
6  ! Calculer et imprimer l'eau totale. A utiliser pour verifier
7  ! la conservation de l'eau
8
9  INTEGER, INTENT(IN) :: klon, klev
10  REAL, INTENT(IN) :: paprs(klon, klev + 1), q(klon, klev), ql(klon, klev)
11  REAL, INTENT(IN) :: aire(klon)
12  REAL qtotal, zx, qcheck
13  INTEGER i, k
14
15  zx = 0.0
16  DO i = 1, klon
17    zx = zx + aire(i)
18  ENDDO
19  qtotal = 0.0
20  DO k = 1, klev
21    DO i = 1, klon
22      qtotal = qtotal + (q(i, k) + ql(i, k)) * aire(i) &
23              * (paprs(i, k) - paprs(i, k + 1)) / RG
24    ENDDO
25  ENDDO
26
27  qcheck = qtotal / zx
28
29END FUNCTION qcheck
30
Note: See TracBrowser for help on using the repository browser.