--- trunk/libf/phylmd/clqh.f90	2011/08/24 11:43:14	49
+++ trunk/libf/phylmd/clqh.f90	2013/06/24 15:39:52	70
@@ -13,7 +13,7 @@
        z0_new, flux_t, flux_q, dflux_s, dflux_l, fqcalving, ffonte, &
        run_off_lic_0, flux_o, flux_g, tslab, seaice)
 
-    ! Author: Z.X. Li (LMD/CNRS)!
+    ! Author: Z. X. Li (LMD/CNRS)
     ! Date: 1993/08/18
     ! Objet : diffusion verticale de "q" et de "h"
 
@@ -22,19 +22,21 @@
     USE dimphy, ONLY : klev, klon, zmasq
     USE dimsoil, ONLY : nsoilmx
     USE indicesol, ONLY : is_ter, nbsrf
-    USE interface_surf, ONLY : interfsurf_hq
+    USE interfsurf_hq_m, ONLY : interfsurf_hq
     USE suphec_m, ONLY : rcpd, rd, rg, rkappa
 
     ! Arguments:
-    INTEGER knon
+    INTEGER, intent(in):: knon
     REAL, intent(in):: dtime              ! intervalle du temps (s)
-    real date0
+    real, intent(in):: date0
     REAL u1lay(klon)        ! vitesse u de la 1ere couche (m/s)
     REAL v1lay(klon)        ! vitesse v de la 1ere couche (m/s)
-    REAL coef(klon, klev)    ! le coefficient d'echange (m**2/s)
-    !                               multiplie par le cisaillement du 
-    !                               vent (dV/dz); la premiere valeur
-    !                               indique la valeur de Cdrag (sans unite)
+
+    REAL, intent(in):: coef(:, :) ! (knon, klev)
+    ! Le coefficient d'echange (m**2/s) multiplie par le cisaillement
+    ! du vent (dV/dz). La premiere valeur indique la valeur de Cdrag
+    ! (sans unite).
+
     REAL t(klon, klev)       ! temperature (K)
     REAL q(klon, klev)       ! humidite specifique (kg/kg)
     REAL ts(klon)           ! temperature du sol (K)
@@ -51,11 +53,11 @@
     REAL agesno(klon)
     REAL rugoro(klon)
     REAL run_off_lic_0(klon)! runof glacier au pas de temps precedent
-    integer jour            ! jour de l'annee en cours
-    real rmu0(klon)         ! cosinus de l'angle solaire zenithal
+    integer, intent(in):: jour            ! jour de l'annee en cours
+    real, intent(in):: rmu0(klon)         ! cosinus de l'angle solaire zenithal
     real rugos(klon)        ! rugosite
     integer knindex(klon)
-    real pctsrf(klon, nbsrf)
+    real, intent(in):: pctsrf(klon, nbsrf)
     real, intent(in):: rlon(klon), rlat(klon)
     real cufi(klon), cvfi(klon)
     logical ok_veget 
@@ -268,28 +270,22 @@
     peqBcoef = 0.
     p1lay =0.
 
-    !      do i = 1, knon
     petAcoef(1:knon) = zx_ch(1:knon, 1)
     peqAcoef(1:knon) = zx_cq(1:knon, 1)
     petBcoef(1:knon) =  zx_dh(1:knon, 1)
     peqBcoef(1:knon) = zx_dq(1:knon, 1)
-    tq_cdrag(1:knon) =coef(1:knon, 1)
+    tq_cdrag(1:knon) =coef(:knon, 1)
     temp_air(1:knon) =t(1:knon, 1)
     epot_air(1:knon) =local_h(1:knon, 1)
     spechum(1:knon)=q(1:knon, 1)
     p1lay(1:knon) = pplay(1:knon, 1)
     zlev1(1:knon) = delp(1:knon, 1)
-    !        swnet = swdown * (1. - albedo)
 
-    !IM swdown=flux SW incident sur terres
-    !IM swdown=flux SW net sur les autres surfaces 
-    !IM     swdown(1:knon) = swnet(1:knon)
     if(nisurf.eq.is_ter) THEN 
        swdown(1:knon) = swnet(1:knon)/(1-albedo(1:knon)) 
     else 
        swdown(1:knon) = swnet(1:knon)
     endif
-    !      enddo
     ccanopy = co2_ppm
 
     CALL interfsurf_hq(itime, dtime, date0, jour, rmu0, &