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

Timestamp:

Apr 29, 2019, 10:07:47 AM (6 years ago)

Author:

aboissinot

Message:

Fix in convadj.F
New version of the thermal plume model (new entrainment and detrainment formulae, alimentation is removed)

File:

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

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

@@ -6 +6 @@
 ! Flags for computations
                                                                   !  default
-      INTEGER,PARAMETER :: iflag_thermals_optflux     = 1         !  0        !
-      INTEGER,PARAMETER :: iflag_thermals_closure     = 2         !  2        !
-      INTEGER,PARAMETER :: iflag_thermals             = 18        !  18       !
-      
-! Flags for (terrestrial) diagnoses
-      
-      LOGICAL,PARAMETER :: sorties                    = .false.   !  false
-      INTEGER,PARAMETER :: iflag_trig_bl              = 1         !  1
-      INTEGER,PARAMETER :: iflag_clos_bl              = 1         !  1
-      INTEGER,PARAMETER :: iflag_coupl                = 5         !  5
+      INTEGER,PARAMETER :: iflag_thermals_optflux     = 1         !  0        -
+      INTEGER,PARAMETER :: dqimpl                     = 1         !  1        flag for thermcell_dq version (1 : implicit scheme || 0 : explicit scheme)
       
 ! Physical parameters
       
-      REAL,PARAMETER :: fact_thermals_ed_dz           = 0.007     !  0.007    !
       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                      = 0.9       !  0.9      !
-      REAL,PARAMETER :: afact                         = 2./3.     !  2./3.    !
-      REAL,PARAMETER :: fact_epsilon                  = 0.000     !  0.002    !
-      REAL,PARAMETER :: detr_q_power                  = 0.5       !  0.5      !
-      REAL,PARAMETER :: detr_q_coef                   = 0.012     !  0.012    !
-      REAL,PARAMETER :: mix0                          = 0.        !  0.       !
-      REAL,PARAMETER :: detr_min                      = 1.d-5     !  1.e-5    Minimal detrainment value
-      REAL,PARAMETER :: entr_min                      = 1.d-5     !  1.e-5    Maximal detrainment value
-      REAL,PARAMETER :: alphamax                      = 0.7       !           Maximal permitted updraft fraction
-      REAL,PARAMETER :: fomass_max                    = 0.5       !           Maximal permitted outgoing layer mass fraction
-      REAL,PARAMETER :: pres_limit                    = 1.e5      !  1.e5     !
+      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 :: nu                            = 0.000     !           Geometrical contributions to entrainment and detrainment
+      REAL,PARAMETER :: alphamax                      = 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     -
       
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -43 +30 @@
 !      ngrid. 
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-      INTEGER,PARAMETER :: linf                       = 2         !     1
+      INTEGER,PARAMETER :: linf                       = 1         !     1
       
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -50 +37 @@
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
       REAL,PARAMETER :: d_temp                        = 0.        !     0.
-      
-! Parameters for diagnoses
-      
-      REAL,PARAMETER :: alp_bl_k                      = 0.5       !     0.5
       
 ! Physical constants
@@ -64 +47 @@
       
 !$OMP THREADPRIVATE(RTT, RG, RKAPPA, RPI, RD)
-      
-!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-! AB : Parameters needed only for a loop in thermcell_alp (diagnoses).
-!      Maybe to be removed.
-!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-      INTEGER,PARAMETER :: nbsrf                      = 1