Index: trunk/LMDZ.COMMON/libf/evolution/pem.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/pem.F90	(revision 3123)
+++ trunk/LMDZ.COMMON/libf/evolution/pem.F90	(revision 3129)
@@ -258,5 +258,5 @@
 
     ! Dummy variables to use the subroutine 'init_testphys1d'
-    logical                             :: startfiles_1D, therestart1D, therestartfi
+    logical                             :: therestart1D, therestartfi
     integer                             :: ndt, day0
     real                                :: ptif, pks, day, gru, grv, atm_wat_profile, atm_wat_tau
@@ -308,6 +308,20 @@
     allocate(q(1,llm,nqtot))
     allocate(longitude(1),latitude(1),cell_area(1))
-    call init_testphys1d(.true.,ngrid,nlayer,610.,nq,q,time_0,ps(1),ucov,vcov,teta,startfiles_1D,therestart1D, &
-                         therestartfi,ndt,ptif,pks,dtphys,zqsat,dq,dqdyn,day0,day,gru,grv,w,                   &
+
+    therestart1D = .false.
+    inquire(file = 'start1D_evol.txt',exist = therestart1D)
+    if (.not. therestart1D) then
+        write(*,*) 'There is no "start1D_evol.txt" file!'
+        error stop 'Initialization cannot be done for the 1D PEM.'
+    endif
+    therestartfi = .false.
+    inquire(file = 'startfi_evol.nc',exist = therestartfi)
+    if (.not. therestartfi) then
+        write(*,*) 'There is no "startfi_evol.nc" file!'
+        error stop 'Initialization cannot be done for the 1D PEM.'
+    endif
+
+    call init_testphys1d('start1D_evol.txt','startfi_evol.nc',.true.,therestart1D,therestartfi,ngrid,nlayer,610., &
+                         nq,q,time_0,ps(1),ucov,vcov,teta,ndt,ptif,pks,dtphys,zqsat,dq,dqdyn,day0,day,gru,grv,w,  &
                          play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau)
     ps(2) = ps(1)