[253] | 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 = .false. |
---|
| 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 = .true. |
---|
| 11 | # Tidal resonance ratio ? ratio T_orbit to T_rotation |
---|
| 12 | nres = 10 |
---|
| 13 | # Write some more output on the screen ? |
---|
| 14 | lwrite = .false. |
---|
| 15 | # Save statistics in file "stats.nc" ? |
---|
| 16 | callstats = .false. |
---|
| 17 | # Test energy conservation of model physics ? |
---|
| 18 | enertest = .true. |
---|
[539] | 19 | # Check to see if cpp, mugaz values used match gas mixture defined in gases.def (recommended) ? |
---|
| 20 | check_cpp_match=.true. |
---|
[253] | 21 | |
---|
| 22 | ## Radiative transfer options |
---|
| 23 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 24 | # call radiative transfer? |
---|
| 25 | callrad = .true. |
---|
| 26 | # the rad. transfer is computed every "iradia" physical timestep |
---|
| 27 | iradia = 4 |
---|
| 28 | # call multilayer correlated-k radiative transfer ? |
---|
| 29 | corrk = .true. |
---|
[539] | 30 | # Include continuum absorption in radiative transfer (note CO2 is treated separately) ? |
---|
| 31 | continuum = .true. |
---|
[253] | 32 | # folder in which correlated-k data is stored ? |
---|
| 33 | corrkdir = CO2_H2Ovar |
---|
| 34 | # call visible gaseous absorption in radiative transfer ? |
---|
| 35 | callgasvis = .true. |
---|
| 36 | # Include Rayleigh scattering in the visible ? |
---|
| 37 | rayleigh = .true. |
---|
| 38 | # Characteristic planetary equilibrium (black body) temperature |
---|
| 39 | # This is used only in the aerosol radiative transfer setup. (see aerave.F) |
---|
| 40 | tplanet = 215. |
---|
| 41 | # Output spectral OLR in 1D/3D? |
---|
| 42 | specOLR = .false. |
---|
| 43 | # Output global radiative balance in file 'rad_bal.out' - slow for 1D!! |
---|
| 44 | meanOLR = .true. |
---|
| 45 | # Variable gas species: Radiatively active ? |
---|
| 46 | varactive = .true. |
---|
| 47 | # Variable gas species: Fixed vertical distribution ? |
---|
| 48 | varfixed = .false. |
---|
| 49 | # Variable gas species: Saturation percentage value at ground ? |
---|
| 50 | satval = 0.0 |
---|
| 51 | |
---|
| 52 | ## Star type |
---|
| 53 | ## ~~~~~~~~~ |
---|
| 54 | startype = 2 |
---|
| 55 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 56 | # The choices are: |
---|
| 57 | # |
---|
| 58 | # startype = 1 Sol (G2V-class main sequence) |
---|
| 59 | # startype = 2 Ad Leo (M-class, synthetic) |
---|
| 60 | # startype = 3 GJ644 |
---|
| 61 | # startype = 4 HD128167 |
---|
[1784] | 62 | # startype = 9 TRAPPIST-1 |
---|
| 63 | # startype = 10 Proxima Centauri |
---|
[253] | 64 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 65 | # Stellar flux at 1 AU. Examples: |
---|
| 66 | # 1366.0 W m-2 Sol today |
---|
| 67 | # 1024.5 W m-2 Sol today x 0.75 = weak early Sun |
---|
| 68 | # 18.462 W m-2 The feeble Gl581 |
---|
| 69 | # 19.960 W m-2 Gl581 with e=0.38 orbital average |
---|
| 70 | Fat1AU = 18.46 |
---|
| 71 | |
---|
| 72 | ## Tracer and aerosol options |
---|
| 73 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 74 | # Gravitational sedimentation of tracers (KEEP FALSE FOR NOW) ? |
---|
| 75 | sedimentation = .false. |
---|
| 76 | |
---|
| 77 | ## Other physics options |
---|
| 78 | ## ~~~~~~~~~~~~~~~~~~~~~ |
---|
| 79 | # call turbulent vertical diffusion ? |
---|
| 80 | calldifv = .true. |
---|
[727] | 81 | # use turbdiff instead of vdifc ? |
---|
| 82 | UseTurbDiff = .true. |
---|
[253] | 83 | # call convective adjustment ? |
---|
| 84 | calladj = .true. |
---|
| 85 | # call thermal conduction in the soil ? |
---|
| 86 | callsoil = .true. |
---|
| 87 | |
---|
| 88 | ######################################################################### |
---|
| 89 | ## extra non-standard definitions for Gliese 581d |
---|
| 90 | ######################################################################### |
---|
| 91 | |
---|
| 92 | ## Radiative transfer options |
---|
| 93 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 94 | # Parameterized Earth-like ozone absorption ? |
---|
| 95 | ozone = .false. |
---|
| 96 | |
---|
| 97 | ## Tracer and aerosol options |
---|
| 98 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
[727] | 99 | # Radiatively active CO2 aerosol? |
---|
| 100 | aeroco2 = .true. |
---|
| 101 | # Fixed CO2 aerosol distribution? |
---|
| 102 | aerofixco2 = .true. |
---|
| 103 | # Radiatively active water aerosol? |
---|
| 104 | aeroh2o = .false. |
---|
| 105 | # Fixed water aerosol distribution? |
---|
| 106 | aerofixh2o = .false. |
---|
| 107 | # basic dust opacity |
---|
| 108 | dusttau = 0.0 |
---|
[253] | 109 | # Varying H2O cloud fraction? |
---|
| 110 | CLFvarying = .false. |
---|
| 111 | # H2O cloud fraction? |
---|
| 112 | CLFfixval = 1.0 |
---|
| 113 | # number mixing ratio of CO2 ice particles |
---|
| 114 | Nmix_co2 = 100000. |
---|
| 115 | # number mixing ratio of water ice particles |
---|
| 116 | Nmix_h2o = 100000. |
---|
| 117 | |
---|
| 118 | ## Water options |
---|
| 119 | ## ~~~~~~~~~~~~~ |
---|
| 120 | # Model water cycle |
---|
| 121 | water = .true. |
---|
| 122 | # Model water cloud formation |
---|
| 123 | watercond = .true. |
---|
| 124 | # Model water precipitation (including coagulation etc.) |
---|
| 125 | waterrain = .true. |
---|
| 126 | # WATER: Precipitation threshold (simple scheme only) ? |
---|
[539] | 127 | rainthreshold = 0.001 |
---|
[253] | 128 | # Include hydrology ? |
---|
| 129 | hydrology = .true. |
---|
[1498] | 130 | # Spectral Dependant Albedo ? |
---|
| 131 | albedo_spectral_mode=.false. |
---|
| 132 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 133 | # If albedo_spectral_mode=.true., albedosnow becomes the 0.5 micron snow albedo. |
---|
| 134 | # |
---|
| 135 | # albedosnow = 0.95 (0.73 Sun-integrated) for fresh snow. |
---|
| 136 | # = 0.50 (0.39 Sun-integrated) for dirty snow. |
---|
| 137 | # = 0.645 (0.50 Sun-integrated) for 'realistic' snow. |
---|
| 138 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
[253] | 139 | # H2O snow (and ice) albedo ? |
---|
[1498] | 140 | albedosnow = 0.50 |
---|
[253] | 141 | # Maximum sea ice thickness ? |
---|
| 142 | maxicethick = 0.05 |
---|
| 143 | # Freezing point of seawater (degrees C) ? |
---|
| 144 | Tsaldiff = 0.0 |
---|
| 145 | |
---|
| 146 | ## CO2 options |
---|
| 147 | ## ~~~~~~~~~~~ |
---|
[1498] | 148 | # Co2 ice albedo ? |
---|
| 149 | albedoco2ice = 0.5 |
---|
[253] | 150 | # gas is non-ideal CO2 ? |
---|
| 151 | nonideal = .false. |
---|
| 152 | # call CO2 condensation ? |
---|
| 153 | co2cond = .true. |
---|
| 154 | # Set initial temperature profile to 1 K above CO2 condensation everywhere? |
---|
| 155 | nearco2cond = .false. |
---|
| 156 | |
---|