/[lmdze]/trunk/Sources/phylmd/qcheck.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/qcheck.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC revision 98 by guez, Tue May 13 17:23:16 2014 UTC
# Line 4  module qcheck_m Line 4  module qcheck_m
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 SUPHEC_M      ! 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    
Legend:
Removed from v.82  
changed lines
  Added in v.98
  ViewVC Help
Powered by ViewVC 1.1.21