libmpdata++ 0.1: a library of parallel MPDATA solvers for systems of generalised transport equations A. Jaruga, S. Arabas, D. Jarecka, H. Pawlowska, P. K. Smolarkiewicz, and M. Waruszewski Institute of Geophysics, Faculty of Physics, University of Warsaw, Warsaw, Poland National Center for Atmospheric Research, Boulder, Colorado, USA European Centre for Medium-Range Weather Forecasts, Reading, ...

Multidimensional Positive Definite Advection Transport Algorithm (MPDATA) and Monotonic Upstream Schemes for Conservation Laws (MUSCL) are compared and evaluated in a biogeochemical-like tracer model. Based on the MUSCL scheme, an adapted scheme for r-coordinate models (MUSCL_VD) is presented. These schemes are first implemented in idealized experiments configured for a closed channel with and ...

Abstract Advection scheme is one of the core challenges in computational fluid dynamics, which restricts capacity model performance many research areas including oceanographic modelling studies. Here this study, we compared newly developed algorithm HSIMT (advection with High-order Spatial Interpolation at Middle Temporal level) well-known MPDATA (Multidimensional Positive-Definite Advective Tr...

To better accommodate the highly disparate length scales encountered in geophysical flows, we have extended EULAG with a solution-adaptive moving mesh capability [1]. The development builds on [2], where the authors set forth a time-dependent curvilinear coordinate formulation of the governing PDEs to enable dynamic mesh adaptivity in EULAG. Here, the anelastic flux-form dynamical core of EULAG...

Load balancing is a widely accepted technique for performance optimization of scientific applications on parallel architectures. Indeed, balanced applications do not waste processor cycles on waiting at points of synchronization and data exchange, maximizing this way the utilization of processors. In this paper, we challenge the universality of the load-balancing approach to optimization of the...

We present a new finite difference code for modeling three-dimensional thermal convection in a spherical shell using the “cubed sphere” method of Ronchi et al. [1]. The equation of motion is solved using a poloidal potential formulation of the equation of motion for an iso-viscous, infinite Prandtl number fluid on a finite difference grid and advective transport is implemented using the 2-order...

Twenty years ago, the multi-dimensional, positive definite, advection transport algorithm was introduced by Smolarkiewicz. Over the two decades since, it has been applied countless times to numerous problems, however almost always on rectilinear grids. One of the few exceptions is the Operational Multiscale Environment model with Grid Adaptivity (OMEGA), an atmospheric simulation system origina...