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

Timestamp:

Feb 16, 2026, 11:41:37 AM (13 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/pem.F90 (modified) (7 diffs)

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

@@ -28 +28 @@
 use config,             only: read_rundef
 use display,            only: print_ini, print_end, print_msg
-use evolution,          only: n_yr_run, n_yr_sim, nmax_yr_sim, nmax_yr_run, dt, idt, r_plnt2earth_yr, pem_ini_date
+use evolution,          only: n_yr_run, n_yr_sim, ntot_yr_sim, nmax_yr_run, dt, idt, r_plnt2earth_yr, pem_ini_date
 use geometry,           only: ngrid, nslope, nday, nsoil_PCM, nsoil, cell_area, total_surface, nlayer
 use glaciers,           only: h2oice_flow, co2ice_flow, flow_co2glaciers, flow_h2oglaciers
 use ice_table,          only: icetable_equilibrium, icetable_dynamic, icetable_depth, icetable_thickness, ice_porefilling, evolve_ice_table
-use info,               only: read_info, update_info
 use layered_deposits,   only: layering, do_layering, del_layering, evolve_layering, ptrarray, layering2surfice, surfice2layering
 use maths,              only: pi
@@ -51 +50 @@
 use tracers,            only: adapt_tracers2pressure, build4PCM_tracers, nq
 use utility,            only: real2str
+use workflow_status,    only: read_workflow_status, update_workflow_status
 use xios_data,          only: load_xios_data
 
@@ -163 +163 @@
 call ini_allocation()
 
-! Read the duration information of the simulation
-call read_info()
+! Read the duration information of the workflow
+call read_workflow_status()
 
 ! Read the PEM parameters
@@ -283 +283 @@
 n_yr_run = 0
 idt = 0
-do while (n_yr_run < min(nmax_yr_runorb,nmax_yr_run) .and. n_yr_sim < nmax_yr_sim)
-    call print_msg('**** Iteration of the PEM leg (Planetary years): '//real2str(n_yr_run + dt))
+do while (n_yr_run < min(nmax_yr_runorb,nmax_yr_run) .and. n_yr_sim < ntot_yr_sim)
+    call print_msg('**** Iteration of the PEM run [Planetary years]: '//real2str(n_yr_run + dt))
     ! Evolve global surface pressure
     call evolve_pressure(d_co2ice,delta_co2_ads,do_sorption,ps_avg_glob_old,ps_avg_glob,ps_avg)
@@ -445 +445 @@
     if (stopcrit%stop_code() == 0 .and. n_yr_run >= nmax_yr_run) stopcrit%nmax_yr_run_reached = .true.
     if (stopcrit%stop_code() == 0 .and. n_yr_run >= nmax_yr_runorb) stopcrit%nmax_yr_runorb_reached = .true.
-    if (stopcrit%stop_code() == 0 .and. n_yr_sim >= nmax_yr_sim) stopcrit%nmax_yr_sim_reached = .true.
+    if (stopcrit%stop_code() == 0 .and. n_yr_sim >= ntot_yr_sim) stopcrit%nmax_yr_sim_reached = .true.
     if (stopcrit%stop_code() == 0 .and. timewall .and. real(c2 - c1,dp)/real(cr,dp) >= timelimit - antetime) stopcrit%time_limit_reached = .true.
     if (stopcrit%is_any_set()) then
@@ -452 +452 @@
     else
         call print_msg('**** This run has achieved '//real2str(n_yr_run)//' Planetary years.')
-        call print_msg('**** The chained simulation has achieved '//real2str(n_yr_sim)//' Planetary years.')
+        call print_msg('**** The workflow has achieved '//real2str(n_yr_sim)//' Planetary years.')
         call print_msg('**** This run is continuing!')
         call print_msg('****')
@@ -501 +501 @@
 call write_restart(ps4PCM,pa4PCM,preff4PCM,q4PCM,teta4PCM,air_mass4PCM)
 
-! Update the duration information of the simulation
-call update_info(n_yr_run,stopcrit%stop_code(),n_yr_sim,nmax_yr_sim)
+! Update the duration information of the workflow
+call update_workflow_status(n_yr_run,stopcrit%stop_code(),n_yr_sim,ntot_yr_sim)
 
 ! Deallocation