--- trunk/phylmd/Interface_surf/interfsurf_hq.f	2014/09/04 10:05:52	104
+++ trunk/phylmd/Interface_surf/interfsurf_hq.f	2014/12/04 16:35:03	116
@@ -39,7 +39,7 @@
     integer, intent(IN):: nisurf ! index de la surface a traiter
     integer, intent(IN):: knon ! nombre de points de la surface a traiter
 
-    integer, intent(in):: knindex(klon)
+    integer, intent(in):: knindex(:) ! (knon)
     ! index des points de la surface a traiter
 
     real, intent(IN):: pctsrf(klon, nbsrf)
@@ -90,8 +90,7 @@
     ! ps pression au sol
     REAL, DIMENSION(klon), INTENT(INOUT):: radsol
     ! radsol rayonnement net aus sol (LW + SW)
-    real, dimension(klon), intent(INOUT):: evap
-    ! evap evaporation totale
+    real, intent(INOUT):: evap(klon) ! evaporation totale
     real, dimension(klon), intent(OUT):: fluxsens, fluxlat
     ! fluxsens flux de chaleur sensible
     ! fluxlat flux de chaleur latente
@@ -207,8 +206,8 @@
 
        ! calcul albedo: lecture albedo fichier boundary conditions
        ! puis ajout albedo neige
-       call interfsur_lim(itime, dtime, jour, nisurf, knon, knindex, &
-            debut, alb_new, z0_new)
+       call interfsur_lim(itime, dtime, jour, nisurf, knindex, debut, &
+            alb_new, z0_new)
 
        ! calcul snow et qsurf, hydrol adapté
        CALL calbeta(nisurf, snow(:knon), qsol(:knon), beta(:knon), &
@@ -221,16 +220,19 @@
        ELSE
           cal = RCPD * capsol
        ENDIF
-       CALL calcul_fluxs(nisurf, dtime, tsurf, p1lay(:knon), cal(:knon), beta(:knon), &
-            tq_cdrag(:knon), ps(:knon), qsurf(:knon), radsol(:knon), &
-            dif_grnd(:knon), temp_air(:knon), spechum(:knon), u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), &
-            petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap(:knon), fluxlat(:knon), fluxsens(:knon), dflux_s(:knon), &
-            dflux_l(:knon))
+       CALL calcul_fluxs(nisurf, dtime, tsurf, p1lay(:knon), cal(:knon), &
+            beta(:knon), tq_cdrag(:knon), ps(:knon), qsurf(:knon), &
+            radsol(:knon), dif_grnd(:knon), temp_air(:knon), spechum(:knon), &
+            u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), &
+            petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap(:knon), &
+            fluxlat(:knon), fluxsens(:knon), dflux_s(:knon), dflux_l(:knon))
 
        CALL fonte_neige(nisurf, dtime, tsurf, p1lay(:knon), beta(:knon), &
-            tq_cdrag(:knon), ps(:knon), precip_rain(:knon), precip_snow(:knon), snow(:knon), qsol(:knon), &
-            temp_air(:knon), spechum(:knon), u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), &
-            peqBcoef(:knon), tsurf_new, evap(:knon), fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon))
+            tq_cdrag(:knon), ps(:knon), precip_rain(:knon), &
+            precip_snow(:knon), snow(:knon), qsol(:knon), temp_air(:knon), &
+            spechum(:knon), u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), &
+            peqAcoef(:knon), petBcoef(:knon), peqBcoef(:knon), tsurf_new, &
+            evap(:knon), fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon))
 
        call albsno(klon, knon, dtime, agesno, alb_neig, precip_snow)
        where (snow(1 : knon) < 0.0001) agesno(1 : knon) = 0.
@@ -245,8 +247,8 @@
     case (is_oce)
        ! Surface "ocean" appel à l'interface avec l'océan
        ! lecture conditions limites
-       call interfoce_lim(itime, dtime, jour, klon, knon, knindex, debut, &
-            tsurf_temp, pctsrf_new)
+       call interfoce_lim(itime, dtime, jour, knindex, debut, tsurf_temp, &
+            pctsrf_new)
 
        cal = 0.
        beta = 1.
@@ -287,8 +289,8 @@
        ! Surface "glace de mer" appel a l'interface avec l'ocean
 
        ! ! lecture conditions limites
