Changeset 3498 for trunk/LMDZ.COMMON/libf/evolution/time_evol_mod.F90

Timestamp:

Nov 7, 2024, 2:48:08 PM (3 weeks ago)

Author:

jbclement

Message:

PEM:

• Correction of the variable name for the ice table depth in "pemetat0.F90". So it is now got as intended from the "startpem.nc" file;
• Renaming of the tendencies in the PEM with the prefix 'd_' instead of 'tend_';
• Modification of the PEM time step type from integer to real. As a consequence, all time variables are now of real type. This change adds the possibility to consider fractions of year as time step.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/time_evol_mod.F90 (modified) (1 diff)

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

@@ -3 +3 @@
 implicit none
 
-integer :: year_bp_ini    ! Initial year (in Planetary years) of the simulation of the PEM defined in run.def (in Earth years)
-integer :: dt_pem         ! Time step used by the PEM in Planetary years
+real    :: year_bp_ini    ! Initial year (in Planetary years) of the simulation of the PEM defined in run.def (in Earth years)
+real    :: dt             ! Time step used by the PEM in Planetary years
 real    :: convert_years  ! Conversion ratio from Planetary years to Earth years
 real    :: h2o_ice_crit   ! Percentage of change of the surface of h2o ice sublimating before stopping the PEM
 real    :: co2_ice_crit   ! Percentage of change of the surface of co2 ice sublimating before stopping the PEM
 real    :: ps_criterion   ! Percentage of change of averaged surface pressure before stopping the PEM
-integer :: Max_iter_pem   ! Maximal number of iteration when converging to a steady state, read in evol.def
+integer :: Max_iter_pem   ! Maximal number of iterations when converging to a steady state, read in evol.def
 logical :: evol_orbit_pem ! True if we want to follow the orbital parameters of obl_ecc_lsp.asc, read in evol.def
 logical :: var_obl        ! True if we want the PEM to follow obl_ecc_lsp.asc parameters for obliquity