--- trunk/libf/phylmd/clmain.f	2008/02/27 13:16:39	3
+++ trunk/libf/phylmd/clmain.f	2008/07/21 16:05:07	12
@@ -116,14 +116,14 @@
 c======================================================================
 c$$$ PB ajout pour soil
 c
-      REAL dtime
+      REAL, intent(in):: dtime
       real date0
-      integer itap
+      integer, intent(in):: itap
       REAL t(klon,klev), q(klon,klev)
       REAL u(klon,klev), v(klon,klev)
 cIM 230604 BAD  REAL radsol(klon) ??? 
       REAL, intent(in):: paprs(klon,klev+1)
-      real pplay(klon,klev)
+      real, intent(in):: pplay(klon,klev)
       REAL, intent(in):: rlon(klon), rlat(klon)
       real cufi(klon), cvfi(klon)
       REAL d_t(klon, klev), d_q(klon, klev)
@@ -149,7 +149,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 +176,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
@@ -356,7 +356,6 @@
       ENDDO
       ENDDO
       DO i = 1, klon  ! vent de la premiere couche
-ccc         zx_alf1 = (paprs(i,1)-pplay(i,2))/(pplay(i,1)-pplay(i,2))
          zx_alf1 = 1.0
          zx_alf2 = 1.0 - zx_alf1
          u1lay(i) = u(i,1)*zx_alf1 + u(i,2)*zx_alf2
@@ -658,7 +657,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
@@ -912,7 +911,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
@@ -936,7 +934,6 @@
          ENDIF
         ENDDO
        ENDIF !nsrf.EQ.is_sic
-      ENDIF !OCEAN
 c
       IF(OCEAN.EQ.'slab  ') THEN
        IF(nsrf.EQ.is_oce) then