A Predictor-Corrector Method for the Numerical Simulation of Sea Pollution

The behaviour of pollutants in the sea is described by parabolicequations with dominant advection terms, mixed derivatives andnon-linear interaction terms; these equations are solved bya special predictor-corrector scheme, consisting of 3 distinctsteps. The mixed derivative is split up into second derivativesalong the grid axes and an auxiliary axis. The 2 first steps(implicit and locally one dimensional) form the predictor. Thethird one is the corrector and yields a consistent discreteapproximation of the time derivative. The accuracy is O(²+h+h²)and a stability condition is derived.