Index: trunk/LMDZ.COMMON/libf/evolution/pem.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/pem.F90	(revision 3188)
+++ trunk/LMDZ.COMMON/libf/evolution/pem.F90	(revision 3189)
@@ -202,4 +202,5 @@
 real, dimension(:,:),     allocatable :: tsurf_ave                          ! Physic x SLOPE field: Averaged Surface Temperature [K]
 real, dimension(:,:,:),   allocatable :: tsoil_ave                          ! Physic x SOIL x SLOPE field: Averaged Soil Temperature [K]
+real, dimension(:,:,:),   allocatable :: tsoil_anom                         ! Amplitude between instataneous and yearly average soil temperature [K]
 real, dimension(:,:,:),   allocatable :: tsurf_PCM_timeseries               ! ngrid x SLOPE XTULES field: Surface Temperature in timeseries [K]
 real, dimension(:,:,:,:), allocatable :: tsoil_phys_PEM_timeseries          ! IG x SLOPE XTULES field: Non averaged Soil Temperature [K]
@@ -476,4 +477,5 @@
 
 allocate(tsoil_ave(ngrid,nsoilmx,nslope))
+allocate(tsoil_anom(ngrid,nsoilmx,nslope))
 allocate(watersoil_density_PEM_ave(ngrid,nsoilmx_PEM,nslope))
 allocate(vmr_co2_PCM(ngrid,timelen))
@@ -523,4 +525,5 @@
 call ini_ice_table_porefilling(ngrid,nslope)
 
+tsoil_anom(:,:,:) = tsoil(:,:,:) - tsoil_ave(:,:,:) ! compute anomaly between Tsoil(t) in the startfi - <Tsoil> to recompute properly tsoil in the restart
 if (soil_pem) then
     call soil_settings_PEM(ngrid,nslope,nsoilmx_PEM,nsoilmx,inertiesoil,TI_PEM)
@@ -985,5 +988,5 @@
 if (soil_pem) then
     call interpol_TI_PEM2PCM(ngrid,nslope,nsoilmx_PEM,nsoilmx,TI_PEM,inertiesoil)
-    tsoil = tsoil_phys_PEM_timeseries(:,1:nsoilmx,:,timelen)
+    tsoil = tsoil_PEM(:,1:nsoilmx,:) + tsoil_anom
 #ifndef CPP_STD
     flux_geo = fluxgeo