Index: /LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 3838)
+++ /LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 3839)
@@ -1876,5 +1876,6 @@
                yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, &
                yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), &
-               yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :))
+               yt2m, yq2m, yt10m, yq10m, yu10m, yustar, &
+               yn2mout(:, nsrf, :))
           ELSE 
           CALL stdlevvar(klon, knon, is_ter, zxli, &
@@ -2712,5 +2713,6 @@
             uzon, vmer, tair1, qair1, zgeo1, &
             tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, &
-            yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :))
+            yt2m, yq2m, yt10m, yq10m, yu10m, yustar, &
+            yn2mout(:, nsrf, :))
         ELSE
         CALL stdlevvar(klon, knon, nsrf, zxli, &
@@ -2724,9 +2726,11 @@
             uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, &
             tairsol_x, qairsol, yz0m, yz0h_oupas, psfce, patm, &
-            yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yn2mout_x(:, nsrf, :))
+            yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yustar_x, &
+            yn2mout_x(:, nsrf, :))
          CALL stdlevvarn(klon, knon, nsrf, zxli, &
             uzon_w, vmer_w, tair1_w, qair1_w, zgeo1_w, &
             tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, &
-            yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yn2mout_w(:, nsrf, :))
+            yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w, &
+            yn2mout_w(:, nsrf, :))
         ELSE
         CALL stdlevvar(klon, knon, nsrf, zxli, &
Index: /LMDZ6/trunk/libf/phylmd/stdlevvar_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/stdlevvar_mod.F90	(revision 3838)
+++ /LMDZ6/trunk/libf/phylmd/stdlevvar_mod.F90	(revision 3839)
@@ -311,5 +311,6 @@
                            u1, v1, t1, q1, z1, &
                            ts1, qsurf, z0m, z0h, psol, pat1, &
-                           t_2m, q_2m, t_10m, q_10m, u_10m, n2mout)
+                           t_2m, q_2m, t_10m, q_10m, u_10m, ustar, &
+                           n2mout)
 !
       USE ioipsl_getin_p_mod, ONLY : getin_p
@@ -348,4 +349,5 @@
 ! u_2m--output-R- vitesse du vent a 2m
 ! u_10m--output-R- vitesse du vent a 10m
+! ustar--output-R- u*
 !AM
 ! t_10m--output-R- temperature de l'air a 10m
@@ -358,5 +360,5 @@
       REAL, dimension(klon), intent(in) :: psol, pat1
 !
-      REAL, dimension(klon), intent(out) :: t_2m, q_2m
+      REAL, dimension(klon), intent(out) :: t_2m, q_2m, ustar
       REAL, dimension(klon), intent(out) :: u_10m, t_10m, q_10m
       INTEGER, dimension(klon, 6), intent(out) :: n2mout
@@ -449,4 +451,5 @@
         tpot(i) = t1(i)* (psol(i)/pat1(i))**RKAPPA
         zdu2 = MAX(CEPDUE*CEPDUE, speed(i)**2)
+        ustar(i) = sqrt(cdram(i) * zdu2)
 !
       ENDDO