--- trunk/libf/dyn3d/ce0l.f90	2012/11/14 16:59:30	68
+++ trunk/dyn3d/ce0l.f	2018/12/10 15:54:30	313
@@ -2,21 +2,36 @@
 
   ! This program sets the initial and boundary values.
 
+  use comdissnew, only: read_comdissnew
   use conf_gcm_m, only: conf_gcm
-  use etat0_mod, only: etat0
-  use jumble, only: new_unit
+  use dimensions, only: iim, jjm
+  use dimphy, only: klon
+  use dynetat0_chosen_m, only: read_serre
+  use etat0_m, only: etat0
+  use grilles_gcm_netcdf_sub_m, only: grilles_gcm_netcdf_sub
+  use indicesol, only: nbsrf
   use limit_mod, only: limit
-  use unit_nml_m, only: unit_nml
+  use unit_nml_m, only: unit_nml, set_unit_nml
 
   implicit none
 
+  REAL phis(iim + 1, jjm + 1) ! surface geopotential, in m2 s-2
+
+  REAL pctsrf(klon, nbsrf)
+  ! ("pctsrf(i, :)" is the composition of the surface at horizontal
+  ! position "i")
+
   !-------------------------------------
 
-  call new_unit(unit_nml)
-  open(unit_nml, file="used_namelists", status="replace", action="write")
+  call set_unit_nml
+  open(unit_nml, file="used_namelists.txt", status="replace", action="write")
   CALL conf_gcm
-  CALL etat0
-  CALL limit
+  call read_comdissnew
+  call read_serre
+  CALL etat0(phis, pctsrf)
+  CALL limit(pctsrf)
+  CALL grilles_gcm_netcdf_sub(phis)
   close(unit_nml)
+  print *, "ce0l: done"
 
 END PROGRAM ce0l