[987] | 1 | {\footnotesize |
---|
| 2 | \begin{verbatim} |
---|
| 3 | |
---|
| 4 | ##General options |
---|
| 5 | ##~~~~~~~~~~~~~~~ |
---|
| 6 | #Run with or without tracer transport ? |
---|
| 7 | tracer=.false. |
---|
| 8 | |
---|
| 9 | #Diurnal cycle ? if diurnal=False, diurnal averaged solar heating |
---|
| 10 | diurnal=.true. |
---|
| 11 | |
---|
| 12 | #Seasonal cycle ? if season=False, Ls stays constant, to value set in "start" |
---|
| 13 | season = .true. |
---|
| 14 | |
---|
| 15 | #write some more output on the screen ? |
---|
| 16 | lwrite = .false. |
---|
| 17 | |
---|
| 18 | #Save statistics in file "stats.nc" ? |
---|
| 19 | stats =.true. |
---|
| 20 | |
---|
| 21 | #Save EOF profiles in file "profiles" for Climate Database? |
---|
| 22 | calleofdump = .false. |
---|
| 23 | |
---|
| 24 | ## Dust scenario. Used if the dust is prescribed (i.e. if tracer=F or active=F) |
---|
| 25 | ## ~~~~~~~~~~~~~ |
---|
| 26 | # =1 Dust opt.deph read in startfi; =2 Viking scenario; =3 MGS scenario, |
---|
| 27 | # =4 Mars Year 24 from TES assimilation (same as =24 for now) |
---|
| 28 | # =24 Mars Year 24 from TES assimilation (ie: MCD reference case) |
---|
| 29 | # =25 Mars Year 25 from TES assimilation (ie: a year with a global dust storm) |
---|
| 30 | # =26 Mars Year 26 from TES assimilation |
---|
| 31 | iaervar = 24 |
---|
| 32 | # Dust vertical distribution: |
---|
| 33 | # (=0: old distrib. (Pollack90), =1: top set by "topdustref", |
---|
| 34 | # =2: Viking scenario; =3 MGS scenario) |
---|
| 35 | iddist = 3 |
---|
| 36 | # Dust top altitude (km). (Matters only if iddist=1) |
---|
| 37 | topdustref = 55. |
---|
| 38 | |
---|
| 39 | ## Physical Parameterizations : |
---|
| 40 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 41 | # call radiative transfer ? |
---|
| 42 | callrad = .true. |
---|
| 43 | # call NLTE radiative schemes ? matters only if callrad=T |
---|
| 44 | callnlte = .true. |
---|
| 45 | # call CO2 NIR absorption ? matters only if callrad=T |
---|
| 46 | callnirco2 = .true. |
---|
| 47 | # call turbulent vertical diffusion ? |
---|
| 48 | calldifv = .true. |
---|
| 49 | # call convective adjustment ? |
---|
| 50 | calladj = .true. |
---|
| 51 | # call CO2 condensation ? |
---|
| 52 | callcond =.true. |
---|
| 53 | # call thermal conduction in the soil ? |
---|
| 54 | callsoil = .true. |
---|
| 55 | # call Lott's gravity wave/subgrid topography scheme ? |
---|
| 56 | calllott = .true. |
---|
| 57 | |
---|
| 58 | ## Radiative transfer options : |
---|
| 59 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 60 | # the rad.transfer is computed every "iradia" physical timestep |
---|
| 61 | iradia = 1 |
---|
| 62 | # Output of the exchange coefficient mattrix ? for diagnostic only |
---|
| 63 | callg2d = .false. |
---|
| 64 | # Rayleigh scattering : (should be .false. for now) |
---|
| 65 | rayleigh = .false. |
---|
| 66 | |
---|
| 67 | ## Tracer (dust water, ice and/or chemical species) options (used if tracer=T): |
---|
| 68 | ## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 69 | # DUST: Transported dust ? (if >0, use 'dustbin' dust bins) |
---|
| 70 | dustbin = 0 |
---|
| 71 | # DUST: Radiatively active dust ? (matters if dustbin>0) |
---|
| 72 | active = .false. |
---|
| 73 | # DUST: use mass and number mixing ratios to predict dust size ? |
---|
| 74 | # (must also have dustbin=1) |
---|
| 75 | doubleq = .false. |
---|
| 76 | # DUST: lifted by GCM surface winds ? |
---|
| 77 | lifting = .false. |
---|
| 78 | # DUST: lifted by dust devils ? |
---|
| 79 | callddevil = .false. |
---|
| 80 | # DUST: Scavenging by CO2 snowfall ? |
---|
| 81 | scavenging = .false. |
---|
| 82 | # DUST/WATERICE: Gravitationnal sedimentation ? |
---|
| 83 | sedimentation = .false. |
---|
| 84 | # WATERICE: Radiatively active transported atmospheric water ice ? |
---|
| 85 | activice = .false. |
---|
| 86 | # WATER: Compute water cycle |
---|
| 87 | water = .false. |
---|
| 88 | # WATER: current permanent caps at both poles. True IS RECOMMENDED |
---|
| 89 | # (with .true., North cap is a source of water and South pole |
---|
| 90 | # is a cold trap) |
---|
| 91 | caps = .true. |
---|
| 92 | # PHOTOCHEMISTRY: include chemical species |
---|
| 93 | photochem = .false. |
---|
| 94 | |
---|
| 95 | ## Thermospheric options (relevant if tracer=T) : |
---|
| 96 | ##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 97 | # call thermosphere ? |
---|
| 98 | callthermos = .false. |
---|
| 99 | # WATER: included without cycle (only if water=.false.) |
---|
| 100 | thermoswater = .false. |
---|
| 101 | # call thermal conduction ? (only if callthermos=.true.) |
---|
| 102 | callconduct = .false. |
---|
| 103 | # call EUV heating ? (only if callthermos=.true.) |
---|
| 104 | calleuv=.false. |
---|
| 105 | # call molecular viscosity ? (only if callthermos=.true.) |
---|
| 106 | callmolvis = .false. |
---|
| 107 | # call molecular diffusion ? (only if callthermos=.true.) |
---|
| 108 | callmoldiff = .false. |
---|
| 109 | # call thermospheric photochemistry ? (only if callthermos=.true.) |
---|
| 110 | thermochem = .false. |
---|
| 111 | # date for solar flux calculation: (1985 < date < 2002) |
---|
| 112 | ## (Solar min=1996.4 ave=1993.4 max=1990.6) |
---|
| 113 | solarcondate = 1993.4 |
---|
| 114 | \end{verbatim} |
---|
| 115 | } |
---|