Index: trunk/LMDZ.COMMON/libf/evolution/changelog.txt
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/changelog.txt	(revision 4074)
+++ trunk/LMDZ.COMMON/libf/evolution/changelog.txt	(revision 4075)
@@ -878,4 +878,7 @@
 == 17/02/2026 == JBC
 - Correct management of H2O ice tendency in 1D when there is not enough ice anymore.
-- Clean initialization of allocatable module arrays (especially needed when no slope)
+- Clean initialization of allocatable module arrays (especially needed when no slope).
 - One more renaming for consistency + few small updates thoughout the code.
+
+== 17/02/2026 == JBC
+Addition of runs ID for the current cycle in "pem_workflow.sts".
Index: trunk/LMDZ.COMMON/libf/evolution/display.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/display.F90	(revision 4074)
+++ trunk/LMDZ.COMMON/libf/evolution/display.F90	(revision 4075)
@@ -123,5 +123,5 @@
 write(stdout,*)
 write(stdout,*) '****** PEM final information *******'
-write(*,'(a,f16.4,a)')         ' + This PEM run '//int2str(i_pem_run)//' achieved ', n_yr_run, ' Planetary years in '//format_duration(dur_secs)//' (computing time).'
+write(*,'(a,f16.4,a)')         ' + The run PEM #'//int2str(i_pem_run)//' achieved ', n_yr_run, ' Planetary years, completed in '//format_duration(dur_secs)//'.'
 write(*,'(a,f16.4,a,f16.4,a)') ' + The workflow has achieved ', n_yr_sim, ' Planetary years =', n_yr_sim*r_plnt2earth_yr, ' Earth years.'
 write(*,'(a,f16.4,a)')         ' + The reached date is now ', (pem_ini_date + n_yr_sim)*r_plnt2earth_yr, ' Earth years.'
Index: trunk/LMDZ.COMMON/libf/evolution/workflow_status.F90
===================================================================
--- trunk/LMDZ.COMMON/libf/evolution/workflow_status.F90	(revision 4074)
+++ trunk/LMDZ.COMMON/libf/evolution/workflow_status.F90	(revision 4075)
@@ -121,7 +121,8 @@
 ! LOCAL VARIABLES
 ! ---------------
-logical(k4)   :: here
-integer(di)   :: cstat, ierr, funit
-character(20) :: fch1, fch2, fch3
+logical(k4)               :: here
+integer(di)               :: i, cstat, ierr, funit, id_PCM_1, n_id_PCM 
+character(20)             :: header1, header2, header3
+character(:), allocatable :: headline
 
 ! CODE
@@ -132,8 +133,8 @@
 
 ! Modify the header (first line)
-write(fch1,'(f'//int2str(nb_digits(n_yr_sim) + 5)//'.4)') n_yr_sim
-write(fch2,'(f'//int2str(nb_digits(ntot_yr_sim) + 5)//'.4)') ntot_yr_sim
-write(fch3,'(f6.4)') r_plnt2earth_yr ! 4 digits to the right of the decimal point to respect the precision of year in "pem_workflow_lib.sh"
-call execute_command_line('sed -i "1s/.*/'//trim(fch1)//' '//trim(fch2)//' '//trim(fch3)//' '//int2str(i_pcm_run)//' '//int2str(i_pem_run + 1)//' '//int2str(n_pcm_runs)//' '//int2str(n_pcm_runs_ini)//'/" pem_workflow.sts',cmdstat = cstat)
+write(header1,'(f'//int2str(nb_digits(n_yr_sim) + 5)//'.4)') n_yr_sim
+write(header2,'(f'//int2str(nb_digits(ntot_yr_sim) + 5)//'.4)') ntot_yr_sim
+write(header3,'(f6.4)') r_plnt2earth_yr ! 4 digits to the right of the decimal point to respect the precision of year in "pem_workflow_lib.sh"
+call execute_command_line('sed -i "1s/.*/'//trim(header1)//' '//trim(header2)//' '//trim(header3)//' '//int2str(i_pcm_run)//' '//int2str(i_pem_run + 1)//' '//int2str(n_pcm_runs)//' '//int2str(n_pcm_runs_ini)//'/" pem_workflow.sts',cmdstat = cstat)
 if (cstat > 0) then
     call stop_clean(__FILE__,__LINE__,'command execution failed!',1)
@@ -143,8 +144,20 @@
 
 ! Add the information of current PEM run at the end of file
+! Headline (cycle infomation), Date, Number of years done by the chained simulation, Number of years done by the PEM run, Code of the stopping criterion
+if (i_pcm_run - 1 <= n_pcm_runs_ini) then
+    id_PCM_1 = i_pcm_run - n_pcm_runs_ini
+    n_id_PCM = n_pcm_runs_ini
+else
+    id_PCM_1 = i_pcm_run - n_pcm_runs
+    n_id_PCM = n_pcm_runs
+end if
+headline = 'cycle=PCM('//int2str(id_PCM_1)
+do i = 1,n_id_PCM - 1
+    headline = headline//'+'//int2str(id_PCM_1 + i)
+end do
+headline = headline//')+PEM('//int2str(i_pem_run)//')'
 open(newunit = funit,file = statusfile_name,status = "old",position = "append",action = "write",iostat = ierr)
 if (ierr /= 0) call stop_clean(__FILE__,__LINE__,'error opening file "'//statusfile_name//'"!',ierr)
-! Date, Number of years done by the PEM run, Number of years done by the chained simulation, Code of the stopping criterion
-write(funit,'(f20.4,f20.4,f20.4,i20)') pem_ini_date + n_yr_sim, n_yr_run, n_yr_sim, stopPEM
+write(funit,'(a,f20.4,f20.4,f20.4,i4)') headline, pem_ini_date + n_yr_sim, n_yr_sim, n_yr_run, stopPEM
 close(funit)