Index: LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90	(revision 4881)
+++ LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90	(revision 4882)
@@ -187,7 +187,4 @@
   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 [-]
Index: LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90	(revision 4881)
+++ LMDZ6/trunk/libf/phylmd/lmdz_lscp_poprecip.F90	(revision 4882)
@@ -263,5 +263,5 @@
                           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,5 +344,5 @@
 
 !--Melting
-REAL :: dqsmelt_max
+REAL :: dqsmelt_max, air_thermal_conduct
 REAL :: nb_snowflake_clr, nb_snowflake_cld
 REAL :: capa_snowflake, temp_wetbulb
@@ -715,4 +715,5 @@
       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 &