American Geophysical Union Become an AGU Member
Subscribe to AGU Journals		 AGU Home AGU Publications

Subscriber Access to Full Article (Nonsubscribers may purchase for $9.00, Includes print PDF, file size: 689655 bytes)

JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113, C05015, doi:10.1029/2007JC004482, 2008

Finite element ocean circulation model based on triangular prismatic elements, with application in studying the effect of topography representation

Q. Wang

Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, Germany


S. Danilov

Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, Germany


J. Schröter

Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, Germany


Abstract

Topographic steering of waves and currents is one of the most important issues of the general ocean circulation. Its representation in ocean models depends on the spatial discretization, which heavily influences the ocean circulation even on the global scale. Models capable of working on unstructured grids allow to represent topography and coastlines continuously and to refine mesh locally. Finite Element Ocean circulation Model (FEOM) belongs to this class of ocean general circulation models. Its discretization is based on unstructured triangular meshes on the surface and prismatic elements in the volume. The model uses continuous linear representation for the horizontal velocity, surface elevation, temperature and salinity, and solves the standard set of hydrostatic primitive equations. The method to stabilize the numerical pressure modes and the implementation of hybrid vertical grids are discussed. The model supports hybrid vertical grids including z and σ grids and their combination within the same numerical kernel. A method for calculating pressure gradient forces is introduced to reduce pressure gradient errors on inclined computational surfaces. Its performance is assessed via numerical experiments. The influence of different vertical discretization is illustrated using a seamount configuration. The model is MPI parallelized and has the functionality required for ocean general circulation modeling. It can be used as an efficient and versatile numerical tool for regional and global oceanographic applications and climate studies.

Received 3 August 2007; accepted 31 December 2007; published 10 May 2008.

Keywords: Finite element; ocean modeling; vertical discretization.

Index Terms: 4255 Oceanography: General: Numerical modeling (0545, 0560); 0550 Computational Geophysics: Model verification and validation; 4562 Oceanography: Physical: Topographic/bathymetric interactions.

Subscriber Access to Full Article (Nonsubscribers may purchase for $9.00, Includes print PDF, file size: 689655 bytes)

Citation: Wang, Q., S. Danilov, and J. Schröter (2008), Finite element ocean circulation model based on triangular prismatic elements, with application in studying the effect of topography representation, J. Geophys. Res., 113, C05015, doi:10.1029/2007JC004482.