The work reported is in pursuit of these goals: high-quality unstructured, non-uniform Voronoi and Delaunay grids; improved finite element and finite volume discretization schemes; and improved finite element and finite volume discretization schemes. These are sought for application to spherical and three-dimensional applications suitable for ocean, atmosphere, ice-sheet, and other climate modeling applications.