A global shallow‐water model using an expanded spherical cube: Gnomonic versus conformal coordinates

A model using shallow-water equations with an Arakawa-type scheme for momentum terms is tested on a quasi-uniform geometry on the sphere, derived by a spherical expansion of the inscribed cube based on the gnomonic projection. Thereby, a quasi-homogeneous distribution of grid points is achieved, and a global finite-difference model is designed which does not require Fourier filtering or suffer from the burden of redundant computational points at high latitudes. Difficulties resulting from the directional discontinuity of the coordinate lines crossing the edges of the expanded cube are almost completely eliminated by using the Arakawa B-grid, so that only scalar points are placed along the edges. An alternative approach is developed based on numerical orthogonalization of the grid whereby, inter alia, the directional discontinuity at the edges is avoided at the cost of some accumulation of points in the vicinity of the vertices of the cube. In the customary Rossby-Haurwitz wave-4 tests, both approaches are shown to converge to a visually indistinguishable solution as the resolution is increased. However, with the orthogonalized, conformal grid, convergence towards the asymptotic solution was substantially faster.

[1]  A. Arakawa Computational design for long-term numerical integration of the equations of fluid motion: two-dimen , 1997 .

[2]  Akio Arakawa,et al.  Integration of the Nondivergent Barotropic Vorticity Equation with AN Icosahedral-Hexagonal Grid for the SPHERE1 , 1968 .

[3]  David L. Williamson,et al.  Integration of the barotropic vorticity equation on a spherical geodesic grid , 1968 .

[4]  David L. Williamson,et al.  INTEGRATION OF THE PRIMITIVE BAROTROPIC MODEL OVER A SPHERICAL GEODESIC GRID , 1970 .

[5]  Richard Asselin,et al.  Frequency Filter for Time Integrations , 1972 .

[6]  R. Sadourny Conservative Finite-Difference Approximations of the Primitive Equations on Quasi-Uniform Spherical Grids , 1972 .

[7]  Fedor Mesinger,et al.  A method for construction of second-order accuracy difference schemes permitting no false two-grid-interval wave in the height field , 1973 .

[8]  Zavisa Janjic,et al.  A Stable Centered Difference Scheme Free of Two-Grid-Interval Noise , 1974 .

[9]  M. J. P. Cullen,et al.  Integrations of the primitive equations on a sphere using the finite element method , 1974 .

[10]  Brian J. Hoskins,et al.  A multi-layer spectral model and the semi-implicit method , 1975 .

[11]  Joseph Sela,et al.  Spectral Modeling at the National Meteorological Center , 1980 .

[12]  Zavisa Janjic,et al.  Nonlinear Advection Schemes and Energy Cascade on Semi-Staggered Grids , 1984 .

[13]  Lawrence L. Takacs,et al.  A Two-Step Scheme for the Advection Equation with Minimized Dissipation and Dispersion Errors , 1985 .

[14]  Fourth-Order Horizontal Advection Schemes on the Semi-staggered Grid , 1988 .

[15]  Slobodan Nickovic,et al.  The Step-Mountain Coordinate: Model Description and Performance for Cases of Alpine Lee Cyclogenesis and for a Case of an Appalachian Redevelopment , 1988 .

[16]  Saulo R. M. Barros,et al.  Integration of the shallow water equations on the sphere using a vector semi-Lagrangian scheme with a multigrid solver , 1990 .

[17]  Zavisa Janjic,et al.  The Step-Mountain Coordinate: Physical Package , 1990 .

[18]  Jean Côté,et al.  An Accurate and Efficient Finite-Element Global Model of the Shallow-Water Equations , 1990 .

[19]  P. Swarztrauber,et al.  A standard test set for numerical approximations to the shallow water equations in spherical geometry , 1992 .

[20]  R. W. Higgins,et al.  A Global Multilevel Atmospheric Model Using a Vector Semi-Lagrangian Finite-Difference Scheme. Part I: Adiabatic Formulation , 1993 .

[21]  Extension and Application of a Local, Minimum Aliasing Method to Multidimensional Problems in Limited-Area Domains , 1993 .

[22]  R. Heikes,et al.  Numerical Integration of the Shallow-Water Equations on a Twisted Icosahedral Grid , 1995 .

[23]  Z. Janjic The Step-Mountain Eta Coordinate Model: Further Developments of the Convection, Viscous Sublayer, and Turbulence Closure Schemes , 1994 .

[24]  T. Black The new NMC mesoscale Eta Model: description and forecast examples , 1994 .

[25]  Pier Stanislao Paolucci,et al.  Finite difference approximation to the shallow water equations on a quasi-uniform spherical grid , 1995, HPCN Europe.

[26]  William A. Gallus,et al.  A non-hydrostatic version of the NMC's regional Eta model , 1996 .

[27]  John L. McGregor Semi-Lagrangian Advection on a Cubic Gnomonic Projection of the Sphere , 1997 .