--- trunk/phylmd/Interface_surf/soil.f	2014/09/11 15:09:15	107
+++ trunk/phylmd/Interface_surf/soil.f	2014/09/16 14:00:41	108
@@ -99,9 +99,8 @@
     ! Depthts:
     ! --------
 
-    REAL fz, rk, fz1, rk1, rk2
+    REAL rk, fz1, rk1, rk2
 
-    fz(rk) = fz1*(dalph_soil**rk-1.)/(dalph_soil-1.)
     pfluxgrd(:) = 0.
     ! calcul de l'inertie thermique a partir de la variable rnat.
     ! on initialise a iice meme au-dessus d'un point de mer au cas
@@ -248,6 +247,13 @@
             lambda*zc(ig,1,indice)-ptsrf(ig))/ptimestep
     END DO
 
+  contains
+
+    real function fz(rk)
+      real rk
+      fz = fz1*(dalph_soil**rk-1.)/(dalph_soil-1.)
+    end function fz
+
   END SUBROUTINE soil
 
 end module soil_m