Pith. sign in

REVIEW 4 major objections 6 minor 64 references

HOPE: An Automatically Differentiable High-Order Non-Oscillatory Finite-Volume Shallow-Water Dynamic Core

T0 review · 4 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A cubed-sphere shallow-water core reaches 11th-order convergence, pushing smooth-flow errors to round-off on 1° grids.

desk verdict HOPE has a clever autograd-based ghost interpolation and a convincing TPP convergence table, but the paper never verifies WENO order, so the headline combined claim remains unsupported. read the letter →

arxiv 2411.19185 v3 pith:JOEZVOJS submitted 2024-11-28 physics.flu-dyn physics.geo-ph

Lilong Zhou This is my paper · ORCID
classification physics.flu-dynphysics.geo-ph MSC 65M0876M1286A10
keywords shallowwaterequationscubedspheregridfinitevolumemethodtensorproductpolynomialWENOautomaticdifferentiationhigh-orderaccuracydynamiccore
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper constructs HOPE, a finite-volume shallow-water dynamic core on a cubed-sphere grid whose spatial reconstruction is built from tensor-product polynomials and a genuine two-dimensional WENO limiter. It claims that, with a new central-stencil ghost-cell interpolation across panel boundaries, the scheme attains $k$th-order spatial accuracy from a $k$-cell-wide stencil, removing the 7th-order ceiling that earlier cubed-sphere finite-volume approaches hit. In the steady-state geostrophic test the 11th-order version reduces 12-day errors to roughly $10^{-16}$ on a 1° grid, near double-precision round-off, and the 9th-order version keeps a Rossby-Haurwitz wavenumber-4 wave intact through day 100. The implementation is written so that gradients flow through every operation, which the paper intends as the basis for coupling a high-order dynamic core with neural-network physical parameterization.

What carries the argument

The load-bearing object is the tensor-product polynomial (TPP) reconstruction matrix: a $k\times k$ stencil is filled with $k^2$ monomial terms $x^{i-1}y^{j-1}$, so evaluating the reconstruction at quadrature points is a fixed matrix-vector product $R=C A^{-1}$ applied to cell averages. The same pattern is used for the genuine 2D WENO sub-stencils, whose optimal linear weights are found by least squares and split into positive and negative parts to stabilize the nonlinear blend. Around panel edges, the paper replaces one-sided interpolation with a coupled interpolation across the boundary: ghost-cell quadrature values solve $(I-C)v = \bar A q$, so ghost values are a linear map $G = B D^{-1} \bar A$ of the known cells. The paper computes $G$ once per mesh and order, using automatic differentiation of the fixed-point iteration rather than forming the inverse directly. Fluxes are determined by the low-Mach LMARS Riemann solver and integrated by Gaussian quadrature, and the whole spatial operator is written as convolutions and matrix products.

What would settle it

Run the same steady-state geostrophic test with a 13th- or 15th-order stencil at C90 and C180 and monitor per-panel corner errors; if the convergence rate drops below the formal order or the ghost-cell fixed-point iteration fails to converge within 10 loops for any panel-corner configuration, the arbitrary-order boundary claim is refuted. More directly, compute the spectral radius of the iteration matrix $D=I-C$ from the appendix; if it is at least 1 for any stencil width, the interpolation mapping does not exist as claimed.

Watch

Extended reading notes

Core claim

The central claim is that a single construction—tensor-product polynomial reconstruction on square stencils, a genuine 2D WENO blend with least-squares optimal weights, and a ghost-cell interpolation that couples both sides of each panel edge through a linear mapping—gives a conservative, non-oscillatory, GPU-friendly shallow-water model whose convergence order can be pushed essentially arbitrarily. The paper demonstrates 3rd, 5th, 7th, 9th, and 11th order schemes on the standard steady-state geostrophic test, with measured convergence rates matching formal order at every level; the 11th-order scheme reaches $L^1$ error around $7\times 10^{-16}$ on a 1-degree grid, and at C180 the error drops to roughly $3.3\times 10^{-19}$. The same core handles topography without added viscosity, sustains a 100-day Rossby-Haurwitz wave at 9th order in the C90 test, and shows clearly better isotropy in a cylindrical dam-break test with genuine 2D WENO than with dimension-by-dimension WENO.

Load-bearing premise

The method's accuracy depends on the assumption that the panel-boundary ghost-cell equations have one unique answer and that repeating the interpolation about ten times always reaches it; the paper shows this works in the tested cases but does not prove it must work for every stencil size or panel corner.

Editorial extensions

