Index: LMDZ6/trunk/libf/phylmd/yamada_c.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/yamada_c.F90	(revision 5961)
+++ LMDZ6/trunk/libf/phylmd/yamada_c.F90	(revision 5962)
@@ -99,14 +99,4 @@
       integer nlay,nlev
 
-!ym      logical first
-!ym      integer ipas
-!ym      save first,ipas
-!FH/IM     data first,ipas/.true.,0/
-!ym      data first,ipas/.false.,0/
-!ym!$OMP THREADPRIVATE( first,ipas)
-!ym       INTEGER, SAVE :: iflag_tke_diff=0
-!ym!$OMP THREADPRIVATE(iflag_tke_diff)
-
-
       integer ig,k
 
@@ -122,13 +112,7 @@
       real leff(klon,klev+1)
       real l0(klon)
-!ym      real,allocatable,save :: l0(:)
-!ym!$OMP THREADPRIVATE(l0)      
       real sq(klon),sqz(klon),zz(klon,klev+1)
       integer iter
 
-!ym      real ric,rifc,b1,kap
-!ym      save ric,rifc,b1,kap
-!ym      data ric,rifc,b1,kap/0.195,0.191,16.6,0.4/
-!ym!$OMP THREADPRIVATE(ric,rifc,b1,kap)
       real, parameter :: ric=0.195,rifc=0.191,b1=16.6,kap=0.4
       real frif,falpha,fsm
@@ -138,6 +122,5 @@
       real lyam(klon,klev),knyam(klon,klev)
       real w2yam(klon,klev),t2yam(klon,klev)
-!ym      logical,save :: firstcall=.true.
-!ym!$OMP THREADPRIVATE(firstcall)       
+
       CHARACTER(len=20),PARAMETER :: modname="yamada_c"
 REAL, DIMENSION(klon,klev+1) :: fluxu,fluxv,fluxt
@@ -160,9 +143,4 @@
 
       okiophys=klon==1
-!ym      if (firstcall) then
-!ym        CALL getin_p('iflag_tke_diff',iflag_tke_diff)
-!ym        allocate(l0(klon))
-!ym        firstcall=.false.
-!ym      endif
 
    IF (ngrid<=0) RETURN ! Bizarre : on n a pas ce probeleme pour coef_diff_turb
@@ -246,5 +224,4 @@
 
 
-!ym      ipas=ipas+1
 
 
@@ -338,6 +315,6 @@
       l0(ig)=0.2*sqz(ig)/sq(ig)
                                                           enddo
-      l(:,1) = 0.
-      l(:,klev+1) = 0.
+      l(:,1) = 0.         !ym <- fix unitialized level
+      l(:,klev+1) = 0.    !ym <- fix unitialized level
       do k=2,klev
                                                           do ig=1,ngrid
@@ -354,6 +331,6 @@
 
           l0(:)=150.
-          l(:,1) = 0.
-          l(:,klev+1) = 0.
+          l(:,1) = 0.       !ym <- fix unitialized level
+          l(:,klev+1) = 0.  !ym <- fix unitialized level
           do k=2,klev
                                                           do ig=1,ngrid