Index: /LMDZ6/trunk/libf/phylmd/ocean_forced_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/ocean_forced_mod.F90	(revision 6179)
+++ /LMDZ6/trunk/libf/phylmd/ocean_forced_mod.F90	(revision 6180)
@@ -584,7 +584,7 @@
                f_cond_s = 0.
                ! 1 / ice-layer heat capacity
-               cal(i) = 2.*RCPD/(hice(ki)*ice_den*ice_cap)
+               cal(i) = 2.*RCPD/(max(hice(ki),h_ice_min)*ice_den*ice_cap)
                ! adjustment time-scale of conductive flux
-               dif_grnd(i) = cal(i)*ice_cond/(hice(ki)*ice_den)/RCPD
+               dif_grnd(i) = cal(i)*ice_cond/(max(hice(ki),h_ice_min)*ice_den)/RCPD
                ! penetrative shortwave flux...
                f_swpen = swnet(i)*pen_frac*exp(-pen_ext*hice(ki))
@@ -645,6 +645,6 @@
                f_cond_s = sno_cond*(tice(ki) - tsurf_new(i))/snow(i)
                ! 1 / heat capacity and conductive timescale
-               uscap = 2./ice_cap/(snow(i) + hice(ki)*ice_den)
-               ustau = uscap*ice_cond/(hice(ki)*ice_den)
+               uscap = 2./ice_cap/(snow(i) + max(hice(ki),h_ice_min)*ice_den)
+               ustau = uscap*ice_cond/(max(hice(ki),h_ice_min)*ice_den)
                ! update ice temp
                tice(ki) = (tice(ki) + dtime*(ustau*t_freeze - uscap*f_cond_s)) &
@@ -654,5 +654,5 @@
             END IF
             ! ice conductive flux (pos up)
-            f_cond_i = ice_cond*(t_freeze - tice(ki))/(hice(ki)*ice_den)
+            f_cond_i = ice_cond*(t_freeze - tice(ki))/(max(hice(ki),h_ice_min)*ice_den)
             f_bot(i) = fsea - f_cond_i
             fcdi(ki) = f_cond_i - fsea
@@ -746,4 +746,5 @@
 !               fsic(i) = pctsrf(ki, is_sic)
 !               hice(ki) = (0.166 + 0.911 * (fsic(i)**2))*(1. + 2.*(fsic(i)-0.351))
+!               hice(ki) = max(h_ice_min , hice(ki))
 !            END DO
           END IF