trunk/phylmd/qcheck.f

module qcheck_m

  IMPLICIT none

contains

  FUNCTION qcheck(klon, klev, paprs, q, ql, aire)

    ! 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 verifier
    ! la conservation de l'eau

    INTEGER klon, klev
    REAL, intent(in):: paprs(klon, klev+1)
    real q(klon, klev), ql(klon, klev)
    REAL aire(klon)
    REAL qtotal, zx, qcheck
    INTEGER i, k

    zx = 0.0
    DO i = 1, klon
       zx = zx + aire(i)
    ENDDO
    qtotal = 0.0
    DO k = 1, klev
       DO i = 1, klon
          qtotal = qtotal + (q(i, k)+ql(i, k)) * aire(i) &
               *(paprs(i, k)-paprs(i, k+1))/RG
       ENDDO
    ENDDO

    qcheck = qtotal/zx

  END FUNCTION qcheck

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