The performance of implicit ocean models on B- and C-grids

Fully-implicit primitive equation ocean models are useful to study the sensitivity of steady ocean flows to parameters, to determine bifurcations of these flows associated with instabilities and to use relatively large time steps in transient flow computations. This paper addresses a problem related to the origin of wiggles occurring in fully-implicit C-grid models. The situation considered is the computation of three-dimensional thermally-driven steady flows in a midlatitude spherical sector. We determine the reason why in a coarse resolution C-grid implicit model, the values of the lateral friction coefficients are restricted to far higher values than for the same B-grid model. The analysis also reveals why the B-grid discretization is superior for the computation of this type of flows.

[1]  W. Weijer,et al.  Stability of the global ocean circulation: The connection of equilibria within a hierarchy of models , 2003 .

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

[3]  Fred Wubs,et al.  A fully implicit model of the three-dimensional thermohaline ocean circulation: 685 , 2001 .

[4]  J. Willebrand,et al.  Numerical Ocean Circulation Modelling: Present Status and Future Directions , 2001 .

[5]  Fred Wubs,et al.  A fully-implicit model of the global ocean circulation , 2003 .

[6]  R. Sadourny.  The Dynamics of Finite-Difference Models of the Shallow-Water Equations , 1975 .

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

[8]  A. Arakawa.  Design of the UCLA general circulation model , 1972 .

[9]  H. Dijkstra,et al.  Instability of the Thermohaline Ocean Circulation on Interdecadal Timescales , 2002 .

[10]  R. Seydel.  Practical bifurcation and stability analysis : from equilibrium to chaos , 1994 .

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

[12]  H. Dijkstra,et al.  Identification of the Mechanism of Interdecadal Variability in the North Atlantic Ocean , 2004 .

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

[14]  Alberto M. Mestas-Nuñez,et al.  The Atlantic Multidecadal Oscillation and its relation to rainfall and river flows in the continental U.S. , 2001 .

[15]  J. Ozer,et al.  Numerical boundary layers and spurious residual flows , 1986 .

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

[17]  S. Manabe,et al.  A Global Ocean-Atmosphere Climate Model. Part II. The Oceanic Circulation , 1975 .

[18]  P. Gent,et al.  Parameterizing eddy-induced tracer transports in ocean circulation models , 1995 .

[19]  E. F. F. Botta,et al.  Matrix Renumbering ILU: An Effective Algebraic Multilevel ILU Preconditioner for Sparse Matrices , 1999, SIAM J. Matrix Anal. Appl..

[20]  Alistair Adcroft,et al.  A New Treatment of the Coriolis Terms in C-Grid Models at Both High and Low Resolutions , 1999 .

[21]  A. Arakawa,et al.  A Potential Enstrophy and Energy Conserving Scheme for the Shallow Water Equations , 1981 .

[22]  W. Weijer,et al.  A systematic approach to determine thresholds of the ocean’s thermohaline circulation , 2004 .