Difference between revisions of "The traceur.def Input File"
From Planets
(Created page with "== What is it for? == Tracers need to be listed in ''traceur.def''. Characteristics of the tracers can be also specified in ''traceur.def''. == What is a tracer? == Tracer...") |
|||
| Line 38: | Line 38: | ||
</pre> | </pre> | ||
| − | An update version named '''ModernTrac-v1''' need to be used | + | An update version named '''ModernTrac-v1''' need to be used with chemistry in order to set up properly the chemical tracers. This gives also the opportunity to specify for each tracers specific characteristics. |
| + | |||
| + | <pre> | ||
| + | #ModernTrac-v1 | ||
| + | # | ||
| + | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| + | # !! README !! | ||
| + | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| + | # Welcome, this is the modern traceur.def ! | ||
| + | # | ||
| + | # ---> DON'T TOUCH/REMOVE THE FIRST LINE OF HEADER #ModernTrac-v1 <--- | ||
| + | # | ||
| + | # * You can use this area to write as many comments as you want starting lines with #. | ||
| + | # * You cannot have comment lines in the middle nor at the end, only at the top. | ||
| + | # | ||
| + | # | ||
| + | # 0. Philosophy This modern version of traceur.def gather in a single file what was | ||
| + | # ~~~~~~~~~~~~~ previously in traceur.def AND gases.def plus some hard-coded things. | ||
| + | # | ||
| + | # 1. Structure After the header, first line is the total number of modern tracers in | ||
| + | # ~~~~~~~~~~~~ your file, ie not only tracers sent to dynamics, but also species that | ||
| + | # takes part in the composition, radiative transfer, chemistry, etc. | ||
| + | # | ||
| + | # Then, each line needs to start with the specie name ... and that's the | ||
| + | # only mandatory information ! But you can specify all the options you | ||
| + | # want in separate blocks following 'option=value', assuming this option | ||
| + | # is defined in traceur_h.F90 & initracer.F for physics, infotrac.F90 | ||
| + | # for dynamics or chimiedata_h.F90 & calchim_asis.F90 for chemistry. | ||
| + | # Indeed this file is read once by dynamics, once by physics and once by | ||
| + | # chemistry who keep only the information needed. | ||
| + | # | ||
| + | # Note that by default a tracer listed below will be sent to dynamics | ||
| + | # except if you specify is_dyn=0. If nothing is given, then is_dyn=1. | ||
| + | # Not yet fully implemented. | ||
| + | # | ||
| + | # 3. Options. Implemented options listed below. | ||
| + | # ~~~~~~~~~~~~ For dynamic see "infotrac.F90". | ||
| + | # For physic see "initracer.F90". | ||
| + | # For chemistry see "calchim_asis.F90". | ||
| + | # init see "inichim_1D.F90" and "inichim_newstart.F90" | ||
| + | # | ||
| + | # Dynamic: vadv ! index of vertical trasport schema | ||
| + | # hadv ! index of horizontal trasport schema | ||
| + | # tnom_transp ! transporting fluid short name: CRisi | ||
| + | # | ||
| + | # Physic: mmol ! mole mass of tracer [g.mol-1] | ||
| + | # aki ! coefficient of thermal concduction | ||
| + | # cpi ! heat capacity [J.kg-1.K-1] | ||
| + | # is_chim ! 1 if tracer used in chemistry, else 0 | ||
| + | # is_rad ! 1 if tracer used in radiative transfer, else 0 | ||
| + | # is_recomb ! 1 if tracer used in recombining scheme, else 0 | ||
| + | # (if 1, must have is_rad=1) | ||
| + | # is_recomb_qset ! 1 if tracer k-corr assume predefined vmr, else 0 | ||
| + | # (if 1, must have is_recomb=1) | ||
| + | # is_recomb_qotf ! 1 if tracer recombination is done on-the-fly, else 0 | ||
| + | # (if 1, must have is_recomb_qset=0) | ||
| + | # | ||
| + | # Chemistry: SF_mode ! 1 if surface set up value, else 2 sedimentation velocity | ||
| + | # SF_value ! [vmr] if SF_mode=1, else [cm.s-1] | ||
| + | # prod_rate ! if SF_mode=2 production flux [molecules.m-2.s-1] | ||
| + | # | ||
| + | # init: qx ! value that initialize constant profile [vmr] | ||
| + | # qxf ! file that initialize profile [Pa,vmr] (1 line header) | ||
| + | # | ||
| + | # | ||
| + | # | ||
| + | # Random Ex : | ||
| + | # 4 | ||
| + | # n2 is_dyn=0 | ||
| + | # co2 cpi=450.0 is_dyn=1 mmol=44. | ||
| + | # hdo hadv=14 vadv=10 tnom_transp=air | ||
| + | # h2s | ||
| + | # | ||
| + | # Insert your tracers list and options below ! | ||
| + | # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | ||
| + | 29 | ||
| + | h2o_ice mmol=18 | ||
| + | co2 mmol=44 is_chim=1 SF_mode=1 SF_value=4.0e-4 is_rad=1 | ||
| + | co mmol=28 is_chim=1 | ||
| + | o mmol=16 is_chim=1 | ||
| + | o1d mmol=16 is_chim=1 | ||
| + | o2 mmol=32 is_chim=1 SF_mode=1 SF_value=2.1e-1 | ||
| + | o3 mmol=48 is_chim=1 is_rad=1 is_recomb=1 is_recomb_qotf=1 | ||
| + | h mmol=1 is_chim=1 | ||
| + | h2 mmol=2 is_chim=1 | ||
| + | oh mmol=17 is_chim=1 | ||
| + | ho2 mmol=33 is_chim=1 | ||
| + | h2o2 mmol=34 is_chim=1 | ||
| + | h2o_vap mmol=18 is_chim=1 is_rad=1 | ||
| + | ch4 mmol=16 is_chim=1 SF_mode=1 SF_value=1.8e-6 is_rad=1 | ||
| + | ch3 mmol=15 is_chim=1 | ||
| + | cho mmol=29 is_chim=1 | ||
| + | ch2o mmol=30 is_chim=1 | ||
| + | ch3o mmol=31 is_chim=1 | ||
| + | ch3o2 mmol=47 is_chim=1 | ||
| + | ch3o2h mmol=48 is_chim=1 | ||
| + | n mmol=14 is_chim=1 | ||
| + | no mmol=30 is_chim=1 | ||
| + | no2 mmol=46 is_chim=1 | ||
| + | n2 mmol=28 is_chim=1 is_rad=1 | ||
| + | hno3 mmol=63 is_chim=1 | ||
| + | hno4 mmol=79 is_chim=1 | ||
| + | n2o mmol=44 is_chim=1 SF_mode=1 SF_value=3.3e-7 | ||
| + | no3 mmol=62 is_chim=1 | ||
| + | n2o5 mmol=108 is_chim=1 | ||
| + | </pre> | ||
Revision as of 09:43, 11 May 2022
What is it for?
Tracers need to be listed in traceur.def.
Characteristics of the tracers can be also specified in traceur.def.
What is a tracer?
Tracers are elements which will be defined at each grid point and tracked in the model. They may interact with the atmosphere in different ways.
The model may include different types of tracers:
- condensed species (e.g., CO2, H2O, dust)
- chemically active species
- radiatively active gases (e.g., water vapor)
In the code, all tracers are stored in one three-dimensional array q, the third index of which corresponds to each individual tracer. In input (“start.nc”, “startfi.nc”) [1] and output files tracers are stored separately using their individual names. Loading specific tracers requires that the approriate tracer names are set in the traceur.def file, and specific computations for given tracers (e.g. computing the water or CO2 cycles) is controlled by setting the corresponding options in the callphys.def file [2].
Structure
Number of tracers need to be specified, followed by tracers names.
Example
The following example will set up water tracers to compute water cloud and water cycle in general.
2 h2o_vap h2o_ice
An update version named ModernTrac-v1 need to be used with chemistry in order to set up properly the chemical tracers. This gives also the opportunity to specify for each tracers specific characteristics.
#ModernTrac-v1 # # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # !! README !! # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ # Welcome, this is the modern traceur.def ! # # ---> DON'T TOUCH/REMOVE THE FIRST LINE OF HEADER #ModernTrac-v1 <--- # # * You can use this area to write as many comments as you want starting lines with #. # * You cannot have comment lines in the middle nor at the end, only at the top. # # # 0. Philosophy This modern version of traceur.def gather in a single file what was # ~~~~~~~~~~~~~ previously in traceur.def AND gases.def plus some hard-coded things. # # 1. Structure After the header, first line is the total number of modern tracers in # ~~~~~~~~~~~~ your file, ie not only tracers sent to dynamics, but also species that # takes part in the composition, radiative transfer, chemistry, etc. # # Then, each line needs to start with the specie name ... and that's the # only mandatory information ! But you can specify all the options you # want in separate blocks following 'option=value', assuming this option # is defined in traceur_h.F90 & initracer.F for physics, infotrac.F90 # for dynamics or chimiedata_h.F90 & calchim_asis.F90 for chemistry. # Indeed this file is read once by dynamics, once by physics and once by # chemistry who keep only the information needed. # # Note that by default a tracer listed below will be sent to dynamics # except if you specify is_dyn=0. If nothing is given, then is_dyn=1. # Not yet fully implemented. # # 3. Options. Implemented options listed below. # ~~~~~~~~~~~~ For dynamic see "infotrac.F90". # For physic see "initracer.F90". # For chemistry see "calchim_asis.F90". # init see "inichim_1D.F90" and "inichim_newstart.F90" # # Dynamic: vadv ! index of vertical trasport schema # hadv ! index of horizontal trasport schema # tnom_transp ! transporting fluid short name: CRisi # # Physic: mmol ! mole mass of tracer [g.mol-1] # aki ! coefficient of thermal concduction # cpi ! heat capacity [J.kg-1.K-1] # is_chim ! 1 if tracer used in chemistry, else 0 # is_rad ! 1 if tracer used in radiative transfer, else 0 # is_recomb ! 1 if tracer used in recombining scheme, else 0 # (if 1, must have is_rad=1) # is_recomb_qset ! 1 if tracer k-corr assume predefined vmr, else 0 # (if 1, must have is_recomb=1) # is_recomb_qotf ! 1 if tracer recombination is done on-the-fly, else 0 # (if 1, must have is_recomb_qset=0) # # Chemistry: SF_mode ! 1 if surface set up value, else 2 sedimentation velocity # SF_value ! [vmr] if SF_mode=1, else [cm.s-1] # prod_rate ! if SF_mode=2 production flux [molecules.m-2.s-1] # # init: qx ! value that initialize constant profile [vmr] # qxf ! file that initialize profile [Pa,vmr] (1 line header) # # # # Random Ex : # 4 # n2 is_dyn=0 # co2 cpi=450.0 is_dyn=1 mmol=44. # hdo hadv=14 vadv=10 tnom_transp=air # h2s # # Insert your tracers list and options below ! # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 29 h2o_ice mmol=18 co2 mmol=44 is_chim=1 SF_mode=1 SF_value=4.0e-4 is_rad=1 co mmol=28 is_chim=1 o mmol=16 is_chim=1 o1d mmol=16 is_chim=1 o2 mmol=32 is_chim=1 SF_mode=1 SF_value=2.1e-1 o3 mmol=48 is_chim=1 is_rad=1 is_recomb=1 is_recomb_qotf=1 h mmol=1 is_chim=1 h2 mmol=2 is_chim=1 oh mmol=17 is_chim=1 ho2 mmol=33 is_chim=1 h2o2 mmol=34 is_chim=1 h2o_vap mmol=18 is_chim=1 is_rad=1 ch4 mmol=16 is_chim=1 SF_mode=1 SF_value=1.8e-6 is_rad=1 ch3 mmol=15 is_chim=1 cho mmol=29 is_chim=1 ch2o mmol=30 is_chim=1 ch3o mmol=31 is_chim=1 ch3o2 mmol=47 is_chim=1 ch3o2h mmol=48 is_chim=1 n mmol=14 is_chim=1 no mmol=30 is_chim=1 no2 mmol=46 is_chim=1 n2 mmol=28 is_chim=1 is_rad=1 hno3 mmol=63 is_chim=1 hno4 mmol=79 is_chim=1 n2o mmol=44 is_chim=1 SF_mode=1 SF_value=3.3e-7 no3 mmol=62 is_chim=1 n2o5 mmol=108 is_chim=1
