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

Diff of /trunk/Sources/phylmd/concvl.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 182 by guez, Tue Mar 15 17:07:47 2016 UTC revision 183 by guez, Wed Mar 16 14:42:58 2016 UTC
# Line 5  module concvl_m Line 5  module concvl_m
5  contains  contains
6    
7    SUBROUTINE concvl(dtime, paprs, play, t, q, u, v, sig1, w01, d_t, d_q, d_u, &    SUBROUTINE concvl(dtime, paprs, play, t, q, u, v, sig1, w01, d_t, d_q, d_u, &
8         d_v, rain, snow_con, kbas, itop_con, upwd, dnwd, dnwd0, ma, cape, &         d_v, rain, kbas, itop_con, upwd, dnwd, dnwd0, ma, cape, iflag, qcondc, &
9         iflag, qcondc, wd, pmflxr, da, phi, mp)         wd, pmflxr, da, phi, mp)
10    
11      ! From phylmd/concvl.F, version 1.3, 2005/04/15 12:36:17      ! From phylmd/concvl.F, version 1.3, 2005/04/15 12:36:17
12      ! Author: Z. X. Li (LMD/CNRS)      ! Author: Z. X. Li (LMD/CNRS)
13      ! Date: 1993 August 18      ! Date: 1993 August 18
14      ! Objet : schéma de convection d'Emanuel (1991), interface      ! Objet : schéma de convection d'Emanuel (1991), interface
     ! (driver commun aux versions 3 et 4)  
15    
16      use cv_driver_m, only: cv_driver      use cv_driver_m, only: cv_driver
17      USE dimphy, ONLY: klev, klon      USE dimphy, ONLY: klev, klon
# Line 31  contains Line 30  contains
30      REAL, intent(out):: d_q(klon, klev) ! increment de la vapeur d'eau      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)      REAL, intent(out):: d_u(klon, klev), d_v(klon, klev)
32      REAL, intent(out):: rain(klon) ! pluie (mm / s)      REAL, intent(out):: rain(klon) ! pluie (mm / s)
     REAL, intent(out):: snow_con(klon) ! neige (mm / s)  
33      INTEGER, intent(out):: kbas(klon)      INTEGER, intent(out):: kbas(klon)
34      integer itop_con(klon)      integer, intent(inout):: itop_con(klon)
35    
36      REAL, intent(out):: upwd(klon, klev)      REAL, intent(out):: upwd(klon, klev)
37      ! saturated updraft mass flux (kg / m2 / s)      ! saturated updraft mass flux (kg / m2 / s)
# Line 56  contains Line 54  contains
54      REAL zx_qs, cor      REAL zx_qs, cor
55      INTEGER i, k      INTEGER i, k
56      REAL qs(klon, klev)      REAL qs(klon, klev)
     REAL, save:: cbmf(klon)  
     INTEGER:: ifrst = 0  
57    
58      !-----------------------------------------------------------------      !-----------------------------------------------------------------
59    
     snow_con = 0.  
   
     IF (ifrst == 0) THEN  
        ifrst = 1  
        cbmf = 0.  
     END IF  
   
60      DO k = 1, klev      DO k = 1, klev
61         DO i = 1, klon         DO i = 1, klon
62            zx_qs = min(0.5, r2es * foeew(t(i, k), rtt >= t(i, k)) / play(i, k))            zx_qs = min(0.5, r2es * foeew(t(i, k), rtt >= t(i, k)) / play(i, k))
# Line 77  contains Line 66  contains
66      END DO      END DO
67    
68      CALL cv_driver(t, q, qs, u, v, play / 100., paprs / 100., iflag, d_t, &      CALL cv_driver(t, q, qs, u, v, play / 100., paprs / 100., iflag, d_t, &
69           d_q, d_u, d_v, rain, pmflxr, cbmf, sig1, w01, kbas, itop_con, dtime, &           d_q, d_u, d_v, rain, pmflxr, sig1, w01, kbas, itop_con, dtime, ma, &
70           ma, upwd, dnwd, dnwd0, qcondc, wd, cape, da, phi, mp)           upwd, dnwd, dnwd0, qcondc, wd, cape, da, phi, mp)
71    
72      rain = rain / 86400.      rain = rain / 86400.
73      d_t = dtime * d_t      d_t = dtime * d_t
Legend:
Removed from v.182  
changed lines
  Added in v.183
  ViewVC Help
Powered by ViewVC 1.1.21