Index: trunk/LMDZ.COMMON/libf/evolution/dyn_ss_ice_m_wrapper.f90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/dyn_ss_ice_m_wrapper.f90	(revision 3487)
+++ trunk/LMDZ.COMMON/libf/evolution/dyn_ss_ice_m_wrapper.f90	(revision 3487)
@@ -0,0 +1,27 @@
+subroutine dyn_ss_ice_m_wrapper(ngrid,nsoil,tHIn,p0,pfrost,T_in,ssi_depth_in,porefill_in,porefill,ssi_depth)
+
+  use allinterfaces
+  implicit none
+  integer               ngrid,nsoil
+  integer               i
+  real(8), dimension(ngrid) ::  p0, pfrost
+  real(8), dimension(ngrid) ::  thIn, ssi_depth_in,T1,ssi_depth
+  real(8), dimension(nsoil,ngrid) :: T_in
+  real(8), dimension(nsoil,ngrid) :: porefill, porefill_in
+
+  do i=1,ngrid
+  !tHIn(i)=250
+  !p0(i)=611
+  !pfrost(i)=0.31
+  !T_in(1:nsoil,i)=180
+  !ssi_depth_in(i)=0.3*i
+  porefill(1:nsoil,i) = -9999
+  ssi_depth(i)=-9999
+  !porefill_in(1:nsoil,i) =0.01*i
+  T1(i)=T_in(1,i)
+  call  dyn_ss_ice_m(ssi_depth_in(i),T1(i),T_in(:,i),nsoil,thIn(i),p0(i),pfrost(i),porefill_in(:,i),porefill(:,i),ssi_depth(i))
+
+  !print *,'new porefill =',porefill(:,i)
+  enddo
+
+end subroutine dyn_ss_ice_m_wrapper