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

Annotation of /trunk/phylmd/concvl.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 189 - (hide annotations)
Tue Mar 29 15:20:23 2016 UTC (8 years, 1 month ago) by guez
Original Path: trunk/Sources/phylmd/concvl.f
File size: 2506 byte(s) 
There was a function gr_phy_write_3d in dyn3d and a function
gr_phy_write_2d in module grid_change. Moved them into a new module
gr_phy_write_m under a generic interface gr_phy_write. Replaced calls
to gr_fi_ecrit by calls to gr_phy_write.

Removed arguments len, nloc and nd of cv30_compress.

Removed arguments wd and wd1 of cv30_uncompress, wd of cv30_yield, wd
of concvl, wd1 of cv_driver. Was just filled with 0. Removed option
ok_gust in physiq, never used.

In cv30_unsat, cv30_yield and cv_driver, we only need to define b to
level nl - 1.
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 97 SUBROUTINE concvl(dtime, paprs, play, t, q, u, v, sig1, w01, d_t, d_q, d_u, &
8 guez 183 d_v, rain, kbas, itop_con, upwd, dnwd, dnwd0, ma, cape, iflag, qcondc, &
9 guez 189 pmflxr, da, 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 52 use cv_driver_m, only: cv_driver
17 guez 69 USE dimphy, ONLY: klev, klon
18     USE fcttre, ONLY: foeew
19     USE suphec_m, ONLY: retv, rtt
20     USE yoethf_m, ONLY: r2es
21 guez 13
22 guez 69 REAL, INTENT (IN):: dtime ! pas d'integration (s)
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     real, intent(out):: dnwd0(klon, klev)
43 guez 180 ! unsaturated downdraft mass flux (kg / m2 / s)
44 guez 62
45 guez 180 REAL ma(klon, klev)
46     real cape(klon) ! output (J / kg)
47 guez 47 INTEGER iflag(klon)
48     REAL qcondc(klon, klev)
49 guez 180 REAL pmflxr(klon, klev + 1)
50 guez 99 REAL, intent(inout):: da(klon, klev), phi(klon, klev, klev), mp(klon, klev)
51 guez 13
52 guez 97 ! Local:
53 guez 180 REAL zx_qs, cor
54 guez 97 INTEGER i, k
55 guez 47 REAL qs(klon, klev)
56 guez 13
57 guez 47 !-----------------------------------------------------------------
58 guez 13
59 guez 47 DO k = 1, klev
60     DO i = 1, klon
61 guez 180 zx_qs = min(0.5, r2es * foeew(t(i, k), rtt >= t(i, k)) / play(i, k))
62     cor = 1. / (1. - retv * zx_qs)
63     qs(i, k) = zx_qs * cor
64 guez 47 END DO
65     END DO
66 guez 3
67 guez 180 CALL cv_driver(t, q, qs, u, v, play / 100., paprs / 100., iflag, d_t, &
68 guez 183 d_q, d_u, d_v, rain, pmflxr, sig1, w01, kbas, itop_con, dtime, ma, &
69 guez 189 upwd, dnwd, dnwd0, qcondc, cape, da, phi, mp)
70 guez 180 rain = rain / 86400.
71     d_t = dtime * d_t
72     d_q = dtime * d_q
73     d_u = dtime * d_u
74     d_v = dtime * d_v
75 guez 3
76 guez 47 END SUBROUTINE concvl
77    
78     end module concvl_m
  ViewVC Help
Powered by ViewVC 1.1.21