Many large mathematical models, not only models arising and used in environmental studies, are described by systems of partial differential equations. The discretization of the spatial derivatives in such models leads to the solution of very large systems of ordinary differential equations. These systems contain many millions of equations and have to be handled over large time intervals by applying many time-steps (up to several hundred thousand time-steps)...