| 1 | ## Orbit / general options |
|---|
| 2 | ## ~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 3 | # Run with or without tracer transport ? |
|---|
| 4 | tracer = .false. |
|---|
| 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 = .false. |
|---|
| 17 | # Test energy conservation of model physics ? |
|---|
| 18 | enertest = .true. |
|---|
| 19 | # check if cpp matches gases.def (should be false in 1d) |
|---|
| 20 | check_cpp_match = .false. |
|---|
| 21 | # force cpp and mugaz values from gases.def (should be false in 1d) |
|---|
| 22 | force_cpp = .false. |
|---|
| 23 | |
|---|
| 24 | ## Directory where external input files are |
|---|
| 25 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 26 | datadir = ../../datagcm |
|---|
| 27 | |
|---|
| 28 | ## Radiative transfer options |
|---|
| 29 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 30 | # call radiative transfer? |
|---|
| 31 | callrad = .true. |
|---|
| 32 | # the rad. transfer is computed every "iradia" physical timestep |
|---|
| 33 | iradia = 5 |
|---|
| 34 | # Use blackbody for stellar spectrum ? |
|---|
| 35 | stelbbody = .false. |
|---|
| 36 | # Stellar blackbody temperature ? |
|---|
| 37 | stelTbb = 2000.000 |
|---|
| 38 | # call multilayer correlated-k radiative transfer ? |
|---|
| 39 | corrk = .true. |
|---|
| 40 | # folder in which correlated-k data is stored ? |
|---|
| 41 | corrkdir = 23x23 |
|---|
| 42 | # call visible gaseous absorption in radiative transfer ? |
|---|
| 43 | callgasvis = .true. |
|---|
| 44 | # call continuum in radiative transfer ? |
|---|
| 45 | Continuum = .true. |
|---|
| 46 | # Allow temperatures above kmatrix tmax |
|---|
| 47 | strictboundcorrk = .false. |
|---|
| 48 | # Include Rayleigh scattering in the visible ? |
|---|
| 49 | rayleigh = .true. |
|---|
| 50 | # Uniform absorption coefficient in radiative transfer? |
|---|
| 51 | graybody = .false. |
|---|
| 52 | # Constant absorption coefficient in visible |
|---|
| 53 | # (in m^2/kg; only if graybody=true): |
|---|
| 54 | # tau_surf= kappa*P/g |
|---|
| 55 | kappa_VI = 0. |
|---|
| 56 | # Constant absorption coefficient in IR |
|---|
| 57 | # (in m^2/kg; only if graybody=true): |
|---|
| 58 | kappa_IR = 0. |
|---|
| 59 | # Use Newtonian cooling in place of radiative transfer ? |
|---|
| 60 | newtonian = .false. |
|---|
| 61 | # Radiative timescale for Newtonian cooling ? [only if newtonian = T] |
|---|
| 62 | tau_relax = 30.00000 |
|---|
| 63 | # Test physics timescale in 1D ? |
|---|
| 64 | testradtimes = .false. |
|---|
| 65 | # Characteristic planetary equilibrium (black body) temperature |
|---|
| 66 | # This is used only in the aerosol radiative transfer setup. (see aerave.F) |
|---|
| 67 | tplanet = 94. |
|---|
| 68 | # Output spectral OLR in 1D/3D? |
|---|
| 69 | specOLR = .false. |
|---|
| 70 | # Output global radiative balance in file 'rad_bal.out' - slow for 1D!! |
|---|
| 71 | meanOLR = .false. |
|---|
| 72 | # Remove lower boundary (e.g. for gas giant sims) |
|---|
| 73 | noradsurf = .false. |
|---|
| 74 | |
|---|
| 75 | |
|---|
| 76 | # Compute effective variations of gravity field with altitude in physics ? |
|---|
| 77 | ## Better be true if you run with chemistry |
|---|
| 78 | eff_gz = .true. |
|---|
| 79 | |
|---|
| 80 | ## Chemistry |
|---|
| 81 | ## ~~~~~~~~~ |
|---|
| 82 | # call chemistry ? (must have tracer=true) |
|---|
| 83 | callchim = .false. |
|---|
| 84 | # Activate zonal mean ? |
|---|
| 85 | moyzon_ch = .false. |
|---|
| 86 | # chemistry is computed every "ichim" phys. timestep |
|---|
| 87 | ichim = 200 |
|---|
| 88 | |
|---|
| 89 | ## Microphysics |
|---|
| 90 | ## ~~~~~~~~~~~~ |
|---|
| 91 | # call microphysics ? (must have tracer=true) |
|---|
| 92 | callmufi = .false. |
|---|
| 93 | ## If yes, compute clouds ? |
|---|
| 94 | callclouds = .false. |
|---|
| 95 | ### If yes, number of ices ? (must be compatible with traceur.def AND microphysical model ) |
|---|
| 96 | nices = 3 |
|---|
| 97 | # Desactivate radiative coupling of microphysics ( i.e. use mean opacity profile ) |
|---|
| 98 | uncoupl_optic_haze = .true. |
|---|
| 99 | # ~~~~~ Parameters for microphysics ~~~~~ |
|---|
| 100 | # Path to microphys. config file ? |
|---|
| 101 | config_mufi = ../../datagcm/microphysics/config.ne15.cfg |
|---|
| 102 | # Pressure level of aer. production (Pa) ? |
|---|
| 103 | p_prod = 1.0 |
|---|
| 104 | # Aer. production rate (kg.m-2.s-1) ? |
|---|
| 105 | tx_prod = 3.5e-13 |
|---|
| 106 | # Equivalent radius production (m) ? |
|---|
| 107 | rc_prod = 2.0e-8 |
|---|
| 108 | # Radius of air (nitrogen) molecule (m) ? |
|---|
| 109 | air_rad = 1.75e-10 |
|---|
| 110 | |
|---|
| 111 | |
|---|
| 112 | ## Star type |
|---|
| 113 | ## ~~~~~~~~~ |
|---|
| 114 | startype = 1 |
|---|
| 115 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 116 | # The choices are: |
|---|
| 117 | # |
|---|
| 118 | # startype = 1 Sol (G2V-class main sequence) |
|---|
| 119 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 120 | # Stellar flux at 1 AU. Examples: |
|---|
| 121 | # 1366.0 W m-2 Sol today |
|---|
| 122 | # 1024.5 W m-2 Sol today x 0.75 = weak early Sun |
|---|
| 123 | Fat1AU = 1366.0 |
|---|
| 124 | |
|---|
| 125 | ## Other physics options |
|---|
| 126 | ## ~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 127 | # call turbulent vertical diffusion ? |
|---|
| 128 | calldifv = .true. |
|---|
| 129 | # use turbdiff instead of vdifc ? |
|---|
| 130 | UseTurbDiff = .true. |
|---|
| 131 | # call convective adjustment ? |
|---|
| 132 | calladj = .true. |
|---|
| 133 | # call thermal conduction in the soil ? |
|---|
| 134 | callsoil = .true. |
|---|