/[lmdze]/trunk/libf/phylmd/qcheck.f90
ViewVC logotype

Contents of /trunk/libf/phylmd/qcheck.f90

Parent Directory Parent Directory | Revision Log Revision Log

Revision 38 - (show annotations)
Thu Jan 6 17:52:19 2011 UTC (13 years, 4 months ago) by guez
File size: 760 byte(s) 
Extracted ASCII art from "inigeom" into a separate text file in the
documentation.

"test_disvert" now creates a separate file for layer thicknesses.

Moved variables from module "yomcst" to module "suphec_m" because this
is where those variables are defined. Kept in "yomcst" only parameters
of Earth orbit. Gave the attribute "parameter" to some variables of
module "suphec_m".

Variables of module "yoethf" were defined in procedure "suphec". Moved
these definitions to a new procedure "yoethf" in module "yoethf_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
  ViewVC Help
Powered by ViewVC 1.1.21