An algorithm for ODEs from atmospheric dispersion problems Original Research Article

The solution of large systems of ordinary differential equations (ODEs), arising from atmospheric dispersion problems is considered. An algorithm using a method due to Klopfenstein is adopted as the main method and combined with an approximate Jacobian factorization and a Gauss-Seidel iteration to provide an efficient solver. The approach is contrasted with that of using implicit-explicit multistep methods. Numerical experiments are presented to illustrate the performance of the method.