[3308] | 1 | #------------------------------------# |
---|
| 2 | # Run control parameters for the PEM # |
---|
| 3 | #------------------------------------# |
---|
[2963] | 4 | |
---|
[3308] | 5 | #---------- Planet type ----------# |
---|
| 6 | planet_type=mars |
---|
| 7 | |
---|
[3082] | 8 | #---------- Output parameters ----------# |
---|
| 9 | # Frequency of outputs for the PEM? Default = 1 |
---|
| 10 | ecritpem=1 |
---|
| 11 | |
---|
[3038] | 12 | #---------- Orbital parameters ----------# |
---|
| 13 | # Do you want to follow an orbital forcing read in "obl_ecc_lsp.asc"? Default = .false. |
---|
[3327] | 14 | # evol_orbit_pem=.false. |
---|
[2963] | 15 | |
---|
[3038] | 16 | # If evol_orbit_pem=.true., number of Earth years before present to start the PEM run? Default = 0 |
---|
[3327] | 17 | # year_earth_bp_ini=0 |
---|
[2963] | 18 | |
---|
[3038] | 19 | # Do you want to vary the obliquity when following "obl_ecc_lsp.asc"? Default = .true. |
---|
[3327] | 20 | # var_obl=.false. |
---|
[2963] | 21 | |
---|
[3038] | 22 | # Do you want to vary the eccentricity when following "obl_ecc_lsp.asc"? Default = .true. |
---|
[3327] | 23 | # var_ecc=.false. |
---|
[2963] | 24 | |
---|
[3038] | 25 | # Do you want to vary the ls perihelie when following "obl_ecc_lsp.asc"? Default = .true. |
---|
[3327] | 26 | # var_lsp=.false. |
---|
[2963] | 27 | |
---|
[3038] | 28 | #---------- Stopping criteria parameters ----------# |
---|
| 29 | # If evol_orbit_pem=.false., maximal number of iterations if no stopping criterion is reached? Default=100000000 |
---|
[3327] | 30 | # Max_iter_pem=10 |
---|
[2963] | 31 | |
---|
[3161] | 32 | # Acceptance rate of sublimating H2O ice surface change? Default = 0.2 |
---|
[3159] | 33 | # h2o_ice_crit=0.2 |
---|
[2963] | 34 | |
---|
[3038] | 35 | # Acceptance rate of sublimating CO2 ice surface change? Default = 0.2 |
---|
[3159] | 36 | # co2_ice_crit=0.2 |
---|
[2963] | 37 | |
---|
[3038] | 38 | # Acceptance rate of pressure surface change? Default = 0.15 |
---|
[2963] | 39 | # ps_criterion = 0.15 |
---|
| 40 | |
---|
[3498] | 41 | # Time step length of the PEM in Martian years? Default = 1. |
---|
| 42 | # dt=1 |
---|
[2963] | 43 | |
---|
[3038] | 44 | #---------- Subsurface parameters ----------# |
---|
| 45 | # Do you want to run with subsurface physical processes in the PEM? Default = .true. |
---|
[2963] | 46 | # soil_pem=.true. |
---|
| 47 | |
---|
[3327] | 48 | # Do you want to run with adsoprtion in the PEM? Default = .false. |
---|
[3038] | 49 | # adsorption_pem=.true. |
---|
[2963] | 50 | |
---|
[3038] | 51 | # Do you want to modify the soil thermal properties with the pressure? Default = .false. |
---|
| 52 | # reg_thprop_dependp=.false. |
---|
[2963] | 53 | |
---|
[3161] | 54 | # Value of the geothermal flux? Default = 0. |
---|
[3038] | 55 | # fluxgeo=0. |
---|
[2963] | 56 | |
---|
[3038] | 57 | # Depth at which the breccia layer begins? Default = 10 m |
---|
| 58 | # depth_breccia=10. |
---|
[2963] | 59 | |
---|
[3038] | 60 | # Depth at which the bedrock layer begins? Default = 1000 m |
---|
| 61 | # depth_bedrock=1000. |
---|
[2963] | 62 | |
---|
[3038] | 63 | # Is the ice table computed at equilibrium? Default = .true. |
---|
[3339] | 64 | # icetable_equilibrium=.false. |
---|
[2963] | 65 | |
---|
[3038] | 66 | # Is the ice table computed with the dynamic method? Default = .false. |
---|
| 67 | # icetable_dynamic=.false. |
---|
[2963] | 68 | |
---|
[3159] | 69 | #---------- Ice management parameters ----------# |
---|
[3256] | 70 | # Amount of H2O ice to initialize the huge reservoir if the variable is not present in "startfi_PEM.nc"? Default = 9200. kg.m-2 (= 10 m) |
---|
[3161] | 71 | # ini_huge_h2oice=1.e4 |
---|
[3136] | 72 | |
---|
[3256] | 73 | # Threshold to consider the amount of H2O ice as an infinite reservoir? Default = 460. kg.m-2 (= 0.5 m) |
---|
[3159] | 74 | # inf_h2oice_threshold=2.e3 |
---|
| 75 | |
---|
[3161] | 76 | # Do you want H2O frost to transform into perennial H2O ice? Default = .false. |
---|
| 77 | # metam_h2oice=.false. |
---|
| 78 | |
---|
[3256] | 79 | # Threshold to consider frost is becoming perennial H2O ice? Default = 460. kg.m-2 (= 0.5 m) |
---|
[3159] | 80 | # metam_h2oice_threshold=5.e-2 |
---|
| 81 | |
---|
[3161] | 82 | # Do you want the H2O ice to flow along subslope inside a cell? Default = .true. |
---|
| 83 | # h2oice_flow=.true. |
---|
| 84 | |
---|
| 85 | # Do you want CO2 frost to transform into perennial CO2 ice? Default = .false. |
---|
| 86 | # metam_co2ice=.false. |
---|
| 87 | |
---|
[3256] | 88 | # Threshold to consider frost is becoming perennial CO2 ice? Default = 16500. kg.m-2 (= 10 m) |
---|
[3159] | 89 | # metam_co2ice_threshold=16.e3 |
---|
| 90 | |
---|
[3161] | 91 | # Do you want the CO2 ice to flow along subslope inside a cell? Default = .true. |
---|
[3327] | 92 | # co2ice_flow=.false. |
---|
[3161] | 93 | |
---|
[3320] | 94 | #---------- Layering parameters ----------# |
---|
[3319] | 95 | # Do you want to run with the layering algorithm? Default = .false. |
---|
| 96 | # layering=.true. |
---|
| 97 | |
---|
[3136] | 98 | # Some definitions for the physics, in file 'callphys.def' |
---|
| 99 | INCLUDEDEF=callphys.def |
---|