Abstract

The finite-difference scheme for the Laplace and flux-divergence operators described in the companion paper (Part I) is consistent when applied on a grid consisting of perfect hexagons. The authors describe a necessary and sufficient condition for this finite-difference scheme to be consistent when applied on a grid consisting of imperfect hexagons and pentagons, and present an algorithm for generating a spherical geodesic grid on a sphere that guarantees that this condition is satisfied. Also, the authors qualitatively describe the error associated with the operators and estimate their order of accuracy when applied on the new grid.

This content is only available as a PDF.