Index: trunk/LMDZ.COMMON/libf/evolution/clim_state_init.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/clim_state_init.F90	(revision 4157)
+++ trunk/LMDZ.COMMON/libf/evolution/clim_state_init.F90	(revision 4170)
@@ -205,4 +205,5 @@
 use surf_ice,  only: set_is_h2o_perice_PCM, set_co2_perice_PCM
 use soil,      only: set_TI_PCM, set_inertiedat_PCM
+use ice_table, only: set_coef_ssdif_PCM
 use display,   only: print_msg, LVL_NFO
 
@@ -276,4 +277,7 @@
 call set_tsoil_PCM(tmp3d)
 
+call get_var_nc('coef_ssdif',tmp2d)
+call set_coef_ssdif_PCM(tmp2d)
+
 deallocate(tmp2d); allocate(tmp2d(ngrid,nsoil_PCM))
 call get_var_nc('inertiedat',tmp2d)
@@ -282,9 +286,4 @@
 call get_var_nc('inertiesoil',tmp3d)
 call set_TI_PCM(tmp3d)
-
-! To do?
-!   h2oice_depth
-!   d_coef
-!   lag_co2_ice
 
 ! Close/Deallocate
@@ -398,7 +397,7 @@
                     !!! transition
                     delta = depth_breccia
-                    TI(i,index_breccia + 1,islope) = sqrt((layer(index_breccia + 1) - layer(index_breccia))/                       &
-                                                               (((delta - layer(index_breccia))/(TI(i,index_breccia,islope)**2)) + &
-                                                               ((layer(index_breccia + 1) - delta)/(TI_breccia**2))))
+                    TI(i,index_breccia + 1,islope) = sqrt((layer(index_breccia + 1) - layer(index_breccia))/                  &
+                                                          (((delta - layer(index_breccia))/(TI(i,index_breccia,islope)**2)) + &
+                                                          ((layer(index_breccia + 1) - delta)/(TI_breccia**2))))
                     do k = index_breccia + 2,index_bedrock
                         TI(i,k,islope) = TI_breccia
@@ -411,7 +410,7 @@
                 !! transition
                 delta = depth_bedrock
-                TI(i,index_bedrock + 1,islope) = sqrt((layer(index_bedrock + 1) - layer(index_bedrock))/                   &
-                                                       (((delta - layer(index_bedrock))/(TI(i,index_bedrock,islope)**2)) + &
-                                                       ((layer(index_bedrock + 1) - delta)/(TI_breccia**2))))
+                TI(i,index_bedrock + 1,islope) = sqrt((layer(index_bedrock + 1) - layer(index_bedrock))/                  &
+                                                      (((delta - layer(index_bedrock))/(TI(i,index_bedrock,islope)**2)) + &
+                                                      ((layer(index_bedrock + 1) - delta)/(TI_breccia**2))))
                 do k = index_bedrock + 2,nsoil
                     TI(i,k,islope) = TI_bedrock
@@ -427,7 +426,7 @@
         do i = 1,ngrid
             if (inertiedat(i,index_breccia) < TI_breccia) then
-                inertiedat(i,index_breccia + 1) = sqrt((layer(index_breccia + 1) - layer(index_breccia))/                    &
-                                                        (((delta - layer(index_breccia))/(inertiedat(i,index_breccia)**2)) + &
-                                                        ((layer(index_breccia + 1) - delta)/(TI_breccia**2))))
+                inertiedat(i,index_breccia + 1) = sqrt((layer(index_breccia + 1) - layer(index_breccia))/                   &
+                                                       (((delta - layer(index_breccia))/(inertiedat(i,index_breccia)**2)) + &
+                                                       ((layer(index_breccia + 1) - delta)/(TI_breccia**2))))
                 do k = index_breccia + 2,index_bedrock
                     inertiedat(i,k) = TI_breccia
@@ -442,7 +441,7 @@
         delta = depth_bedrock
         do i = 1,ngrid
-            inertiedat(i,index_bedrock + 1) = sqrt((layer(index_bedrock + 1) - layer(index_bedrock))/                    &
-                                                    (((delta - layer(index_bedrock))/(inertiedat(i,index_bedrock)**2)) + &
-                                                    ((layer(index_bedrock + 1) - delta)/(TI_bedrock**2))))
+            inertiedat(i,index_bedrock + 1) = sqrt((layer(index_bedrock + 1) - layer(index_bedrock))/                   &
+                                                   (((delta - layer(index_bedrock))/(inertiedat(i,index_bedrock)**2)) + &
+                                                   ((layer(index_bedrock + 1) - delta)/(TI_bedrock**2))))
         end do