Index: LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 3401)
+++ LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90	(revision 3402)
@@ -3227,17 +3227,22 @@
                 ! F. Codron sensible default values for ocean and sea ice
                 IF (nsrf.EQ.is_oce) THEN
-                    tsurf(i,nsrf) = 271.35 ! Freezing sea water
-                    DO k=1,nsw
+                   tsurf(i,nsrf) = 271.35
+                   ! (temperature of sea water under sea ice, so that
+                   ! is also the temperature of appearing sea water)
+                   DO k=1,nsw
                       alb_dir(i,k,nsrf) = 0.06 ! typical Ocean albedo
                       alb_dif(i,k,nsrf) = 0.06
-                    ENDDO
+                   ENDDO
                 ELSE IF (nsrf.EQ.is_sic) THEN
-                    tsurf(i,nsrf) = 273.15 ! Melting ice 
-                    DO k=1,nsw
+                   tsurf(i,nsrf) = 271.35
+                   ! (Temperature at base of sea ice. Surface
+                   ! temperature could be higher, up to 0 Celsius
+                   ! degrees. We set it to -1.8 Celsius degrees for
+                   ! consistency with the ocean slab model.)
+                   DO k=1,nsw
                       alb_dir(i,k,nsrf) = 0.3 ! thin ice
                       alb_dif(i,k,nsrf) = 0.3
-                    ENDDO
+                   ENDDO
                 ENDIF
-                ! F. Codron
              ELSE
                 ! The continents have changed. The new fraction receives the mean sum of the existent fractions