Sources/phylmd/qcheck.f

module qcheck_m

  IMPLICIT none

contains

  pure FUNCTION qcheck(paprs, q, ql)

    ! From phylmd/physiq.F, v 1.22 2006/02/20 09:38:28

    ! Calculer et imprimer l'eau totale. A utiliser pour vérifier
    ! la conservation de l'eau.

    use comgeomphy, only: airephy
    use dimphy, only: klon, klev
    use SUPHEC_M, ONLY: rg

    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 + airephy(i)
    ENDDO
    qtotal = 0.0
    DO k = 1, klev
       DO i = 1, klon
          qtotal = qtotal + (q(i, k)+ql(i, k)) * airephy(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	pure FUNCTION qcheck(paprs, q, ql)
8
9	! From phylmd/physiq.F, v 1.22 2006/02/20 09:38:28
10
11	! Calculer et imprimer l'eau totale. A utiliser pour vérifier
12	! la conservation de l'eau.
13
14	use comgeomphy, only: airephy
15	use dimphy, only: klon, klev
16	use SUPHEC_M, ONLY: rg
17
18	REAL, intent(in):: paprs(:, :) ! (klon, klev + 1)
19	real, intent(in):: q(:, :), ql(:, :) ! (klon, klev)
20
21	! Local:
22	REAL qtotal, zx, qcheck
23	INTEGER i, k
24
25	!---------------------------------------------------------
26
27	zx = 0.0
28	DO i = 1, klon
29	zx = zx + airephy(i)
30	ENDDO
31	qtotal = 0.0
32	DO k = 1, klev
33	DO i = 1, klon
34	qtotal = qtotal + (q(i, k)+ql(i, k)) * airephy(i) &
35	*(paprs(i, k)-paprs(i, k+1))/RG
36	ENDDO
37	ENDDO
38
39	qcheck = qtotal / zx
40
41	END FUNCTION qcheck
42
43	end module qcheck_m