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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 17 - (hide annotations)
Tue Aug 5 13:31:32 2008 UTC (15 years, 9 months ago) by guez
Original Path: trunk/libf/phylmd/qcheck.f90
File size: 758 byte(s) 
Created rule for "compare_sampl_*" files in
"Documentation/Manuel_LMDZE.texfol/Graphiques/GNUmakefile".

Extracted "qcheck", "radiornpb", "minmaxqfi" into separate files.

Read pressure coordinate of ozone coefficients once per run instead of
every day.

Added some "intent" attributes.

Added argument "nq" to "ini_histday". Replaced calls to "gr_fi_ecrit"
by calls to "gr_phy_write_2d". "Sigma_O3_Royer" is written to
"histday.nc" only if "nq >= 4". Moved "ini_histrac" to module
"ini_hist".

Compute "zmasse" in "physiq", pass it to "phytrac".

Removed computations of "pftsol*" and "ppsrf*" in "phytrac".

Do not use variable "rg" from module "YOMCST" in "TLIFT".
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     use YOMCST
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