--- trunk/libf/phylmd/clmain.f	2008/04/18 14:45:53	10
+++ trunk/libf/phylmd/clmain.f	2008/07/25 19:59:34	13
@@ -116,7 +116,7 @@
 c======================================================================
 c$$$ PB ajout pour soil
 c
-      REAL dtime
+      REAL, intent(in):: dtime
       real date0
       integer, intent(in):: itap
       REAL t(klon,klev), q(klon,klev)
@@ -141,7 +141,8 @@
       REAL run_off_lic_0(klon), y_run_off_lic_0(klon)
 
       REAL flux_u(klon,klev, nbsrf), flux_v(klon,klev, nbsrf)
-      REAL rugmer(klon), agesno(klon,nbsrf),rugoro(klon)
+      REAL rugmer(klon), agesno(klon,nbsrf)
+      real, intent(in):: rugoro(klon)
       REAL cdragh(klon), cdragm(klon)
       integer jour            ! jour de l'annee en cours
       real rmu0(klon)         ! cosinus de l'angle solaire zenithal
@@ -149,7 +150,7 @@
       LOGICAL, intent(in):: debut
       logical, intent(in):: lafin
       logical ok_veget
-      character*6 ocean
+      character(len=*), intent(IN):: ocean
       integer npas, nexca
 c
       REAL pctsrf(klon,nbsrf)
@@ -176,7 +177,7 @@
       REAL zv1(klon)
 cAA
 c$$$ PB ajout pour soil
-      LOGICAL soil_model
+      LOGICAL, intent(in):: soil_model
 cIM ajout seuils cdrm, cdrh
       REAL cdmmax, cdhmax
 cIM: 261103
@@ -253,9 +254,7 @@
       LOGICAL first_appel
       SAVE first_appel
       DATA first_appel/.true./
-      LOGICAL debugindex
-      SAVE debugindex
-      DATA debugindex/.false./
+      LOGICAL:: debugindex = .false.
       integer idayref
       REAL t2m(klon,nbsrf), q2m(klon,nbsrf)
       REAL u10m(klon,nbsrf), v10m(klon,nbsrf)
@@ -657,7 +656,7 @@
          call ustarhb(knon,yu,yv,y_cd_m, yustar)
 
         if (prt_level > 9) THEN
-          WRITE(lunout,*)'USTAR = ',yustar
+          print *,'USTAR = ',yustar
         ENDIF
 
 c   iflag_pbl peut etre utilise comme longuer de melange
@@ -911,7 +910,6 @@
          enddo
       enddo
 cIM "slab" ocean 
-      IF(OCEAN.EQ.'slab  '.OR.OCEAN.EQ.'force ') THEN
        IF (nsrf.EQ.is_oce) THEN
         DO j = 1, knon
 c on projette sur la grille globale
@@ -935,7 +933,6 @@
          ENDIF
         ENDDO
        ENDIF !nsrf.EQ.is_sic
-      ENDIF !OCEAN
 c
       IF(OCEAN.EQ.'slab  ') THEN
        IF(nsrf.EQ.is_oce) then