Index: LMDZ6/trunk/libf/phylmd/yamada4.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/yamada4.F90	(revision 4017)
+++ LMDZ6/trunk/libf/phylmd/yamada4.F90	(revision 4019)
@@ -986,5 +986,5 @@
  INTEGER            ngrid              ! Nombre de points concern??s sur l'horizontal
  INTEGER            iflag_pbl          ! Choix du sch??ma de turbulence
- REAL            pbl_lmixmin_alpha  ! on active ou non le calcul de la longueur de melange minimum
+ REAL               pbl_lmixmin_alpha  ! on active ou non le calcul de la longueur de melange minimum en fonction du relief
  REAL               lmixmin            ! Minimum absolu de la longueur de m??lange
  REAL               zlay(klon, klev)   ! altitude du centre de la couche
@@ -1101,13 +1101,14 @@
   END IF
 
-!=================================================================================
-!  CALCUL d'une longueur de melange en fonctions de la topographie sous maille: l2
-! si plb_lmixmin_alpha=TRUE et si on se trouve sur de la terre ( pas actif sur les 
+!===========================================================================================
+!  CALCUL d'une longueur de melange minimum en fonctions de la topographie sous maille: l2
+! si pbl_lmixmin_alpha=TRUE et si on se trouve sur de la terre ( pas actif sur les 
 ! glacier, la glace de mer et les oc??ans)
-!=================================================================================
+!===========================================================================================
 
    l2(1:ngrid,:)=0.0
    l_mixmin(1:ngrid,:,nsrf)=0.
    l_mix(1:ngrid,:,nsrf)=0.
+   hlim(1:ngrid)=0.
 
    IF (nsrf .EQ. 1) THEN
@@ -1167,5 +1168,5 @@
         ENDDO
      ENDDO
-   ENDIF                                                                        ! pbl_lmixmin_alpha
+   ENDIF                                                                           ! pbl_lmixmin_alpha
 
 !==================================================================================
@@ -1183,15 +1184,11 @@
 ! Diagnostics
 
- DO k=2,klev
+ DO k=1,klev+1
     DO ig=1,ngrid
        jg=ni(ig)
        l_mix(jg,k,nsrf)=lmix(ig,k)
-       l_mixmin(jg,k,nsrf)=l2(ig,k)
+       l_mixmin(jg,k,nsrf)=MAX(l2(ig,k),lmixmin)
     ENDDO
  ENDDO
- DO ig=1,ngrid
-    jg=ni(ig)
-    l_mix(jg,1,nsrf)=hlim(ig)
- ENDDO