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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 38 - (hide 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 guez 17 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 guez 38 use SUPHEC_M
12 guez 17
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