Numerical integration of the barotropic vorticity equation