/[lmdze]/trunk/phylmd/concvl.f
ViewVC logotype

Annotation of /trunk/phylmd/concvl.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 206 - (hide annotations)
Tue Aug 30 12:52:46 2016 UTC (7 years, 8 months ago) by guez
Original Path: trunk/Sources/phylmd/concvl.f
File size: 2575 byte(s) 
Removed dimension klev of flux_[tquv] and y_flux_[tquv] in
clmain. Removed dimension klev of flux_[tquv] in physiq. Removed
dimension klev of flux_[tq] in hbtm. Removed dimension klev of
flux_[tq] in clqh and computations for layers other than the surface
layer. Removed dimension klev of flux_v in clvent and computations for
layers other than the surface layer. Values for layers other than the
surface layer were not used nor output (not even in LMDZ).

Removed argument dnwd0 of concvl. Simply write - mp in physiq
(following LMDZ).

Removed useless intermediary variables zxflux[tquv] in physiq.
1 guez 47 module concvl_m
2 guez 3
3 guez 47 IMPLICIT NONE
4 guez 3
5 guez 47 contains
6 guez 13
7 guez 206 SUBROUTINE concvl(paprs, play, t, q, u, v, sig1, w01, d_t, d_q, d_u, d_v, &
8     rain, kbas, itop_con, upwd, dnwd, ma, cape, iflag, qcondc, pmflxr, da, &
9     phi, mp)
10 guez 13
11 guez 180 ! From phylmd/concvl.F, version 1.3, 2005/04/15 12:36:17
12 guez 62 ! Author: Z. X. Li (LMD/CNRS)
13 guez 103 ! Date: 1993 August 18
14 guez 62 ! Objet : schéma de convection d'Emanuel (1991), interface
15 guez 13
16 guez 195 use comconst, only: dtphys
17 guez 52 use cv_driver_m, only: cv_driver
18 guez 69 USE dimphy, ONLY: klev, klon
19     USE fcttre, ONLY: foeew
20     USE suphec_m, ONLY: retv, rtt
21     USE yoethf_m, ONLY: r2es
22 guez 13
23 guez 180 REAL, INTENT (IN):: paprs(klon, klev + 1)
24 guez 72 REAL, INTENT (IN):: play(klon, klev)
25 guez 189 REAL, intent(in):: t(klon, klev) ! temperature (K)
26     real, intent(in):: q(klon, klev) ! fraction massique de vapeur d'eau
27 guez 91 real, INTENT (IN):: u(klon, klev), v(klon, klev)
28 guez 72 REAL, intent(inout):: sig1(klon, klev), w01(klon, klev)
29 guez 103 REAL, intent(out):: d_t(klon, klev)
30     REAL, intent(out):: d_q(klon, klev) ! increment de la vapeur d'eau
31     REAL, intent(out):: d_u(klon, klev), d_v(klon, klev)
32 guez 180 REAL, intent(out):: rain(klon) ! pluie (mm / s)
33     INTEGER, intent(out):: kbas(klon)
34 guez 183 integer, intent(inout):: itop_con(klon)
35 guez 62
36     REAL, intent(out):: upwd(klon, klev)
37 guez 180 ! saturated updraft mass flux (kg / m2 / s)
38 guez 62
39     real, intent(out):: dnwd(klon, klev)
40 guez 180 ! saturated downdraft mass flux (kg / m2 / s)
41 guez 62
42 guez 180 REAL ma(klon, klev)
43     real cape(klon) ! output (J / kg)
44 guez 195 INTEGER, intent(out):: iflag(klon)
45 guez 47 REAL qcondc(klon, klev)
46 guez 180 REAL pmflxr(klon, klev + 1)
47 guez 205 REAL, intent(out):: da(:, :) ! (klon, klev)
48     REAL, intent(out):: phi(:, :, :) ! (klon, klev, klev)
49 guez 13
50 guez 205 REAL, intent(out):: mp(:, :) ! (klon, klev) Mass flux of the
51     ! unsaturated downdraft, defined positive downward, in kg m-2
52     ! s-1. M_p in Emanuel (1991 928).
53    
54 guez 97 ! Local:
55 guez 180 REAL zx_qs, cor
56 guez 97 INTEGER i, k
57 guez 47 REAL qs(klon, klev)
58 guez 13
59 guez 47 !-----------------------------------------------------------------
60 guez 13
61 guez 47 DO k = 1, klev
62     DO i = 1, klon
63 guez 180 zx_qs = min(0.5, r2es * foeew(t(i, k), rtt >= t(i, k)) / play(i, k))
64     cor = 1. / (1. - retv * zx_qs)
65     qs(i, k) = zx_qs * cor
66 guez 47 END DO
67     END DO
68 guez 3
69 guez 180 CALL cv_driver(t, q, qs, u, v, play / 100., paprs / 100., iflag, d_t, &
70 guez 195 d_q, d_u, d_v, rain, pmflxr, sig1, w01, kbas, itop_con, ma, upwd, &
71 guez 205 dnwd, qcondc, cape, da, phi, mp)
72 guez 180 rain = rain / 86400.
73 guez 195 d_t = dtphys * d_t
74     d_q = dtphys * d_q
75     d_u = dtphys * d_u
76     d_v = dtphys * d_v
77 guez 3
78 guez 47 END SUBROUTINE concvl
79    
80     end module concvl_m
  ViewVC Help
Powered by ViewVC 1.1.21