Changeset 4072 for trunk/LMDZ.COMMON/libf/evolution/workflow_status.F90

Timestamp:

Feb 16, 2026, 11:41:37 AM (12 days ago)

Author:

jbclement

Message:

PEM:
Modify the PEM semantics to standardize file/functions/variable names and avoid ambiguity: workflow = repetion of cycles; cycle = PCM phase + PEM phase; PCM phase = 2+ PCM runs of 1 year each, PEM phase = 1 PEM run of 'x' years.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/workflow_status.F90 (moved) (moved from trunk/LMDZ.COMMON/libf/evolution/info.F90) (8 diffs)

trunk/LMDZ.COMMON/libf/evolution/workflow_status.F90

@@ -1 @@
-MODULE info
+MODULE workflow_status
 !-----------------------------------------------------------------------
 ! NAME
-!     info
+!     workflow_status
 !
 ! DESCRIPTION
 !     Handles counters for PCM/PEM coupled runs and duration information
-!     of the simulation metadata.
+!     of the simulation.
 !
 ! AUTHORS & DATE
@@ -25 +25 @@
 ! PARAMETERS
 ! ----------
-character(14), parameter :: infofile_name = 'launchPEM.info'
-
-! VARIABLES
-! ---------
-integer(di) :: iPCM, iPEM, nPCM, nPCM_ini ! Data about the chained simulation of PCM/PEM runs
+character(16), parameter :: statusfile_name = 'pem_workflow.sts'
+integer(di),   protected :: i_pcm_run, i_pem_run, n_pcm_runs, n_pcm_runs_ini ! Data about the chained simulation of PCM/PEM runs
 
 contains
@@ -35 +32 @@
 
 !=======================================================================
-SUBROUTINE read_info()
+SUBROUTINE read_workflow_status()
 !-----------------------------------------------------------------------
 ! NAME
-!    read_info
+!    read_workflow_status
 !
 ! DESCRIPTION
-!    Read the file "launchPEM.info" to get the number of simulated
+!    Read the file "pem_workflow.sts" to get the number of simulated
 !    years.
 !
@@ -54 +51 @@
 ! ------------
 use stoppage,  only: stop_clean
-use evolution, only: n_yr_sim, nmax_yr_sim, set_r_plnt2earth_yr
+use evolution, only: n_yr_sim, ntot_yr_sim, set_r_plnt2earth_yr
 use display,   only: print_msg
+use utility,   only: int2str
 
 ! DECLARATION
@@ -69 +67 @@
 ! CODE
 ! ----
-inquire(file = infofile_name,exist = here)
-if (.not. here) call stop_clean(__FILE__,__LINE__,'cannot find required file "'//infofile_name//'"! It should be created by the launching script.',1)
-call print_msg('> Reading "'//infofile_name//'"')
-open(newunit = funit,file = infofile_name,status = 'old',form = 'formatted',action = 'read',iostat = ierr)
-if (ierr /= 0) call stop_clean(__FILE__,__LINE__,'error opening file "'//infofile_name//'"!',ierr)
-read(funit,*) n_yr_sim, nmax_yr_sim, tmp, iPCM, iPEM, nPCM, nPCM_ini
+inquire(file = statusfile_name,exist = here)
+if (.not. here) call stop_clean(__FILE__,__LINE__,'cannot find required file "'//statusfile_name//'"! It should be created by the PEM workflow script.',1)
+call print_msg('> Reading "'//statusfile_name//'"')
+open(newunit = funit,file = statusfile_name,status = 'old',form = 'formatted',action = 'read',iostat = ierr)
+if (ierr /= 0) call stop_clean(__FILE__,__LINE__,'error opening file "'//statusfile_name//'"!',ierr)
+read(funit,*) n_yr_sim, ntot_yr_sim, tmp, i_pcm_run, i_pem_run, n_pcm_runs, n_pcm_runs_ini
 call set_r_plnt2earth_yr(tmp)
 close(funit)
 
-END SUBROUTINE read_info
+call print_msg('Current PEM run: '//int2str(i_pem_run))
+call print_msg('Input PCM runs : '//int2str(i_pcm_run - 1)//' and '//int2str(i_pcm_run))
+
+END SUBROUTINE read_workflow_status
 !=======================================================================
 
 !=======================================================================
-SUBROUTINE update_info(n_yr_run,stopPEM,n_yr_sim,nmax_yr_sim)
+SUBROUTINE update_workflow_status(n_yr_run,stopPEM,n_yr_sim,ntot_yr_sim)
 !-----------------------------------------------------------------------
 ! NAME
-!    update_info
+!    update_workflow_status
 !
 ! DESCRIPTION
-!    Update the first line of "launchPEM.info" to count the number of
-!    simulated years. Write in "launchPEM.info" the reason why the PEM
+!    Update the first line of "pem_workflow.sts" to count the number of
+!    simulated years. Write in "pem_workflow.sts" the reason why the PEM
 !    stopped and the number of simulated years.
 !
@@ -116 +117 @@
 real(dp),    intent(in) :: n_yr_run    ! # of years
 real(dp),    intent(in) :: n_yr_sim    ! Current simulated year
-real(dp),    intent(in) :: nmax_yr_sim ! Maximum number of years to be simulated
+real(dp),    intent(in) :: ntot_yr_sim ! Maximum number of years to be simulated
 
 ! LOCAL VARIABLES
@@ -126 +127 @@
 ! CODE
 ! ----
-call print_msg('> Updating "'//infofile_name//'"')
-inquire(file = infofile_name,exist = here)
-if (.not. here) call stop_clean(__FILE__,__LINE__,'cannot find required file "'//infofile_name//'"! It should be created by the launching script.',1)
+call print_msg('> Updating "'//statusfile_name//'"')
+inquire(file = statusfile_name,exist = here)
+if (.not. here) call stop_clean(__FILE__,__LINE__,'cannot find required file "'//statusfile_name//'"! It should be created by the PEM workflow script.',1)
 
 ! 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(nmax_yr_sim) + 5)//'.4)') nmax_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 "launch_pem.sh"
-call execute_command_line('sed -i "1s/.*/'//trim(fch1)//' '//trim(fch2)//' '//trim(fch3)//' '//int2str(iPCM)//' '//int2str(iPEM + 1)//' '//int2str(nPCM)//' '//int2str(nPCM_ini)//'/" launchPEM.info',cmdstat = cstat)
+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)
 if (cstat > 0) then
     call stop_clean(__FILE__,__LINE__,'command execution failed!',1)
@@ -142 +143 @@
 
 ! Add the information of current PEM run at the end of file
-open(newunit = funit,file = infofile_name,status = "old",position = "append",action = "write",iostat = ierr)
-if (ierr /= 0) call stop_clean(__FILE__,__LINE__,'error opening file "'//infofile_name//'"!',ierr)
-! Date, Number of years done by the PEM run, Number of years done by the chainded simulation, Code of the stopping criterion
+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
 close(funit)
 
-END SUBROUTINE update_info
+END SUBROUTINE update_workflow_status
 !=======================================================================
 
-END MODULE info
+END MODULE workflow_status