A comparison of numerical methods for solving diffusion-reaction equations in air quality models