Changeset 4882

Timestamp:

Mar 28, 2024, 9:36:08 PM (5 weeks ago)

Author:

evignon

Message:

nouvelle formule pour la conductivite thermique de l'air pour lscp poprecip

Location:

LMDZ6/trunk/libf/phylmd

Files:

• lmdz_lscp_ini.F90 (modified) (1 diff)
• lmdz_lscp_poprecip.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90

@@ -187 +187 @@
   REAL, SAVE, PROTECTED :: eps=1.E-10                       ! A COMMENTER TODO [-]
   !$OMP THREADPRIVATE(eps)
-
-  REAL, SAVE, PROTECTED :: air_thermal_conduct=2.4e-2      ! A COMMENTER TODO [-]
-  !$OMP THREADPRIVATE(air_thermal_conduct)
 
   REAL, SAVE, PROTECTED :: coef_ventil=1.                   ! A COMMENTER TODO [-]

LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90

@@ -263 +263 @@
                           rho_rain, rho_snow, r_rain, r_snow, rho_ice, r_ice,  &
                           tau_auto_snow_min, tau_auto_snow_max,                &
-                          thresh_precip_frac, eps, air_thermal_conduct,        &
+                          thresh_precip_frac, eps,                             &
                           coef_ventil, alpha_freez, beta_freez, temp_nowater,  &
                           iflag_cloudth_vert, iflag_rain_incloud_vol,          &
@@ -344 +344 @@
 
 !--Melting
-REAL :: dqsmelt_max
+REAL :: dqsmelt_max, air_thermal_conduct
 REAL :: nb_snowflake_clr, nb_snowflake_cld
 REAL :: capa_snowflake, temp_wetbulb
@@ -715 +715 @@
       nb_snowflake_clr = snowclr(i) / precipfracclr(i) / snow_fallspeed_clr &
                        / ( 4. / 3. * RPI * r_snow**3. * rho_snow )
+      air_thermal_conduct=(5.69+0.017*(temp(i)-RTT))*1.e-3*4.184                ! thermal conductivity of the air, SI
       dqsclrmelt = - nb_snowflake_clr * 4. * RPI * air_thermal_conduct &
                    * capa_snowflake / RLMLT * coef_ventil &