Abstract
A spectral finite-volume (SFV) method is proposed for the numerical solution of the shallow water equations. This is the first phase in the development of a layered (isopycnal) ocean model. Its target applications include, in particular, the simulation of the wind-driven oceanic circulation in geometrically complex basins where layer outcropping and/or isopycnal-bathymetry intersection must be handled explicitly. The present formulation is geometrically flexible and can extend accuracy to arbitrary high order with no change to the basic algorithm. A flux-corrected transport (FCT) algorithm ensures the stability of the computations in regions of vanishing layer thickness and in areas where the flow features are underresolved. The spatial discretization is based on a two-level grid: a globally unstructured elemental grid and a locally structured grid consisting of N x N quadrilateral cells within each element. The numerical solution is continuous within each element but discontinuous across elements; the discontinuity is resolved by upwinding along characteristics. The accuracy and convergence rate of the SFV method are verified on two linearized problems amenable to analytical solution; the SFV solution exhibits a convergence order of N + 1 for smooth solutions. The FCT portion of the model is tested by simulating the formation of an oblique hydraulic jump in a supercritical channel flow. The model is then applied to simulate, in reduced-gravity mode, the double-gyre and wind-driven upper-ocean circulations in a square basin. Finally, the previous experiment is repeated in the North Atlantic basin to illustrate the application of the model in a realistic geometry.