Index: LMDZ4/branches/LMDZ4-dev/libf/phylmd/cpl_mod.F90
===================================================================
--- LMDZ4/branches/LMDZ4-dev/libf/phylmd/cpl_mod.F90	(revision 1124)
+++ LMDZ4/branches/LMDZ4-dev/libf/phylmd/cpl_mod.F90	(revision 1125)
@@ -359,7 +359,5 @@
                read_u0(:,:), read_v0(:,:), tmp_r0(:,:))
 !$OMP END MASTER
-          CALL WriteField('read_u0',read_u0)
-          CALL WriteField('read_v0',read_v0)
-          CALL WriteField('read_r0',tmp_r0)
+
        ELSE
           read_u0(:,:) = 0.
@@ -442,5 +440,5 @@
 
   SUBROUTINE cpl_receive_seaice_fields(knon, knindex, &
-       tsurf_new, alb_new)
+       tsurf_new, alb_new, u0_new, v0_new)
 !
 ! This routine returns the fields for the seaice that have been read from the coupler
@@ -459,4 +457,6 @@
     REAL, DIMENSION(klon), INTENT(OUT)      :: tsurf_new
     REAL, DIMENSION(klon), INTENT(OUT)      :: alb_new
+    REAL, DIMENSION(klon), INTENT(OUT)      :: u0_new
+    REAL, DIMENSION(klon), INTENT(OUT)      :: v0_new
 
 ! Local variables
@@ -472,4 +472,6 @@
     CALL cpl2gath(read_alb_sic, alb_new, knon, knindex)
     CALL cpl2gath(read_sic, sic_new, knon, knindex)
+    CALL cpl2gath(read_u0, u0_new, knon, knindex)
+    CALL cpl2gath(read_v0, v0_new, knon, knindex)
 
 !*************************************************************************************
Index: LMDZ4/branches/LMDZ4-dev/libf/phylmd/ocean_cpl_mod.F90
===================================================================
--- LMDZ4/branches/LMDZ4-dev/libf/phylmd/ocean_cpl_mod.F90	(revision 1124)
+++ LMDZ4/branches/LMDZ4-dev/libf/phylmd/ocean_cpl_mod.F90	(revision 1125)
@@ -259,5 +259,5 @@
 
     CALL cpl_receive_seaice_fields(knon, knindex, &
-         tsurf_cpl, alb_cpl)
+         tsurf_cpl, alb_cpl, u0, v0)
 
     alb1_new(1:knon) = alb_cpl(1:knon)
@@ -273,9 +273,8 @@
     beta = 1.0
     
-! Suppose zero surface speed
-    u0(:)=0.0
-    v0(:)=0.0
-    u1_lay(:) = u1(:) - u0(:)
-    v1_lay(:) = v1(:) - v0(:)
+    DO i = 1, knon
+       u1_lay(i) = u1(i) - u0(i)
+       v1_lay(i) = v1(i) - v0(i)
+    END DO
 
     CALL calcul_fluxs(knon, is_sic, dtime, &