--- trunk/Sources/dyn3d/etat0.f	2016/02/05 16:02:34	175
+++ trunk/Sources/dyn3d/etat0.f	2017/03/22 13:40:27	214
@@ -20,7 +20,7 @@
     use caldyn0_m, only: caldyn0
     use comconst, only: cpp, kappa, iniconst
     use comgeom, only: aire_2d, apoln, apols, cu_2d, cv_2d, inigeom
-    use conf_gcm_m, only: nday, day_step, iphysiq
+    use conf_gcm_m, only: nday
     use dimens_m, only: iim, jjm, llm, nqmx
     use dimphy, only: zmasq
     use dimsoil, only: nsoilmx
@@ -33,7 +33,7 @@
     use fxhyp_m, only: fxhyp
     use fyhyp_m, only: fyhyp
     use geopot_m, only: geopot
-    use grid_atob, only: grille_m
+    use grille_m_m, only: grille_m
     use grid_change, only: init_dyn_phy, dyn_phy
     use indicesol, only: is_oce, is_sic, is_ter, is_lic, epsfra
     use iniadvtrac_m, only: iniadvtrac
@@ -43,8 +43,7 @@
     use netcdf95, only: nf95_close, nf95_get_var, nf95_gw_var, nf95_put_var, &
          nf95_inq_varid, nf95_open
     use nr_util, only: pi, assert
-    use paramet_m, only: ip1jm, ip1jmp1
-    use phyetat0_m, only: rlat, rlon
+    use phyetat0_m, only: rlat, rlon, itau_phy
     use phyredem0_m, only: phyredem0, ncid_restartphy
     use phyredem_m, only: phyredem
     use q_sat_m, only: q_sat
@@ -270,7 +269,7 @@
     ! Ad\'equation avec le maque terre/mer
     WHERE (pctsrf(:, is_lic) < EPSFRA) pctsrf(:, is_lic) = 0.
     WHERE (zmasq < EPSFRA) pctsrf(:, is_lic) = 0.
-    pctsrf(:, is_ter) = zmasq
+    where (zmasq <= EPSFRA) pctsrf(:, is_ter) = zmasq
     where (zmasq > EPSFRA)
        where (pctsrf(:, is_lic) >= zmasq)
           pctsrf(:, is_lic) = zmasq
@@ -336,7 +335,8 @@
     w01 = 0.
 
     nday = 0
-    call phyredem0(lmt_pas = day_step / iphysiq, itau_phy = 0)
+    itau_phy = 0 ! side effect
+    call phyredem0
 
     call nf95_inq_varid(ncid_restartphy, "trs", varid)
     call nf95_put_var(ncid_restartphy, varid, null_array)