Integration of reaction equations containing extremely different reaction rates -application to realistic active tracers in the troposphere
Abstract
If reaction rates differ by 18 orders of magnitude or more - as it happens with chemical reaction equations for active tracers in the atmosphere - codes just using double precision are bound to produce erroneous results. I apply a novel scaling procedure that allows application of i) double precision routines ii) fast extrapolation schemes(1) when integrating the reaction equations. (1) P. Deuflhard, F.A. Bornemann, Numerische Mathematik II, de Gruyter, Berlin 1994
- Publication:
-
EGS - AGU - EUG Joint Assembly
- Pub Date:
- April 2003
- Bibcode:
- 2003EAEJA....13499L