If this is right

  • At the same stencil width, HOPE gains one order of accuracy over earlier cubed-sphere finite-volume reconstructions that used one-sided ghost cells, so a $k$-cell stencil is enough for $k$th-order convergence.
  • Because the ghost-cell mapping is one sparse matrix multiplication computed once per mesh, ultra-high order (9th and 11th) becomes practical on the cubed sphere rather than being stopped by corner oscillations.
  • The same algorithm stack—TPP reconstruction as convolution, quadrature as matrix multiplication, no loops or branches—means the dynamic core maps directly onto GPU execution and retains gradients through every Runge-Kutta stage.
  • In smooth flow tests the error at 1° resolution is already near machine round-off, so further resolution increases are unnecessary for smooth dynamics; the practical error budget shifts to non-smooth regions and physical parameterizations.
  • WENO2D's substantially better isotropy in the dam-break test implies that genuine 2D reconstruction, despite extra cost, is the right choice for cubed-sphere cores rather than dimension-by-dimension application of 1D WENO.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: the linear boundary mapping suggests a concrete convergence-stability criterion—check the spectral radius of $I-C$ for each stencil width; if it stays below 1 for all panel corners, the arbitrary-order claim can be extended from empirical to provable.
  • Editorial inference: the automatic-differentiation design implies a natural next step the paper does not run: backpropagating forecast loss through the full model to train neural subgrid parameterizations in a coupled rather than offline way, and adjoint-free 4D-Var data assimilation.
  • Editorial inference: the isotropy advantage of genuine 2D WENO in the dam-break test suggests the reconstruction would also reduce grid imprinting in tracer transport on other non-orthogonal structured grids, though the paper does not test those cases.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. The manuscript presents HOPE, a finite-volume shallow-water dynamic core on a cubed sphere, implemented in PyTorch for GPU acceleration and automatic differentiation. The core spatial reconstruction is a tensor-product polynomial (TPP) scheme of odd orders 3–11, optionally switched to a genuine two-dimensional WENO scheme for non-oscillatory behavior, with a newly proposed central-stencil ghost-cell interpolation across panel boundaries. The paper reports four claimed properties: arbitrary high-order accuracy, essential non-oscillation, mass conservation, and scalability. Numerical experiments include the Williamson case 2 steady-state geostrophic flow, case 5 flow over an isolated mountain, case 6 Rossby–Haurwitz wave, the Galewsky perturbed jet, and a designed cylinder dam-break shock test. The central quantitative evidence is Table 5.1, which shows clean convergence rates for TPP3, TPP5, TPP7, TPP9, and TPP11 on the geostrophic flow; the WENO scheme is only exercised in the qualitative dam-break experiment. The paper also briefly discusses a PyTorch implementation with convolution and matrix-multiplication operators and reports a GPU timing comparison of two einsum/permute strategies.

Significance. If the claims are fully substantiated, HOPE would be a valuable contribution: an automatically differentiable, high-order finite-volume dynamic core is directly useful for hybrid ML–physics modeling, and the central-stencil ghost-cell interpolation could remove a known barrier to very high order on the cubed sphere. The paper's convergence table for the linear TPP schemes is strong and quantitatively convincing, and the framing around GPU-friendly, branch-free kernels is timely. However, the evidence currently supports only the high-order property for the linear TPP reconstruction; the non-oscillatory, mass-conservation, and scalability properties are asserted but not quantitatively demonstrated. The WENO2D scheme, which is advertised as essential to the combined 'high-order non-oscillatory' claim, is never subjected to a smooth-flow order test, and the dam-break comparison is qualitative. Thus the central composite claim is defensible but not yet established to the standard required by the paper's own framing.

