| 1 | # |
|---|
| 2 | #------------------------------------- |
|---|
| 3 | # Run control parameters for the PEM |
|---|
| 4 | #------------------------------------- |
|---|
| 5 | |
|---|
| 6 | #---------- Output parameters ----------# |
|---|
| 7 | # Frequency of outputs for the PEM? Default = 1 |
|---|
| 8 | ecritpem=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' |
|---|
| 76 | INCLUDEDEF=callphys.def |
|---|