|
Last change
on this file since 3152 was
3149,
checked in by jbclement, 19 months ago
|
|
PEM:
- Simplification of the algorithm managing the stopping criteria;
- Complete rework of the ice management in the PEM (H2O & CO2);
Subroutines to evolve the H2O and CO2 ice are now in the same module "evol_ice_mod.F90".
Tendencies are computed from the variation of "ice + frost" between the 2 PCM runs.
Evolving ice in the PEM is now called 'h2o_ice' or 'co2_ice' (not anymore in 'qsurf' and free of 'water_reservoir').
Default value 'ini_h2o_bigreservoir' (= 10 m) initializes the H2O ice of the first PEM run where there is 'watercap'. For the next PEM runs, initialization is done with the value kept in "startpem.nc". CO2 ice is taken from 'perennial_co2ice' of the PCM (paleoclimate flag must be true).
Simplification of the condition to compute the surface ice cover needed for the stopping criteria.
Frost ('qsurf') is not evolved by the PEM and given back to the PCM.
New default threshold value 'inf_h2oice_threshold' (= 2 m) to decide at the end of the PEM run if the H2O ice should be 'watercap' or not for the next PCM runs. If H2O ice cannot be 'watercap', then the remaining H2O ice is transferred to the frost ('qsurf').
- Renaming of variables/subroutines for clarity;
- Some cleanings throughout the code;
- Small updates in files of the deftank.
JBC
|
|
File size:
1.2 KB
|
| Line | |
|---|
| 1 | MODULE time_evol_mod |
|---|
| 2 | |
|---|
| 3 | implicit none |
|---|
| 4 | |
|---|
| 5 | integer :: year_bp_ini ! Initial year (in Planetary years) of the simulation of the PEM defined in run.def (in Earth years) |
|---|
| 6 | integer :: dt_pem ! Time step used by the PEM in Planetary years |
|---|
| 7 | real :: convert_years ! Conversion ratio from Planetary years to Earth years |
|---|
| 8 | real :: water_ice_criterion ! Percentage of change of the surface of water ice sublimating before stopping the PEM |
|---|
| 9 | real :: co2_ice_criterion ! Percentage of change of the surface of co2 ice sublimating before stopping the PEM |
|---|
| 10 | real :: ps_criterion ! Percentage of change of averaged surface pressure before stopping the PEM |
|---|
| 11 | integer :: Max_iter_pem ! Maximal number of iteration when converging to a steady state, read in evol.def |
|---|
| 12 | logical :: evol_orbit_pem ! True if we want to follow the orbital parameters of obl_ecc_lsp.asc, read in evol.def |
|---|
| 13 | logical :: var_obl ! True if we want the PEM to follow obl_ecc_lsp.asc parameters for obliquity |
|---|
| 14 | logical :: var_ecc ! True if we want the PEM to follow obl_ecc_lsp.asc parameters for eccenticity |
|---|
| 15 | logical :: var_lsp ! True if we want the PEM to follow obl_ecc_lsp.asc parameters for ls perihelie |
|---|
| 16 | |
|---|
| 17 | END MODULE time_evol_mod |
|---|
Note: See
TracBrowser
for help on using the repository browser.
