| 1 | ## Orbit / general options |
|---|
| 2 | ## ~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 3 | # Run with or without tracer transport ? |
|---|
| 4 | tracer = .true. |
|---|
| 5 | # Diurnal cycle ? if diurnal=false, diurnally averaged solar heating |
|---|
| 6 | diurnal = .true. |
|---|
| 7 | # Seasonal cycle ? if season=false, Ls stays constant, to value set in "start" |
|---|
| 8 | season = .true. |
|---|
| 9 | # Tidally resonant orbit ? must have diurnal=false, correct rotation rate in newstart |
|---|
| 10 | tlocked = .false. |
|---|
| 11 | # Tidal resonance ratio ? ratio T_orbit to T_rotation |
|---|
| 12 | nres = 1 |
|---|
| 13 | # Write some more output on the screen ? |
|---|
| 14 | lwrite = .false. |
|---|
| 15 | # Save statistics in file "stats.nc" ? |
|---|
| 16 | callstats = .true. |
|---|
| 17 | # Test energy conservation of model physics ? |
|---|
| 18 | enertest = .false. |
|---|
| 19 | |
|---|
| 20 | ## Radiative transfer options |
|---|
| 21 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 22 | # call radiative transfer? |
|---|
| 23 | callrad = .true. |
|---|
| 24 | # the rad. transfer is computed every "iradia" physical timestep |
|---|
| 25 | iradia = 4 |
|---|
| 26 | # call multilayer correlated-k radiative transfer ? |
|---|
| 27 | corrk = .true. |
|---|
| 28 | # folder in which correlated-k data is stored ? |
|---|
| 29 | corrkdir = earlymars |
|---|
| 30 | # call visible gaseous absorption in radiative transfer ? |
|---|
| 31 | callgasvis = .true. |
|---|
| 32 | # Include Rayleigh scattering in the visible ? |
|---|
| 33 | rayleigh = .true. |
|---|
| 34 | # Characteristic planetary equilibrium (black body) temperature |
|---|
| 35 | # This is used only in the aerosol radiative transfer setup. (see aerave.F) |
|---|
| 36 | tplanet = 215. |
|---|
| 37 | # Output spectral OLR in 1D/3D? |
|---|
| 38 | specOLR = .false. |
|---|
| 39 | # Output global radiative balance in file 'rad_bal.out' - slow for 1D!! |
|---|
| 40 | meanOLR = .true. |
|---|
| 41 | # Variable gas species: Radiatively active ? |
|---|
| 42 | varactive = .true. |
|---|
| 43 | # Variable gas species: Fixed vertical distribution ? |
|---|
| 44 | varfixed = .false. |
|---|
| 45 | # Variable gas species: Saturation percentage value at ground ? |
|---|
| 46 | satval = 0.0 |
|---|
| 47 | |
|---|
| 48 | ## Star type |
|---|
| 49 | ## ~~~~~~~~~ |
|---|
| 50 | startype = 1 |
|---|
| 51 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 52 | # The choices are: |
|---|
| 53 | # |
|---|
| 54 | # startype = 1 Sol (G2V-class main sequence) |
|---|
| 55 | # startype = 2 Ad Leo (M-class, synthetic) |
|---|
| 56 | # startype = 3 GJ644 |
|---|
| 57 | # startype = 4 HD128167 |
|---|
| 58 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 59 | # Stellar flux at 1 AU. Examples: |
|---|
| 60 | # 1366.0 W m-2 Sol today |
|---|
| 61 | # 1024.5 W m-2 Sol today x 0.75 = weak Sun 3.8 GYa |
|---|
| 62 | # 18.462 W m-2 The feeble GJ581 |
|---|
| 63 | # 19.960 W m-2 GJ581 with e=0.38 orbital average |
|---|
| 64 | Fat1AU = 1024.5 |
|---|
| 65 | |
|---|
| 66 | ## Tracer and aerosol options |
|---|
| 67 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 68 | # Gravitational sedimentation of tracers (KEEP FALSE FOR NOW) ? |
|---|
| 69 | sedimentation = .false. |
|---|
| 70 | |
|---|
| 71 | ## Other physics options |
|---|
| 72 | ## ~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 73 | # call turbulent vertical diffusion ? |
|---|
| 74 | calldifv = .true. |
|---|
| 75 | # call convective adjustment ? |
|---|
| 76 | calladj = .true. |
|---|
| 77 | # call thermal conduction in the soil ? |
|---|
| 78 | callsoil = .true. |
|---|
| 79 | |
|---|
| 80 | ######################################################################### |
|---|
| 81 | ## extra non-standard definitions for Gliese 581d |
|---|
| 82 | ######################################################################### |
|---|
| 83 | |
|---|
| 84 | ## Tracer and aerosol options |
|---|
| 85 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 86 | # Fixed aerosol distributions? |
|---|
| 87 | aerofixed = .false. |
|---|
| 88 | # Varying H2O cloud fraction? |
|---|
| 89 | CLFvarying = .false. |
|---|
| 90 | # H2O cloud fraction? |
|---|
| 91 | CLFfixval = 0.5 |
|---|
| 92 | # number mixing ratio of CO2 ice particles |
|---|
| 93 | Nmix_co2 = 100000. |
|---|
| 94 | # number mixing ratio of water ice particles |
|---|
| 95 | Nmix_h2o = 100000. |
|---|
| 96 | |
|---|
| 97 | ## Water options |
|---|
| 98 | ## ~~~~~~~~~~~~~ |
|---|
| 99 | # Model water cycle |
|---|
| 100 | water = .true. |
|---|
| 101 | # Model water cloud formation |
|---|
| 102 | watercond = .true. |
|---|
| 103 | # Model water precipitation (including coagulation etc.) |
|---|
| 104 | waterrain = .true. |
|---|
| 105 | # WATER: Precipitation threshold (simple scheme only) ? |
|---|
| 106 | rainthreshold = 0.011 |
|---|
| 107 | # WATER: hydrology ? |
|---|
| 108 | hydrology = .true. |
|---|
| 109 | |
|---|
| 110 | ## CO2 options |
|---|
| 111 | ## ~~~~~~~~~~~ |
|---|
| 112 | # gas is non-ideal CO2 ? |
|---|
| 113 | nonideal = .false. |
|---|
| 114 | # call CO2 condensation ? |
|---|
| 115 | co2cond = .true. |
|---|
| 116 | # Set initial temperature profile to 1 K above CO2 condensation everywhere? |
|---|
| 117 | nearco2cond = .false. |
|---|
| 118 | |
|---|
| 119 | ## N2 options |
|---|
| 120 | ## ~~~~~~~~~~ |
|---|
| 121 | # include N2 in cp and Rayleigh scattering calculation ? |
|---|
| 122 | addn2 = .false. |
|---|
| 123 | # N2 mixing ratio |
|---|
| 124 | n2mixratio = 0.0 |
|---|
| 125 | |
|---|
| 126 | ## H2 options |
|---|
| 127 | ## ~~~~~~~~~~ |
|---|
| 128 | # include H2 in cp and Rayleigh scattering calculation ? |
|---|
| 129 | # for now we assume H2 dominant if included |
|---|
| 130 | addh2 = .false. |
|---|
