source: trunk/LMDZ.COMMON/libf/evolution/deftank/run_PEM.def @ 3154

Last change on this file since 3154 was 3149, checked in by jbclement, 2 years 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: 2.5 KB
Line 
1#
2#-------------------------------------
3# Run control parameters for the PEM
4#-------------------------------------
5
6#---------- Output parameters ----------#
7# Frequency of outputs for the PEM? Default = 1
8ecritpem=1
9
10#---------- Orbital parameters ----------#
11# Do you want to follow an orbital forcing read in "obl_ecc_lsp.asc"? Default = .false.
12  evol_orbit_pem=.true.
13
14# If evol_orbit_pem=.true., number of Earth years before present to start the PEM run? Default = 0
15  year_earth_bp_ini=-1000000
16
17# Do you want to vary the obliquity when following "obl_ecc_lsp.asc"? Default = .true.
18  var_obl=.true.
19
20# Do you want to vary the eccentricity when following "obl_ecc_lsp.asc"? Default = .true.
21  var_ecc=.true.
22
23# Do you want to vary the ls perihelie when following "obl_ecc_lsp.asc"? Default = .true.
24  var_lsp=.true.
25
26#---------- Stopping criteria parameters ----------#
27# If evol_orbit_pem=.false., maximal number of iterations if no stopping criterion is reached? Default=100000000
28#  Max_iter_pem=100000000
29
30# Acceptance rate of sublimating water ice surface change? Default = 0.2
31# water_ice_criterion=0.2
32
33# Acceptance rate of sublimating CO2 ice surface change? Default = 0.2
34# co2_ice_criterion=0.2
35
36# Acceptance rate of pressure surface change? Default = 0.15
37# ps_criterion = 0.15
38
39# Time step length of the PEM in Martian years? Default = 1
40# dt_pem=1
41
42#---------- Subsurface parameters ----------#
43# Do you want to run with subsurface physical processes in the PEM? Default = .true.
44# soil_pem=.true.
45
46# Do you want to run with adsoprtion in the PEM? Default = .true.
47# adsorption_pem=.true.
48
49# Do you want the glaciers to flow along subslope inside a mesh? Default = .true.
50# co2glaciersflow=.true.
51
52# Do you want to modify the soil thermal properties with the pressure? Default = .false.
53# reg_thprop_dependp=.false.
54
55#---------- Layering parameters ----------#
56# Value of the geothermal flux? Default= 0.
57# fluxgeo=0.
58
59# Depth at which the breccia layer begins? Default = 10 m
60# depth_breccia=10.
61
62# Depth at which the bedrock layer begins? Default = 1000 m
63# depth_bedrock=1000.
64
65# Is the ice table computed at equilibrium? Default = .true.
66# icetable_equilibrium=.true.
67
68# Is the ice table computed with the dynamic method? Default = .false.
69# icetable_dynamic=.false.
70
71#---------- H2O reservoir parameters ----------#
72# Amount of h2o ice to initialize the big reservoir if the variable is not present in startfi_PEM.nc? Default = 1.e4
73# ini_h2o_bigreservoir=1e4
74
75# Some definitions for the physics, in file 'callphys.def'
76INCLUDEDEF=callphys.def
Note: See TracBrowser for help on using the repository browser.