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.

   e.g. $INPUTDIR_MESSY/offlem/regions4ozone_T42_7regions.nc

  (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!

- Namelists:

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;  

                                    NML=./offlem/regions4ozone_origin.nml',
b) ./offlem/regions4ozone_origin.nml

     description of region-definition file

c) o3orig.nml

! -*- f90 -*-!
&CTRL
! NOTES:
! - 1: Euler backward (default)
! - 2: Euler forward
!i_integrate = 1 ! op_pj_20100922
!
! number of tracers (=regions+2errors, if l_err=t else =regions) ! op_pj_20100922
i_trac_orig = 14

! names of regions

sn_o3orig = 'O3ONHTS','O3OTRTS','O3OSHTS','O3ONPLS','O3ONMLS','O3OTRLS','O3OTRMS','O3OSMLS',

'O3OSPLS','O3ONPUS','O3ONMUS','O3OTRUS','O3OSMUS','O3OSPUS'

! ,'O3OERRP','O3OERRN'
l_err = .false. ! Switch on error diagnostics
/
&CPL
! NOTES:
! - WMO : 'tropop', 'tp_p_wmo'
! - PV : 'tropop', 'tp_p_pv'
! - combined : 'tropop', 'tp_p' ! op_pj_20100922
! - average : 'offlem', 'tpmean' ! op_pj_20100922
c_tropop = 'tropop', 'tp_p_wmo'
!

- 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!)

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.

- Namelist

  ! Name of input file

        infile ="gas_mim1.eqn",   

  ! Name of output file, which can be icluded in the diagtrac.tex file (outfile.tex)
        outfile ="diagtrac.tex",           !

  ! Species, which define the family (max. 30)
        fam = "O3", "O3P", "O1D", "ClO" , "HOCl", "Cl2O2", "OClO",
          "NO2", "HNO3", "NO3" , "N2O5", "HNO4", "PAN", "NACA" ,
          "MPAN","IC3H7NO3","LC4H9NO3", "ISON", "ClNO3", "BrO","HOBr",
          "BrNO3"

   ! weitghting factor for each species:
       fac = 1.,1.,1.,1.,1.,2.,1.,
          1.,1.,2.,3.,1.,1.,1.,
          1.,1.,1.,1.,2.,1.,1.,
          2.,
  ! Name of family. Loss and Prod terms are then called Loss$cfam Prod$cfam

      cfam ="O3"

 ! 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!

Compile prodloss.f90    

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