KPP solvers FAQ
Welcome to the "KPP and GEOS-Chem - Frequently Asked Questions" page.
- 1 What is KPP?
- 2 What are the pros of using KPP?
- 3 What are the cons of using KPP?
- 4 Can I still use the KPP-generated solver if I modify GEOS-Chem?
- 5 Which GC simulation can use the KPP solvers?
- 6 Can I still use SMVGEAR?
- 7 Which KPP solvers can I use?
- 8 How do I switch solver?
- 9 How do I choose the absolute and relative tolerance?
- 10 How does the default Rodas-3 Rosenbrock integrator compare with the SMVGEAR ?
- 11 References
What is KPP?
The Kinetic PreProcessor is a software that automatically generates F90 code that solves the chemistry (defined by input files) in a box model.
The KPP input files have been generated from the globchem.dat GEOS-Chem file with a Perl script, and KPP was run. The Kpp-generated code has been interfaced with GEOS-Chem.
KPP has been run twice for the standard NOx-Ox-VOC GEOS-Chem simulation: once for the 43 transported tracers model, and once for the 54 tracers version (ie with in-line secondary aerosols), since the two models have a slightly different globchem.dat.
What are the pros of using KPP?
In one word, faster and better. In many points:
- It gives you a choice of different solvers that can have relative advantages depending on the level of desired accuracy [Sandu and Sander, 2006].
- One of the solver (radau5) is also recommended for "accurate reference solutions" [Sandu and Sander, 2006].
- By default, GEOS-Chem is compiled with the Rosenbrock Rodas3 solver:
- Tests reported in Henze et al.  and Eller et al.  show that the Rosenbrock solver is twice faster than SMVGEAR for the same accuracy. However this accuracy is with respect to a reference output obtained with the same solver.
- "At the very least, [using the Rosenbrock solver] results in an improvement in the numerical solution [...] for slightly less computational cost." [Henze et al., 2007]
- At low NOx regimes, like paleo simulations, the Rosenbrock solver is found more stable than SMVGEAR (L. Murray, personal communication, GEOS-Chem v8-01-04, 2009).
To get a measure of (3.3), a detailed comparison between SMVGEAR, Rosenbrock (Rodas3), and a reference solution [see pdf] was performed. It shows that Rodas3 can be almost 10% faster for a more accurate solution. SMVGEAR is particularly bad at reproducing isoprene and the inorganic sulfur nitrates.
Here is the full list of transported species that show an RMS error larger than 10% for SMVGEAR and/or Rodas3. The numbers are RMS of relative error in %, from a comparison with a reference run, obtained with a Runge-Kutta order 5 integrator and very tight tolerances. All runs are for 13 days without transport on 4x5, the comparison is done with concentrations averaged over the 13th day.
SMVGEAR Rodas3 ==================== ISOP 222.22 NIT 87.36 11.82 MVK 33.26 MACR 30.40 PMN 23.91 NH3 19.10 DMS 11.43 ALK4 10.72
What are the cons of using KPP?
If you are changing the chemistry (species and/or reactions): you have to run KPP each time you modify globchem.dat to generate the gckpp*.f90 files. Then these files have to be adapted for GC, and replace the previous ones in GC code source.
Generally, if you are developing GC, be aware that the interface with KPP variables is done through both CSPEC and RRATE arrays (see physproc.f, calcrate.f, chemdr.f, and chemistry_mod.f).
Can I still use the KPP-generated solver if I modify GEOS-Chem?
Yes. See "What are the cons of using KPP?"
Which GC simulation can use the KPP solvers?
For now, only the standard NOx-Ox-VOC simulations, with in- or off-line secondary aerosols. As of v8-02-03, the difference is in the number of transported tracers: 43 and 54.
By default, the code is compiled for 43 tracers. For the 54 tracers simulation, you can use the NTRAC flag when calling make:
When switching between 43 and 54, you need to call:
Can I still use SMVGEAR?
yes. There is a switch in the chemistry menu of the input.geos to choose between SMVGEAR and the KPP-derived solver at runtime.
Which KPP solvers can I use?
As of GEOS-Chem v8-02-03, you can choose between several stiff numerical solvers that depend on either the backward differentiation formula (BDF), Runge Kutta, or Rosenbrock methods.
Runge-Kutta (Fully Implicit 3-stage Runge-Kutta methods), based on several quadratures (*): Radau-2A quadrature (order 5)  (default, stiffly accurate, robust, expensive) Radau-1A quadrature (order 5)  Lobatto-3C quadrature (order 4)  Gauss quadrature (order 6) 
LSODES (BDF, Livermore ODE solver, sparse array version)
Rosenbrock methods (*): Ros-2  Ros-3  Ros-4  Rodas-3  (default) Rodas-4 
(*) For the Rosenbrock and Runge-Kutta methods, the default method can be overwritten by setting (before compilation, in chemistry_mod.f) the variable ictrl(3) to the number in square brackets.
For details on those methods and their implementation, see Sandu and Sander  and its electronic supplement, and references therein.
How do I switch solver?
Solvers are chosen at three steps: before compilation, at compilation, and at runtime.
At runtime, you can switch between SMVGEAR and the KPP-derived solver with a single flag in the input.geos file. See the GEOS-Chem manual.
The KKP solver is chosen at compilation. When calling make, you can overwrite the default rosenbrock solver, with the KPPSOLVER flag:
make KPPSOLVER=<solver> where <solver> must be one of: radau5, lsodes, runge_kutta, rosenbrock (radau5 is a standalone Runge-Kutta solver of order 5 with Radau-2A quadrature)
You can further fine tune your choice of integrator by setting ictrl(3) in chemistry_mod.f before compilation. See "Which KPP solvers can I use?" above.
How do I choose the absolute and relative tolerance?
The tolerances used by solvers to check convergence are set differently for SMVGEAR and the KPP solvers.
For SMVGEAR, the range of absolute tolerance (ATOL) is defined in the mglob.dat file by YLOWU and YHIU. The relative tolerance (RTOL) is also set in mglob.dat with ERRMAXU, but it is overwritten in the code if you set it above 1d-3 (this is its max value, in other words).
For KPP-derived solvers, the tolerance are set in gckpp_driver routine (in chemistry_mod.f). The code is distributed with the suggested value of: RTOL=2d-1 and ATOL=1d-2.
For the Rodas-3 Rosenbrock solver, time and accuracy of the solution for various set of tolerances are described in this pdf. Some results are reproduced in <How does the default Rodas-3 Rosenbrock integrator compare with the SMVGEAR ?>
If you switch to Livermore solver (LSODE), you do need to set ATOL to 1d6 or around.
For an accurate solution with Radau5, you need small tolerances: RTOL=1d-8 for example.
How does the default Rodas-3 Rosenbrock integrator compare with the SMVGEAR ?
Faster by 10% for a more accurate solution. See this pdf for details. Here we summarize some of the findings:
Absolute tolerance ATOL has little to no effect in the [1e-2,1e+2] examined range. Using ATOL=1d-2, here is a comparison between SMVGEAR and Rodas-3 with several relative tolerances RTOL. The RMS of relative error in % of species with an error larger than 10%, and the runtime with respect to SMVGEAR are:
Rodas3 SMVGEAR Rodas3 Rodas3 RTol=0.3 RTol=0.2 RTol=0.1 ======================================================================= ISOP 222.22 NIT 14.75 87.36 16.33 11.82 MVK 33.26 MACR 30.40 PMN 23.91 NH3 19.10 DMS 11.43 ALK4 10.72 H2O2 2.23e9 HNO4 26.12 PAN 23.43 PPN 11.47 runtime 87% 100% 92% 93% notes recommended used in the adjoint-code
Eller et al., Geosci. Model Dev., 2, 89-96, 2009
Henze et al., ACP, 7 ,2413-2433, 2007
Sandu and Sander, ACP, 6, 187-195, 2006
--phs 16:54, 17 September 2009 (EDT)