DYNAMICO, an icosahedral hydrostatic dynamical core designed for consistency and versatility

Dubos, Thomas; Dubey, S.; Tort, Marine; Mittal, R. C.; Meurdesoif, Yann; Hourdin, F.

doi:10.5194/gmdd-8-1749-2015

Cited by 8 publications

(12 citation statements)

References 63 publications

(69 reference statements)

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Under the first work package, higher-level scientific concepts are used to develop a dialect for each of three icosahedral models: DYNAMICO [17], ICON [54], and NICAM [47]. A domain-specific language (DSL) is the result of finding commonalities in the three dialects.…”

Section: The Aimes Projectmentioning

confidence: 99%

AIMES: Advanced Computation and I/O Methods for Earth-System Simulations

Kunkel

Jumah

Novikova

et al. 2020

Lecture Notes in Computational Science and Engineering

View full text Add to dashboard Cite

Dealing with extreme scale earth system models is challenging from the computer science perspective, as the required computing power and storage capacity are steadily increasing. Scientists perform runs with growing resolution or aggregate results from many similar smaller-scale runs with slightly different initial conditions (the so-called ensemble runs). In the fifth Coupled Model Intercomparison Project (CMIP5), the produced datasets require more than three Petabytes of storage and the compute and storage requirements are increasing significantly for CMIP6. Climate scientists across the globe are developing next-generation models based on improved numerical formulation leading to grids that are discretized in alternative forms such as an icosahedral (geodesic) grid. The developers of these models face similar problems in scaling, maintaining and optimizing code. Performance portability and the maintainability of code are key concerns of scientists as, compared to industry projects, model code is continuously revised and extended to incorporate further levels of detail. This leads to a rapidly growing code base that is rarely refactored. However, code modernization is important to maintain productivity

show abstract

Section: The Aimes Projectmentioning

confidence: 99%

AIMES: Advanced Computation and I/O Methods for Earth-System Simulations

Kunkel

Jumah

Novikova

et al. 2020

Lecture Notes in Computational Science and Engineering

View full text Add to dashboard Cite

show abstract

“…To develop the language extensions, we worked together with the domain scientists, each of whom is an expert with one of the three icosahedral models DYNAMICO [1], ICON [2], and NICAM [3], in a co-design approach:…”

Section: Collaborative Extension Developmentmentioning

confidence: 99%

“…In GGDML, we provide support for climate model development, and focus on support for icosahedral models, which have not been widely considered in software engineering. GGDML is developed under the AIMES project to offer higher level code design for the three icosahedral models DYNAMICO [1], ICON [2], and NICAM [3].…”

Section: Introductionmentioning

confidence: 99%

GGDML: Icosahedral Models Language Extensions

Jumah¹

2017

J. Computer Sci. Techno. Updates

View full text Add to dashboard Cite

The optimization opportunities of a codebase are not completely exploited by compilers. In fact, there are optimizations that must be done within the source code. Hence, if the code developers skip some details, some performance is lost. Thus, the use of a general-purpose language to develop a performance-demanding software-e.g. climate models-needs more care from the developers. They should take into account hardware details of the target machine. Besides, writing a high-performance code for one machine will have a lower performance on another one. The developers usually write multiple optimized sections or even code versions for the different target machines. Such codes are complex and hard to maintain. In this article we introduce a higher-level code development approach, where we develop a set of extensions to the language that is used to write a model's code. Our extensions form a domainspecific language (DSL) that abstracts domain concepts and leaves the lower level details to a configurable source-to-source translation process. The purpose of the developed extensions is to support the icosahedral climate/atmospheric model development. We have started with the three icosahedral models: DYNAMICO, ICON, and NICAM. The collaboration with the scientists from the weather/climate sciences enabled agreedupon extensions. When we have suggested an extension we kept in mind that it represents a higher-level domain-based concept, and that it carries no lower-level details. The introduced DSL (GGDML-General Grid Definition and Manipulation Language) hides optimization details like memory layout. It reduces code size of a model to less than one third its original size in terms of lines of code. The development costs of a model with GGDML are therefore reduced significantly.

show abstract

“…In the spirit of Molemaker et al (2005) and Dubos et al (2015) we then deduce the discrete Poisson equation directly from the discretized kinetic energy and the discretized divergence operator. The price of this approach is to prognose either the covariant components or the contravariant ones.…”

Section: A C C E P T E D Mmentioning

confidence: 99%

Compact symmetric Poisson equation discretization for non-hydrostatic sigma coordinates ocean model

et al. 2017

View full text Add to dashboard Cite

In anticipation of relaxing the hydrostatic assumption in a sigma coordinates primitive equations ocean model, we show how a projection method can be designed with the use of a compact symmetric 15point stencil for the Poisson equation. This is achieved by recognizing that, owing the non-orthogonality of the grid, the velocity has a contravariant and covariant set of components. The two sets play a different role in the primitive equations: the contravariant components enter the definition of the model fluxes, whereas the covariant components experience the forces and in particular the pressure gradient. By treating these two sets separately, the discretized gradient and divergence operators are simple finite differences. The two sets of components are related via a linear transformation, the metric tensor, which is entirely determined by the kinetic energy. We show how the spatial discretization of the kinetic energy fully controls the Poisson equation discretization, including its boundary conditions. The discretization of the Poisson equation is shown to converge at second order and to behave as well or better than alternative methods. This approach is a prerequisite to implement later an efficient Poisson solver, such as a multigrid algorithm. Highlights► A new formulation of the non-hydrostatic primitive equations in sigma coordinates is proposed. ► Momentum and flux are two different objects related via a linear transformation. ► The Poisson equation for pressure is discretized with a compact symmetric stencil. ► The approach is a prerequisite to implement a fast Poisson solver for a forthcoming non-hydrostatic version of ROMS.

show abstract

DYNAMICO, an icosahedral hydrostatic dynamical core designed for consistency and versatility

Cited by 8 publications

References 63 publications

AIMES: Advanced Computation and I/O Methods for Earth-System Simulations

AIMES: Advanced Computation and I/O Methods for Earth-System Simulations

GGDML: Icosahedral Models Language Extensions

Compact symmetric Poisson equation discretization for non-hydrostatic sigma coordinates ocean model

Contact Info

Product

Resources

About