-       CALL interfoce_lim(itime, dtime, jour, klon, knon, knindex, &
-            debut, tsurf_new, pctsrf_new)
+       CALL interfoce_lim(itime, dtime, jour, knindex, debut, tsurf_new, &
+            pctsrf_new)
 
        DO ii = 1, knon
           tsurf_new(ii) = tsurf(ii)
@@ -318,9 +320,10 @@
 
        CALL calcul_fluxs(nisurf, dtime, tsurf_temp, p1lay(:knon), cal(:knon), &
             beta(:knon), tq_cdrag(:knon), ps(:knon), qsurf(:knon), &
-            radsol(:knon), dif_grnd(:knon), temp_air(:knon), spechum(:knon), u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), &
-            peqAcoef(:knon), petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap(:knon), fluxlat(:knon), fluxsens(:knon), &
-            dflux_s(:knon), dflux_l(:knon))
+            radsol(:knon), dif_grnd(:knon), temp_air(:knon), spechum(:knon), &
+            u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), &
+            petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap(:knon), &
+            fluxlat(:knon), fluxsens(:knon), dflux_s(:knon), dflux_l(:knon))
 
        !IM: flux entre l'ocean et la glace de mer pour le "slab" ocean
        DO i = 1, knon
@@ -330,9 +333,11 @@
        ENDDO
 
        CALL fonte_neige(nisurf, dtime, tsurf_temp, p1lay(:knon), beta(:knon), &
-            tq_cdrag(:knon), ps(:knon), precip_rain(:knon), precip_snow(:knon), snow(:knon), qsol(:knon), &
-            temp_air(:knon), spechum(:knon), u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), &
-            peqBcoef(:knon), tsurf_new, evap(:knon), fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon))
+            tq_cdrag(:knon), ps(:knon), precip_rain(:knon), &
+            precip_snow(:knon), snow(:knon), qsol(:knon), temp_air(:knon), &
+            spechum(:knon), u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), &
+            peqAcoef(:knon), petBcoef(:knon), peqBcoef(:knon), tsurf_new, &
+            evap(:knon), fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon))
 
        ! calcul albedo
 
@@ -372,16 +377,19 @@
        beta = 1.0
        dif_grnd = 0.0
 
-       call calcul_fluxs(nisurf, dtime, tsurf, p1lay(:knon), cal(:knon), beta(:knon), &
-            tq_cdrag(:knon), ps(:knon), qsurf(:knon), radsol(:knon), &
-            dif_grnd(:knon), temp_air(:knon), spechum(:knon), u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), &
-            petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap(:knon), fluxlat(:knon), fluxsens(:knon), dflux_s(:knon), &
-            dflux_l(:knon))
+       call calcul_fluxs(nisurf, dtime, tsurf, p1lay(:knon), cal(:knon), &
+            beta(:knon), tq_cdrag(:knon), ps(:knon), qsurf(:knon), &
+            radsol(:knon), dif_grnd(:knon), temp_air(:knon), spechum(:knon), &
+            u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), &
+            petBcoef(:knon), peqBcoef(:knon), tsurf_new, evap(:knon), &
+            fluxlat(:knon), fluxsens(:knon), dflux_s(:knon), dflux_l(:knon))
 
        call fonte_neige(nisurf, dtime, tsurf, p1lay(:knon), beta(:knon), &
-            tq_cdrag(:knon), ps(:knon), precip_rain(:knon), precip_snow(:knon), snow(:knon), qsol(:knon), &
-            temp_air(:knon), spechum(:knon), u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), peqAcoef(:knon), petBcoef(:knon), &
-            peqBcoef(:knon), tsurf_new, evap(:knon), fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon))
+            tq_cdrag(:knon), ps(:knon), precip_rain(:knon), &
+            precip_snow(:knon), snow(:knon), qsol(:knon), temp_air(:knon), &
+            spechum(:knon), u1_lay(:knon), v1_lay(:knon), petAcoef(:knon), &
+            peqAcoef(:knon), petBcoef(:knon), peqBcoef(:knon), tsurf_new, &
+            evap(:knon), fqcalving(:knon), ffonte(:knon), run_off_lic_0(:knon))
 
        ! calcul albedo
        CALL albsno(klon, knon, dtime, agesno, alb_neig, precip_snow)