Skip to content

Research at St Andrews

Assessing the numerical accuracy of complex spherical shallow-water flows

Research output: Contribution to journalArticle



Ali R. Mohebalhojeh, David G. Dritschel

School/Research organisations


The representation of nonlinear shallow-water flows poses severe challenges for numerical modeling. The use of contour advection with contour surgery for potential vorticity (PV) within the contour-advective semi-Lagrangian (CASL) algorithm makes it possible to handle near-discontinuous distributions of PV with an accuracy beyond what is accessible to conventional algorithms used in numerical weather and climate prediction. The emergence of complex distributions of the materially conserved quantity PV, in the absence of forcing and dissipation, results from large-scale shearing and deformation and is a common feature of high Reynolds number flows in the atmosphere and oceans away from boundary layers. The near-discontinuous PV in CASL sets a limit on the actual numerical accuracy of the Eulerian, grid-based part of CASL. For the spherical shallow-water equations, the limit is studied by comparing the accuracy of CASL algorithms with second-order-centered, fourth-order-compact, and sixth-order-super compact finite differencing in latitude in conjunction with a spectral treatment in longitude. The comparison is carried out on an unstable midlatitude jet at order one Rossby number and low Froude number that evolves into complex vortical structures with sharp gradients of PV. Quantitative measures of global conservation of energy and angular momentum, and of imbalance as diagnosed using PV inversion by means of Bolin-Charney balance, indicate that fourth-order differencing attains the highest numerical accuracy achievable for such nonlinear, advectively dominated flows.



Original languageEnglish
Pages (from-to)3876-3894
Number of pages19
JournalMonthly Weather Review
Issue number11
Publication statusPublished - Nov 2007

    Research areas


Discover related content
Find related publications, people, projects and more using interactive charts.

View graph of relations

Related by author

  1. Comparison of the Moist Parcel-In-Cell (MPIC) model with large-eddy simulation for an idealized cloud

    Böing, S. J., Dritschel, D. G., Parker, D. J. & Blyth, A. M., 29 Apr 2019, In : Quarterly Journal of the Royal Meteorological Society. In press, 17 p.

    Research output: Contribution to journalArticle

  2. On the regularity of the Green-Naghdi equations for a rotating shallow fluid layer

    Dritschel, D. G. & Jalali, M. R., 25 Apr 2019, In : Journal of Fluid Mechanics. 865, p. 100-136

    Research output: Contribution to journalArticle

  3. Scale-invariant singularity of the surface quasigeostrophic patch

    Scott, R. K. & Dritschel, D. G., 25 Mar 2019, In : Journal of Fluid Mechanics. 863, 12 p., R2.

    Research output: Contribution to journalArticle

Related by journal

  1. Toward a PV-based algorithm for the dynamical core of hydrostatic global models

    Mohebalhojeh, A. R., Joghataei, M. & Dritschel, D. G., Jun 2016, In : Monthly Weather Review. 144, 6, p. 2481-2502 22 p.

    Research output: Contribution to journalArticle

  2. On the Spectral Convergence of the Supercompact Finite-Difference Schemes for the f-Plane Shallow-Water Equations

    Ghader, S., Mohebalhojeh, A. R. & Esfahanian, V., Jul 2009, In : Monthly Weather Review. 137, 7, p. 2393-2406 14 p.

    Research output: Contribution to journalArticle

  3. The Diabatic Contour-Advective Semi-Lagrangian Algorithms for the Spherical Shallow Water Equations

    Mohebalhojeh, A. R. & Dritschel, D. G., 29 Sep 2009, In : Monthly Weather Review. 137, 9, p. 2979-2994 16 p.

    Research output: Contribution to journalArticle

  4. The diabatic contour advective semi-Lagrangian algorithm

    Dritschel, D. G. & Ambaum, M. HP., Sep 2006, In : Monthly Weather Review. 134, 9, p. 2503-2514 12 p.

    Research output: Contribution to journalArticle

ID: 627848