Index: trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F
===================================================================
--- trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F	(revision 1533)
+++ trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F	(revision 1536)
@@ -13,5 +13,5 @@
       USE comgeomfi_h, only: lati, long, area
       use phyredem, only: physdem0,physdem1
-      use comgeomphy, only: initcomgeomphy
+      use comgeomphy, only: initcomgeomphy, airephy
       use slab_ice_h, only: noceanmx
       use planete_mod, only: apoastr,periastr,year_day,peri_day,
@@ -140,5 +140,6 @@
 c=======================================================================
 ! initialize "serial/parallel" related stuff
-      CALL init_phys_lmdz(iim,jjm+1,llm,1,(/(jjm-1)*iim+2/))
+!      CALL init_phys_lmdz(iim,jjm+1,llm,1,(/(jjm-1)*iim+2/))
+      CALL init_phys_lmdz(1,1,llm,1,(/1/))
       call initcomgeomphy
 
@@ -237,5 +238,5 @@
 ! OK. now that run.def has been read once -- any variable is in memory.
 ! so we can dump the dummy run.def
-      call system("rm -rf run.def")
+!      call system("rm -rf run.def") ! Ehouarn: delay this to after inifis
 
 ! while we're at it, check if there is a 'traceur.def' file
@@ -330,5 +331,5 @@
      !!! AREA. useless in 1D
       area(1)=1.E+0
-      aire(1)=area(1) !JL+EM to have access to the area in the diagfi.nc files. area in comgeomfi.h and aire in comgeom.h
+      airephy(1)=area(1) !JL+EM to have access to the area in the diagfi.nc files. area in comgeomfi.h and aire in comgeom.h
      !!! GEOPOTENTIAL. useless in 1D because control by surface pressure
       phisfi(1)=0.E+0
@@ -516,4 +517,9 @@
       CALL inifis(1,llm,nq,day0,daysec,nday,dtphys,
      .            latitude,longitude,area,rad,g,r,cpp)
+
+! At this point, both getin() and getin_p() functions have been used,
+! and the run.def file can be removed.
+      call system("rm -rf run.def")
+
 !!! We check everything is OK.
       PRINT *,"CHECK"
@@ -635,4 +641,7 @@
       write(*,*) " v = ",grv
 
+! To be clean, also set vertical winds to zero
+      w(1:nlayer)=0
+
 c     Initialisation des vents  au premier pas de temps
       DO ilayer=1,nlayer