major comments (4)
  1. [§3.2 and §5.5] The high-order convergence claim is demonstrated only for the linear TPP reconstruction in Table 5.1, not for the WENO2D scheme that is the advertised non-oscillatory component. The nonlinear weights in Eqs. (3.22)–(3.23) are normalized separately over the plus and minus branches, so the scheme does not automatically reduce to the designed linear weights in smooth regions; whether the resulting perturbation preserves the formal order is nontrivial and is not analyzed. The only WENO experiment, the dam-break case in §5.5, is discontinuous and contains no error norms or convergence rates. I request a smooth-flow convergence test (e.g., Williamson case 2) run with WENO2D active, together with a discussion or verification that the smooth-region approximation retains the designed order. Without this, the central combined claim of 'high order' plus 'non-oscillation' is not supported.
  2. [§5.5] The dam-break experiment is only qualitative: Figure 5.6 compares WENO1D and WENO2D by visual inspection of circular symmetry. The claim of 'essential non-oscillation' needs a quantitative measure, such as total variation, overshoot/undershoot bounds, or a grid-convergence study of the shock propagation, and the comparison between the 1D and 2D WENO variants would be much stronger with such diagnostics. As written, the non-oscillatory property is asserted rather than demonstrated.
  3. [§3.3.1 and Appendix Eq. (7.10)–(7.12)] The novel ghost-cell interpolation rests on the assumption that the linear system (I − C)v = Ā q in Eq. (7.10) is nonsingular and that the fixed-point iteration initialized at zero converges within 10 loops. The paper states only that 'in practical applications, we have observed that the iteration typically converges within fewer than 10 loops' and gives no proof of the spectral radius or an error bound. Since the claimed arbitrary-order accuracy and stability depend on this mapping, the manuscript should provide either a proof (for the specific TPP and quadrature construction) or a careful numerical conditioning study, including the largest stencil (TPP11/13) and panel-corner configurations. Without this, the boundary treatment is an unverified load-bearing component.
  4. [§1, §6, and §5] Mass conservation is listed as one of the four advantageous properties, but no numerical diagnostic of global mass conservation is presented in any of the test cases. Although the finite-volume formulation is conservative in exact arithmetic, the ghost-cell interpolation, the panel-boundary field conversions in §3.3.2, and the metric source terms can introduce conservation errors in practice. I request a global mass error time series for at least the steady-state geostrophic flow, the Rossby–Haurwitz wave, and the dam-break case, with a clear statement of whether conservation holds to machine precision or only to truncation error. Without such evidence, one of the four headline properties is unverified.
minor comments (6)
  1. [§3, Eq. (3.1)] Equation (3.1) appears to contain a typo: the two flux terms for F are both written as F_{i+1/2,j}; the second should be F_{i−1/2,j}. This should be corrected for clarity.
  2. [§3.2] The sentence beginning 'we introduce three two types of reconstruction methods' contains a redundant phrase; please revise to 'two types'.
  3. [§5.3] The Rossby–Haurwitz discussion compares HOPE's 100-day performance with TRiSK (25 days), Li et al. (89 days), and MCORE (90 days), but the text does not state which HOPE configuration (order, resolution) achieved the 100-day result in the caption or text consistently; please clarify whether it is the 9th-order C90 simulation from Figure 5.3.
  4. [References] Reference [17] is attributed to 'Battaglia' but appears to be a duplicate of the GraphCast paper by Lam et al. (currently reference [16]). Please verify and merge/remove the duplicate citation.
  5. [§6] The sentence 'thereby inheriting its automatic differentiation capability' is repeated in the same paragraph; please remove the duplicate statement.
  6. [§5.1] The statement that 7th/9th/11th-order errors 'have surpassed the limits expressible by double-precision numbers' should be reconciled with the reported error values in Table 5.1 (e.g., TPP7 at C360 is 2.31E−15), since relative errors near 1e−15 are at the level of double-precision round-off. A short sentence on the use of quadruple precision in the runs would help avoid confusion.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: convergence rates are measured against an external analytical benchmark, and the sole self-citation is a non-load-bearing comparison.

full rationale

The paper's central convergence claims are verified against the external Williamson et al. (1992) analytical solution (case 2), so the reported rates in Table 5.1 are not fitted to the data being predicted. The TPP reconstruction matrix (3.15) and the ghost-cell mapping G (7.12) are derived from polynomial interpolation and Gaussian quadrature independently of the test cases; no parameter is tuned to reduce the measured error. The only self-citation, [47], is used as a comparison for Rossby-Haurwitz wave longevity and does not support the arbitrary-order or non-oscillation claims. The skeptic's point that WENO2D lacks a smooth-flow convergence test is a missing verification or correctness risk, not a circular reduction: the WENO construction (3.19)-(3.23) does not take the measured error as an input. Therefore no claim in the paper reduces by construction to its inputs.

Assumptions & free parameters 3 free parameters · 3 assumptions · 0 invented entities

The scheme introduces no new physical entities. The free parameters are standard WENO constants and implementation tolerances. The key unproven axioms are the invertibility and convergence of the reconstruction and ghost-cell linear systems, which are load-bearing for the arbitrary-order claim.

free parameters (3)
  • WENO epsilon = not specified
    Used in Eq. (3.23) to avoid division by zero in the nonlinear weights. The paper does not state its value; it is a standard WENO parameter but can affect the non-oscillation behavior.
  • theta = 3
    Constant in the negative-weight split Eq. (3.20), chosen by hand following Shi et al. [33], not fitted to test data.
  • Ghost interpolation iteration count and tolerance = 10 iterations, eps 1e-14 double / 1e-5 single
    The ghost cell mapping G is computed by fixed-point iteration; the number of iterations and tolerance are chosen empirically and affect the accuracy of the boundary treatment.
