pith. sign in

arxiv: 2604.07103 · v1 · submitted 2026-04-08 · 🧮 math.NA · cs.NA· physics.flu-dyn

A new high-order finite-volume advection scheme on spherical Voronoi grids and a comparative study in a mimetic finite-volume moist shallow-water model

Luan F. Santos , Jeferson B. Granjeiro , Pedro S. Peixoto This is my paper

Pith reviewed 2026-05-10 17:28 UTC · model grok-4.3

classification 🧮 math.NA cs.NAphysics.flu-dyn
keywords high-order advectionspherical Voronoi gridsk-exact reconstructionfinite-volume methodmoist shallow-water modeltracer transportnumerical weather predictiongrid robustness
0
0 comments X

The pith

k-exact reconstruction extends high-order advection schemes to irregular spherical Voronoi grids while preserving accuracy.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops new finite-volume advection schemes on spherical centroidal Voronoi tessellations by extending the k-exact reconstruction technique from flat domains to the sphere. It demonstrates through standard advection tests and a mimetic moist shallow-water model that the schemes reach their target high-order accuracy, remain largely insensitive to grid irregularity, and perform comparably to methods already used in models such as MPAS. The work matters because these grids permit local refinement for regional detail without changing the overall discretization, and higher-order advection reduces artificial smoothing of fine tracer features such as moisture. If the results hold, the limiting factor for robust variable-resolution atmospheric modeling becomes the shallow-water equations discretization rather than the choice of advection scheme.

Core claim

The proposed schemes based on k-exact reconstruction achieve high-order accuracy in classical spherical advection tests, exhibit little sensitivity to grid distortion on locally refined SCVTs, and produce results comparable to existing advection schemes when transporting moisture in a mimetic finite-volume moist shallow-water model. As a result, overall grid robustness with respect to irregularity is governed by the sensitivity of the shallow-water model discretization itself, independent of the advection scheme.

What carries the argument

The k-exact reconstruction approach, which builds high-degree polynomial approximations from cell averages on each spherical Voronoi cell to evaluate accurate interface fluxes for finite-volume advection.

If this is right

  • High-order accuracy holds on irregular and locally refined spherical grids, supporting flexible resolution in atmospheric models.
  • Reduced numerical diffusion improves preservation of fine-scale moisture structures during transport.
  • Advection-scheme choice does not determine the overall robustness to grid irregularity in the moist shallow-water model.
  • The new schemes yield performance similar to current MPAS methods, making them interchangeable for variable-resolution forecasts.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The approach may enable adaptive mesh refinement in global spherical models without requiring separate advection code for refined patches.
  • Extending the same reconstruction to full three-dimensional nonhydrostatic dynamical cores could be tested next to check behavior under realistic vertical motions.
  • Because the advection step itself is no longer the bottleneck, future gains in model robustness would come from improving the mimetic finite-volume treatment of the shallow-water equations.

Load-bearing premise

That the k-exact reconstruction approach can be extended to spherical Voronoi tessellations while preserving high-order accuracy and stability on highly irregular and locally refined grids.

What would settle it

An advection test on a highly distorted SCVT grid in which the measured convergence rate falls below the design order or the scheme produces instability.

Figures

Figures reproduced from arXiv: 2604.07103 by Jeferson B. Granjeiro, Luan F. Santos, Pedro S. Peixoto.

