--- trunk/phylmd/clqh.f	2014/05/13 17:23:16	98
+++ trunk/phylmd/clqh.f	2014/09/09 12:54:30	106
@@ -5,12 +5,11 @@
 contains
 
   SUBROUTINE clqh(dtime, itime, jour, debut, rlat, knon, nisurf, knindex, &
-       pctsrf, soil_model, tsoil, qsol, ok_veget, ocean, rmu0, co2_ppm, &
-       rugos, rugoro, u1lay, v1lay, coef, t, q, ts, paprs, pplay, delp, &
-       radsol, albedo, alblw, snow, qsurf, precip_rain, precip_snow, fder, &
-       swnet, fluxlat, pctsrf_new, agesno, d_t, d_q, d_ts, z0_new, flux_t, &
-       flux_q, dflux_s, dflux_l, fqcalving, ffonte, run_off_lic_0, flux_o, &
-       flux_g, tslab, seaice)
+       pctsrf, tsoil, qsol, rmu0, co2_ppm, rugos, rugoro, u1lay, &
+       v1lay, coef, t, q, ts, paprs, pplay, delp, radsol, albedo, alblw, &
+       snow, qsurf, precip_rain, precip_snow, fder, swnet, fluxlat, &
+       pctsrf_new, agesno, d_t, d_q, d_ts, z0_new, flux_t, flux_q, dflux_s, &
+       dflux_l, fqcalving, ffonte, run_off_lic_0, flux_o, flux_g)
 
     ! Author: Z. X. Li (LMD/CNRS)
     ! Date: 1993/08/18
@@ -37,7 +36,7 @@
 
     REAL t(klon, klev)       ! temperature (K)
     REAL q(klon, klev)       ! humidite specifique (kg/kg)
-    REAL ts(klon)           ! temperature du sol (K)
+    REAL, intent(in):: ts(klon) ! temperature du sol (K)
     REAL evap(klon)         ! evaporation au sol
     REAL paprs(klon, klev+1) ! pression a inter-couche (Pa)
     REAL pplay(klon, klev)   ! pression au milieu de couche (Pa)
@@ -47,23 +46,27 @@
     REAL alblw(klon)
     REAL snow(klon)         ! hauteur de neige
     REAL qsurf(klon)         ! humidite de l'air au dessus de la surface
-    real precip_rain(klon), precip_snow(klon)
+
+    real, intent(in):: precip_rain(klon)
+    ! liquid water mass flux (kg/m2/s), positive down
+
+    real, intent(in):: precip_snow(klon)
+    ! solid water mass flux (kg/m2/s), positive down
+
     REAL agesno(klon)
     REAL rugoro(klon)
     REAL run_off_lic_0(klon)! runof glacier au pas de temps precedent
     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)
+    integer, intent(in):: knindex(:) ! (knon)
     real, intent(in):: pctsrf(klon, nbsrf)
     real, intent(in):: rlat(klon)
-    logical ok_veget 
-    REAL co2_ppm            ! taux CO2 atmosphere
-    character(len=*), intent(in):: ocean
+    REAL, intent(in):: co2_ppm            ! taux CO2 atmosphere
 
     REAL d_t(klon, klev)     ! incrementation de "t"
     REAL d_q(klon, klev)     ! incrementation de "q"
-    REAL d_ts(klon)         ! incrementation de "ts"
+    REAL, intent(out):: d_ts(:) ! (knon) incrementation de "ts"
     REAL flux_t(klon, klev)  ! (diagnostic) flux de la chaleur
     !                               sensible, flux de Cp*T, positif vers
     !                               le bas: j/(m**2 s) c.a.d.: W/m2
@@ -76,11 +79,13 @@
     ! Flux d'eau "perdue" par la surface et nécessaire pour que limiter la
     ! hauteur de neige, en kg/m2/s
     REAL fqcalving(klon)
+
     !IM "slab" ocean
-    REAL tslab(klon)  !temperature du slab ocean (K) (OCEAN='slab  ')
-    REAL seaice(klon) ! glace de mer en kg/m2
-    REAL flux_o(klon) ! flux entre l'ocean et l'atmosphere W/m2
-    REAL flux_g(klon) ! flux entre l'ocean et la glace de mer W/m2
+
+    REAL, intent(out):: flux_o(klon) ! flux entre l'ocean et l'atmosphere W/m2
+
+    REAL, intent(out):: flux_g(klon) 
+    ! flux entre l'ocean et la glace de mer W/m2
 
     INTEGER i, k
     REAL zx_cq(klon, klev)
@@ -116,13 +121,14 @@
     real swnet(klon), swdown(klon)
     real p1lay(klon)
     !$$$C PB ajout pour soil
-    LOGICAL, intent(in):: soil_model
     REAL tsoil(klon, nsoilmx)
-    REAL qsol(klon)
+
+    REAL, intent(inout):: qsol(klon)
+    ! column-density of water in soil, in kg m-2
 
     ! Parametres de sortie
     real fluxsens(klon), fluxlat(klon)
-    real tsurf_new(klon), alb_new(klon)
+    real tsurf_new(knon), alb_new(klon)
     real z0_new(klon)
     real pctsrf_new(klon, nbsrf)
     ! JLD
@@ -273,16 +279,13 @@
     endif
     ccanopy = co2_ppm
 
-    CALL interfsurf_hq(itime, dtime, jour, rmu0, iim, jjm, &
-         nisurf, knon, knindex, pctsrf, rlat, debut, &
-         ok_veget, soil_model, nsoilmx, tsoil, qsol,  u1lay, v1lay, &
-         temp_air, spechum, tq_cdrag, petAcoef, peqAcoef, &
-         petBcoef, peqBcoef, precip_rain, precip_snow, &
-         fder, rugos, rugoro, &
-         snow, qsurf, ts, p1lay, psref, radsol, ocean, &
-         evap, fluxsens, fluxlat, dflux_l, dflux_s, tsurf_new, &
-         alb_new, alblw, z0_new, pctsrf_new, agesno, fqcalving, &
-         ffonte, run_off_lic_0, flux_o, flux_g, tslab, seaice)
+    CALL interfsurf_hq(itime, dtime, jour, rmu0, nisurf, knon, knindex, &
+         pctsrf, rlat, debut, nsoilmx, tsoil, qsol, u1lay, v1lay, temp_air, &
+         spechum, tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef, &
+         precip_rain, precip_snow, fder, rugos, rugoro, snow, qsurf, &
+         ts(:knon), p1lay, psref, radsol, evap, fluxsens, fluxlat, dflux_l, &
+         dflux_s, tsurf_new, alb_new, alblw, z0_new, pctsrf_new, agesno, &
+         fqcalving, ffonte, run_off_lic_0, flux_o, flux_g)
 
     do i = 1, knon
        flux_t(i, 1) = fluxsens(i)