assumptions (3)
  • domain assumption The TPP cell-average reconstruction matrix A in Eq. (3.13) is invertible for all stencils up to order 11 and beyond.
    The paper solves a = A^{-1} q_bar without proving non-singularity; convergence tests suggest it holds for the tested orders, but there is no general proof for arbitrary order.
  • ad hoc to paper The linear system (I-C)v = A_bar q in Eq. (7.10) is nonsingular and the fixed-point iteration initialized at zero converges within 10 iterations.
    The entire panel-boundary treatment relies on this; the paper provides empirical evidence but no convergence proof. This is the weakest assumption.
  • standard math The metric and governing equations for the cubed sphere (Eq. 2.1-2.8) are the correct standard formulation and the C-property reconstruction of sqrt(G)phi_t eliminates topography oscillations.
    The equations are standard from prior literature; the C-property fix follows [11], but the paper does not demonstrate the oscillation suppression with a test that would expose failure.

how reviews work

0 comments
Cite this review

Pith. "Pith review of HOPE: An Automatically Differentiable High-Order Non-Oscillatory Finite-Volume Shallow-Water Dynamic Core." pith.science (2026). https://pith.science/paper/JOEZVOJS

@misc{pith2026241119185,
  author       = {Pith},
  title        = {Pith review of: HOPE: An Automatically Differentiable High-Order Non-Oscillatory Finite-Volume Shallow-Water Dynamic Core},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JOEZVOJS}},
  note         = {Machine review of arXiv:2411.19185}
}
read the original abstract

An automatically differentiable, high-order non-oscillatory finite volume shallow water dynamic core has been constructed on a cubed sphere grid. This dynamic core has four advantageous properties: high order accuracy, essential non-oscillation, mass conservation, and scalability. Besides, the code development is based on PyTorch, enabling the model to run seamlessly on both CPU and GPU, while naturally possessing the capability of automatic differentiation. We named the new dynamic core as High Order Prediction Environment (HOPE). The spatial reconstruction is based on the two-dimensional tensor product polynomial (TPP) and the genuine two-dimensional Weighted Essentially Non-Oscillatory scheme. A novel panel boundary approach ensures that the accuracy can reach arbitrary order. These algorithms have very high degree of compatibility with GPU architecture, allowing the computational overhead to be mitigated through the utilization of GPU. The Low Mach number Approximate Riemann Solver scheme is adopted as Riemann solvers to determine fluxes on the Gaussian points on edges. Flux across the interface between each cell edge is computed using Gaussian quadrature, and the tendencies of prognostic variables are obtained by integration all the source terms and the fluxes across the cell boundaries. This shallow water dynamic core exhibits outstanding performance in ideal shallow water test cases. In the steady-state geostrophic flow, the 11th order scheme reduces errors to nearly double precision round-off error even on coarse grids. Furthermore, HOPE maintains the Rossby-Haurwitz wave over 100 days without collapse. To test the non-oscillation property, we designed a cylinder dam break case, the WENO approach effectively suppresses non-physical oscillation, and the genuine two-dimensional reconstruction exhibits significantly better isotropy than the dimension-by-dimension scheme.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

64 extracted references · 61 canonical work pages

  1. [1]

    Tsinghua University, Beijing, China

    work page

  2. [2]

    Abstract An automatically differentiable, high-order non -oscillatory finite volume shallow water dynamic core has been constructed on a cubed sphere grid

    CMA Earth System Modeling and Prediction Centre (CEMC), Beijing, China Corresponding author: Zhou Lilong (zhoull@cma.gov.cn) Key words: Automatic differentiation; High-order accuracy; Non-Oscillation; Finite-volume methods; Shallow-water Equations; Dynamic core. Abstract An automatically differentiable, high-order non -oscillatory finite volume shallow wa...

    work page

  3. [3]

    Introduction Numerical weather prediction (NWP) bec ame a foundation of weather forecasting in past decades. The horizontal resolution of operational models in each NWP center is now finer than 10km, European Centre for Medium-Range Weather Forecasts (ECMWF) enabled real -time global weather forecasting with about 9km horizontal grid since 2016 [24], and ...

    work page 2022

  4. [4]

    We develop HOPE, a shallow water model has four advantageous properties: arbitrary high order accuracy, essential non-oscillation, mass conservation, and scalability

    work page

  5. [5]

    We desire a novel high order ghost cell interpolation scheme for cubed sphere grid, it needs only one sparse matrix multiplication to reach arbitrary odd convergence order

    work page

  6. [6]

    We implement genuine two -dimensional reconstruction on cubed sphere, comparing to the dimension -by- dimension scheme, the genuine two-dimensional provides less dimension split error

    work page

  7. [7]

    HOPE is developed on PyTorch, the auto-differentiate capability is naturally obtained, it’s easy to couple with any neural-network (NN) based functions, such as sub-grid physical parameterization

    work page

  8. [8]

    HOPE demonstrates excellent computational efficiency on GPU platforms

    The algorithms of the HOPE model primarily involve convolution and matrix multiplications, which are widely used in the artificial neural networks and highly compatible with GPUs. HOPE demonstrates excellent computational efficiency on GPU platforms. In the following part of the introduction, we introduce the relevant work on constructing the HOPE model, ...

    work page 2010

