--- trunk/Sources/phylmd/Interface_surf/soil.f	2015/07/16 17:39:10	156
+++ trunk/Sources/phylmd/Interface_surf/soil.f	2015/07/20 16:01:49	157
@@ -85,12 +85,11 @@
     REAL zc(klon, nsoilmx, nbsrf), zd(klon, nsoilmx, nbsrf)
     REAL lambda
     SAVE dz1, dz2, zc, zd, lambda
-    LOGICAL firstcall, firstsurf(nbsrf)
-    SAVE firstcall, firstsurf
+    LOGICAL firstsurf(nbsrf)
+    SAVE firstsurf
     REAL isol, isno, iice
     SAVE isol, isno, iice
 
-    DATA firstcall/.TRUE./
     DATA firstsurf/.TRUE., .TRUE., .TRUE., .TRUE./
 
     DATA isol, isno, iice/2000., 2000., 2000./
@@ -131,8 +130,6 @@
        STOP 1
     END IF
 
-
-    ! $$$      IF (firstcall) THEN
     IF (firstsurf(indice)) THEN
 
        ! -----------------------------------------------------------------------
@@ -175,12 +172,11 @@
        END DO
        ! PB
        firstsurf(indice) = .FALSE.
-       ! $$$         firstcall =.false.
 
        ! Initialisations:
        ! ----------------
 
-    ELSE !--not firstcall
+    ELSE
        ! -----------------------------------------------------------------------
        ! Computation of the soil temperatures using the Cgrd and Dgrd
        ! coefficient computed at the previous time-step:
@@ -200,7 +196,7 @@
           END DO
        END DO
 
-    END IF !--not firstcall
+    END IF
     ! -----------------------------------------------------------------------
     ! Computation of the Cgrd and Dgrd coefficient for the next step:
     ! ---------------------------------------------------------------