Changeset 4895

Timestamp:

Apr 9, 2024, 8:29:30 AM (4 weeks ago)

Author:

evignon

Message:

changement de nom d'un paramètre de proprecip. coef_ventil -> gamma_melt

Location:

LMDZ6/trunk/libf/phylmd

Files:

• lmdz_lscp_ini.F90 (modified) (3 diffs)
• lmdz_lscp_poprecip.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90

@@ -188 +188 @@
   !$OMP THREADPRIVATE(eps)
 
-  REAL, SAVE, PROTECTED :: coef_ventil=1.                   ! A COMMENTER TODO [-]
-  !$OMP THREADPRIVATE(coef_ventil)
+  REAL, SAVE, PROTECTED :: gamma_melt=1.                   ! A COMMENTER TODO [-]
+  !$OMP THREADPRIVATE(gamma_melt)
 
   REAL, SAVE, PROTECTED :: alpha_freez=4.                 ! A COMMENTER TODO [-]
@@ -306 +306 @@
     CALL getin_p('gamma_rim',gamma_rim)
     CALL getin_p('gamma_freez',gamma_freez)
+    CALL getin_p('gamma_melt',gamma_melt)
     CALL getin_p('r_snow',r_snow)
     CALL getin_p('rain_fallspeed',rain_fallspeed)
@@ -363 +364 @@
     WRITE(lunout,*) 'lscp_ini, gamma_rim:', gamma_rim
     WRITE(lunout,*) 'lscp_ini, gamma_freez:', gamma_freez
+    WRITE(lunout,*) 'lscp_ini, gamma_melt:', gamma_melt
     WRITE(lunout,*) 'lscp_ini, r_snow:', r_snow
     WRITE(lunout,*) 'lscp_ini, rain_fallspeed_clr:', rain_fallspeed_clr

LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90

@@ -264 +264 @@
                           tau_auto_snow_min, tau_auto_snow_max,                &
                           thresh_precip_frac, eps,                             &
-                          coef_ventil, alpha_freez, beta_freez, temp_nowater,  &
+                          gamma_melt, alpha_freez, beta_freez, temp_nowater,   &
                           iflag_cloudth_vert, iflag_rain_incloud_vol,          &
                           cld_lc_lsc_snow, cld_lc_con_snow, gamma_freez,       &
@@ -689 +689 @@
   !--snow_fallspeed: snow fall velocity (in clear/cloudy sky) [m/s]
   !--air_thermal_conduct: thermal conductivity of the air [J/m/K/s]
-  !--coef_ventil: ventilation coefficient [-]
+  !--gamma_melt: melting tuning parameter [-]
   !--nb_snowflake: number of snowflakes (in clear/cloudy air) [-]
 
@@ -718 +718 @@
                        / ( 4. / 3. * RPI * r_snow**3. * rho_snow )
       dqsclrmelt = - nb_snowflake_clr * 4. * RPI * air_thermal_conduct &
-                   * capa_snowflake / RLMLT * coef_ventil &
+                   * capa_snowflake / RLMLT * gamma_melt &
                    * MAX(0., temp_wetbulb - RTT) * dtime
       
@@ -733 +733 @@
                        / ( 4. / 3. * RPI * r_snow**3. * rho_snow )
       dqscldmelt = - nb_snowflake_cld * 4. * RPI * air_thermal_conduct &
-                   * capa_snowflake / RLMLT * coef_ventil &
+                   * capa_snowflake / RLMLT * gamma_melt &
                    * MAX(0., temp_wetbulb - RTT) * dtime