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). Furthermore, many scenarios are as a rule to be run. This explains the fact that the computational tasks in this situation are enormous. Therefore, it is necessary to select fast numerical methods; to develop parallel codes and, what is most important when the problems solved are very large to organize the computational process in a proper way...