source: LMDZ6/branches/WETDEP_DECOUPLE/libf/phylmd/qcheck.F90 @ 5407

Last change on this file since 5407 was 2419, checked in by Ehouarn Millour, 9 years ago

Extracted "qcheck" function from physiq_mod.F90 module; it is self-standing and general, and moreover its being in the physiq module causes compilation issues with ifort.
EM

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