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