Changeset 2835 for trunk/LMDZ.COMMON/libf/evolution/temps_mod_evol.F90

Timestamp:

Nov 30, 2022, 11:29:29 AM (2 years ago)

Author:

romain.vande

Message:

Mars PEM:
Introduction of the possibility to follow an orbital forcing.
Introduction of new control parameters.
Cleaning of the PEM (removing unused files, add comments and new files)

A file named run_PEM.def can be added to the run.def. It contains the following variables:

_ evol_orbit_pem: Boolean. Do you want to follow an orbital forcing predefined (read in ob_ex_lsp.asc for example)? (default=false)
_ year_bp_ini: Integer. Number of year before present to start the pem run if evol_orbit_pem=.true. , default=0
_ Max_iter_pem: Integer. Maximal number of iteration if none of the stopping criterion is reached and if evol_orbit_pem=.false., default=99999999
_ dt_pem: Integer. Time step of the PEM in year, default=1
_ alpha_criterion: Real. Acceptance rate of sublimating ice surface change, default=0.2
_ soil_pem: Boolean. Do you want to run with subsurface physical processes in the PEM? default=.true.

RV

File:

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

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

@@ -3 +3 @@
 IMPLICIT NONE  
 
-  INTEGER   nyear           !     nyear   : Maximun number of year over which the PEM can interpolate
+  INTEGER   year_bp_ini     !     year_bp_ini : Initial year of the simulation of the PEM (in evol.def)
   INTEGER   dt_pem          !     dt_pem  : in years, the time step used by the PEM    
   REAL      alpha_criterion !     alpha_criterion : percentage of change before stopping the PEM 
-
-!$OMP THREADPRIVATE(nyear)        
-
-!WARNING: when adding a threadprivate variable in this module
-!        do not forget to add it to the copyin clause when opening an OpenMP
-!        parallel section. e.g. in gcm before call leapfrog_loc and/or
-!        possibly in iniphysiq
+  INTEGER   year_PEM        !     year written in startfiPEM.nc
+  INTEGER   Max_iter_pem    !     Maximal number of iteration when converging to a steady state, read in evol.def
+  LOGICAL   evol_orbit_pem  !     True if we want to follow the orbital parameters of ob_ex_lsp.asc, read in evol.def
 
 END MODULE temps_mod_evol