--- trunk/Sources/phylmd/yamada4.f	2017/11/03 12:38:47	228
+++ trunk/Sources/phylmd/yamada4.f	2017/11/15 13:56:45	246
@@ -8,13 +8,13 @@
 
 contains
 
-  SUBROUTINE yamada4(dt, g, zlev, zlay, u, v, teta, cd, q2, km, kn, kq, &
-       ustar, iflag_pbl)
+  SUBROUTINE yamada4(dt, g, zlev, zlay, u, v, teta, q2, km, kn, ustar)
 
     ! From LMDZ4/libf/phylmd/yamada4.F, version 1.1 2004/06/22 11:45:36
 
-    use nr_util, only: assert, assert_eq
+    USE conf_phys_m, ONLY: iflag_pbl
     USE dimphy, ONLY: klev
+    use nr_util, only: assert, assert_eq
 
     REAL, intent(in):: dt ! pas de temps
     real, intent(in):: g
@@ -34,35 +34,25 @@
     ! temp\'erature potentielle au centre de chaque couche (en entr\'ee :
     ! la valeur au d\'ebut du pas de temps)
 
-    REAL, intent(in):: cd(:) ! (knon) cdrag, valeur au d\'ebut du pas de temps
-
     REAL, intent(inout):: q2(:, :) ! (knon, klev + 1)
     ! $q^2$ au bas de chaque couche 
     ! En entr\'ee : la valeur au d\'ebut du pas de temps ; en sortie : la
     ! valeur \`a la fin du pas de temps.
 
-    REAL km(:, :) ! (knon, klev + 1)
+    REAL, intent(out):: km(:, 2:) ! (knon, 2:klev)
     ! diffusivit\'e turbulente de quantit\'e de mouvement (au bas de
     ! chaque couche) (en sortie : la valeur \`a la fin du pas de temps)
 
-    REAL kn(:, :) ! (knon, klev + 1)
+    REAL, intent(out):: kn(:, 2:) ! (knon, 2:klev)
     ! diffusivit\'e turbulente des scalaires (au bas de chaque couche)
     ! (en sortie : la valeur \`a la fin du pas de temps)
 
-    REAL kq(:, :) ! (knon, klev + 1)
     real, intent(in):: ustar(:) ! (knon)
 
-    integer, intent(in):: iflag_pbl
-    ! iflag_pbl doit valoir 6, 8 ou 9
-    ! l = 6, on prend syst\'ematiquement une longueur d'\'equilibre
-    ! iflag_pbl = 6 : Mellor and Yamada 2.0
-    ! iflag_pbl = 8 : Mellor and Yamada 2.5
-    ! iflag_pbl = 9 : un test ?
-
     ! Local:
     integer knon
     real kmin, qmin
-    real pblhmin(size(cd)), coriol(size(cd)) ! (knon)
+    real pblhmin(size(ustar)), coriol(size(ustar)) ! (knon)
     real qpre
     REAL unsdz(size(zlay, 1), size(zlay, 2)) ! (knon, klev)
     REAL unsdzdec(size(zlev, 1), size(zlev, 2)) ! (knon, klev + 1)
@@ -82,8 +72,8 @@
     real zq
     real dtetadz(size(zlev, 1), size(zlev, 2)) ! (knon, klev + 1)
     real l(size(zlev, 1), size(zlev, 2)) ! (knon, klev + 1)
-    real l0(size(cd)) ! (knon)
-    real sq(size(cd)), sqz(size(cd)) ! (knon)
+    real l0(size(ustar)) ! (knon)
+    real sq(size(ustar)), sqz(size(ustar)) ! (knon)
     real zz(size(zlev, 1), size(zlev, 2)) ! (knon, klev + 1)
     integer iter
     real:: ric = 0.195, rifc = 0.191, b1 = 16.6
@@ -92,11 +82,11 @@
 
     call assert(any(iflag_pbl == [6, 8, 9]), "yamada4 iflag_pbl")
     knon = assert_eq([size(zlev, 1), size(zlay, 1), size(u, 1), size(v, 1), &
-         size(teta, 1), size(cd), size(q2, 1), size(km, 1), size(kn, 1), &
-         size(kq, 1)], "yamada4 knon")
+         size(teta, 1), size(ustar), size(q2, 1), size(km, 1), size(kn, 1)], &
+         "yamada4 knon")
     call assert(klev == [size(zlev, 2) - 1, size(zlay, 2), size(u, 2), &
-         size(v, 2), size(teta, 2), size(q2, 2) - 1, size(km, 2) - 1, &
-         size(kn, 2) - 1, size(kq, 2) - 1], "yamada4 klev")
+         size(v, 2), size(teta, 2), size(q2, 2) - 1, size(km, 2) + 1, &
+         size(kn, 2) + 1], "yamada4 klev")
 
     ipas = ipas + 1
 
@@ -258,7 +248,6 @@
           zq = sqrt(q2(ig, k))
           km(ig, k) = l(ig, k)*zq*sm(ig, k)
           kn(ig, k) = km(ig, k)*alpha(ig, k)
-          kq(ig, k) = l(ig, k)*zq*0.2
        enddo
     enddo
 
@@ -284,7 +273,6 @@
           if (kn(ig, k).lt.kmin.or.km(ig, k).lt.kmin) then
              kn(ig, k) = kmin
              km(ig, k) = kmin
-             kq(ig, k) = kmin
              ! la longueur de melange est suposee etre l = kap z
              ! K = l q Sm d'ou q2 = (K/l Sm)**2
              q2(ig, k) = (qmin/sm(ig, k))**2