Volker Grewe's Submodels and Tools for EMAC
Status: June, 16 2018
|O3ORIG||1.0||ready to use||Diagnostic to trace the geographic origin of ozone|
|TAGGING||1.1||ready to use||Diagnostic to trace the contribution of emission categories to the concentration of species such as NOy, O3, HOx|
|AirTraf||1.0||test phase||Emissions form air traffic (please contact H. Yamashita - DLR)|
|ProdLoss||1.0||operational||Diagnoses equation-file and identifies all production and loss terms of a family|
Back to homepage
last change: May, 4 2011
Description follows MESSy2
Diagnostic for a source-receptor relationship for ozone.
Ozone source regions are defined and a tracer to each allocated, which shows how far this ozone is transported through the atmosphere.
(- Correct implementation of the prod and loss terms; done available in MESSY2.5)
- More flexible region definition
- More flexible wrt. species "O3ORIG" -> "SpecieORIG"
set diagtracfile = o3orig in messy/mbm/caaba/mecca/batch/<yourfile>.bat
or choose o3orig.bat
This guaranties that ozone pro and loss terms are available for the submodel.
- File describing the source regions.
(If you want to change the regions, you have to change the input file + submodel O3ORIG
- Chemical equations defining prod and loss of ozone
a) use ../messy/mbm/caaba/mecca/diagtrac/o3orig.tex:
activate respective equations ProdO3 and LossO3
(Warning: Correct only for MESSY>2.5)
b) No changes in gas.spc and tracdef.tex necessary!
a) import.nml :
activate lines to initiate import of region definition
!!$ RG_TRIG(75) = 1,'years', 'first',0, 'O3ORIG', 1,1,1,1, 'GP=0;
description of region-definition file
! -*- f90 -*-!
! names of regions
sn_o3orig = 'O3ONHTS','O3OTRTS','O3OSHTS','O3ONPLS','O3ONMLS','O3OTRLS','O3OTRMS','O3OSMLS',
- Additional ozone tracers in *tracer_gp.nc for each region.
- Total ozone loss and prod. rates in *tracer_gp.nc
- Additional diagnostics in *o3orig*.nc
Needs submoduls: MECCA1, ++?
- Switch on Module-switch USE_O3ORIG in switch.nml
- Create region definition file, according to your requirements and adapt the submodel
(initialisation and integration)
- Adapt namelists (see above)
- Error tracers not yet tested.
- hopefully not
- Use ProdLoss for the definition of ozone loss and production terms.
Stefanie Meul, Sophie Oberländer, FU Berlin,
Hella Garny, DLR-Oberpfaffenhofen.
Grewe, 2006, The origin of ozone, ACP 6, 2006.
Garny et al., Attribution of ozone changes to dynamical and chemical processes in CCMs and CTMs, Geosci. Model Dev., 4, 271-286, 2011. (See appendix!)
last change: April, 14 2011
last change: April, 14 2011
last change: April, 26 2011
Diagnoses a file with chemical reactions and identifies all production and loss terms of a family, based on Crutzen and Schmailzl, 1983.
Crutzen and Schmailzl, Chemical Budgets of the Stratosphere, Planet. .Space Sci., Vol. 31, No. 9, pp. 1009-1032, 1983.
- KPP equation file: mecca1.eqn, gas_mim1.eqn, etc.
! Name of input file
! Name of output file, which can be icluded in
the diagtrac.tex file (outfile.tex)
! Species, which define the family (max. 30)
! weitghting factor for each species:
! Version of MECCA =1 or =2
mecca = 1
- $outfile.tex: List of reactions, which affect the family
- $outfile.log: List of all reactions and description of why they are (not) regarded.
This file includes at the end also warnings!
prodloss.x < namelist
- Warning: species not found in equations
Check whether the species is actually in the file
Check whether species is correctly written: HOCl and ClOH are different!
Check whether blanks directly before and behind
the species are missing in the equation file.
- Warning: Error found in reaction
A line is identified as an equation but it is not in the right format
- The equation file needs a line "#EQUATION" other wise no equations are identified
- Program is tested, but not fool-proof, please check results carefully
- Don't forget to include the variables in gas.spc and tracdef.tex
- Variables can be easily separated into classes, e.g. ozone loss from HOx, etc.
Rolf Sander, MPI-Mainz; Patrick Jöckel, DLR-Oberpfaffenhofen
download the zip-file, which contains prodloss.f90 and prodloss.nml: prodloss
or see MESSY/messy_2.41/messy/tools/prodloss