[3193] | 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" |
---|
[3195] | 8 | season = .true. |
---|
[3193] | 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 | |
---|
[3195] | 84 | ## Tracer and aerosol options |
---|
[3193] | 85 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 86 | # Gravitational sedimentation of tracers (just H2O ice for now) ? |
---|
| 87 | sedimentation = .true. |
---|
| 88 | |
---|
[3197] | 89 | naerkind = 1 |
---|
[3329] | 90 | hazeprop_file = optprop_tholins_fractal100nm.dat |
---|
[3193] | 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 | ######################################################################### |
---|
[3195] | 107 | |
---|
| 108 | ## Tracer and aerosol options |
---|
[3193] | 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 | |
---|
[3195] | 125 | ## Condensation |
---|
[3193] | 126 | ## ~~~~~~~~~~~ |
---|
| 127 | # call N2 condensation ? |
---|
| 128 | n2cond = .false. |
---|
| 129 | # Set initial temperature profile to 1 K above N2 condensation everywhere? |
---|
| 130 | #nearn2cond = .false. |
---|
| 131 | |
---|