--- trunk/phylmd/Interface_surf/interfsurf_hq.f	2018/07/20 16:46:48	282
+++ trunk/phylmd/Interface_surf/interfsurf_hq.f	2018/07/20 17:08:44	283
@@ -67,7 +67,7 @@
     real, intent(IN):: rugos(:) ! (knon) rugosite
     real, intent(IN):: rugoro(:) ! (knon) rugosite orographique
     real, intent(INOUT):: snow(:) ! (knon)
-    real, intent(INOUT):: qsurf(klon)
+    real, intent(OUT):: qsurf(:) ! (knon)
     real, intent(IN):: ts(:) ! (knon) temp\'erature de surface
     real, intent(IN):: p1lay(:) ! (knon) pression 1er niveau (milieu de couche)
     real, intent(IN):: ps(:) ! (knon) pression au sol
@@ -163,10 +163,10 @@
        CALL soil(dtime, is_ter, snow, ts, tsoil, soilcap, soilflux)
        cal = RCPD / soilcap
 
-       CALL calcul_fluxs(dtime, ts, p1lay, cal, beta, tq_cdrag, ps, &
-            qsurf(:knon), radsol + soilflux, dif_grnd(:knon), temp_air, &
-            spechum, u1_lay, v1_lay, petAcoef, peqAcoef, petBcoef, &
-            peqBcoef, tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l)
+       CALL calcul_fluxs(dtime, ts, p1lay, cal, beta, tq_cdrag, ps, qsurf, &
+            radsol + soilflux, dif_grnd(:knon), temp_air, spechum, u1_lay, &
+            v1_lay, petAcoef, peqAcoef, petBcoef, peqBcoef, tsurf_new, evap, &
+            fluxlat, flux_t, dflux_s, dflux_l)
        CALL fonte_neige(is_ter, dtime, precip_rain(:knon), precip_snow(:knon), &
             snow, qsol, tsurf_new, evap, fqcalving, ffonte(:knon), &
             run_off_lic_0(:knon))
@@ -183,10 +183,10 @@
        cal = 0.
        beta = 1.
        dif_grnd = 0.
-       call calcul_fluxs(dtime, tsurf, p1lay, cal, beta, tq_cdrag, ps, &
-            qsurf(:knon), radsol, dif_grnd(:knon), temp_air, spechum, u1_lay, &
-            v1_lay, petAcoef, peqAcoef, petBcoef, peqBcoef, tsurf_new, evap, &
-            fluxlat, flux_t, dflux_s, dflux_l)
+       call calcul_fluxs(dtime, tsurf, p1lay, cal, beta, tq_cdrag, ps, qsurf, &
+            radsol, dif_grnd(:knon), temp_air, spechum, u1_lay, v1_lay, &
+            petAcoef, peqAcoef, petBcoef, peqBcoef, tsurf_new, evap, fluxlat, &
+            flux_t, dflux_s, dflux_l)
        agesno = 0.
        albedo = alboc_cd(rmu0(knindex)) * fmagic
        z0_new = sqrt(rugos**2 + rugoro**2)
@@ -210,10 +210,10 @@
        tsurf = tsurf_new
        beta = 1.
 
-       CALL calcul_fluxs(dtime, tsurf, p1lay, cal, beta, tq_cdrag, ps, &
-            qsurf(:knon), radsol + soilflux, dif_grnd(:knon), temp_air, &
-            spechum, u1_lay, v1_lay, petAcoef, peqAcoef, petBcoef, &
-            peqBcoef, tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l)
+       CALL calcul_fluxs(dtime, tsurf, p1lay, cal, beta, tq_cdrag, ps, qsurf, &
+            radsol + soilflux, dif_grnd(:knon), temp_air, spechum, u1_lay, &
+            v1_lay, petAcoef, peqAcoef, petBcoef, peqBcoef, tsurf_new, evap, &
+            fluxlat, flux_t, dflux_s, dflux_l)
        CALL fonte_neige(is_sic, dtime, precip_rain(:knon), &
             precip_snow(:knon), snow, qsol, tsurf_new, evap, &
             fqcalving, ffonte(:knon), run_off_lic_0(:knon))
@@ -234,10 +234,10 @@
        beta = 1.
        dif_grnd = 0.
 
-       call calcul_fluxs(dtime, ts, p1lay, cal, beta, tq_cdrag, ps, &
-            qsurf(:knon), radsol + soilflux, dif_grnd(:knon), temp_air, &
-            spechum, u1_lay, v1_lay, petAcoef, peqAcoef, petBcoef, &
-            peqBcoef, tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l)
+       call calcul_fluxs(dtime, ts, p1lay, cal, beta, tq_cdrag, ps, qsurf, &
+            radsol + soilflux, dif_grnd(:knon), temp_air, spechum, u1_lay, &
+            v1_lay, petAcoef, peqAcoef, petBcoef, peqBcoef, tsurf_new, evap, &
+            fluxlat, flux_t, dflux_s, dflux_l)
        call fonte_neige(is_lic, dtime, precip_rain(:knon), &
             precip_snow(:knon), snow, qsol, tsurf_new, evap, &
             fqcalving, ffonte(:knon), run_off_lic_0(:knon))