Figure 1
Figure 1. Figure 1: Illustration of SCVT grids used in this work. (a) Spherical grid composed of [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Illustration of the stencils using first-level neighbors (a) and first- and second￾level neighbors (b), used to construct the reconstruction polynomials. The reconstruction polynomial is then written as ϕ R i (x, y) = c0 + c1x + c2y + c3x 2 + c4xy + c5y 2 , (14) with coefficients {cm} determined by enforcing ϕ R i (xk, yk) = ϕk , k ∈ Si , (15) where (xk, yk) are the projected coordinates of the control vol… view at source ↗
Figure 3
Figure 3. Figure 3: Initial condition for the zonal wind test case: a Gaussian hill centered at [PITH_FULL_IMAGE:figures/full_fig_p020_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: L∞ errors (a) and convergence rates (b) for the advection of a Gaussian hill under a zonal wind, computed with OG (solid) and SG (dashed) schemes on quasi-uniform SCVT grids (levels 2-7) without flux limiting. Line colors indicate the nominal order of accuracy. The Gaussian hill is initially centered at latitude and longitude (0, 0) (Figure 3a) [PITH_FULL_IMAGE:figures/full_fig_p021_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: L2 errors (a) and convergence rates (b) for the advection of a Gaussian hill under a zonal wind, computed with OG (solid) and SG (dashed) schemes on quasi-uniform SCVT grids (levels 2-7) with flux limiting. Line colors indicate the nominal order of accuracy. The Gaussian hill is initially centered at latitude and longitude (0, 0) (Figure 3a). The L2 errors and convergence rates on the quasi-uniform SCVT gr… view at source ↗
Figure 6
Figure 6. Figure 6: L2 error convergence without (a) and with flux limiting (b) for the advection of a Gaussian hill under a zonal wind, computed with OG (solid) and SG (dashed) schemes on variable-resolution SCVT grids (levels 2-7). Line colors indicate the nominal order of accuracy. The Gaussian hill is initially centered at latitude and longitude 7π 18 , − π 12  , located over the Andes (Figure 3b). Next, we repeat the te… view at source ↗
Figure 7
Figure 7. Figure 7: Time evolution of the L∞ error for the advection of a Gaussian hill under a zonal wind over 5 time units, without (a) and with (b) flux limiting, computed using OG (solid) and SG (dashed) schemes on level 7 variable-resolution SCVT grid. Line colors indicate the nominal order of accuracy. The Gaussian hill is initially centered at latitude and longitude 7π 18 , − π 12  , located over the Andes (Figure 3b)… view at source ↗
Figure 8
Figure 8. Figure 8: Initial condition for the deformational flow test case: two Gaussian hills on (a) [PITH_FULL_IMAGE:figures/full_fig_p024_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: L∞ errors (a) and convergence rates (b) for the advection of two Gaussian hills under a deformational flow velocity field, computed with OG (solid) and SG (dashed) schemes on quasi-uniform SCVT grids (levels 2-7) without flux limiting. Line colors indi￾cate the nominal order of accuracy. The initial condition is shown in Figure 8a [PITH_FULL_IMAGE:figures/full_fig_p025_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: L2 errors (a) and convergence rates (b) for the advection of two Gaussian hills under a deformational flow velocity field, computed with OG (solid) and SG (dashed) schemes on variable-resolution SCVT grids (levels 2-7) with flux limiting. Line colors indicate the nominal order of accuracy. The initial condition is shown in Figure 8a. this behavior is not observed here [PITH_FULL_IMAGE:figures/full_fig_p0… view at source ↗
Figure 11
Figure 11. Figure 11: L∞ error convergence without (a) and with flux limiting (b) for the advection of a Gaussian hill under a deformational flow velocity field, computed with OG (solid) and SG (dashed) schemes on variable-resolution SCVT grids (levels 2-7). Line colors indicate the nominal order of accuracy. The initial condition is shown in Figure 8b. We again use variable-resolution grids to assess the sensitivity of each s… view at source ↗
Figure 12
Figure 12. Figure 12: Slotted cylinder advection under the deformational flow velocity field on a [PITH_FULL_IMAGE:figures/full_fig_p027_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Slotted cylinder advection under the deformational flow velocity field on a [PITH_FULL_IMAGE:figures/full_fig_p028_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Slotted cylinder advection under the deformational flow velocity field on a [PITH_FULL_IMAGE:figures/full_fig_p029_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: L2 error convergence for the moist shallow-water model in the moist steady geostrophic test case from Zerroukat and Allen [37], considering 12 days of simulation. Results are shown for simulations with flux limiting applied to the fluid depth field, using the OG (solid lines) and SG (dashed lines) schemes on quasi-uniform SCVT grids (levels 2-7). Line colors indicate the nominal order of accuracy. From [… view at source ↗
Figure 16
Figure 16. Figure 16: Spatial distribution of the errors for the cloud water field at the uniform SCVT [PITH_FULL_IMAGE:figures/full_fig_p033_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: L2 error convergence for the moist shallow-water model in the moist steady geostrophic test case from Zerroukat and Allen [37], considering 12 days of simulation. Results are shown for simulations with flux limiting applied to the fluid depth field, using the OG (solid lines) and SG (dashed lines) schemes on variable-resolution SCVT grids (levels 2-7). Line colors indicate the nominal order of accuracy. t… view at source ↗
Figure 18
Figure 18. Figure 18: Cloud field in the Southern Hemisphere at day 7 for the moist shallow-water [PITH_FULL_IMAGE:figures/full_fig_p035_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: As in Fig. 18, but for the Andes topography–based refined SCVT grid. [PITH_FULL_IMAGE:figures/full_fig_p036_19.png] view at source ↗
read the original abstract

Spherical centroidal Voronoi tessellations (SCVTs), currently used in numerical weather forecasting models such as the Model for Prediction Across Scales (MPAS), are a type of spherical grid that is highly flexible, allowing the construction of locally refined regions with higher resolution without requiring modifications to the numerical discretization or its implementation. However, the irregularity of SCVT grids makes the construction of robust high-order schemes challenging. In particular, in atmospheric modeling, high-order advection schemes are desirable since they reduce numerical diffusion and improve the representation of fine-scale tracer structures. Therefore, in this work, we propose a new class of high-order advection schemes on the sphere based on the $k$-exact reconstruction approach, extending their successful use on planar domains to the spherical surface. We assess the performance of the proposed method and compare it with existing advection schemes for SCVT grids used in MPAS. The evaluation includes classical advection test cases on the sphere as well as simulations with a mimetic finite-volume moist shallow-water model, in which the advection scheme is applied to the transport of moisture tracers. Grid-related robustness was investigated using locally refined spherical grids with a local focus on the Andes topography. Our results show that the proposed schemes achieve high-order accuracy in the advection tests, exhibit little sensitivity to grid distortion, and produce comparable results to existing schemes in the moist shallow-water model. Overall, grid robustness is therefore limited to the sensitivity of the discretization of the shallow-water model, irrespective of the advection scheme.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper claims to introduce a new class of high-order advection schemes on spherical centroidal Voronoi tessellations (SCVTs) by extending the k-exact reconstruction approach from planar domains. These schemes are tested on classical spherical advection cases and embedded in a mimetic finite-volume moist shallow-water model for moisture transport. Using locally refined SCVT grids with an Andes topography focus, the authors report that the schemes attain high-order accuracy, exhibit little sensitivity to grid distortion, and produce results comparable to existing MPAS schemes. The final claim is that overall grid robustness is governed by the shallow-water discretization rather than the advection scheme choice.

Significance. If the numerical results hold under scrutiny, the work is significant for numerical weather prediction and geophysical fluid dynamics. SCVT grids enable flexible local refinement in models like MPAS without altering the core discretization, but their irregularity has hindered high-order advection. Demonstrating a robust k-exact extension with low distortion sensitivity and practical performance in moist tests could reduce numerical diffusion for fine-scale tracers. The comparative study and topography-focused robustness tests provide actionable insights for model developers.

major comments (3)
  1. [Abstract] Abstract: The assertion that the schemes 'achieve high-order accuracy in the advection tests' is not accompanied by any quantitative convergence rates, error norms (L1/L2/L∞), or observed orders versus resolution. This quantitative support is load-bearing for the central claim and must appear in the results section with explicit tables or figures.
  2. [Method description] Method (reconstruction procedure): No truncation-error analysis or proof is supplied showing that the polynomial reproduction property of k-exact reconstruction is preserved under the spherical metric and non-planar Voronoi cell geometry. If the moment-fitting step implicitly assumes local planarity, formal high-order accuracy may degrade on highly irregular or refined SCVTs; the observed 'little sensitivity to grid distortion' then risks being test-specific rather than general.
  3. [Numerical results] Moist shallow-water experiments: The statement that 'grid robustness is therefore limited to the sensitivity of the discretization of the shallow-water model, irrespective of the advection scheme' requires explicit quantification (e.g., norm differences between schemes versus model variations) to substantiate independence; comparable results alone do not establish the claim.
minor comments (2)
  1. [Method] Notation for the spherical reconstruction basis and weight computation should be defined more explicitly, perhaps with an equation for the local polynomial fit on the sphere.
  2. [Introduction] Add references to prior k-exact implementations on unstructured planar grids to contextualize the spherical extension.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thorough and constructive review of our manuscript. We have carefully addressed each major comment below and revised the manuscript to incorporate the suggested improvements where possible.

read point-by-point responses

  1. Referee: The assertion that the schemes 'achieve high-order accuracy in the advection tests' is not accompanied by any quantitative convergence rates, error norms (L1/L2/L∞), or observed orders versus resolution. This quantitative support is load-bearing for the central claim and must appear in the results section with explicit tables or figures.

    Authors: We agree that explicit quantitative support strengthens the central claim. Although convergence behavior was demonstrated via plots in the original results section, we have added a dedicated table in the revised results section that reports L1, L2, and L∞ error norms together with the observed orders of accuracy across successive resolutions for the advection test cases. The abstract has also been updated to reference these observed orders. revision: yes

  2. Referee: No truncation-error analysis or proof is supplied showing that the polynomial reproduction property of k-exact reconstruction is preserved under the spherical metric and non-planar Voronoi cell geometry. If the moment-fitting step implicitly assumes local planarity, formal high-order accuracy may degrade on highly irregular or refined SCVTs; the observed 'little sensitivity to grid distortion' then risks being test-specific rather than general.

    Authors: The referee correctly identifies the absence of a formal truncation-error analysis. The reconstruction is performed in local tangent planes with geodesic distances used to compute moments, thereby incorporating the spherical metric. We have added a new paragraph in the methods section that derives the leading-order error terms under the local-planarity approximation and shows that the polynomial reproduction property is retained to the design order. While a complete a priori analysis for arbitrary SCVTs lies outside the present scope, the added discussion and the consistent high-order convergence observed on both uniform and highly distorted refined grids provide evidence that the behavior is not test-specific. revision: partial

  3. Referee: The statement that 'grid robustness is therefore limited to the sensitivity of the discretization of the shallow-water model, irrespective of the advection scheme' requires explicit quantification (e.g., norm differences between schemes versus model variations) to substantiate independence; comparable results alone do not establish the claim.

    Authors: We accept that the original claim rested on qualitative comparison. In the revised manuscript we have inserted a quantitative comparison: a table of L2-norm differences in the moisture fields obtained when swapping advection schemes versus the differences obtained when varying the shallow-water discretization parameters (time-step size and explicit diffusion coefficient). The tabulated values confirm that scheme-to-scheme differences remain smaller than those induced by changes to the underlying shallow-water discretization, thereby supporting the stated conclusion. revision: yes

Circularity Check

0 steps flagged

No significant circularity; claims rest on independent numerical validation

full rationale

The paper extends the established k-exact reconstruction technique from planar domains to spherical Voronoi grids and demonstrates high-order accuracy plus low distortion sensitivity solely through numerical advection tests and moist shallow-water simulations on independent benchmark cases. No load-bearing step reduces by construction to a fitted parameter, self-citation chain, or renamed input; the central results are externally falsifiable via the reported test outcomes and comparisons with existing MPAS schemes.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The claim relies on the mathematical extension of k-exact reconstruction to spherical surfaces and empirical validation through numerical experiments on standard test cases and the moist shallow-water model.

free parameters (1)
  • reconstruction order k
    The order of the polynomial reconstruction is a user-chosen parameter determining the accuracy level of the scheme.
axioms (1)
  • domain assumption Spherical geometry and Voronoi cell metrics can be incorporated into the k-exact reconstruction without loss of formal order of accuracy
    This is the key extension assumed when moving from planar to spherical domains.

pith-pipeline@v0.9.0 · 5599 in / 1468 out tokens · 106149 ms · 2026-05-10T17:28:28.687600+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

61 extracted references · 61 canonical work pages

  1. [1]

    Staniforth, J

    A. Staniforth, J. Thuburn, Horizontal grids for global weather and climate prediction models: A review, Q. J. Roy. Meteor. Soc. 138 (2012) 1 – 26

    work page 2012

  2. [2]

    D. L. Williamson, The Evolution of Dynamical Cores for Global Atmo- spheric Models, Journal of the Meteorological Society of Japan. Ser. II 85B (2007) 241–269

    work page 2007

  3. [3]

    Mouallem, L

    J. Mouallem, L. Harris, X. Chen, Implementation of the novel duo-grid in GFDL’s fv3 dynamical core, Journal of Advances in Modeling Earth Systems 15 (2023) e2023MS003712

    work page 2023

  4. [4]

    P. S. Peixoto, S. R. Barros, Analysis of grid imprinting on geodesic sphericalicosahedralgrids, JournalofComputationalPhysics237(2013) 61–78

    work page 2013

  5. [5]

    Weller, J

    H. Weller, J. Thuburn, C. J. Cotter, Computational modes and grid imprinting on five quasi-uniform spherical C grids, Monthly Weather Review 140 (2012) 2734 – 2755. 41

    work page 2012

  6. [6]

    P. S. Peixoto, Accuracy analysis of mimetic finite volume operators on geodesic grids and a consistent alternative, Journal of Computational Physics 310 (2016) 127–160

    work page 2016

  7. [7]

    Tomita, M

    H. Tomita, M. Satoh, K. Goto, An optimization of the icosahedral grid modified by spring dynamics, Journal of Computational Physics 183 (2002) 307–331

    work page 2002

  8. [8]

    Miura, M

    H. Miura, M. Kimoto, A comparison of grid quality of optimized spheri- cal hexagonal–pentagonal geodesic grids, Monthly Weather Review 133 (2005) 2817 – 2833

    work page 2005

  9. [9]

    R. P. Heikes, D. A. Randall, C. S. Konor, Optimized icosahedral grids: Performance of finite-difference operators and multigrid solver, Monthly Weather Review 141 (2013) 4450 – 4469

    work page 2013

  10. [10]

    Q. Du, V. Faber, M. Gunzburger, Centroidal Voronoi tessellations: Applications and algorithms, SIAM Review 41 (1999) 637–676

    work page 1999

  11. [11]

    Q. Du, M. D. Gunzburger, L. Ju, Constrained centroidal Voronoi tes- sellations for surfaces, SIAM Journal on Scientific Computing 24 (2003) 1488–1506

    work page 2003

  12. [12]

    D. Engwirda, JIGSAW-GEO (1.0): locally orthogonal staggered un- structured grid generation for general circulation modelling on the sphere, Geoscientific Model Development 10 (2017) 2117–2140

    work page 2017

  13. [13]

    W. C. Skamarock, J. B. Klemp, M. G. Duda, L. D. Fowler, S.-H. Park, T. D. Ringler, A multiscale nonhydrostatic atmospheric model using centroidal Voronoi tesselations and C-grid staggering, Monthly Weather Review 140 (2012) 3090 – 3105

    work page 2012

  14. [14]

    Ringler, M

    T. Ringler, M. Petersen, R. L. Higdon, D. Jacobsen, P. W. Jones, M. Maltrud, A multi-resolution approach to global ocean modeling, Ocean Modelling 69 (2013) 211–232

    work page 2013

  15. [15]

    Thuburn, T

    J. Thuburn, T. Ringler, W. Skamarock, J. Klemp, Numerical represen- tation of geostrophic modes on arbitrarily structured C-grids, Journal of Computational Physics 228 (2009) 8321–8335. 42

    work page 2009

  16. [16]

    Ringler, J

    T. Ringler, J. Thuburn, J. Klemp, W. Skamarock, A unified approach to energy conservation and potential vorticity dynamics for arbitrarily- structured C-grids, Journal of Computational Physics 229 (2010) 3065– 3090

    work page 2010

  17. [17]

    S. Freitas, Model for Ocean-Land-Atmosphere Prediction (MONAN): Avaliação dos Candidatos ao Núcleo Dinâmico do Componente At- mosférico (MONAN-ATM) - Recomendação ao Comitê Científico do MONAN, Technical Report, Instituto Nacional de Pesquisas Espaci- ais (INPE), São José dos Campos, 2023. URL:http://urlib.net/ 8JMKD3MGP3W34T/49MKML8

    work page 2023

  18. [18]

    N. K.-R. Kevlahan, T. Dubos, Wavetrisk-1.0: an adaptive wavelet hy- drostatic dynamical core, Geoscientific Model Development 12 (2019) 4901–4921

    work page 2019

  19. [19]

    Dubos, S

    T. Dubos, S. Dubey, M. Tort, R. Mittal, Y. Meurdesoif, F. Hour- din, DYNAMICO-1.0, an icosahedral hydrostatic dynamical core de- signed for consistency and versatility, Geoscientific Model Development 8 (2015) 3131–3150

    work page 2015

  20. [20]

    L. F. Santos, P. S. Peixoto, Topography-based local spherical Voronoi grid refinement on classical and moist shallow-water finite-volume mod- els, Geoscientific Model Development 14 (2021) 6919–6944

    work page 2021

  21. [21]

    W. C. Skamarock, A. Gassmann, Conservative transport schemes for spherical geodesic grids: High-order flux operators for ode-based time integration, Monthly Weather Review 139 (2011) 2962–2975

    work page 2011

  22. [22]

    Miura, An upwind-biased conservative advection scheme for spherical hexagonal–pentagonal grids, Monthly Weather Review 135 (2007) 4038– 4044

    H. Miura, An upwind-biased conservative advection scheme for spherical hexagonal–pentagonal grids, Monthly Weather Review 135 (2007) 4038– 4044

    work page 2007

  23. [23]

    W. C. Skamarock, M. Menchaca, Conservative transport schemes for spherical geodesic grids: High-order reconstructions for forward-in-time schemes, Monthly Weather Review 138 (2010) 4497–4508

    work page 2010

  24. [24]

    Miura, W

    H. Miura, W. C. Skamarock, An upwind-biased transport scheme us- ing a quadratic reconstruction on spherical icosahedral grids, Monthly Weather Review 141 (2013) 832–847. 43

    work page 2013

  25. [25]

    Y. Zhang, Extending high-order flux operators on spherical icosahedral grids and their applications in the framework of a shallow water model, Journal of Advances in Modeling Earth Systems 10 (2018) 145–164

    work page 2018

  26. [26]

    C. J. Subich, Higher-order finite volume differential operators with se- lective upwinding on the icosahedral spherical grid, Journal of Compu- tational Physics 368 (2018) 21–46

    work page 2018

  27. [27]

    Ollivier-Gooch, High-order ENO schemes for unstructured meshes basedonleast-squaresreconstruction, in: 35thAerospaceSciencesMeet- ing and Exhibit, 1997, p

    C. Ollivier-Gooch, High-order ENO schemes for unstructured meshes basedonleast-squaresreconstruction, in: 35thAerospaceSciencesMeet- ing and Exhibit, 1997, p. 540

    work page 1997

  28. [28]

    Ollivier-Gooch, M

    C. Ollivier-Gooch, M. Van Altena, A high-order-accurate unstructured mesh finite-volume scheme for the advection–diffusion equation, Journal of Computational Physics 181 (2002) 729–752

    work page 2002

  29. [29]

    Ollivier-Gooch, A

    C. Ollivier-Gooch, A. Nejat, K. Michalak, Obtaining and verifying high- order unstructured finite volume solutions to the Euler equations, AIAA journal 47 (2009) 2105

    work page 2009

  30. [30]

    Michalak, C

    C. Michalak, C. Ollivier-Gooch, Accuracy preserving limiter for the high-order accurate solution of the euler equations, Journal of Compu- tational Physics 228 (2009) 8693–8711

    work page 2009

  31. [31]

    Nejat, C

    A. Nejat, C. Ollivier-Gooch, A high-order accurate unstructured finite volume Newton–Krylov algorithm for inviscid compressible flows, Jour- nal of Computational Physics 227 (2008) 2582–2609

    work page 2008

  32. [32]

    Jalali, M

    A. Jalali, M. Sharbatdar, C. Ollivier-Gooch, Accuracy analysis of un- structured finite volume discretization schemes for diffusive fluxes, Com- puters & Fluids 101 (2014) 220–232

    work page 2014

  33. [33]

    Hoshyari, C

    S. Hoshyari, C. F. Ollivier-Gooch, A higher-order unstructured finite volume solver for three-dimensional compressible flows, in: 2018 AIAA Aerospace Sciences Meeting, AIAA, 2018, pp. 2018–1306

    work page 2018

  34. [34]

    Barth, P

    T. Barth, P. Frederickson, Higher order solution of the euler equations on unstructured grids using quadratic reconstruction, in: 28th AIAA Aerospace Sciences Meeting, 1990. AIAA 1990-13. 44

    work page 1990

  35. [35]

    Barth, Recent developments in high-order k-exact reconstruction on unstructured meshes, in: 31st AIAA Aerospace Sciences Meeting, 1993

    T. Barth, Recent developments in high-order k-exact reconstruction on unstructured meshes, in: 31st AIAA Aerospace Sciences Meeting, 1993. AIAA 1993-668

    work page 1993

  36. [36]

    R. D. Nair, P. H. Lauritzen, A class of deformational flow test cases for linear transport problems on the sphere, Journal of Computational Physics 229 (2010) 8868–8887

    work page 2010

  37. [37]

    Zerroukat, T

    M. Zerroukat, T. Allen, A moist Boussinesq shallow water equations set for testing atmospheric models, Journal of Computational Physics 290 (2015) 55–72

    work page 2015

  38. [38]

    L. Ju, T. Ringler, M. Gunzburger, Voronoi Tessellations and Their Ap- plication to Climate and Global Modeling, Springer Berlin Heidelberg, Berlin, Heidelberg, 2011, pp. 313–342

    work page 2011

  39. [39]

    S. T. Zalesak, Fully multidimensional flux-corrected transport algo- rithms for fluids, Journal of computational physics 31 (1979) 335–362

    work page 1979

  40. [40]

    Colella, P

    P. Colella, P. R. Woodward, The Piecewise Parabolic Method (PPM) for gas-dynamical simulations, Journal of Computational Physics 54 (1984) 174–201

    work page 1984

  41. [41]

    W. M. Putman, S.-J. Lin, Finite-volume transport on various cubed- sphere grids, Journal of Computational Physics 227 (2007) 55–78

    work page 2007

  42. [42]

    L. F. Santos, J. Mouallem, P. S. Peixoto, Analysis of finite-volume transport schemes on cubed-sphere grids and an accurate scheme for divergent winds, Journal of Computational Physics 522 (2025) 113618

    work page 2025

  43. [43]

    T.M.Bendall, J.Kent, SWIFT:Amonotonic, flux-formsemi-lagrangian tracer transport scheme for flow with large courant numbers, Monthly Weather Review 153 (2025) 565 – 587

    work page 2025

  44. [44]

    L. J. Wicker, W. C. Skamarock, Time-splitting methods for elastic mod- els using forward time schemes, Monthly Weather Review 130 (2002) 2088 – 2097

    work page 2002

  45. [45]

    Godunov, A difference method for numerical calculation of discon- tinuous solutions of the equations of hydrodynamics, Mat

    S. Godunov, A difference method for numerical calculation of discon- tinuous solutions of the equations of hydrodynamics, Mat. Sb. 47(89):3 (1959) 271–306. 45

    work page 1959

  46. [46]

    Van Leer, Towards the ultimate conservative difference scheme

    B. Van Leer, Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection, Journal of Computational Physics 23 (1977) 276–299

    work page 1977

  47. [47]

    P. H. Lauritzen, P. A. Ullrich, R. D. Nair, Atmospheric Transport Schemes: Desirable Properties and a Semi-Lagrangian View on Finite- Volume Discretizations, Springer Berlin Heidelberg, Berlin, Heidelberg, 2011, pp. 185–250

    work page 2011

  48. [48]

    Gross, H

    M. Gross, H. Wan, P. J. Rasch, P. M. Caldwell, D. L. Williamson, D. Klocke, C. Jablonowski, D. R. Thatcher, N. Wood, M. Cullen, B. Beare, M. Willett, F. Lemarié, E. Blayo, S. Malardel, P. Termonia, A. Gassmann, P. H. Lauritzen, H. Johansen, C. M. Zarzycki, K. Sak- aguchi, R. Leung, Physics–dynamics coupling in weather, climate, and earth system models: Ch...

    work page 2018

  49. [49]

    J. O. Ferguson, C. Jablonowski, H. Johansen, Assessing adaptive mesh refinement (AMR) in a forced shallow-water model with moisture, Monthly Weather Review 147 (2019) 3673 – 3692

    work page 2019

  50. [50]

    Hartney, T

    N. Hartney, T. M. Bendall, J. Shipton, Exploring forms of the moist shallow-water equations using a new compatible finite-element discreti- sation, Quarterly Journal of the Royal Meteorological Society 152 (2026) e70018

    work page 2026

  51. [51]

    P. S. Peixoto, S. R. Barros, On vector field reconstructions for semi- lagrangian transport methods on geodesic staggered grids, Journal of Computational Physics 273 (2014) 185–211

    work page 2014

  52. [52]

    Arakawa, V

    A. Arakawa, V. R. Lamb, A potential enstrophy and energy conserving scheme for the shallow water equations, Monthly Weather Review 109 (1981) 18 – 36

    work page 1981

  53. [53]

    Thuburn, C

    J. Thuburn, C. J. Cotter, T. Dubos, A mimetic, semi-implicit, forward- in-time, finite volume shallow water model: comparison of hexago- nal–icosahedral and cubed-sphere grids, Geoscientific Model Develop- ment 7 (2014) 909–929

    work page 2014

  54. [54]

    Eldred, D

    C. Eldred, D. Randall, Total energy and potential enstrophy conserving schemes for the shallow water equations using hamiltonian methods – 46 part 1: Derivation and properties, Geoscientific Model Development 10 (2017) 791–810

    work page 2017

  55. [55]

    An interpretation of trisk-type schemes from a discrete exterior calculus perspective

    C. Eldred, W. Bauer, An interpretation of trisk-type schemes from a discrete exterior calculus perspective, arXiv preprint arXiv:2210.07476 (2022)

    work page arXiv 2022

  56. [56]

    Weller, Controlling the computational modes of the arbitrarily struc- tured C grid, Monthly Weather Review 140 (2012) 3220 – 3234

    H. Weller, Controlling the computational modes of the arbitrarily struc- tured C grid, Monthly Weather Review 140 (2012) 3220 – 3234

    work page 2012

  57. [57]

    Weller, Non-orthogonal version of the arbitrary polygonal C-grid and a new diamond grid, Geoscientific Model Development 7 (2014) 779–797

    H. Weller, Non-orthogonal version of the arbitrary polygonal C-grid and a new diamond grid, Geoscientific Model Development 7 (2014) 779–797

    work page 2014

  58. [58]

    D. L. Williamson, J. B. Drake, J. J. Hack, R. Jakob, P. N. Swarztrauber, A Standard Test Set for Numerical Approximations to the Shallow Wa- ter Equations in Spherical Geometry, J. Comput. Phys. 102 (1992) 211–224

    work page 1992

  59. [59]

    Galewsky, R

    J. Galewsky, R. K. Scott, L. M. Polvani, An initial-value problem for testing numerical models of the global shallow-water equations, Tellus A: Dynamic Meteorology and Oceanography 56 (2004) 429–440

    work page 2004

  60. [60]

    Preprint

    A.Gassmann, Third-order momentumadvection onthe quasi-hexagonal C-grid on the sphere, ESS Open Archive (2022). Preprint

    work page 2022

  61. [61]

    H. Wang, W. C. Skamarock, G. Feingold, Evaluation of scalar advection schemes in the advanced research WRF model using large-eddy simula- tions of aerosol–cloud interactions, Monthly Weather Review 137 (2009) 2547 – 2558. 47

    work page 2009