Changeset 4910 for LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90

Timestamp:

Apr 18, 2024, 6:16:25 PM (4 weeks ago)

Author:

evignon

Message:

debut du chantier de refonte de cloudth_mpc

File:

• LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90

@@ -12 +12 @@
   !$OMP THREADPRIVATE(ok_bug_fonte_lscp)
 
-  REAL, SAVE, PROTECTED :: seuil_neb=0.001      ! cloud fraction threshold: a cloud really exists when exceeded
+  REAL, SAVE, PROTECTED :: seuil_neb=0.001      ! cloud fraction threshold: a cloud can precipitate when exceeded
   !$OMP THREADPRIVATE(seuil_neb)
+
+  REAL, SAVE, PROTECTED :: min_neb_th=1e-10      ! a cloud produced by bi-gaussian really exists when exceeded
+  !$OMP THREADPRIVATE(min_neb_th)
+
+  REAL, SAVE, PROTECTED :: min_frac_thermals=1.e-10   ! minimum thermals fraction for use of bigaussian distribution
+  !$OMP THREADPRIVATE(min_frac_thermals)
 
   INTEGER, SAVE :: lunout, prt_level            ! Logical unit number and level for standard output
@@ -43 +49 @@
   !$OMP THREADPRIVATE(a_tr_sca)
   
-  INTEGER, SAVE, PROTECTED ::  iflag_mpc_bl=0   ! flag to activate boundary layer mixed phase cloud param
-  !$OMP THREADPRIVATE(iflag_mpc_bl)
-  
+  REAL, SAVE, PROTECTED ::  min_frac_th_cld=1.e-10   ! minimum thermal fraction to compute a thermal cloud fraction
+  !$OMP THREADPRIVATE(min_frac_th_cld)
+
   LOGICAL, SAVE, PROTECTED :: ok_radocond_snow=.false. ! take into account the mass of ice precip in the cloud ice content seen by radiation
   !$OMP THREADPRIVATE(ok_radocond_snow)
@@ -260 +266 @@
     CALL getin_p('iflag_evap_prec',iflag_evap_prec)
     CALL getin_p('seuil_neb',seuil_neb)
-    CALL getin_p('iflag_mpc_bl',iflag_mpc_bl)
     CALL getin_p('ok_radocond_snow',ok_radocond_snow) 
     CALL getin_p('t_glace_max',t_glace_max)
@@ -320 +325 @@
     WRITE(lunout,*) 'lscp_ini, iflag_evap_prec:', iflag_evap_prec
     WRITE(lunout,*) 'lscp_ini, seuil_neb:', seuil_neb
-    WRITE(lunout,*) 'lscp_ini, iflag_mpc_bl:', iflag_mpc_bl
     WRITE(lunout,*) 'lscp_ini, ok_radocond_snow:', ok_radocond_snow
     WRITE(lunout,*) 'lscp_ini, t_glace_max:', t_glace_max