--- trunk/Sources/phylmd/phyetat0.f	2016/09/01 10:30:53	207
+++ trunk/phylmd/phyetat0.f	2018/07/12 14:49:20	276
@@ -9,6 +9,7 @@
   ! by a simple index, in degrees
 
   integer, save:: itau_phy
+  REAL, save:: zmasq(KLON) ! fraction of land
 
   private klon
 
@@ -25,7 +26,7 @@
     ! Objet : lecture de l'état initial pour la physique
 
     USE conf_gcm_m, ONLY: raz_date
-    use dimphy, only: zmasq, klev
+    use dimphy, only: klev
     USE dimsoil, ONLY : nsoilmx
     USE indicesol, ONLY : epsfra, is_lic, is_oce, is_sic, is_ter, nbsrf
     use netcdf, only: nf90_global, nf90_inq_varid, NF90_NOERR, NF90_NOWRITE
@@ -36,7 +37,10 @@
     REAL, intent(out):: ftsol(klon, nbsrf)
     REAL, intent(out):: ftsoil(klon, nsoilmx, nbsrf)
     REAL, intent(out):: qsurf(klon, nbsrf)
-    REAL, intent(out):: qsol(:) ! (klon)
+
+    REAL, intent(out):: qsol(:)
+    ! (klon) column-density of water in soil, in kg m-2
+
     REAL, intent(out):: snow(klon, nbsrf)
     REAL, intent(out):: albe(klon, nbsrf)
     REAL, intent(out):: evap(klon, nbsrf)
@@ -179,8 +183,6 @@
     call NF95_INQ_VARID(ncid_startphy, "QS", varid)
     call nf95_get_var(ncid_startphy, varid, qsurf)
 
-    ! Eau dans le sol (pour le modele de sol "bucket")
-
     ierr = NF90_INQ_VARID(ncid_startphy, "QSOL", varid)
     IF (ierr == NF90_NOERR) THEN
        call nf95_get_var(ncid_startphy, varid, qsol)