| [3175] | 1 | ## Datagcm path |
|---|
| 2 | datadir=/home/tanguy/SIMUS/datagcm |
|---|
| 3 | ## Orbit / general options |
|---|
| 4 | ## ~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 5 | # Run with or without tracer transport ? |
|---|
| 6 | tracer = .true. |
|---|
| 7 | # Diurnal cycle ? if diurnal=false, diurnally averaged solar heating |
|---|
| 8 | diurnal = .false. |
|---|
| 9 | # Seasonal cycle ? if season=false, Ls stays constant, to value set in "start" |
|---|
| 10 | season = .false. |
|---|
| 11 | # Run 1D or 3D version (false) / fast version volatile transport (true) |
|---|
| 12 | fast = .false. |
|---|
| 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 | # run with triton orbit ? |
|---|
| 20 | triton = .false. |
|---|
| 21 | # specific paleo run ? |
|---|
| 22 | paleo = .true. |
|---|
| 23 | # nb Earth years as paleoclimate step ? (if paleo = true) |
|---|
| 24 | paleoyears = 1. |
|---|
| 25 | # conservation of nitrogen |
|---|
| 26 | conservn2 = .false. |
|---|
| 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 = 1 |
|---|
| 34 | # call multilayer correlated-k radiative transfer ? |
|---|
| 35 | corrk = .true. |
|---|
| 36 | ch4fix = .false. |
|---|
| 37 | vmrch4fix = 0.5 |
|---|
| 38 | nlte=.false. |
|---|
| 39 | cooling=.false. |
|---|
| 40 | # folder in which correlated-k data is stored ? |
|---|
| 41 | corrkdir = pluton_ch4var_COfix |
|---|
| 42 | # call visible gaseous absorption in radiative transfer ? |
|---|
| 43 | callgasvis = .true. |
|---|
| 44 | # Output spectral OLR in 3D? |
|---|
| 45 | specOLR = .false. |
|---|
| 46 | # Output global radiative balance in file 'olr.txt' - slow for 1D!! |
|---|
| 47 | meanOLR = .false. |
|---|
| 48 | |
|---|
| 49 | # Star type (inherited from generic model) |
|---|
| 50 | startype = 1 |
|---|
| 51 | Fat1AU = 1366.0 |
|---|
| 52 | |
|---|
| 53 | ## Tracer and aerosol options |
|---|
| 54 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 55 | # Call N2 condensation ? |
|---|
| 56 | N2cond = .false. |
|---|
| 57 | # N2 condensation subtimestep ? |
|---|
| 58 | nbcond = 20 |
|---|
| 59 | # Gravitational sedimentation of tracers ? |
|---|
| 60 | sedimentation = .false. |
|---|
| 61 | # Compute methane cycle ? |
|---|
| 62 | methane = .true. |
|---|
| 63 | condmetsurf=.true. |
|---|
| 64 | # Compute CO cycle ? |
|---|
| 65 | carbox = .true. |
|---|
| 66 | condcosurf=.true. |
|---|
| 67 | # Compute methane clouds ? |
|---|
| 68 | metcloud = .false. |
|---|
| 69 | # Compute CO clouds ? |
|---|
| 70 | monoxcloud = .false. |
|---|
| 71 | # Haze production ? |
|---|
| 72 | haze = .true. |
|---|
| 73 | hazeconservch4 = .true. |
|---|
| 74 | # Fast haze production (with option fast=.true.) ? |
|---|
| 75 | fasthaze = .true. |
|---|
| 76 | # number mixing ratio of haze particles for ch4 clouds |
|---|
| 77 | Nmix_ch4 = 100000. |
|---|
| 78 | # number mixing ratio of haze particles for co clouds |
|---|
| 79 | Nmix_co = 100000. |
|---|
| 80 | # Radiatively active haze ? |
|---|
| 81 | aerohaze=.false. |
|---|
| 82 | proffix = .true. |
|---|
| 83 | # fractal particle ? |
|---|
| 84 | fractal=.false. |
|---|
| 85 | # geothermal flux at the bottom layer |
|---|
| 86 | fluxgeo = 0. |
|---|
| 87 | patchflux=0 ! Warm patch of flux in soil.F |
|---|
| 88 | |
|---|
| 89 | ## Surface properties |
|---|
| 90 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 91 | # options to change albedo with time |
|---|
| 92 | mode_n2=0 ! 0=fixed, 1=change w thickness, 2=special for SP |
|---|
| 93 | deltab=0. ! if mod_n2=1, difference of albedo depending of thickness |
|---|
| 94 | metlateq=35 ! latitude separating dark/bright methane |
|---|
| 95 | tholateq=25 ! latitude separating eq/pole tholins |
|---|
| 96 | mode_ch4=1 |
|---|
| 97 | specalb=.true. |
|---|
| 98 | |
|---|
| 99 | alb_n2b=0.7 |
|---|
| 100 | alb_n2a=0.7 |
|---|
| 101 | alb_ch4=0.5 |
|---|
| 102 | alb_ch4_eq=0.5 |
|---|
| 103 | alb_ch4_s=0.5 |
|---|
| 104 | alb_co=0.4 |
|---|
| 105 | alb_tho=0.2 |
|---|
| 106 | alb_tho_eq=0.1 |
|---|
| 107 | emis_n2b=0.8 |
|---|
| 108 | emis_n2a=0.8 |
|---|
| 109 | emis_ch4=0.8 |
|---|
| 110 | emis_co=0.8 |
|---|
| 111 | emis_tho=1.0 |
|---|
| 112 | emis_tho_eq=1.0 |
|---|
| 113 | |
|---|
| 114 | # Change TI with time |
|---|
| 115 | changeti=.true. |
|---|
| 116 | ITCH4=800. |
|---|
| 117 | ITN2=800 |
|---|
| 118 | ITH2O=2000. |
|---|
| 119 | |
|---|
| 120 | ## Other physics options |
|---|
| 121 | ## ~~~~~~~~~~~~~~~~~~~~ |
|---|
| 122 | # call turbulent vertical diffusion ? |
|---|
| 123 | calldifv = .true. |
|---|
| 124 | # call convective adjustment ? |
|---|
| 125 | calladj = .true. |
|---|
| 126 | # call thermal conduction in the soil ? |
|---|
| 127 | callsoil = .true. |
|---|
| 128 | # call thermal conduction |
|---|
| 129 | callconduct=.true. |
|---|
| 130 | # call molecular viscosity 3D only |
|---|
| 131 | callmolvis=.false. |
|---|
| 132 | # call molecular diffusion |
|---|
| 133 | callmoldiff=.false. |
|---|
| 134 | # flux qui arrive au-dessus de latmosphere |
|---|
| 135 | phitop=0. |
|---|
| 136 | |
|---|
