Index: trunk/LMDZ.COMMON/libf/evolution/changelog.txt
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/changelog.txt	(revision 3711)
+++ trunk/LMDZ.COMMON/libf/evolution/changelog.txt	(revision 3738)
@@ -625,16 +625,16 @@
 Updating the average pressure at each iteration.
 
-== 1/4/2025 == EV
+== 01/04/2025 == EV
 updating paramaters to work with the layering
 
-== 2/4/2025 == EV
+== 02/04/2025 == EV
 updating the h2o and co2 ice values that will be used by the PCM according to the top layer in the stratum
 updating h2o ice recomp tendency to reduce the ice table loss rate due increase of depth and a decrease in the humidity around the ice.
 
-== 3/4/2025 == EV
+== 03/04/2025 == EV
 adding the growth of the lag layer
 changing the stopping critiria to take into account the top startum layer
 removing the case where if we have 0 ice the PEM stops.
 
-
-
+== 25/04/2025 == JBC
+Resulting change due to revision r3737 in the 1D Mars PCM.
Index: trunk/LMDZ.COMMON/libf/evolution/pem.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/pem.F90	(revision 3711)
+++ trunk/LMDZ.COMMON/libf/evolution/pem.F90	(revision 3738)
@@ -111,5 +111,4 @@
     use comvert_mod,              only: ap, bp
     use writerestart1D_mod,       only: writerestart1D
-   ! use layering_mod,             only: layering, array2stratif, nb_str_max, layering_algo
 #endif
 
@@ -287,5 +286,5 @@
     real, dimension(:),     allocatable :: zqsat
     real, dimension(:,:,:), allocatable :: dq, dqdyn
-    real, dimension(nlayer)             :: play, w
+    real, dimension(nlayer)             :: play, w, q_def_h2o_vap
     real, dimension(nlayer + 1)         :: plev
 #else
@@ -323,5 +322,5 @@
         endif
     enddo
-    
+
     if (num_str) then ! This is a numeric sting so we considerer this is the job id
            ! Execute the system command
@@ -436,5 +435,5 @@
     call init_testphys1d('start1D.txt','startfi.nc',therestart1D,therestartfi,ngrid,nlayer,610.,nq,q,         &
                          time_0,pdyn(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)
+                         play,plev,latitude,longitude,cell_area,atm_wat_profile,atm_wat_tau,q_def_h2o_vap)
     nsplit_phys = 1
     day_step = steps_per_sol
@@ -1302,5 +1301,5 @@
 
 write(*,*)
-write(*,*) '********* PEM finalization *********'
+write(*,*) '****** PEM final information *******'
 write(*,'(a,f10.2,a)') " + The PEM leg has run for ", i_myear_leg, " Martian years."
 write(*,'(a,f10.2,a,f10.2,a)') " + The chained simulation has run for ", i_myear, " Martian years =", i_myear*convert_years, " Earth years."