[3175] | 1 | datadir=/home1/tbertran/MODELES/KBO/datagcm/ |
---|
| 2 | ## Orbit / general options |
---|
| 3 | ## ~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 4 | # Run with or without tracer transport ? |
---|
| 5 | tracer = .true. |
---|
| 6 | # Diurnal cycle ? if diurnal=false, diurnally averaged solar heating |
---|
| 7 | diurnal = .true. |
---|
| 8 | # Seasonal cycle ? if season=false, Ls stays constant, to value set in "start" |
---|
| 9 | season = .true. |
---|
| 10 | # Run 1D or 3D version (false) / fast version volatile transport (true) |
---|
| 11 | fast = .false. |
---|
| 12 | # Write some more output on the screen ? |
---|
| 13 | lwrite = .false. |
---|
| 14 | # Save statistics in file "stats.nc" ? |
---|
| 15 | callstats = .false. |
---|
| 16 | # Test energy conservation of model physics ? |
---|
| 17 | enertest = .false. |
---|
| 18 | # run with triton orbit ? |
---|
| 19 | triton = .false. |
---|
| 20 | |
---|
| 21 | ## Radiative transfer options |
---|
| 22 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 23 | # call radiative transfer? |
---|
| 24 | callrad = .true. |
---|
| 25 | # the rad. transfer is computed every "iradia" physical timestep |
---|
| 26 | iradia = 5 |
---|
| 27 | # call multilayer correlated-k radiative transfer ? |
---|
| 28 | corrk = .true. |
---|
| 29 | ch4fix = .true. |
---|
| 30 | vmrch4fix = 0.5 |
---|
| 31 | nlte=.false. |
---|
| 32 | # folder in which correlated-k data is stored ? |
---|
| 33 | corrkdir = pluton_ch4var_COfix |
---|
| 34 | # call visible gaseous absorption in radiative transfer ? |
---|
| 35 | callgasvis = .true. |
---|
| 36 | # Output spectral OLR in 3D? |
---|
| 37 | specOLR = .false. |
---|
| 38 | # Output global radiative balance in file 'olr.txt' - slow for 1D!! |
---|
| 39 | meanOLR = .false. |
---|
| 40 | |
---|
| 41 | ## Star type (inherited from generic model) |
---|
| 42 | ## ~~~~~~~~~ |
---|
| 43 | startype = 1 |
---|
| 44 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 45 | # The choices are: |
---|
| 46 | # |
---|
| 47 | # startype = 1 Sol (G2V-class main sequence) |
---|
| 48 | # startype = 2 Ad Leo (M-class, synthetic) |
---|
| 49 | # startype = 3 GJ644 |
---|
| 50 | # startype = 4 HD128167 |
---|
| 51 | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 52 | # Stellar flux at 1 AU. Examples: |
---|
| 53 | # 1366.0 W m-2 Sol today |
---|
| 54 | # 1017.0 W m-2 Sol today x 0.75 = weak early Sun |
---|
| 55 | # 18.462 W m-2 The feeble Gl581 |
---|
| 56 | Fat1AU = 1366.0 |
---|
| 57 | |
---|
| 58 | ## Tracer and aerosol options |
---|
| 59 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 60 | # Call N2 condensation ? |
---|
| 61 | N2cond = .true. |
---|
| 62 | # Gravitational sedimentation of tracers ? |
---|
| 63 | sedimentation = .true. |
---|
| 64 | # Compute glacial flow (for fast option only)? |
---|
| 65 | glaflow = .false. |
---|
| 66 | # Compute methane cycle ? |
---|
| 67 | methane = .true. |
---|
| 68 | # Compute CO cycle ? |
---|
| 69 | carbox = .true. |
---|
| 70 | # Compute methane clouds ? |
---|
| 71 | metcloud = .true. |
---|
| 72 | # Compute CO clouds ? |
---|
| 73 | monoxcloud = .true. |
---|
| 74 | # Haze production ? |
---|
| 75 | haze = .true. |
---|
| 76 | # source additional haze ? |
---|
| 77 | source_haze= .false. |
---|
| 78 | # number mixing ratio of haze particles for ch4 clouds |
---|
| 79 | Nmix_ch4 = 100000. |
---|
| 80 | # number mixing ratio of haze particles for co clouds |
---|
| 81 | Nmix_co = 100000. |
---|
| 82 | # Radiatively active haze ? |
---|
| 83 | aerohaze=.false. |
---|
| 84 | # fractal particle ? |
---|
| 85 | fractal=.false. |
---|
| 86 | # Day fraction for limited cold trap in SP? |
---|
| 87 | dayfrac=0. |
---|
| 88 | # geothermal flux at the bottom layer |
---|
| 89 | fluxgeo = 0. |
---|
| 90 | |
---|
| 91 | ## Surface properties |
---|
| 92 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 93 | constalb=.true. |
---|
| 94 | alb_n2b=0.7 |
---|
| 95 | alb_n2a=0.7 |
---|
| 96 | alb_ch4_eq=0.5 |
---|
| 97 | alb_ch4=0.5 |
---|
| 98 | alb_co=0.4 |
---|
| 99 | alb_tho=0.2 |
---|
| 100 | alb_tho_eq=0.1 |
---|
| 101 | emis_n2b=0.8 |
---|
| 102 | emis_n2a=0.8 |
---|
| 103 | emis_ch4=0.8 |
---|
| 104 | emis_co=0.8 |
---|
| 105 | emis_tho=1.0 |
---|
| 106 | emis_tho_eq=1.0 |
---|
| 107 | |
---|
| 108 | # Change TI with time |
---|
| 109 | changeti=.false. |
---|
| 110 | ITCH4=800. |
---|
| 111 | ITN2=800 |
---|
| 112 | ITH2O=800. |
---|
| 113 | |
---|
| 114 | ## Other physics options |
---|
| 115 | ## ~~~~~~~~~~~~~~~~~~~~~ |
---|
| 116 | # call turbulent vertical diffusion ? |
---|
| 117 | calldifv = .true. |
---|
| 118 | # call convective adjustment ? |
---|
| 119 | calladj = .true. |
---|
| 120 | # call thermal conduction in the soil ? |
---|
| 121 | callsoil = .true. |
---|
| 122 | # call thermal conduction |
---|
| 123 | callconduct=.true. |
---|
| 124 | # call molecular viscosity 3D only |
---|
| 125 | callmolvis=.true. |
---|
| 126 | # call molecular diffusion |
---|
| 127 | callmoldiff=.true. |
---|
| 128 | # flux qui arrive au-dessus de latmosphere |
---|
| 129 | phitop=0. |
---|
| 130 | |
---|