Journal metrics

Journal metrics

  • IF value: 4.252 IF 4.252
  • IF 5-year value: 4.890 IF 5-year 4.890
  • CiteScore value: 4.49 CiteScore 4.49
  • SNIP value: 1.539 SNIP 1.539
  • SJR value: 2.404 SJR 2.404
  • IPP value: 4.28 IPP 4.28
  • h5-index value: 40 h5-index 40
  • Scimago H index value: 51 Scimago H index 51
Discussion papers | Copyright
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.

Model description paper 02 Feb 2018

Model description paper | 02 Feb 2018

Review status
This discussion paper is a preprint. It is a manuscript under review for the journal Geoscientific Model Development (GMD).

Thetis coastal ocean model: discontinuous Galerkin discretization for the three-dimensional hydrostatic equations

Tuomas Kärnä1, Stephan C. Kramer3, Lawrence Mitchell2,4, David A. Ham2, Matthew D. Piggott3, and António M. Baptista1 Tuomas Kärnä et al.
  • 1NSF Science and Technology Center for Coastal Margin Observation & Prediction, Oregon Health & Science University, Portland, Oregon, USA
  • 2Department of Mathematics, Imperial College London, London, United Kingdom
  • 3Department of Earth Science and Engineering, Imperial College London, London, United Kingdom
  • 4Department of Computing, Imperial College London, London, United Kingdom

Abstract. Unstructured grid ocean models are advantageous for simulating the coastal ocean and river-estuary-plume systems. However, unstructured grid models tend to be diffusive and/or computationally expensive which limits their applicability to real life problems. In this paper, we describe a novel discontinuous Galerkin (DG) finite element discretization for the hydrostatic equations. The formulation is fully conservative and second-order accurate in space and time. Monotonicity of the advection scheme is ensured by using a strong stability preserving time integration method and slope limiters. Compared to previous DG models advantages include a more accurate mode splitting method, revised viscosity formulation, and new second-order time integration scheme. We demonstrate that the model is capable of simulating baroclinic flows in the eddying regime with a suite of test cases. Numerical dissipation is well-controlled, being comparable or lower than in existing state-of-the-art structured grid models.

Download & links
Tuomas Kärnä et al.
Interactive discussion
Status: open (extended)
Status: open (extended)
AC: Author comment | RC: Referee comment | SC: Short comment | EC: Editor comment
[Subscribe to comment alert] Printer-friendly Version - Printer-friendly version Supplement - Supplement
Tuomas Kärnä et al.
Model code and software

Thetis coastal ocean model Thetis

A Compiler for Fast Expression Evaluation COFFEE

The Finite Element Automated Tabulator FIAT

A smarter library of finite elements FInAT

An automated finite element system Firedrake

Portable, Extensible Toolkit for Scientific Computation PETSc

The Python interface to PETSc petsc4py

Framework for performance-portable parallel computations on unstructured meshes PyOP2

The Two Stage Form Compiler TSFC

The Unified Form Language UFL

Tuomas Kärnä et al.
Total article views: 495 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
359 121 15 495 14 17
  • HTML: 359
  • PDF: 121
  • XML: 15
  • Total: 495
  • BibTeX: 14
  • EndNote: 17
Views and downloads (calculated since 02 Feb 2018)
Cumulative views and downloads (calculated since 02 Feb 2018)
Viewed (geographical distribution)
Total article views: 493 (including HTML, PDF, and XML) Thereof 488 with geography defined and 5 with unknown origin.
Country # Views %
  • 1
No saved metrics found.
No discussed metrics found.
Latest update: 18 Aug 2018
Publications Copernicus
Short summary
Unstructured meshes are attractive for coastal ocean modeling as they allow more accurate representation of complex coastal topography. Unstructured mesh models are, however, often perceived as slow and inaccurate. We present a novel Discontinuous Galerkin ocean model, Thetis. We demonstrate that the model is able to simulate baroclinic ocean flows with high accuracy on a triangular prismatic mesh. This work paves the way for highly accurate and efficient three-dimensional coastal ocean models.
Unstructured meshes are attractive for coastal ocean modeling as they allow more accurate...