A New Treatment of the Coriolis Terms in C-Grid Models at Both High and Low Resolutions

Abstract Numerical models of the ocean typically employ gridpoint techniques in which the dynamical variables defining the state of the ocean are held on a staggered grid. One common arrangement of the variables, known as the Arakawa C-grid, is particularly prone to gridscale noise that is due to spatial averaging of Coriolis terms and that is manifest when the grid resolution is coarse with respect to the deformation radius. Here, the authors analyze the problem in the context of linear inertia–gravity waves and discuss the reason for the prevalence of noise. They suggest a solution to the problem in which the C-grid model variables are augmented with D-grid velocity variables. An analysis of the resulting C–D grid indicates favorable behavior and numerical results are presented to demonstrate this. Finally, they discuss the similarity in nature between the C–D grid and the Z-grid, to explain why the C–D grid works well at both high and low resolution.

[1]  David A. Randall,et al.  Geostrophic Adjustment and the Finite-Difference Shallow-Water Equations , 1994 .

[2]  Akio Arakawa,et al.  Computational Design of the Basic Dynamical Processes of the UCLA General Circulation Model , 1977 .

[3]  A. E. Gill,et al.  Adjustment of the Ocean under Buoyancy Forces. Part I: The Role of Kelvin Waves , 1986 .

[4]  Roxana C. Wajsowicz,et al.  Free Planetary Waves in Finite-Difference Numerical Models , 1986 .

[5]  Rainer Bleck,et al.  A wind-driven isopycnic coordinate model of the north and equatorial Atlantic Ocean 2. The Atlantic Basin Experiments , 1990 .

[6]  A. Arakawa,et al.  Numerical methods used in atmospheric models , 1976 .

[7]  L. Perelman,et al.  A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers , 1997 .

[8]  Bruno Blanke,et al.  Variability of the Tropical Atlantic Ocean Simulated by a General Circulation Model with Two Different Mixed-Layer Physics , 1993 .

[9]  Michael S. Fox-Rabinovitz.  Computational Dispersion Properties of Horizontal Staggered Grids for Atmospheric and Ocean Models , 1991 .

[10]  R. C. Wajsowicz.  Adjustment of the Ocean under Buoyancy Forces. Part II: The Role of Planetary Waves , 1986 .

[11]  Dale B. Haidvogel,et al.  A semi-spectral primitive equation ocean circulation model using vertical sigma and orthogonal curvilinear horizontal coordinates , 1991 .

[12]  K. Bryan.  A Numerical Method for the Study of the Circulation of the World Ocean , 1997 .