--- trunk/phylmd/qcheck.f90 2013/11/15 18:45:49 76 +++ trunk/phylmd/qcheck.f 2014/05/13 17:23:16 98 @@ -4,35 +4,39 @@ contains - FUNCTION qcheck(klon, klev, paprs, q, ql, aire) + pure FUNCTION qcheck(paprs, q, ql) ! From phylmd/physiq.F, v 1.22 2006/02/20 09:38:28 - use SUPHEC_M + ! Calculer et imprimer l'eau totale. A utiliser pour vérifier + ! la conservation de l'eau. - ! Calculer et imprimer l'eau totale. A utiliser pour verifier - ! la conservation de l'eau + use comgeomphy, only: airephy + use dimphy, only: klon, klev + use SUPHEC_M, ONLY: rg - INTEGER klon, klev - REAL, intent(in):: paprs(klon, klev+1) - real q(klon, klev), ql(klon, klev) - REAL aire(klon) + REAL, intent(in):: paprs(:, :) ! (klon, klev + 1) + real, intent(in):: q(:, :), ql(:, :) ! (klon, klev) + + ! Local: REAL qtotal, zx, qcheck INTEGER i, k + !--------------------------------------------------------- + zx = 0.0 DO i = 1, klon - zx = zx + aire(i) + zx = zx + airephy(i) ENDDO qtotal = 0.0 DO k = 1, klev DO i = 1, klon - qtotal = qtotal + (q(i, k)+ql(i, k)) * aire(i) & + qtotal = qtotal + (q(i, k)+ql(i, k)) * airephy(i) & *(paprs(i, k)-paprs(i, k+1))/RG ENDDO ENDDO - qcheck = qtotal/zx + qcheck = qtotal / zx END FUNCTION qcheck