| 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 = .false. |
|---|
| 17 | # Test energy conservation of model physics ? |
|---|
| 18 | enertest = .false. |
|---|
| 19 | # check if cpp and mugaz from start.nc are consistent with values computed by comp_cpp_mugaz with gases.def |
|---|
| 20 | check_cpp_match = .true. |
|---|
| 21 | |
|---|
| 22 | |
|---|
| 23 | ### run with triton orbit ? |
|---|
| 24 | ##triton = .false. |
|---|
| 25 | ### conservation n2 |
|---|
| 26 | ##conservn2=.false. |
|---|
| 27 | |
|---|
| 28 | ## Directory where external input files are |
|---|
| 29 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 30 | datadir = datadir |
|---|
| 31 | |
|---|
| 32 | ## Radiative transfer options |
|---|
| 33 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 34 | # call radiative transfer? |
|---|
| 35 | callrad = .true. |
|---|
| 36 | # the rad. transfer is computed every "iradia" physical timestep |
|---|
| 37 | iradia = 5 |
|---|
| 38 | # call multilayer correlated-k radiative transfer ? |
|---|
| 39 | corrk = .true. |
|---|
| 40 | # folder in which correlated-k data is stored ? |
|---|
| 41 | corrkdir = pluton_ch4var_COfix |
|---|
| 42 | # strictboundcorrk |
|---|
| 43 | strictboundcorrk = .true. |
|---|
| 44 | # call visible gaseous absorption in radiative transfer ? |
|---|
| 45 | callgasvis = .true. |
|---|
| 46 | # call continuum in radiative transfer ? |
|---|
| 47 | continuum = .false. |
|---|
| 48 | # Include Rayleigh scattering in the visible ? |
|---|
| 49 | rayleigh = .true. |
|---|
| 50 | # Test physics timescale in 1D ? |
|---|
| 51 | testradtimes = .false. |
|---|
| 52 | # Output spectral OLR in 1D/3D? |
|---|
| 53 | specOLR = .false. |
|---|
| 54 | # Output global radiative balance in file 'rad_bal.out' - slow for 1D!! |
|---|
| 55 | meanOLR = .false. |
|---|
| 56 | # Variable gas species: Radiatively active ? |
|---|
| 57 | varactive = .true. |
|---|
| 58 | # Variable gas species: Fixed vertical distribution ? |
|---|
| 59 | # (not to be used in time integration mode) |
|---|
| 60 | varfixed = .false. |
|---|
| 61 | # Variable gas species: Saturation percentage value at ground ? |
|---|
| 62 | satval = .8 |
|---|
| 63 | # Use fixed vertical profile, 1 step, no iteration ? |
|---|
| 64 | kastprof = .false. |
|---|
| 65 | |
|---|
| 66 | ## Star type |
|---|
| 67 | ## ~~~~~~~~~ |
|---|
| 68 | startype = 1 |
|---|
| 69 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 70 | # The choices are: |
|---|
| 71 | # |
|---|
| 72 | # startype = 1 Sol (G2V-class main sequence) |
|---|
| 73 | # startype = 2 Ad Leo (M-class, synthetic) |
|---|
| 74 | # startype = 3 GJ644 |
|---|
| 75 | # startype = 4 HD128167 |
|---|
| 76 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 77 | # Stellar flux at 1 AU. Examples: |
|---|
| 78 | # 1366.0 W m-2 Sol today |
|---|
| 79 | # 1024.5 W m-2 Sol today x 0.75 = weak early Sun |
|---|
| 80 | # 18.462 W m-2 The feeble Gl581 |
|---|
| 81 | # 19.960 W m-2 Gl581 with e=0.38 orbital average |
|---|
| 82 | Fat1AU = 1366.0 |
|---|
| 83 | |
|---|
| 84 | ## Tracer and aerosol options |
|---|
| 85 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 86 | # Gravitational sedimentation of tracers (just H2O ice for now) ? |
|---|
| 87 | sedimentation = .true. |
|---|
| 88 | |
|---|
| 89 | naerkind = 1 |
|---|
| 90 | hazeprop_file = optprop_tholins_fractal100nm.dat |
|---|
| 91 | |
|---|
| 92 | ## Other physics options |
|---|
| 93 | ## ~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 94 | # call turbulent vertical diffusion ? |
|---|
| 95 | calldifv = .true. |
|---|
| 96 | # use turbdiff instead of vdifc ? |
|---|
| 97 | UseTurbDiff = .true. |
|---|
| 98 | # call convective adjustment ? |
|---|
| 99 | calladj = .true. |
|---|
| 100 | # call thermal conduction in the soil ? |
|---|
| 101 | callsoil = .true. |
|---|
| 102 | lay1_soil = 1.4e-2 |
|---|
| 103 | |
|---|
| 104 | ######################################################################### |
|---|
| 105 | ## extra non-standard definitions for Earth |
|---|
| 106 | ######################################################################### |
|---|
| 107 | |
|---|
| 108 | ## Tracer and aerosol options |
|---|
| 109 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 110 | # atm mass update due to tracer evaporation/condensation? |
|---|
| 111 | mass_redistrib = .false. |
|---|
| 112 | ## Radiatively active CH4 aerosol? |
|---|
| 113 | #aeroch4 = .false. |
|---|
| 114 | |
|---|
| 115 | # haze production? |
|---|
| 116 | haze = .false. |
|---|
| 117 | # Radiatively active haze aerosol? |
|---|
| 118 | aerohcaer = .false. |
|---|
| 119 | # Fixed haze distribution? |
|---|
| 120 | aerofixhcaer = .false. |
|---|
| 121 | # haze particule mean radius |
|---|
| 122 | #rad_hcaer = 10e-9 |
|---|
| 123 | |
|---|
| 124 | |
|---|
| 125 | ## Condensation |
|---|
| 126 | ## ~~~~~~~~~~~ |
|---|
| 127 | # call N2 condensation ? |
|---|
| 128 | n2cond = .false. |
|---|
| 129 | # Set initial temperature profile to 1 K above N2 condensation everywhere? |
|---|
| 130 | #nearn2cond = .false. |
|---|
| 131 | |
|---|