Show all 64 references
  1. [9]

    C -property

    Governing Equation on Cubed Sphere Cubed sphere grid decomposes sphere to six panels, the computational space is structured and rectangular in each panel, these features make it easy to take high order reconstruction and massive threads parallel, details in Figure 2.1. The ear...

    1972

  2. [10]

    Genuine 2D

    Numerical Discretization The finite volume method evaluates the temporal tendency of cell average by net flux, the flux across cell edges is able to be obtained by gaussian quadrature, we calculate the field value on gaussian quadrature point by spatial reconstruction and dete...

    2002

  3. [11]

    torch.tensor.permute

    High Performance Implementation and Automatic Differentiation The spatial operator and temporal integration of HOPE can be easily implemented using PyTorch. Specifically, the spatial reconstruction given by Eq.(3.14) is analogous to a convolution operation, while the Gaussian ...

  4. [12]

    (1992)[41]

    Numerical Experiments The standard test cases for spherical shallow water model are mentioned by Williamson et al. (1992)[41]. In this article, we test HOPE dynamic core using case number 2, 5, 6 desired in [41], and the case of perturbed jet flow mentioned in [8]. Besides, we...

    1992

  5. [13]

    When the flow field s are sufficiently smooth, the high -order accuracy characteristics of HOPE enable rapid convergence of simulation errors

    Conclusions In this article, we present HOPE, an innovative high -order finite volume model that boasts the capability of achieving arbitrary odd-order convergence. When the flow field s are sufficiently smooth, the high -order accuracy characteristics of HOPE enable rapid con...

  6. [14]

    Appendix In this appendix, we introduce a novel boundary ghost cell interpolation scheme for cubed sphere , which is able to support HOPE to reach the accuracy over 11th order or even higher. There are two types of cells, in-domain and out-domain (also named ghost cell, as sho...

  7. [15]

    Acker F., Borges B. de R. Borges, Costa, B. An improved WENO-Z scheme. Journal of Computational Physics, 313, 726-753, 2016

    2016

  8. [16]

    D., Tufo, H

    Bao, L., Nair, R. D., Tufo, H. M. A Mass and Momentum Flux-Form High-Order Discontinuous Galerkin Shallow Water Model on the Cubed-Sphere. Journal of Computational Physics, 271, 224-243

  9. [17]

    Accurate medium-range global weather forecasting with 3D neural networks

    Bi, K., Xie, L., Zhang, H., Chen, X., Gu, X., Tian, Q. Accurate medium-range global weather forecasting with 3D neural networks. Nature, 619, 533-538, 2023

    2023

  10. [18]

    An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws

    Borges R., Carmona, M., Costa, B., Don, W. An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws. Journal of Computational Physics, 227, 3191-3211, 2008

    2008

  11. [19]

    P., Jablonowski, C., Lin, S.-J

    Chen, X., Andronova, N., V an L., Bram, P., Joyce E., Boyd, J. P., Jablonowski, C., Lin, S.-J.. A Control-V olume Model of the Compressible Euler Equations with a Vertical Lagrangian Coordinate. Monthly Weather Review, 141, 2526-2544, 2013

    2013

  12. [20]

    FengWu: Pushing the Skillful Global Medium-range Weather Forecast beyond 10 Days Lead

    Chen, Kang, Han, Tao, Chen, Xi, Ma, Leiming, Gong, Junchao, Zhang, Tianning, Yang, Xiaokang, Bai, Lei, Ling, Fenghua, Su, Rui, Ci, Yuanzheng, Ouyang, Wanli. FengWu: Pushing the Skillful Global Medium-range Weather Forecast beyond 10 Days Lead. arXiv:2304.02948v1 [cs.AI] 6 Apr 2023

    2023

  13. [21]

    Shallow water model on cubed-sphere by multi-moment finite volume method

    Chen, C., Xiao, F.. Shallow water model on cubed-sphere by multi-moment finite volume method. Journal of Computational Physics, 227, 5019-5044, 2008

    2008

  14. [22]

    An Initial-Value Problem for Testing Numerical Models of the Global Shallow-Water Equations

    Galewsky, Joseph, Scott, Richard K., Polvani, Lorenzo M. An Initial-Value Problem for Testing Numerical Models of the Global Shallow-Water Equations. Tellus, 56A, 429-440, 2004

    2004

  15. [23]

    Weighted Essentially Non-oscillatory Schemes on Triangular Meshes

    Hu, C., Shu, C.-W. Weighted Essentially Non-oscillatory Schemes on Triangular Meshes. Journal of Computational Physics, 150, 97-127, 1999

    1999

  16. [24]

    Main Technical Improvements of GRAPES- Meso V4.0 and Verification (in Chinese)

    Huang, L., Chen, D., Deng, L., Xu, Z., Yu, F., Jiang, Y., Zhou, F.. Main Technical Improvements of GRAPES- Meso V4.0 and Verification (in Chinese). Journal of Applied Meteorological Science, 28, 25-37, 2017

    2017

  17. [25]

    and Xiao F

    Ii, S. and Xiao F. A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid. Journal of Computational Physics, 229, 1774-1796, 2010

    2010

  18. [26]

    Efficient Implementation of Weighted ENO Schemes

    Jiang, G.-S., Shu C.-W. Efficient Implementation of Weighted ENO Schemes. Journal of Computational Physics, 126, 202-228, 1996

    1996

  19. [27]

    K., Nair, R

    Katta, K. K., Nair, R. D. and Kumar, V.. High-order finite volume shallow water model on the cubed-sphere: 1D reconstruction scheme. Applied Mathematics and Computation, 266, 316-327, 2015

    2015

  20. [28]

    K., Nair, R

    Katta, K. K., Nair, R. D. and Kumar, V.. High-Order Finite-V olume Transport on the Cubed Sphere: Comparison between 1D and 2D Reconstruction Schemes. Monthly Weather Review, 143, 2937-2954, 2015

    2015

  21. [29]

    Neural General Circulation Models for Weather and Climate

    Kochkov, Dmitrii, Yuval, Janni, Langmore, Ian, Norgaard, Peter, Smith, Jamie, Mooers, Griffin, Lottes, James, Rasp, Stephan, Duben, Peter, Klower, Milan, Hatfield, Sam, Battaglia, Peter, Sanchez-Gonzalez, Alvaro, Willson, Matthew, Brenner, Michael P ., Hoyer, Stephan. Neural G...

  22. [30]

    Learning Skillful Medium-Range Global Weather Forecasting

    Lam, Remi, Sanchez-Gonzalez, Alvaro, Willson, Matthew, Wirnsberger, Peter, Fortunato, Meire, Alet, Ferran, Ravuri, Suman, Ewalds, Timo, Eaton-Rosen, Zach, Hu, Weihua, Merose, Alexander, Hoyer, Stephan, Holland, George, Vinyals, Oriol, Stott, Jacklynn, Pritzel, Alexander, Moham...

    2023

  23. [31]

    GraphCast: Learning Skillful Medium-Range Global Weather Forecasting

    Battaglia, Peter. GraphCast: Learning Skillful Medium-Range Global Weather Forecasting. arXiv:2212.12794v1 [cs.LG] 24 Dec 2022

    2022 arXiv

  24. [32]

    Analysis of and Solution to the Polar Numerical Noise Within the Shallow‐ Water Model on the Latitude‐Longitude Grid

    Li, Jianghao, Wang, Bin, Dong, Li. Analysis of and Solution to the Polar Numerical Noise Within the Shallow‐ Water Model on the Latitude‐Longitude Grid. Journal of Advances in Modeling Earth Systems, 12, 8, 2020

    2020

  25. [33]

    A Sequel to AUSM: AUSM+

    Liou, M.-S. A Sequel to AUSM: AUSM+. Journal of Computational Physics, 129, 364–382, 1996

    1996

  26. [34]

    A sequel to AUSM, Part II: AUSM+-up for all speeds

    Liou, M.-S. A sequel to AUSM, Part II: AUSM+-up for all speeds. Journal of Computational Physics, 214, 137- 170, 2006

    2006

  27. [35]

    and Chan T

    Liu, X.-D., Osher, S. and Chan T. Weighted Essentially Non-oscillatory Schemes. Journal of Computational Physics, 115, 200-212, 1994

    1994

  28. [36]

    Combination of WENO and Explicit Runge– Kutta Methods for Wind Transport in the Meso-NH Model

    Lunet, T., Lac, C., Auguste, F., Visentin, F., Masson, V ., Escobar, J. Combination of WENO and Explicit Runge– Kutta Methods for Wind Transport in the Meso-NH Model. Monthly Weather Review, 145, 3817-3838, 2017

    2017

  29. [37]

    and Wu, S.-P

    Luo, X. and Wu, S.-P. An improved WENO-Z+ scheme for solving hyperbolic conservation laws, Journal of Computational Physics, 445, 110608, 2021

    2021

  30. [38]

    A new grid for the IFS, Meteorology section of ECMWF Newsletter, 146, 23–28, 2016

    Malardel, S., Wedi, N., Deconinck, W., Diamantakis, M., Kühnlein, C., Mozdzynski, G., Hamrud, M., Smolarkiewicz, P . A new grid for the IFS, Meteorology section of ECMWF Newsletter, 146, 23–28, 2016

    2016

  31. [39]

    Norman, M. R. A high‐order WENO‐limited finite‐volume algorithm for atmospheric flow using the ADER‐ differential transform time discretization. Quarterly Journal of the Royal Meteorological Society, 147, 1661-1690, 2021

    2021

  32. [40]

    D., Thomas, S

    Nair, R. D., Thomas, S. J., Loft, R. D. A Discontinuous Galerkin Transport Scheme on the Cubed Sphere. Monthly Weather Review, 133, 814-828, 2005

    2005

  33. [41]

    D., Thomas, S

    Nair, R. D., Thomas, S. J., Loft, R. D. A Discontinuous Galerkin Global Shallow Water Model. Monthly Weather Review, 133, 876-887, 2005

    2005

  34. [42]

    D., Choi, H

    Nair, R. D., Choi, H. W., Tufo, H. M.. Computational aspects of a scalable high-order discontinuous Galerkin atmospheric dynamical core. Computers & Fluids, 38, 309-319, 2009

    2009

  35. [43]

    D., Bao, L., Toy, M

    Nair, R. D., Bao, L., Toy, M. D., Robert, K. A High-Order Multiscale Global Atmospheric Model. 8th AIAA Atmospheric and Space Environments Conference, 1-13, 2016

    2016

  36. [44]

    A., Brenowitz, N

    Perkins, W. A., Brenowitz, N. D., Bretherton, C. S.. Nugent, J. M. Emulation of Cloud Microphysics in a Climate Model. Journal of Advances in Modeling Earth Systems, 16(4), 2024

    2024

  37. [45]

    G., Kaul, C

    Pressel, K. G., Kaul, C. M., Schneider, T., Tan, Z., Mishra, S. Large-eddy simulation in an anelastic framework with closed water and entropy balances. Journal of Advances in Modeling Earth Systems, 7, 1425-1456, 2015

    2015

  38. [46]

    Conservative Finite-Difference Approximations of the Primitive Equations on Quasi-Uniform Spherical Grids

    Sadourny, R. Conservative Finite-Difference Approximations of the Primitive Equations on Quasi-Uniform Spherical Grids. Monthly Weather Review, 100, 136-144, 1972

    1972

  39. [47]

    A Technique of Treating Negative Weights in WENO Schemes

    Shi, J., Hu, C., Shu, C.-W.. A Technique of Treating Negative Weights in WENO Schemes. Journal of Computational Physics, 175 , 108-127, 2002

    2002

  40. [48]

    M., Wilhelmson, R

    Straka, J. M., Wilhelmson, R. B., Wicker, L. J., Anderson, J. R., Droegemeier, K. K.. Numerical Solutions of A Non-Linear Density Current: A Benchmark Solution and Comparisons. International Journal for Numerical Methods In Engineering, 17, 1-22, 1993

    1993

  41. [49]

    NumericalWeather Prediction (in Chinese), 2nd Edition., China Meteorological Press, Beijing, 487 pp, 2015

    Shen, T., Tian, Y ., Shen, D., Shen, X., Sun, X., Liang, X., Lu, W. NumericalWeather Prediction (in Chinese), 2nd Edition., China Meteorological Press, Beijing, 487 pp, 2015

    2015

  42. [50]

    Numerical Simulations of Rossby-Haurwitz Waves

    Thuburn, J, Li, Yong. Numerical Simulations of Rossby-Haurwitz Waves. Tellus A, 52, 181-189, 2000

    2000

  43. [51]

    A., Jablonowski, C., van Leer, B

    Ullrich, P. A., Jablonowski, C., van Leer, B.. High-Order Finite-V olume Methods for the Shallow-Water Equations on the Sphere. Journal of Computational Physics, 229, 6104-6134, 2010

    2010

  44. [52]

    A., Jablonowski, C

    Ullrich, P. A., Jablonowski, C.. MCore: A Non-hydrostatic Atmospheric Dynamical Core Utilizing High-order Finite-volume Methods. Journal of Computational Physics, 231, 5078-5108, 2012

    2012

  45. [53]

    Stable Climate Simulations Using a Realistic General Circulation Model with Neural Network Parameterizations for Atmospheric Moist Physics and Radiation Processes

    Wang, Xin, Han, Yilun, Xue, Wei, Yang, Guangwen, Zhang, Guang J. Stable Climate Simulations Using a Realistic General Circulation Model with Neural Network Parameterizations for Atmospheric Moist Physics and Radiation Processes. Geoscientific Model Development, 15, 3923-3940, 2022

    2022

  46. [54]

    J., Skamarock, W

    Wicker, L. J., Skamarock, W. C.. Time-Splitting Methods for Elastic Models Using Forward Time Schemes. Monthly Weather Review, 130, 2088-2097, 2002

    2002

  47. [55]

    L., Drake, J

    Williamson, D. L., Drake, J. B., Hack, J. J., Jakob, R., and Swarztrauber, P. N.: A Standard Test Set for Numerical Approximations of the Shallow Water Equations in Spherical Geometry, J. Comput. Phys., 102, 211–224, 1992

    1992

  48. [56]

    Adjusting Wind Speed Prediction of Numerical Weather Forecast Model Based on Maching Learning Methods

    Sun, Q., Jiao, R., Xia, J., Yan, Z., Li, H., Sun, J., Wang, L., Liang, Z.. Adjusting Wind Speed Prediction of Numerical Weather Forecast Model Based on Maching Learning Methods. Meteorological Monthly, 45(3), 426-

  49. [57]

    FengWu-4DV ar: Coupling the Data-driven Weather Forecasting Model with 4D Variational Assimilation

    Xiao, Yi, Bai, Lei, Xue, Wei, Chen, Kang, Han, Tao, Ouyang, Wanli. FengWu-4DV ar: Coupling the Data-driven Weather Forecasting Model with 4D Variational Assimilation. arXiv:2312.12455 [physics.ao-ph], 2023

    2023 arXiv

  50. [58]

    Mathematical and Physical Bases of Numerical Weather Prediction, V olume 1 (in Chinese), Beijing: Science Press, 1979

    Zeng, Q.. Mathematical and Physical Bases of Numerical Weather Prediction, V olume 1 (in Chinese), Beijing: Science Press, 1979

    1979

  51. [59]

    I., Wang, J

    Zhang, Y ., Long, M., Chen, K., Xing, L., Jin, R., Jordan, M. I., Wang, J. Skilful Nowcasting of Extreme Precipitation with NowcastNet. 619, 536-532, 2023

    2023

  52. [60]

    Weighted essentially non-oscillatory scheme on unstructured quadrilateral and triangular meshes for hyperbolic conservation laws

    Zhao, F., Pan, L., Wang S.. Weighted essentially non-oscillatory scheme on unstructured quadrilateral and triangular meshes for hyperbolic conservation laws. Journal of Computational Physics, 374, 605-624, 2018

    2018

  53. [61]

    Extending Square Conservation to Arbitrarily Structured C-grids with Shallow Water Equations

    Zhou, Lilong, Feng, Jinming, Hua, Lijuan, Zhong, Linhao. Extending Square Conservation to Arbitrarily Structured C-grids with Shallow Water Equations. Geoscientific Model Development, 13, 581-595, 2020

    2020

  54. [62]

    A New Type of Multi-resolution WENO Schemes with Increasingly Higher Order of Accuracy

    Zhu, J., Shu, C.-W.. A New Type of Multi-resolution WENO Schemes with Increasingly Higher Order of Accuracy. Journal of Computational Physics, 375, 659–683, 2018

    2018

  55. [63]

    A New Type of Multi-resolution WENO Schemes with Increasingly Higher Order of Accuracy on Triangular Meshes

    Zhu, J., Shu, C.-W.. A New Type of Multi-resolution WENO Schemes with Increasingly Higher Order of Accuracy on Triangular Meshes. Journal of Computational Physics, 392, 19-33, 2019

    2019

  56. [64]

    New Finite V olume Weighted Essentially Nonoscillatory Schemes on Triangular Meshs

    Zhu, J., Qiu, J.. New Finite V olume Weighted Essentially Nonoscillatory Schemes on Triangular Meshs. SIAM Journal on Scientific Computing, 40, A903-A928, 2018

    2018

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.