Difference between revisions of "Sandbox"
From Geos-chem
| (156 intermediate revisions by 3 users not shown) | |||
| Line 1: | Line 1: | ||
| − | == | + | == Instructions for adding tracers to GEOS-Chem == |
| − | + | '''''By [mailto:mpayer@seas.harvard.edu Melissa Payer]''''' | |
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | These instructions describe how to add advected tracers to GEOS-Chem. Advected tracers are carried in the <tt>STT</tt> array (which is declared in <tt>GeosCore/tracer_mod.F</tt>), and are transported by the winds. | |
| − | < | + | |
| − | < | + | |
| − | < | + | |
| − | + | ||
| − | + | NOTE: Some advected tracers are also chemical species, which comprise GEOS-Chem's NOx-Ox-HC-aerosol chemistry mechanism. | |
| − | + | ||
| − | + | ||
| − | < | + | #<p>Add tracers to input.geos under Tracer Menu</p> |
| − | < | + | ##<p>Increase “Number of tracers” accordingly</p> |
| − | < | + | ##<p>Add tracer number, name, molecular weight, and emitted species to end of tracer list</p> |
| + | #<p>Add tracer(s) to restart file – see separate instructions</p> | ||
| + | #<p>Modify the following source code routines:</p> | ||
| + | ##<p><tt>Headers/CMN_SIZE_mod.F</tt></p> | ||
| + | ###<p>Increase <tt>NNPAR</tt> by number of tracers being added</p> | ||
| + | ###<p>If the species have anthropogenic emissions, increase NEMPARA accordingly</p> | ||
| + | ###<p>If the species have biogenic emissions, increase NEMPARB accordingly</p> | ||
| + | ###<p>If the species have biomass burning emissions, increase NBIOMAX accordingly</p> | ||
| + | ##<p><tt>GeosCore/tracerid_mod.F</tt></p> | ||
| + | ###<p>Increase <tt>NNNTRID</tt> by the number of tracers being added NOTE: <tt>NNNTRID</tt> must be equal to <tt>NNPAR</tt> in <tt>CMN_SIZE_mod.F</tt></p> | ||
| + | ###<p>Define tracer IDs as needed:</p> | ||
| + | ####<p>IDxxxx – ID for SMVGEAR species</p> | ||
| + | ####<p>IDTxxxx – GEOS-Chem tracer ID</p> | ||
| + | ####<p>IDExxxx – GEOS-Chem emission ID</p> | ||
| + | ####<p>IDBFxxxx – GEOS-Chem biofuel ID</p> | ||
| + | ####<p>IDBxxxx – GEOS-Chem biomass ID</p> | ||
| + | ###<p>In subroutine <tt>TRACERID</tt>,</p> | ||
| + | ####<p>Under “Assign tracer, biomass, biofuel, and anthro emission ID’s,” add case for added tracers to initialize IDTxxxx, IDBFxxxx, IDBxxxx</p> | ||
| + | ####<p>Under “Initialize the IDExxx flags,” add IF statement for added tracer IDs and initialize IDExxxx</p> | ||
| + | ####<p>Under “Fill IDEMS with appropriate tracer ID #’s,” add IF statement for added tracers</p> | ||
| + | ####<p>Under “Print additional information,” add WRITE statements for added tracers</p> | ||
| + | ###<p>In subroutine SETTRACE (NOTE: Only do the following if species is included in SMVGEAR chemical mechanism<p> | ||
| + | ####<p>Add IF statement, to initialize IDxxxx</p> | ||
| + | ####<p>In subroutine INIT_TRACERID,</p> | ||
| + | ####<p>Zero all IDs added (IDxxxx, IDTxxxx, IDExxxx, etc.)</p> | ||
| + | ##<p>Other modules that may need to be updated include:</p> | ||
| + | ###<p>wetscav_mod.F</p> | ||
| + | ###<p>drydep_mod.F</p> | ||
| + | ###<p>emissions_mod.F</p> | ||
| + | ###<p>biomass_mod.F</p> | ||
| + | ###<p>biofuel_mod.F</p> etc. | ||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | |||
| − | + | Instructions for adding species to GEOS-Chem | |
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | Modify run directory files | |
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | 1) Add tracers to input.geos under Tracer Menu | |
| + | c) Increase “Number of tracers” accordingly | ||
| + | d) Add tracer number, name, molecular weight, and emitted species to end of tracer list | ||
| − | + | 2) Add tracer(s) to restart file – see separate instructions | |
| − | + | ||
| − | + | ||
| − | + | 3) Modify globchem.dat to include additional species, kinetic reactions (including emissions and drydep), or photolysis reactions | |
| − | + | NOTE: If you modify globchem.dat and plan to run GEOS-Chem using the KPP chemical solver, you will need to generate new gckpp*.F90 files. For more information see the following wiki pages: | |
| − | + | http://wiki.seas.harvard.edu/geos-chem/index.php/KPP_solvers_FAQ#What_are_the_cons_of_using_KPP.3F | |
| − | + | http://wiki.seas.harvard.edu/geos-chem/index.php/Interfacing_GEOS-Chem_with_KPP#Generating_KPP_input_files_from_GEOS-Chem_globchem.dat | |
| − | + | ||
| − | + | 4) If necessary, modify: | |
| − | + | a) ratj.d – contains species names and branching ratios for FAST–JX photolysis species | |
| − | + | For more information: http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#5.4.1 | |
| − | + | b) jv_spec.dat – contains cross-sections and quantum yields for FAST–JX photolysis species | |
| − | + | For more information: http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#5.4.3 | |
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
| − | + | ||
Latest revision as of 21:25, 10 May 2023
Instructions for adding tracers to GEOS-Chem
These instructions describe how to add advected tracers to GEOS-Chem. Advected tracers are carried in the STT array (which is declared in GeosCore/tracer_mod.F), and are transported by the winds.
NOTE: Some advected tracers are also chemical species, which comprise GEOS-Chem's NOx-Ox-HC-aerosol chemistry mechanism.
Add tracers to input.geos under Tracer Menu
Increase “Number of tracers” accordingly
Add tracer number, name, molecular weight, and emitted species to end of tracer list
Add tracer(s) to restart file – see separate instructions
Modify the following source code routines:
Headers/CMN_SIZE_mod.F
Increase NNPAR by number of tracers being added
If the species have anthropogenic emissions, increase NEMPARA accordingly
If the species have biogenic emissions, increase NEMPARB accordingly
If the species have biomass burning emissions, increase NBIOMAX accordingly
GeosCore/tracerid_mod.F
Increase NNNTRID by the number of tracers being added NOTE: NNNTRID must be equal to NNPAR in CMN_SIZE_mod.F
Define tracer IDs as needed:
IDxxxx – ID for SMVGEAR species
IDTxxxx – GEOS-Chem tracer ID
IDExxxx – GEOS-Chem emission ID
IDBFxxxx – GEOS-Chem biofuel ID
IDBxxxx – GEOS-Chem biomass ID
In subroutine TRACERID,
Under “Assign tracer, biomass, biofuel, and anthro emission ID’s,” add case for added tracers to initialize IDTxxxx, IDBFxxxx, IDBxxxx
Under “Initialize the IDExxx flags,” add IF statement for added tracer IDs and initialize IDExxxx
Under “Fill IDEMS with appropriate tracer ID #’s,” add IF statement for added tracers
Under “Print additional information,” add WRITE statements for added tracers
In subroutine SETTRACE (NOTE: Only do the following if species is included in SMVGEAR chemical mechanism<p>
- <p>Add IF statement, to initialize IDxxxx
In subroutine INIT_TRACERID,
Zero all IDs added (IDxxxx, IDTxxxx, IDExxxx, etc.)
Other modules that may need to be updated include:
wetscav_mod.F
drydep_mod.F
emissions_mod.F
biomass_mod.F
biofuel_mod.F
etc.
Instructions for adding species to GEOS-Chem
Modify run directory files
1) Add tracers to input.geos under Tracer Menu c) Increase “Number of tracers” accordingly d) Add tracer number, name, molecular weight, and emitted species to end of tracer list
2) Add tracer(s) to restart file – see separate instructions
3) Modify globchem.dat to include additional species, kinetic reactions (including emissions and drydep), or photolysis reactions NOTE: If you modify globchem.dat and plan to run GEOS-Chem using the KPP chemical solver, you will need to generate new gckpp*.F90 files. For more information see the following wiki pages: http://wiki.seas.harvard.edu/geos-chem/index.php/KPP_solvers_FAQ#What_are_the_cons_of_using_KPP.3F http://wiki.seas.harvard.edu/geos-chem/index.php/Interfacing_GEOS-Chem_with_KPP#Generating_KPP_input_files_from_GEOS-Chem_globchem.dat
4) If necessary, modify: a) ratj.d – contains species names and branching ratios for FAST–JX photolysis species For more information: http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#5.4.1 b) jv_spec.dat – contains cross-sections and quantum yields for FAST–JX photolysis species For more information: http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#5.4.3
