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

Timestamp:

Apr 11, 2024, 5:50:38 PM (5 weeks ago)

Author:

evignon

Message:

nettoyage poprecip + changements de valeur par defaut du seuil
pour l'autoconversion solide dans poprecip
Lea, Audran, Etienne

File:

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

LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90

@@ -140 +140 @@
   !$OMP THREADPRIVATE(ok_poprecip)
 
-  REAL, SAVE, PROTECTED :: cld_lc_lsc_snow               ! snow autoconversion coefficient, stratiform rain
+  REAL, SAVE, PROTECTED :: cld_lc_lsc_snow=2.e-5            ! snow autoconversion coefficient, stratiform. default from  Chaboureau and PInty 2006
   !$OMP THREADPRIVATE(cld_lc_lsc_snow)
 
-  REAL, SAVE, PROTECTED :: cld_lc_con_snow                ! snow autoconversion coefficient, convective rain
+  REAL, SAVE, PROTECTED :: cld_lc_con_snow=2.e-5            ! snow autoconversion coefficient, convective 
   !$OMP THREADPRIVATE(cld_lc_con_snow)
 
@@ -166 +166 @@
   REAL, SAVE, PROTECTED :: rho_ice=920.                    ! A COMMENTER TODO [kg/m3]
   !$OMP THREADPRIVATE(rho_ice)
-
-  REAL, SAVE, PROTECTED :: rho_snow                        ! A COMMENTER TODO [kg/m3]
-  !$OMP THREADPRIVATE(rho_snow)
 
   REAL, SAVE, PROTECTED :: r_rain=500.E-6                   ! A COMMENTER TODO [m]
@@ -279 +276 @@
     CALL getin_p('cld_lc_lsc',cld_lc_lsc)
     CALL getin_p('cld_lc_con',cld_lc_con)
-    cld_lc_lsc_snow=cld_lc_lsc
-    cld_lc_con_snow=cld_lc_con
     CALL getin_p('cld_lc_lsc_snow',cld_lc_lsc_snow)
     CALL getin_p('cld_lc_con_snow',cld_lc_con_snow)
@@ -393 +388 @@
 
 
-    !--Initialisations of constants depending on other constants
-    !--rho_snow formula from r_snow (Brandes et al. 2007 - JAMC)
-    rho_snow = 1.e3 * 0.178 * (r_snow * 2. * 1000.)**(-0.922)
-
-
     !AA Temporary initialisation
     a_tr_sca(1) = -0.5