Parent Directory | Revision Log
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 |