--- trunk/phylmd/Interface_surf/coef_diff_turb.f	2018/07/26 16:45:51	298
+++ trunk/phylmd/Interface_surf/coef_diff_turb.f	2018/09/06 14:25:07	303
@@ -12,7 +12,6 @@
     use nr_util, only: assert
 
     USE clesphys, ONLY: ok_kzmin
-    use comconst, only: dtphys
     use coefkz_m, only: coefkz
     use coefkzmin_m, only: coefkzmin
     use coefkz2_m, only: coefkz2
@@ -64,10 +63,10 @@
             * (paprs(:, 1) / pplay(:, k))**rkappa * (1. + 0.61 * q(:, k))
 
        zlev(:, 1) = 0.
-       zlev(:, klev + 1) = 2. * zlay(:, klev) - zlay(:, klev - 1)
        forall (k = 2:klev) zlev(:, k) = 0.5 * (zlay(:, k) + zlay(:, k-1))
+       zlev(:, klev + 1) = 2. * zlay(:, klev) - zlev(:, klev)
 
-       CALL yamada4(dtphys, zlev, zlay, u, v, teta, q2, coefm, coefh, &
+       CALL yamada4(zlev, zlay, u, v, teta, q2, coefm, coefh, &
             ustarhb(u(:, 1), v(:, 1), cdragm))
     else
        CALL coefkz(nsrf, paprs, pplay, ts, u, v, t, q, zgeop, coefm, coefh)