Locally-orthogonal unstructured grid-generation for general circulation modelling on the sphere*
Abstract. An algorithm for the generation of non-uniform, locally-orthogonal staggered unstructured spheroidal grids is described. This technique is designed to generate high-quality staggered Voronoi/Delaunay dual meshes appropriate for general circulation modelling on the sphere, including applications to atmospheric simulation, ocean-modelling and numerical weather prediction. Using a recently developed Frontal-Delaunay refinement technique, a method for the construction of unstructured spheroidal Delaunay triangulations is introduced. A locally-orthogonal polygonal grid, derived from the associated Voronoi diagram, is computed as the staggered dual. It is shown that use of the Delaunay-refinement technique allows for the generation of unstructured grids that satisfy a priori constraints on minimum mesh-quality. The initial staggered Voronoi/Delaunay tessellation is iteratively improved through hill-climbing optimisation techniques. Such an approach is shown to produce grids with very high element quality and smooth grading characteristics, while imposing relatively low computational expense. Initial results are presented for a selection of uniform and non-uniform spheroidal grids appropriate for high-resolution, multi-scale general circulation modelling. The use of user-defined mesh-spacing functions to generate smoothly graded, non-uniform grids for multi-resolution type studies is discussed in detail.