Index: /LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F
===================================================================
--- /LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F	(revision 443)
+++ /LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F	(revision 444)
@@ -267,6 +267,9 @@
       SAVE deltat                 ! ecart avec la SST de reference
 c
-      REAL fqsol(klon,nbsrf)
-      SAVE fqsol                  ! humidite du sol
+      REAL fqsurf(klon,nbsrf)
+      SAVE fqsurf                 ! humidite de l'air au contact de la surface
+c
+      REAL qsol(klon)
+      SAVE qsol                  ! hauteur d'eau dans le sol
 c
       REAL fsnow(klon,nbsrf)
@@ -502,5 +505,5 @@
       REAL cldt(klon),cldq(klon) !nuage total, eau liquide integree
 c
-      REAL zxtsol(klon), zxqsol(klon), zxsnow(klon), zxfluxlat(klon)
+      REAL zxtsol(klon), zxqsurf(klon), zxsnow(klon), zxfluxlat(klon)
 c
       REAL dist, rmu0(klon), fract(klon)
@@ -746,5 +749,5 @@
 c
          CALL phyetat0 ("startphy.nc",dtime,co2_ppm_etat0,solaire_etat0,
-     .       rlat,rlon,pctsrf, ftsol,ftsoil,deltat,fqsol,fsnow,
+     .       rlat,rlon,pctsrf, ftsol,ftsoil,deltat,fqsurf,fqsol,fsnow,
      .       falbe, fevap, rain_fall,snow_fall,solsw, sollwdown,
      .       dlw,radsol,frugs,agesno,clesphy0,
@@ -1099,6 +1102,6 @@
      e            julien, rmu0, 
      e            ok_veget, ocean, npas, nexca, ftsol,
-     $            soil_model,ftsoil,
-     $            paprs,pplay,radsol, fsnow,fqsol,fevap,falbe,falblw,
+     $            soil_model,ftsoil, qsol, 
+     $            paprs,pplay,radsol, fsnow,fqsurf,fevap,falbe,falblw,
      $            fluxlat,
 cIM cf. JLD  e            rain_fall, snow_fall, solsw, sollw, sollwdown, fder,
@@ -1186,5 +1189,5 @@
 cIM cf. JLD
             wfbils(i,nsrf) = ( fsolsw(i,nsrf) + fsollw(i,nsrf)
-     $         + fluxt(i,nsrf) + fluxlat(i,nsrf) ) * pctsrf(i,nsrf)
+     $         + fluxt(i,1,nsrf) + fluxlat(i,nsrf) ) * pctsrf(i,nsrf)
             zxtsol(i) = zxtsol(i) + ftsol(i,nsrf)*pctsrf(i,nsrf)
             zxfluxlat(i) = zxfluxlat(i) + fluxlat(i,nsrf)*pctsrf(i,nsrf)
@@ -1718,13 +1721,13 @@
 c
 c      CALL hydrol(dtime,pctsrf,rain_fall, snow_fall, zxevap,
-c     .            agesno, ftsol,fqsol,fsnow, ruis)
-c
-      DO i = 1, klon
-         zxqsol(i) = 0.0
+c     .            agesno, ftsol,fqsurf,fsnow, ruis)
+c
+      DO i = 1, klon
+         zxqsurf(i) = 0.0
          zxsnow(i) = 0.0
       ENDDO
       DO nsrf = 1, nbsrf
       DO i = 1, klon
-         zxqsol(i) = zxqsol(i) + fqsol(i,nsrf)*pctsrf(i,nsrf)
+         zxqsurf(i) = zxqsurf(i) + fqsurf(i,nsrf)*pctsrf(i,nsrf)
          zxsnow(i) = zxsnow(i) + fsnow(i,nsrf)*pctsrf(i,nsrf)
       ENDDO
@@ -1736,5 +1739,5 @@
 cXXX      DO i = 1, klon
 cXXX         IF (pctsrf(i,nsrf).LT.epsfra) THEN
-cXXX            fqsol(i,nsrf) = zxqsol(i)
+cXXX            fqsurf(i,nsrf) = zxqsurf(i)
 cXXX            fsnow(i,nsrf) = zxsnow(i)
 cXXX         ENDIF
@@ -1993,6 +1996,6 @@
 ccc         IF (ok_oasis) CALL quitcpl
          CALL phyredem ("restartphy.nc",dtime,radpas,
-     .      rlat, rlon, pctsrf, ftsol, ftsoil, deltat, fqsol, fsnow,
-     .      falbe, fevap, rain_fall, snow_fall,
+     .      rlat, rlon, pctsrf, ftsol, ftsoil, deltat, fqsurf, qsol,
+     .      fsnow, falbe, fevap, rain_fall, snow_fall,
      .      solsw, sollwdown,dlw,
      .      radsol,frugs,agesno,