--- trunk/dyn3d/etat0.f	2014/07/02 18:39:15	99
+++ trunk/dyn3d/etat0.f	2014/09/18 19:56:46	113
@@ -13,7 +13,7 @@
 
 contains
 
-  SUBROUTINE etat0
+  SUBROUTINE etat0(phis)
 
     ! From "etat0_netcdf.F", version 1.3, 2005/05/25 13:10:09
 
@@ -46,7 +46,6 @@
     use q_sat_m, only: q_sat
     use regr_lat_time_coefoz_m, only: regr_lat_time_coefoz
     use regr_pr_o3_m, only: regr_pr_o3
-    use serre, only: alphax
     use startdyn, only: start_init_dyn
     USE start_init_orog_m, only: start_init_orog, mask
     use start_init_phys_m, only: start_init_phys
@@ -54,7 +53,10 @@
     use temps, only: itau_phy, annee_ref, day_ref
     use test_disvert_m, only: test_disvert
 
-    ! Variables local to the procedure:
+    REAL, intent(out):: phis(:, :) ! (iim + 1, jjm + 1)
+    ! surface geopotential, in m2 s-2
+
+    ! Local:
 
     REAL latfi(klon), lonfi(klon)
     ! (latitude and longitude of a point of the scalar grid identified
@@ -80,7 +82,6 @@
     real seaice(klon) ! kg m-2
     REAL frugs(klon, nbsrf), agesno(klon, nbsrf)
     REAL rugmer(klon)
-    REAL phis(iim + 1, jjm + 1) ! surface geopotential, in m2 s-2
     real, dimension(iim + 1, jjm + 1):: zmea_2d, zstd_2d, zsig_2d, zgam_2d
     real, dimension(iim + 1, jjm + 1):: zthe_2d, zpic_2d, zval_2d
     real, dimension(iim + 1, jjm + 1):: tsol_2d, qsol_2d, ps
@@ -294,8 +295,6 @@
     ! Calcul interm\'ediaire :
     CALL massdair(p3d, masse)
 
-    print *, 'ALPHAX = ', alphax
-
     forall (l = 1:llm)
        masse(:, 1, l) = SUM(aire_2d(:iim, 1) * masse(:iim, 1, l)) / apoln
        masse(:, jjm + 1, l) = &