Index: LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 3952)
+++ LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 3953)
@@ -1896,5 +1896,6 @@
                          AcoefU_x, AcoefU_w, AcoefV_x, AcoefV_w, &
                          BcoefH_x, BcoefH_w, BcoefQ_x, BcoefQ_w, &
-                         BcoefU_x, BcoefU_w, BcoefV_x, BcoefV_w  &
+                         BcoefU_x, BcoefU_w, BcoefV_x, BcoefV_w, &
+                         Kech_h_x, Kech_h_w, Kech_h  &
                          )
          CALL wx_pbl_prelim_beta(knon, dtime, ywake_s, ybeta,  &
Index: LMDZ6/trunk/libf/phylmd/wx_pbl_var_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/wx_pbl_var_mod.F90	(revision 3952)
+++ LMDZ6/trunk/libf/phylmd/wx_pbl_var_mod.F90	(revision 3953)
@@ -311,5 +311,6 @@
                                  AcoefU_x, AcoefU_w, AcoefV_x, AcoefV_w, &
                                  BcoefT_x, BcoefT_w, BcoefQ_x, BcoefQ_w, &
-                                 BcoefU_x, BcoefU_w, BcoefV_x, BcoefV_w  &
+                                 BcoefU_x, BcoefU_w, BcoefV_x, BcoefV_w, &
+                                 Kech_h_x_out, Kech_h_w_out, Kech_h_out  &
                                  )
 !
@@ -339,4 +340,6 @@
     REAL, DIMENSION(knon),        INTENT(IN)        :: BcoefT_x, BcoefT_w, BcoefQ_x, BcoefQ_w
     REAL, DIMENSION(knon),        INTENT(IN)        :: BcoefU_x, BcoefU_w, BcoefV_x, BcoefV_w
+!
+    REAL, DIMENSION(knon),        INTENT(OUT)       :: Kech_h_x_out, Kech_h_w_out, Kech_h_out
 !
 ! Local variables
@@ -462,4 +465,8 @@
          Kech_q(j) = Kech_q_x(j) + sigw(j)*dd_Kq
          Kech_m(j) = Kech_m_x(j) + sigw(j)*dd_Km
+!
+         Kech_h_x_out(j) = Kech_h_x(j)
+         Kech_h_w_out(j) = Kech_h_w(j)
+         Kech_h_out(j)   = Kech_h(j)
 !
 ! Effective exchange coefficient computation