Precision implicit second order differences for numerical weather simulation

Optimal approximations are derived for locally linearized partial differential equations with matrices of the same form as Galerkin B-spline approximations. The matrix coefficients are chosen for correct frequency and group velocity of long waves and at four mesh lengths. An overall error of 1% for intermediate wavelengths is possible. Splitting methods are found to handle advection terms relatively realistically.