Changeset 2132 for trunk/LMDZ.GENERIC/libf/phystd/thermcell_mod.F90

Timestamp:

Apr 29, 2019, 3:13:04 PM (6 years ago)

Author:

aboissinot

Message:

Alimentation (more precisely lalim) is removed from thermcell_flux as well as corrections which need it.
Consequently, iflag_thermals_optflux flag is removed from thermcell_mod.

File:

• trunk/LMDZ.GENERIC/libf/phystd/thermcell_mod.F90 (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/phystd/thermcell_mod.F90

@@ -6 +6 @@
 ! Flags for computations
                                                                   !  default
-      INTEGER,PARAMETER :: iflag_thermals_optflux     = 1         !  0        -
       INTEGER,PARAMETER :: dqimpl                     = 1         !  1        flag for thermcell_dq version (1 : implicit scheme || 0 : explicit scheme)
       
@@ -13 +12 @@
       REAL,PARAMETER :: r_aspect_thermals             = 1.0       !           Aspect ratio of the thermals (width / height)
       REAL,PARAMETER :: tau_thermals                  = 0.        !  0.       Relaxation time
-      REAL,PARAMETER :: betalpha                      = 1.0       !  0.9      - between 0 (e=d) and 1 (rho*fraca=cst)
-      REAL,PARAMETER :: afact                         = 1.        !  2./3.    - buoyancy contribution, between 0 and 1
-      REAL,PARAMETER :: fact_epsilon                  = 1.e-4     !  2.e-3    - friction
+      REAL,PARAMETER :: betalpha                      = 1.0       !  0.9       - factor between 0 (e=d) and 1 (rho*fraca=cst)
+      REAL,PARAMETER :: afact                         = 1.        !  2./3.    Buoyancy contribution - factor between 0 and 1
+      REAL,PARAMETER :: fact_epsilon                  = 1.e-4     !  2.e-3    Friction at plume borders - exponential decrease
       REAL,PARAMETER :: nu                            = 0.000     !           Geometrical contributions to entrainment and detrainment
-      REAL,PARAMETER :: alphamax                      = 0.7       !  0.7      Maximal permitted updraft fraction
+      REAL,PARAMETER :: alpha_max                     = 0.7       !  0.7      Maximal permitted updraft fraction
       REAL,PARAMETER :: fomass_max                    = 0.5       !  0.5      Maximal permitted outgoing layer mass fraction
-      REAL,PARAMETER :: pres_limit                    = 2.e5      !  1.e5     -
+      REAL,PARAMETER :: pres_limit                    = 2.e5      !  1.e5     Minimal permitted pressure to trigger a thermal plume
       
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~