--- trunk/phylmd/Interface_surf/clqh.f	2018/08/02 15:55:01	300
+++ trunk/phylmd/Interface_surf/clqh.f	2018/08/02 17:23:07	301
@@ -4,11 +4,11 @@
 
 contains
 
-  SUBROUTINE clqh(julien, debut, nisurf, knindex, tsoil, qsol, mu0, rugos, &
-       rugoro, u1lay, v1lay, coef, tq_cdrag, t, q, ts, paprs, pplay, delp, &
-       radsol, albedo, snow, qsurf, precip_rain, precip_snow, fluxlat, &
-       pctsrf_new_sic, agesno, d_t, d_q, d_ts, z0_new, flux_t, flux_q, &
-       dflux_s, dflux_l, fqcalving, ffonte, run_off_lic_0)
+  SUBROUTINE clqh(julien, nisurf, knindex, tsoil, qsol, mu0, rugos, rugoro, &
+       u1lay, v1lay, coef, tq_cdrag, t, q, ts, paprs, pplay, delp, radsol, &
+       albedo, snow, qsurf, precip_rain, precip_snow, fluxlat, pctsrf_new_sic, &
+       agesno, d_t, d_q, d_ts, z0_new, flux_t, flux_q, dflux_s, dflux_l, &
+       fqcalving, ffonte, run_off_lic_0, run_off_lic)
 
     ! Author: Z. X. Li (LMD/CNRS)
     ! Date: 1993 Aug. 18th
@@ -21,7 +21,6 @@
     USE suphec_m, ONLY: rkappa
 
     integer, intent(in):: julien ! jour de l'annee en cours
-    logical, intent(in):: debut
     integer, intent(in):: nisurf
     integer, intent(in):: knindex(:) ! (knon)
     REAL, intent(inout):: tsoil(:, :) ! (knon, nsoilmx)
@@ -98,6 +97,8 @@
     REAL, intent(inout):: run_off_lic_0(:) ! (knon)
     ! run-off glacier au pas de temps precedent
 
+    REAL, intent(OUT):: run_off_lic(:) ! (knon) ruissellement total
+
     ! Local:
 
     INTEGER k
@@ -112,12 +113,12 @@
     ! (La pression de r\'ef\'erence est celle au sol.)
 
     call climb_hq_down(pkf, cq, dq, ch, dh, paprs, pplay, t, coef, delp, q)
-    CALL interfsurf_hq(julien, mu0, nisurf, knindex, debut, tsoil, qsol, &
-         u1lay, v1lay, t(:, 1), q(:, 1), tq_cdrag, ch(:, 1), cq(:, 1), &
-         dh(:, 1), dq(:, 1), precip_rain, precip_snow, rugos, rugoro, snow, &
-         qsurf, ts, pplay(:, 1), paprs(:, 1), radsol, evap, flux_t, fluxlat, &
-         dflux_l, dflux_s, tsurf_new, albedo, z0_new, pctsrf_new_sic, agesno, &
-         fqcalving, ffonte, run_off_lic_0)
+    CALL interfsurf_hq(julien, mu0, nisurf, knindex, tsoil, qsol, u1lay, &
+         v1lay, t(:, 1), q(:, 1), tq_cdrag, ch(:, 1), cq(:, 1), dh(:, 1), &
+         dq(:, 1), precip_rain, precip_snow, rugos, rugoro, snow, qsurf, ts, &
+         pplay(:, 1), paprs(:, 1), radsol, evap, flux_t, fluxlat, dflux_l, &
+         dflux_s, tsurf_new, albedo, z0_new, pctsrf_new_sic, agesno, &
+         fqcalving, ffonte, run_off_lic_0, run_off_lic)
     flux_q = - evap
     d_ts = tsurf_new - ts
     call climb_hq_up(d_t, d_q, cq, dq, ch, dh, flux_t, flux_q, pkf, t, q)