The ICON-1.2 hydrostatic atmospheric dynamical core on triangular grids – Part 1: Formulation and performance of the baseline version
Open Access
- 5 June 2013
- journal article
- Published by Copernicus GmbH in Geoscientific Model Development
- Vol. 6 (3), 735-763
- https://doi.org/10.5194/gmd-6-735-2013
Abstract
As part of a broader effort to develop next-generation models for numerical weather prediction and climate applications, a hydrostatic atmospheric dynamical core is developed as an intermediate step to evaluate a finite-difference discretization of the primitive equations on spherical icosahedral grids. Based on the need for mass-conserving discretizations for multi-resolution modelling as well as scalability and efficiency on massively parallel computing architectures, the dynamical core is built on triangular C-grids using relatively small discretization stencils. This paper presents the formulation and performance of the baseline version of the new dynamical core, focusing on properties of the numerical solutions in the setting of globally uniform resolution. Theoretical analysis reveals that the discrete divergence operator defined on a single triangular cell using the Gauss theorem is only first-order accurate, and introduces grid-scale noise to the discrete model. The noise can be suppressed by fourth-order hyper-diffusion of the horizontal wind field using a time-step and grid-size-dependent diffusion coefficient, at the expense of stronger damping than in the reference spectral model. A series of idealized tests of different complexity are performed. In the deterministic baroclinic wave test, solutions from the new dynamical core show the expected sensitivity to horizontal resolution, and converge to the reference solution at R2B6 (35 km grid spacing). In a dry climate test, the dynamical core correctly reproduces key features of the meridional heat and momentum transport by baroclinic eddies. In the aqua-planet simulations at 140 km resolution, the new model is able to reproduce the same equatorial wave propagation characteristics as in the reference spectral model, including the sensitivity of such characteristics to the meridional sea surface temperature profile. These results suggest that the triangular-C discretization provides a reasonable basis for further development. The main issues that need to be addressed are the grid-scale noise from the divergence operator which requires strong damping, and a phase error of the baroclinic wave at medium and low resolutions.Keywords
This publication has 65 references indexed in Scilit:
- AMIP Simulation with the CAM4 Spectral Element Dynamical CoreJournal of Climate, 2013
- Earth System Modelling - Volume 2Published by Springer Science and Business Media LLC ,2012
- Kernel-based vector field reconstruction in computational fluid dynamic modelsInternational Journal for Numerical Methods in Fluids, 2011
- Inspection of hexagonal and triangular C-grid discretizations of the shallow water equationsJournal of Computational Physics, 2011
- On utility of triangular C-grid type discretization for numerical modeling of large-scale ocean flowsOcean Dynamics, 2010
- Analysis of Discrete Shallow-Water Models on Geodesic Delaunay Grids with C-Type StaggeringMonthly Weather Review, 2005
- Reducing the effect of global communication in GMRES(m) and CG on parallel distributed memory computersApplied Numerical Mathematics, 1995
- The design and application of upwind schemes on unstructured meshesPublished by American Institute of Aeronautics and Astronautics (AIAA) ,1989
- The Piecewise Parabolic Method (PPM) for gas-dynamical simulationsJournal of Computational Physics, 1984
- Frequency Filter for Time IntegrationsMonthly Weather Review, 1972