Solution of Partial Differential Equations on the Sphere (PDEs on the sphere)
Monday 24th September 2012 to Friday 28th September 2012
09:00 to 10:20  Registration  
09:50 to 10:20  Morning Coffee  
10:20 to 10:30  Welcome from the Deputy Director  
10:30 to 10:55 
John Thuburn (University of Exeter) Mimetic Semiimplicit Solution of the Shallow Water Equations on HexagonalIcosahedral and CubedSphere Grids Session: Horizontal Session
A new algorithm is presented for the solution of the shallow water equations on quasiuniform spherical grids. It combines a mimetic spatial discretization with a CrankNicolson time scheme for fast waves and an accurate and conservative forwardintime advection scheme for mass and potential vorticity. The algorithm is tested on two families of grids: hexagonalicosahedral Voronoi grids, and modified equiangular cubesphere grids. For the cubedsphere case, a key ingredient is the development of a suitable discrete Hodge star operator for the nonorthogonal grid.
Results of several test cases will be presented. The results confirm a number of desirable properties for which the scheme was designed: exact mass conservation, very good available energy and potential enstrophy conservation, vanishing curl of grad, steady geostrophic modes, and accurate PV advection. The scheme is stable for large wave Courant numbers and for advective Courant numbers up to about one.
The accuracy of the scheme appears to be limited by the accuracy of the various mimetic spatial operators. On the hexagonal grid there is no evidence for damaging effects of computational Rossby modes, despite attempts to force them explicitly.

INI 1  
10:55 to 11:20 
Thomas Dubos (École Polytechnique) Variational derivation of energyconserving finitedifference schemes for geophysical fluid equations Session: Horizontal Session
At the continuous level, the socalled vectorinvariant form of the equations of fluid motion can be obtained from Hamilton's principle of least action, using the EulerPoincaré formalism. In this form of the equations, the mass flux and Bernoulli function appear as functional derivatives of the Hamiltonian. The integral conservation of energy and the Lagrangian conservation of potential vorticity then follow straightforwardly.
This setting can be imitated at the discrete level to yield energyconserving schemes. Key ingredients include an energyconserving vector product, a discrete rule of integration by parts, a discrete approximation of total energy, and the definition of the discrete mass flux and Bernoulli function from partial derivatives of the discrete energy. For the rotating shallowwater equations, it turns out that schemes by Sadourny (1975) on Cartesian meshes and by Bonaventura & Ringler (2005) on DelaunayVoronoi meshes fit in the above framework.
Furthermore new schemes can be obtained. For instance it is possible to modify the discrete mass flux to match a modification of the discrete energy, as Renner (1981) and Skamarock et al. (2012) did to suppress a numerical instability. Finally, a socalled discrete Hodge star operator is obtained on a nonorthogonal pair of primal and dual meshes. This operator is part of a potentialvorticity conserving scheme of the shallowwater equations on meshes for which an orthogonal dual does not exist, like the equiangular cubed sphere. Potential application to sets of threedimensional equations will be discussed.

INI 1  
11:20 to 11:45 
Colin Cotter (Imperial College London) Finite element exterior calculus framework for geophysical fluid dynamics Session: Horizontal Session
We present a finite element framework that extends the mimetic approach of Ringler, Thuburn, Skamarock and Klemp (2010) to finite element methods. The framework preserves the following benefits of the mimetic approach: local mass conservation, steady geostrophic linear modes on the fplane, consistent conservative potential vorticity advection, but allows for selection of higherorder finite element discretisations and adjustment of the balance of vorticity and mass degrees of freedom. We will explain the foundations of the framework and present some preliminary numerical results which form part of the UK Met Office/NERC/STFC GungHo dynamical core project.

INI 1  
11:45 to 12:10 
Stefan Vater (Freie Universität Berlin) A multilevel time integrator for computing longwave shallow water flows at low Froude numbers Session: Horizontal Session
A new multilevel semiimplicit scheme for the computation of low Froude number shallow water flows is presented. Motivated by the needs of atmospheric flow applications, it aims to minimize dispersion and amplitude errors in the computation of long wave gravity waves. While it correctly balances "slaved" dynamics of shortwave solution components induced by slow forcing, the method eliminates freely propagating compressible shortwave modes, which are under resolved in time. This is achieved through the multilevel approach borrowing ideas from multigrid schemes for elliptic equations. The scheme is second order accurate and admits time steps depending only on the flow velocity. It incorporates a predictor step using a Godunovtype method for hyperbolic conservation laws and two elliptic corrections.
The multilevel method is initially derived for the onedimensional linearized shallow water equations. Scalewise decomposition of the data enables a scale dependent blending of time integrators with different principal features. To guide the selection of these integrators, the discretedispersion relations of some standard secondorder schemes are analyzed, and their response to high wave number low frequency source terms is discussed. The resulting method essentially consists of the solution of a Helmholtz problem on the original fine grid, where the operator and the right hand side incorporate the multiscale information of the discretization. The performance of the method in the linear case is illustrated on a test case with "multiscale" initial data and a problem with a slowly varying high wave number source term.
For the computation of fully nonlinear shallow water flows, a projection method for zero Froude number flows is generalized by incorporating the local time derivatives of the height. This semiimplicit method is combined with the multilevel method for the linearized equations to obtain the above mentioned properties of the scheme. Numerical tests address the scheme's ability to correctly cover the asymptotic flow regime of longwave gravity waves passing over shortrange topography and its balancing properties for the lake at rest.

INI 1  
12:30 to 13:30  Lunch at Wolfson Court  
13:30 to 14:00  Poster Session  
14:00 to 14:25 
Sergey Danilov (AlfredWegenerInstitut für Polar und Meeresforschung (AWI)) Multiscale ocean simulations with FESOM Session: Horizontal Session
FiniteElement Seaice Ocean circulation Model (FESOM) uses unstructured triangular meshes and is employed either as standingalone global ocean model or as part of FESOMECHAM5/6 coupled setup. Several examples pertaining to studies focused on polar oceans will be given to characterize the potential of FESOM and also to indicate typical difficulties related to the strongly variable resolution.
Brief intercomparison of FESOM to a global cellvertex finitevolume setup is presented. While it is more numerically efficient than the finiteelement setup of FESOM, it requires special measures to maintain stable performance because of its too large velocity space. Most effective stabilization is provided by either computing the momentum advection on scalar (mediandual) control volumes, or making use of the vectorinvariant form. Both approaches effectively trim the nonlinear term in the momentum equation to the size of scalar space.

INI 1  
14:25 to 14:50 
Jürgen Steppeler (Deutscher Wetterdienst (DWD)) Uniformly third Order conserving Schems on Polygonal Grids Session: Horizontal Session
Uniformly third Order conserving Schems on Polygonal Grids. The interest in polygonal grids is increasing. They are an alternative to the more commonly used spectral and latitude longitude grids. Among other advantages they offer the possibility of a rather uniform cover of the sphere with grid cells. Other advantages concern the ease of using multiprocessing computers and using special vertical treatments, such as shaved cells. Well known examples of polygonal grids are the cube sphere and the icosahedral grid. After initial research by Sadourny and Williamsson the practicability of this approach was shown by Baumgardner and Steppeler. In particular Baumgardner showed that problems with some approaches can be traced back to the fact that for slightly irregular resolution methods are not uniformly second order. After correcting this problem Baumgardner was able to show that problems arising from irregular grids do not occur. Steppeler generalized this approach to third order. Both Baumgardners and Steppelers approaches were non conser ving. A generalization to conserving schemes will be presented and computational examples given. Another high order approach is the pecral element method, which currently is available for orders 4 an higher only. The approach presented can be considered as a version of third order spectral elements. The advantages of third order schemes over even higher order approaches will be discussed.

INI 1  
14:50 to 15:15 
Ramachandran Nair (National Center for Atmospheric Research) A SemiLagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the CubedSphere Session: Horizontal Session
The discontinuous Galerkin (DG) method combines fine features of highorder accurate finiteelement and finitevolume methods. Because of its geometric flexibility and high parallel efficiency, DG method is becoming increasingly popular in atmospheric and ocean modeling. However, a major drawback of DG method is its stringent CFL stability restriction associated with explicit timestepping. A way to get around this issue is to combine DG method with a Lagrangian approach based on the characteristic LagrangeGalerkin philosophy. Unfortunately, a fully 2D approach combining DG and Lagrangian methods is algorithmically complex and computationally expensive for practical application, particularly for nonorthogonal curvilinear geometry such as the cubedsphere grid system. We adopt a dimensionsplitting approach where a regular semiLagrangian (SL) scheme is combined with the DG method. The resulting SLDG scheme employs a sequence of 1D operations for solving transport equation on the cubedsphere. The SLDG scheme is inherently conservative and has the option to incorporate a local positivitypreserving filter for tracers. A novel feature of the SLDG algorithm is that it can be used for multitracer transport for global models employing spectralelement (structured or unstructured) grids. The SLDG scheme is tested for various benchmark advection testsuites on the sphere and results will be presented in the seminar.

INI 1  
15:15 to 15:45  Afternoon Tea  
15:45 to 16:10 
JeanPierre Croisille (Université de Lorraine) Compact finite difference schemes on the CubedSphere Session: Horizontal Session
The CubedSphere is a spherical grid made of six quasicartesian square like patches. It was originally introduced by Sadourny some forty years ago. We extend to this grid the design of highorder finite difference compact operators. Such discrete operators are used in Computational Fluid Dynamics on structured grids for applications such as Direct Numerical Simulation of turbulent flows, or aeroacoustics problems. We consider in this work the design of a uniformly fourthorder accurate spherical gradient. The main approximation principle consists in defining a network of great circles covering the CubedSphere along which a highorder hermitian gradient can be calculated. This procedure allows a natural treatment at the interface of the six patches. The main interest of the approach is a fully symmetric approximation system on the CubedSphere. We numerically demonstrate the accuracy of the approximate gradient on several test problems, in particular the cosinebell testcase of Williamson et al. for climatology.

INI 1  
16:10 to 16:35 
Nicholas Kevlahan (McMaster University) A conservative adaptive wavelet method for the shallow water equations on staggered grids Session: Horizontal Session
This paper presents the first dynamically adaptive wavelet method for the shallow water equations on a staggered hexagonal Cgrid. Pressure is located at the centres of the primal grid (hexagons) and velocity is located at the edges of the dual grid (triangles). Distinct biorthogonal second generation wavelet transforms are developed for the pressure and the velocity. These wavelet transforms are based on secondorder accurate interpolation and restriction operators. Together with compatible restriction operators for the mass flux and Bernoulli function, they ensure that mass is conserved and that there is no numerical generation of vorticity when solving the shallow water equations. Grid refinement relies on appropriate thresholding of the wavelet coefficients, allowing error control in both the quasigeostrophic and inertiagravity wave regimes. The shallow water equations are discretized on the dynamically adapted multiscale grid using a mass and potentialenstrophy conserving finitedifference scheme. The conservation and error control properties of the method are verified by applying it to a propagating inertiagravity wave packet and to rotating shallow water turbulence. Significant savings in the number of degrees of freedom are achieved even in the case of rotating shallowwater turbulence. The numerical dissipation introduced by the grid adaptation is quantified. The method has been designed so it can be extended easily to the icosahedral subdivision of the sphere. This work provides important building blocks for the development of fully adaptive general circulation models.

INI 1  
16:35 to 17:00 
TaeJin Oh (Korea Institute of Atmospheric Prediction Systems (KIAPS)) Dynamical Core Developments at KIAPS Session: Horizontal Session
Korea Institute of Atmospheric Prediction Systems (KIAPS) is a new organization founded to develop the next generation operational numerical weather prediction (NWP) model for Korea Meteorological Administration (KMA). With the increasing demand for global highresolution simulation, scalability becomes an important issue and highly scalable numerical methods such as spectral element (continuous Galerkin, CG) or discontinuous Galerkin (DG) methods are gaining interest. Although conventional finite difference (FD) or finite volume (FV) methods offer excellent efficiency, it is difficult to formulate high order schemes on nonorthogonal grid structures which is needed to avoid grid singularities. On the other hand, CG/DG methods are not constrained to lower order on unstructured, nonorthogonal grids. We present high order convergence properties of CG/DG methods for advection, shallow water equations on structured and/or unstructured grids in one and/or two dimensions. In order t o match the high order spatial truncation error, an explicit general order mstage m–1 order strong stability preserving (SSP) RungeKutta time integrator is used. With this setup, we can obtain arbitrary order convergence rates.

INI 1  
17:00 to 18:00  Welcome Drinks Reception and Poster Session 
09:00 to 09:25 
Sarvesh Kumar Dubey (Indian Institute of Technology) Slopelimited transport schemes using icosahedral hexagonal grid Session: Advection Session
In this work two simple advection schemes for unstructured meshes are proposed and implemented over spherical icosahedralhexagonal grids. One of them is fully discrete in space and time while the other one is a semi discrete scheme with third order RungeKutta time integration. Both schemes use cellwise linear reconstruction. We therefore also present two possible candidates for consistent gradient discretization over general grids. These gradients are designed in a finite volume sense with an adequate modification to guarantee convergence in the absence of a special grid optimization. Monotonicity of the advection schemes is enforced by a slope limiter, at contrast with the widely used of posterior approach of fluxcorrected transport (FCT). Convergence of the proposed gradient reconstruction operators is verified numerically. It is found that the proposed modification is indeed necessary for convergence on nonoptimized grids. Recently proposed advection test cases are used to evaluate the performance of the slopelimited advection schemes. It is verified that they are convergent and positive. We also compare these schemes to a variant where positivity is enforced by the FCT approach. FCT produces slightly less diffusion but it seems to be at the price of some nonphysical antidiffusion. These results suggest that the proposed slope limited advection schemes are a viable option for icosahedralhexagonal grids over sphere.

INI 1  
09:25 to 09:50 
Oswald Knoth (Leibniz Institute for Tropospheric Research, Leipzig) Numerical Solution of the Advection Equation on Unstructured Spherical Grids with Logarithmic Reconstruction Session: Advection Session
There are numerous approaches for solving hyperbolic differential equations in the context of finite volume methods. One popular approach is the limiter free LocalDoubleLogarithmicReconstruction (LDLR) of Artebrant and Schroll. The aim of this work is to construct a threedimensional reconstructing function based on the LDLR for solving the advection equation on unstructured spherical grids. The new method should preserve the characteristics of the LDLR. That means in particular a reconstruction without use of limiters and with a small stencil of only the nearest neighbors of a particular cell. Also local extrema should be conserved while the local variation of the reconstruction within one cell should be under control.
We propose an ansatz which works on unstructured polyhedral grids. To come up to this, an ansatz function with one logarithmic expression for each face of the polyghedron is constructed. Required gradients at cell face midpoints are determined by use of the MultiPointFluxApproximation (MPFA) method. Further derivative information are obtained with the help of special barycentric coordinates. All necessary integrals of the ansatz functions can be computed exactly. The spatially discretized equations are combined with explicit RungeKutta methods to advance the solution in time.
The new advection procedure is numerically evaluated with standard test cases from the literature on different unstructured spherical grids.

INI 1  
09:50 to 10:15 
Janakiraman Subburathnam (Centre for Development of Advanced Computing, (CDAC)) Linear advection characteristics of a variable resolution global spectral method on the sphere Session: Advection Session
Advection experiments are conducted with a variable resolution global spectral method on the sphere. The variable resolution grid has finer resolution over the tropics and the resolution decreases smoothly as we move towards the poles [Janakiraman et al (2012) ]. An Eulerian formulation of the linear advection is used to for the spectral discretization. Near dispersionfree advection is achieved on the high resolution tropical belt. The transport across homogeneous resolution regions produce very less dispersion errors. Transport over the poles result in severe grid representation errors. It is shown that increasing the resolution over the reatly reduces this error. Transport of a feature from apoint close to poles but not over it does not produce such representation errors. A comparison of timeschemes such as 4th order RungeKutta method, Magazenkov scheme and Leapfrog scheme for the advection experiment is also presented.
Reference(s) : S. Janakiraman, Ravi S Nanjundiah and A.S. Vasudeva Murthy, A novel variable resolution global spectral method on the sphere, Vol. 231, No. 7, pp 27942810, 2012.

INI 1  
10:15 to 10:40 
Kara Peterson (Sandia National Laboratories) Optimizationbased Conservative Transport on the Sphere Session: Advection Session
We present a new optimizationbased conservative transport algorithm for scalar quantities (i.e. mass) that preserves monotonicity without the use of flux limiters. The method is formulated as a singly linearly constrained quadratic program with simple bounds where the net mass updates to the cell are the optimization variables. The objective is to minimize the discrepancy between a target or highorder mass update and a mass update that satisfies physical bounds. In this way, we separate accuracy considerations, handled by the objective functional from the enforcement of physical bounds, handled by the contraints. With this structure mass conservation is incorporated as a constraint and a simple, efficient, and easily parallelizable optimization algorithm is obtained. This algorithm has been extended to a latitudelongitude grid for twodimensional remapping and transport on the sphere. Results for several standard test problems on the sphere will be shown to illustrate the accuracy and robustness of the method.

INI 1  
10:40 to 11:10  Morning Coffee  
11:10 to 11:35 
Colin Zarzycki (University of Michigan) Improving tropical cyclone representation in general circulation models through the use of variable resolution Session: Adaptive Meshing
Modeling of tropical cyclones in General Circulation Models (GCMs) has been a historically difficult task due to issues such as relatively small storm sizes and intense convective processes. However, recent advances in GCM model design coupled with improvements in computing ability now allow for GCM simulations with grid spacings as small as 1530 km. This presentation evaluates the potential of GCMs at these high horizontal resolutions to simulate tropical cyclones. In particular, we explore a novel variableresolution mesh approach that allows for high spatial resolutions in areas of interest, such as lowlatitude ocean basins where tropical cyclones are prevalent. Such GCM designs with variableresolution meshes have the potential to become a future tool in weather forecasting as well as for regional climate assessments.
A staticallynested, variableresolution mesh option has recently been introduced into the National Center for Atmospheric Research (NCAR) Community Atmosphere Model's (CAM) Spectral Element (SE) dynamical core. We present preliminary CAMSE model simulations using an idealized tropical cyclone test case with a variety of grid sizes and refinement scales. We evaluate the evolution of tropical cyclones initialized at various locations in or near grid scale transition regions. Specific focus is centered on factors crucial to storm cyclogenesis and maintenance such as airsea interaction and vertical development. In addition to shortterm, deterministic tests, we also investigate the performance of multiresolution meshes in longerterm climate simulations, including attention paid to the dependance of nonseeded tropical cyclone genesis on spatial resolution and the subsequent implications for regional climate modeling within a global modeling framework. We also discuss pot ential computational consequences of using such a setup in either process or climate studies.

INI 1  
11:35 to 12:00 
Phillip Colella (Lawrence Berkeley National Laboratory) Adaptive Highorder Finite Volume Discretizations on Spherical Thin Shells Session: Adaptive Meshing
We present an adaptive, conservative finite volume approach applicable to solving hyperbolic PDE's on both 2D surface and 3D thin shells on the sphere. The starting point for this method is the equiangular cubedsphere mapping, which maps six rectangular coordinate patches (blocks) onto the sphere. The images of these blocks form a disjoint union covering the sphere, with the mappings of adjacent blocks being continuous, but not differentiable, at block boundaries. Our method uses a fourthorder accurate discretization to compute flux averages on faces, with a higherorder least squaresbased interpolation to compute stencil operations near block boundaries. To suppress oscillations at discontinuities and underresolved gradients, we use a limiter that preserves fourthorder accuracy at smooth extrema, and a redistribution scheme to preserve positivity where appropriate for advected scalars. By using both space and timeadaptive mesh refinement, the solver allocates comp utational effort only where greater accuracy is needed. The resulting method is demonstrated to be fourthorder accurate for advection and shallow water equation model problems, and robust at solution discontinuities. We will also present an approach for the compressible Euler equations on a 3D thin spherical shell. Refinement is performed only in the horizontal directions, The radial direction is treated implicitly (using a fourthorder RK IMEX scheme) to eliminate time step constraints from vertical acoustic waves.

INI 1  
12:00 to 12:25 
Jonathan Lambrechts (Université Catholique de Louvain) Generation of Provably Correct Curvilinear Meshes Session: Adaptive Meshing
The development of highorder numerical technologies for CFD is underway for many years now. For example, Discontinuous Galerkin methods (DGM) have been largely studied in the literature, initially in a quite theoretical context, and now in the application point of view. In many contributions, it is shown that the accuracy of the method strongly depends of the accuracy of the geometrical discretization. In other words, the following question is raised:we have the high order methods, but how do we get the meshes?
This talk focus on the generation of highly curved ocean meshes of polynomial order 2 to 4.
In the first part, we propose a robust procedure that allows to build a curvilinear mesh for which every element is guaranteed to be valid. The technique builds on standard optimization method (BICG) combined with a logbarrier objective function to guarantee the positivity of the elements Jacobian and thus the validity of the elements.
To be valid is not the only requirement for a goodquality mesh. If the temporal discretization is explicit, even a valid element can lead to a very stringent constraint on the stable time step. The second part of the talk is devoted to the optimization of the curvilinear ocean meshes to obtain large stable time steps.

INI 1  
12:30 to 13:30  Lunch at Wolfson Court  
13:30 to 14:00  Poster Session  
14:00 to 14:25 
Donna Calhoun (Boise State University) A logically Cartesian, adaptively refined twopatch sphere grid for modeling transport in the atmosphere Session: Adaptive Meshing
Recently, we demonstrated results using a second order
finite volume scheme on a novel, logically Cartesian, twopatch 2d
sphere grid. The numerical scheme, based on the wavepropagation
algorithms in Clawpack ( R. J. LeVeque, Univ. of Washington), was
easily adapted to the single grid Cartesian layout of our twopatch
sphere grid mapping. Furthermore, we were able to use an existing
adaptive mesh refinement patchbased (AMR) code (Berger, Oliger et
al.) to run computational efficient simulations.
In our current work, we are developing a new AMR code which uses wave
propagation algorithms on nonoverlapping AMR grids stored as leaves
in a forest of quad or octtrees. The underlying treebased code,
p4est (Carsten Burstedde, Univ. of Bonn) manages the multiblock
connectivity in a multiprocessor environment and has been shown to be
highly scalable in applications of interest to geophysicists. This
new AMR code, which we call ForestClaw, will easily handle the
adaptivity for our twopatch sphere grid, as well as the cubedsphere,
and more generally, any multiblock geometry.
We will present preliminary results from our efforts to use ForestClaw
for modeling volcanic ash dispersal in the atmosphere. This is joint
work with Carsten Burstedde (Univ. of Bonn) and researchers at the
Cascade Volcanic Observatory (Vancouver, Washington).

INI 1  
14:25 to 14:50 
Andreas Müller (Naval Postgraduate School) Accuracy of adaptive discontinuous Galerkin simulations Session: Adaptive Meshing
Adaptive mesh refinement generally serves to increase computational efficiency without compromising the accuracy of the numerical solution significantly. However it is an open question in which regions the spatial resolution can actually be coarsened without affecting the accuracy of the result significantly. This question is investigated for a specific example of dry atmospheric convection, namely the simulation of warm air bubbles. For this purpose a novel numerical model is developed that is tailored towards this specific meteorological problem. The compressible Euler equations are solved with a Discontinuous Galerkin method. Time integration is done with a semiimplicit approach and the dynamic grid adaptivity uses space filling curves via the AMATOS function library. The numerical model is validated with a convergence study and five standard test cases.
A method is introduced which allows one to compare the accuracy between different choices of refinement regions even in a case when the exact solution is not known. Essentially this is done by comparing features of the solution that are strongly sensitive to spatial resolution. For a rising warm air bubble the additional error by using adaptivity is smaller than 1% of the total numerical error if the average number of elements used for the adaptive simulation is about 50% smaller than the number used for the simulation with the uniform fineresolution grid. Correspondingly the adaptive simulation is almost two times faster than the uniform simulation.

INI 1  
14:50 to 15:15 
Giovanni Tumolo (Abdus Salam International Centre for Theoretical Physics) A semiimplicit, semiLagrangian, padaptive Discontinuous Galerkin method for the rotating shallow water equations Session: Adaptive Meshing
As a first step towards construction and analysis of a DG based dynamical core for high resolution atmospheric modeling, a semiimplicit and semiLagrangian Discontinuous Galerkin method for the shallow water equations with rotation is proposed and analyzed.
The method is equipped with a simple padaptivity criterion, that allows to adjust dynamically the number of local degrees of freedom employed to the local structure of the solution.
Numerical results in the framework of two dimensional test cases prove that the method captures accurately and effectively the main features of linear gravity and inertial gravity waves. Also the solution of nonlinear Stommel problem is correctly simulated. The effectiveness of the method is also demonstrated by numerical results obtained at high Courant numbers and with automatic choice of the local approximation degree.
The present research has been carried out during the PhD of the autor (supervisors F. Giorgi and L. Bonaventura) with financial support from the {\it Abdus Salam International Center for Theoretical Physics} and in collaboration with L.Bonaventura and M. Restelli {\it MOX  Politecnico di Milano}.
See MOXReport 04/2012 Tumolo, G.; Bonaventura, L.; Restelli, M. A semiimplicit, semiLagrangian, padaptive Discontinuous Galerkin method for the shallow water equations.

INI 1  
15:15 to 15:45  Afternoon Tea  
15:45 to 16:10 
Peter Bosler (University of Michigan) Particle methods for geophysical flow on the sphere Session: Adaptive Meshing
We develop a fluid dynamics solver for spherical domains based on the Lagrangian form of the equations of motion and a treestructured mesh of fluid particles. Initial discretizations are based on the cubed sphere or icosahedral triangles, but these arrangements distort as the particles move with the fluid velocity. A remeshing scheme based on interpolation of the Lagrangian parameter is applied at regular intervals to maintain computational accuracy as the flow evolves. Ongoing work with adaptive mesh refinement is introduced. We present solutions of the advection equation with prescribed velocities from recent test cases (Nair and Lauritzen, 2010), and solutions of the barotropic vorticity equation where velocity is given by an approximate BiotSavart integral and midpoint rule quadrature for test cases including RossbyHaurwitz waves and Gaussian vortices . The discrete equations in the barotropic vorticity equation are those of point vortices on the sphere (e.g. Bogo molov, 1977). For the shallow water equations, we must also include the effects of point sources upon the fluid particles; we discuss the challenges posed by divergent flows in this Lagrangian context and strategies for evaluating nonlinear forcing terms on irregular meshes.

INI 1  
16:10 to 16:35 
John Boyd (University of Michigan) Progress in Radial Basis Function Methods: Adaptive VortexRBF Methods for the Sphere and Other Advances Session: Adaptive Meshing
Vortex methods, in which flow is modelled by overlapping Lagrangian vortex blobs, and radial basis function pseudospectral methods, which have proven their worth for irregular, adaptive grids, are combined into a single algorithm to compute flows on a sphere. The strengths and limitations of vortexRBF methods for vortexdominated flows with and without rotation will be illustrated by case studies including finite amplitude Rossby waves, vortex merger and deformation, and rollup and instability. Other pertinent advances will be described as time permits.

INI 1  
16:35 to 17:00 
Michal Kopera (Naval Postgraduate School) Adaptive mesh refinement for discontinuous Galerkin method on quadrilateral nonconforming grid Session: Adaptive Meshing
In recent years there has been a significant interest in using the discontinuous Galerkin method (DG) for solving fluid dynamics problems. Studies in [1], [2], and [3] have shown that the DG method is indeed a good choice for the construction of future nonhydrostatic numerical weather prediction models. It combines highorder accuracy of the solution with geometric flexibility of unstructured grids and exhibits excellent scaling properties.
In order to increase the scale resolution capabilities of DG, as well as to take better advantage of available computing power, the use of adaptive mesh refinement (AMR) for a quadrilateral nonconforming grid is being investigated. The results of AMR implementation to a 2D version of the Nonhydrostatic Unified Model of the Atmosphere (NUMA) will be presented during the talk. Both static and dynamic hadaptivity (element grid adaptivity) will be considered. Since increased local grid resolution has to impose severe constraints on the timestep of an explicit method, implicitexplicit (IMEX) and multi rate time integration will be discussed as well.
As the latest trend in scientific computing is to employ clusters of GPUs for highperformance computations, the implementation of CUDA kernels into AMR NUMA will be examined.
[1] Giraldo, F. & Restelli, M. (2008). A study of spectral element and discontinuous Galerkin methods for the NavierStokes equations in nonhydrostatic mesoscale atmospheric modelling: Equation sets and test cases. Journal of Computational Physics, 227, 38493877
[2] Restelli, M. & Giraldo, F.X. (2009). A conservative discontinuous Galerkin semiimplicit formulation for the NavierStokes equations in nonhydrostatic mesoscale modelling. SIAM J. Sci. Comp., 31, 22312257
[3] Kelly, J.F. & Giraldo, F.X. (2012). Continuous and discontinuous Galerkin methods for scalable nonhydrostatic atmospheric models: limitedarea mode. Journal of Computational Physics, in review (2012).

INI 1  
17:00 to 17:25 
Evaggelos Kritsikis (Laboratoire des Sciences du Climat et l'Environnement (LSCE)) Secondorder conservative remapping between unstructured spherical meshes
Remapping from one finitedimensional description of a function to another is a common problem in numerical modelling. In particular, information transit between meshes is required for, e.g., model coupling or mesh adaptation. In order to preserve the properties of a numerical scheme such as conservativity, accuracy, positivity, etc., the remapping algorithm must itself possess these properties. A secondorder, conservative remapping between unstructured spherical meshes is presented. Areas are computed exactly by the defect formula and gradients estimated by the Gauss formula. Data is treestructured, so that neighbour search is done in logarithmic time. In addition, the algorithm lends itself well to parallelisation. Numerical tests on various unstructured grids are given.

INI 1 
09:00 to 09:25 
Yiyuan Li (Chinese Academy of Sciences) The orthogonal curvilinear terrainfollowing coordinate for atmospheric models Session: Vertical Session
This study designs a novel terrainfollowing coordinate, called the “Orthogonal curvilinear terrainfollowing coordinate” (Orthogonal sigma coordinate, OScoordinate), which has the ability to handle both the wellknown “high level errors” and the “pressure gradient force (PGF) errors” above the steep terrain in the classic terrainfollowing coordinate (sigmacoordinate).
The OScoordinate is proposed through a way quite different from and against the sigmacoordinate. The basis vectors of OScoordinate which is unit and terrainfollowing are firstly designed, and then the corresponding definitions of every coordinates are solved, so that the scalar equations of OScoordinate are solved. The PGF computational form has only one term in each momentum equation of OScoordinate, thereby avoiding the PGF problem completely.
Finally, a unified framework is proposed to combine the equations of the Cartesian coordinate, sigmacoordinate and OScoordinate together via an “onoff”, and is used to implement idealized advection experiments for testing the OScoordinate. The evaluation on the experimental results show the vertical coordinate surfaces in the OScoordinate can be much smoother than those in the sigmacoordinate, and therefore the OScoordinate outperforms all along the sigmacoordinate. This new coordinate can be used for oceanic models as well, and needs to be tested in both the idealized and realcase experiments.

INI 1  
09:25 to 09:50 
Takeshi Enomoto (Kyoto University) Reexamination of nonhydrostatic formulations using the hydrostaticpressure based coordinates Session: Vertical Session
The rapid increase of computing power is making global nonhydrostatic simulations affordable. A natural approach is to extend the formulation to include the nonhydrostatic effect. The advantage of this approach is that the existing data assimilation and tools require minimal changes. ECMWF and JMA seem to pursue this approach. ECMWF has achieved TL7999 (corresponding to approximately 2.5 km) with a fast Lendre transform using the butterfly algorithm (Nils Wedi, pers. comm.). Hiromasa Yoshimura (MRI/JMA) has built a nonhydrostatic version of JMA GSM using double Fourier series. Their formulations are based on Laprise (1992) that proposes the vertical coordinates based on hydrostatic pressure. Juang (1992, 2000) also adopts hydrostatic sigma–coordinates in the vertical but there are subtle differences. The latter introduces the hydrostatic temperature. In a limitedarea model, such as MSM, the hydrostatic temperature may be given externally. In a GCM, however, the hyd rostatic temperature must be determined internally if is not timeindependent. We investigated the two formations and found the assumption of the hydrostatic state of Laprise (1992) may be used to diagnose the hydrostatic temperature within MSM. Similarly the hydrostatic assumption of Arakawa and Konor (2009) can be used. MSM is found to run stably with any of these diagnosed hydrostatic states. The diagnosed hydrostatic temperature would enable the application of the formulation of MSM to the global domain.

INI 1  
09:50 to 10:15 
Margarida Moragues (Barcelona Supercomputing Center (BSCCNS)) A Variational Multiscale Stabilized Finite Element Method to solve the Euler Equations for Nonhydrostatic Stratified Benchmarks Session: Vertical Session
In this talk we present a Variational Multiscale Stabilization (VMS) for Compressible
Euler Equations applied to the Finite Element (FE) solution of nonhydrostatic stratified
flows. The VMS method was firstly presented by Hughes and coworkers [1] in the context
of incompressible flows. In the present work, recentely presented in [3], we extend
these concepts to Compressible Flows. In the framework of nonhydrostatic atmospheric
dynamics, we test the algorithm for problems at low Mach numbers.
A general version of the current compressible VMS technique was originally devised
for Computational Fluid Dynamics (CFD) of compressible flows without stratification [2].
The present work is justified by the previously observed good performance of VMS and
by the advantages that an elementbased Galerkin formulation offers on massively parallel
architectures, a challange for both CFD and Numerical Weather Prediction (NWP).
Unphysical vertical oscillations that may appear for not wellbalanced approximations
are a relevant problem in NWP, especially in the proximity of steep topography. In that respect, to properly discretize the dominant hydrostatics, a particular interpolation technique is proposed.
To evaluate the performance of the method in this context, some standard test cases of
stratified environments are presented.
References
[1] T. Hughes, Multiscale Phenomena: Green’s functions, the DirichlettoNeumann formulation,
subgrid scale models, bubbles and the origins of stabilized methods, CMAME 127 (1995)
387–401.
[2] M. Moragues, M. V´azquez, G. Houzeaux, R. Aubry, Variational Multiscale Stabilization of Compressible flows in Parallel Architectures, Parallel CFD 2010, Taiwan, May 2010.
[3] S. Marras, M. Moragues, M. V´azquez, O. Jorba, G. Houzeaux, A Variational Multiscale
Stabilized Finite Element Method for the Solution of the Euler Equations of Nonhydrostatic
Stratified Flows, J. Comput. Phys. (submitted 2012)

INI 1  
10:15 to 10:40 
Hiroe Yamazaki (University of Cambridge) 2D and 3D simulations of a nonhydrostatic atmospheric model on a blockstructured Cartesian mesh Session: Vertical Session
Under the rapid development of computing power, this study aims at developing a nextgeneration atmospheric model for ultrahigh resolution simulations at horizontal grid intervals of less than 100 meters. Recently, Cartesian grids are drawing attention as an attractive choice for highresolution atmospheric models that need to handle steep slopes in mountainous areas. They have the advantage of avoiding errors because of the slantwise orientation of grid lines in models using conventional terrainfollowing grids. In our model, a cut cell method is used for representing topography on a Cartesian grid. Also, a blockstructured mesh approach is introduced to achieve computationally efficient Cartesian grid simulations with both high vertical resolution near the ground and reasonable conservation characteristics.
The result of a 2D numerical simulation shows the model successfully reproduces a flow over a semicircular mountain on a locally refined mesh around the mountain. The result agrees well with that using a uniformly fine mesh. Some recent 3D results of our model will also be discussed.

INI 1  
10:40 to 11:10  Morning Coffee  
11:10 to 11:35 
Colm Clancy (Environment Canada) Semiimplicit predictorcorrector methods for atmospheric models Session: Time Session
A class of semiimplicit time integration schemes is proposed for use in atmospheric modelling. Various explicit predictorcorrector methods can be combined with an implicit treatment of the linear terms responsible for fast modes in order to improve computational stability. Linear analysis is used to identify promising algorithms. These are then tested in a finite volume shallow water model on an icosahedral grid. This framework has been developed at Environment Canada and will be extended to a more complete, threedimensional atmospheric model in the future. Experiments with standard test cases show that the proposed time schemes allow for stable integrations with relatively long timesteps, while maintaining a sufficient level of accuracy. In particular, they are more robust than the traditional, singlestage semiimplicit method using the leapfrog discretisation and do not require a time filter to control the computational mode.

INI 1  
11:35 to 12:00 
Peter Lynch (University College Dublin) Laplace transform integration and the slow equations Session: Time Session
We consider the Laplace transform (LT) filtering integration scheme applied to the shallow water equations, and demonstrate how it can be formulated as a finite difference scheme in the time domain by analytical inversion of the transform.
Both Eulerian and semiLagrangian versions of the scheme are analyzed. We show the relationship between the LT scheme and the slow equations. We demonstrate the advantages of the LT scheme by means of numerical integrations.

INI 1  
12:00 to 12:25 
Kevin Viner (Naval Research Laboratory) A stable treatment of conservative thermodynamic variables for semiimplicit semiLagrangian dynamical cores Session: Time Session
Atmospheric motions span a wide array of frequencies, the slowest of which provide us with our day to day weather. In numerical weather prediction it is necessary to narrow the focus to these lower frequencies for efficiency. The way this is typically done in global modeling is to apply an implicit timedifferencing method to terms linked to high frequency motions while applying an explicit method to terms linked to low frequency motions like advection, thereby allowing a larger stable time step. While semiLagrangian schemes further increase efficiency by removing the stability restrictions of the standard CFL condition concerning velocity, they are still held to a slightly different deformation CFL condition concerning the local variation in velocity. As such, it is still necessary to slow these fast wave modes in the semiLagrangian framework to maintain a stable system.
SemiLagrangian systems employing a conservative thermodynamic variable such as potential temperature as a prognostic variable face a unique dilemma in applying this standard method because the term responsible for gravity wave generation is also a vertical advection term which is absorbed into the total derivative. Application of the scheme in the absence of an explicit gravity wave term results in an unstable system since the fast gravity mode frequencies are not properly reduced. Stability can be maintained at the expense of both accuracy and efficiency by way of artificial damping and reduced time steps.
The modification discussed here defines a basic state potential temperature field which is advected in an Eulerian fashion while the semiLagrangian method is applied to the perturbation potential temperature. As the gravitywave term in question is now expressed explicitly through the total derivative of the basic state potential temperature, gravity mode stability is returned to the system; the semiimplicit treatment of the new term manages stability with respect to the associated CFL condition. Tests of the new scheme in the Navy Global Environmental Model (NAVGEM) show that the method allows stable integration at typical semiLagrangian time steps in the absence of artificial damping.

INI 1  
12:25 to 12:50 
Matthew Norman (Oak Ridge National Laboratory) MultiMoment ADERTaylor Methods for Systems of Conservation Laws with Source Terms in One Dimension Session: Time Session
A new integration method combining the ADER time discretization with a multimoment finitevolume framework is introduced. ADER runtime is reduced by performing only one CauchyKowalewski (CK) procedure per cell per time step, and by using the Differential Transform Method for highorder derivatives. Three methods are implemented: (1) singlemoment WENO [WENO], (2) twomoment Hermite WENO [HWENO], and (3) entirely local multimoment [MMLoc]. MMLoc evolves all moments, sharing the locality of Galerkin methods yet with a constant time step during p refinement.
Five experiments validate the methods: (1) linear advection, (2) Burger's equation shock, (3) transient shallowwater (SW) , (4) steadystate SW simulation, and (5) SW shock. WENO and HWENO methods showed expected polynomial h refinement convergence and successfully limited oscillations for shock experiments. MMLoc showed expected polynomial h refinement and exponential p refinement convergence for linear advection and showed subexponential (yet superpolynomial) convergence with p refinement in the SW case.
HWENO accuracy was generally equal to or better than a fivemoment MMLoc scheme. MMLoc was less accurate than RKDG at lower refinements, but with greater h  and p convergence, RKDG accuracy is eventually surpassed. The ADER time integrator of MMLoc also proved more accurate with p refinement at a CFL of unity than a semidiscrete RK analog of MMLoc. Being faster in serial and requiring less frequent internode communication than Galerkin methods, the ADERbased MMLoc and HWENO schemes can be spatially refined and have the same runtime, making them a competitive option for further investigation.

INI 1  
12:50 to 13:30  Lunch at Wolfson Court  
19:30 to 22:00  Conference Dinner at St Catharine's College 
09:00 to 09:25 
Todd Ringler (Los Alamos National Laboratory) A MultiResolution Modeling Approach for Global Ocean Modeling Session: Tests Cases 
INI 1  
09:25 to 09:50 
Jared Whitehead (Los Alamos National Laboratory) Potential Vorticity: A Diagnostic Tool for General Circulation Models Session: Tests Cases
Maintaining correlation between tracers and the dynamical wind and temperature fields with which they interact is a desirable trait of climate and weather models. A systematic, explicit test is developed to measure the consistency between a dynamical core's integration of the momentum equation, and its tracer transport algorithm. Potential vorticity is used as a diagnostic tool allowing direct comparison between the treatment of the dynamics and the integration of passive tracers. Several quantitative and qualitative metrics are suggested to measure this consistency including gridindependent probability density functions. Comparisons between the four primary dynamical cores of the National Center for Atmospheric Research's (NCAR) Community Earth System Model 1.0 (CESM) are presented. It is found that the finite volume (CAMFV) and spectralelement (CAMSE) dynamical cores perform better than the spectraltransform based Eulerian (CAMEUL) and semiLagrangian (CA MSLD) dynamical cores in the presence of a breaking wave.

INI 1  
09:50 to 10:15 
Peter Lauritzen (National Center for Atmospheric Research) A standard test case suite for twodimensional linear transport on the sphere Session: Tests Cases
It is the purpose of this paper to propose a standard test case suite for twodimensional transport schemes on the sphere intended to be used for model development and facilitating scheme intercomparison. The test cases are designed to assess important aspects of accuracy in geophysical fluid dynamics such as numerical order of convergence, “minimal” resolution, the ability of the transport scheme to preserve filaments, transport “rough” distributions, and to preserve preexisting functional relations between species/tracers under challenging flow conditions.
The experiments are designed to be easy to set up. They are specified in terms of two analytical wind fields (one nondivergent and one divergent) and four analytical initial con ditions (varying from smooth to discontinuous). Both conventional error norms as well as novel mixing and filament preservation diagnostics are used that are easy to implement. The experiments pose different challenges for the range of transport approaches from Lagrangian to Eulerian. The mixing and filament preservation diagnostics do not require an analytical/reference solution, which is in contrast to standard error norms where a “true” solution is needed. Results using a dozen stateoftheart transport schemes that participated in the 2011 NCAR workshop on transport are presented.

INI 1  
10:15 to 10:40 
Mike Cullen (Met Office) Evaluating numerical methods by using asymptotic limit solutions Session: Tests Cases
Numerical models of the atmosphere and ocean 'solve' the governing NavierStokes equations, but because of the complexity of the true solutions, can only do so in a grossly averaged sense. Standard numerical analysis theory is then of limited use, because it expects solutions to be close to the truth. The NavierStokes equations are known to have solutions close to computable asymptotic limits in many cases of practical interest. This talk gives examples, including dynamics/boundarylayer interaction, frontal solutions where the asymptotic limit is singular, and longterm behaviour of baroclinic systems, together with numerical demonstrations. The technique should help to validate the new integration schemes being proposed.

INI 1  
10:40 to 11:10  Morning Coffee  
11:10 to 11:35 
Paul Ullrich (University of California, Davis) MCore: Towards a MultiResolution NonHydrostatic FiniteVolume Dynamical Core Session: From Dynamical Cores to Operational Models
This talk will focus on recent efforts towards developing a true multiresolution framework for modeling regional climate in a global domain. In particular, I will present ongoing work on the highorder nonhydrostatic finitevolume MCore atmospheric general circulation model, which has been designed for massively parallel systems and with a focus on using mesh refinement to seamlessly tie together regional and global scales. Further, efforts to develop test cases of intermediate complexity, which include moisture and simplified physics parameterizations, have been undertaken and simulated within MCore. These efforts will be presented, along with a look at some upcoming approaches to numerical modeling using finitevolume methods on the sphere.

INI 1  
11:35 to 12:00 
Jin Lee (National Oceanic and Atmosphere Administration) A ThreeDimensional FiniteVolume Nonhydrostatic Icosahedral Modle (NIM) Session: From Dynamical Cores to Operational Models
The Nonhydrostatic Icosahedral Model (NIM) formulates the latest numerical innovation of the threedimensional finitevolume control volume on the quasiuniform icosahedral grid suitable for ultrahigh resolution simulations. NIM’s modeling goal is to improve numerical accuracy for weather and climate simulations as well as to utilize the stateofart computing architecture such as massive parallel CPUs and GPUs to deliver routine highresolution forecasts in timely manner. NIM uses innovations in model formulation similar to its hydrostatic version of the Flowfollowing Icosahedral Model (FIM) developed by Earth System Research Laboratory (ESRL) which has been tested and accepted for future use by the National Weather Service as part of their operational global prediction ensemble. Innovations from the FIM used in the NIM include:
* A local coordinate system remapped spherical surface to plane for numerical accuracy (Lee and MacDonald, 2009), * Grid points in a tabledriven horizontal loop that allow any horizontal point sequence (A.E. MacDonald, et al., 2010), * FluxCorrected Transport formulated on finitevolume operators to maintain conservative positive definite transport (J.L, Lee, ET. Al., 2010), * All differentials evaluated as finitevolume integrals around the cells, *Icosahedral grid optimization (Wang and Lee, 2011)
NIM extends the twodimensional finitevolume operators used in FIM into the threedimensional finitevolume solvers designed to improve pressure gradient calculation and orographic precipitation over complex terrain. The NIM dynamical core has been successfully verified with various nonhydrostatic benchmark test cases such as warm bubble, density current, internal gravity wave, and mountain waves. Physical parameterizations have been incorporated into the NIM dynamic core and successfully tested with multimonth aquaplanet simulations.

INI 1  
12:00 to 12:25 
Xingliang Li (China Meteorological Administration) A multimoment constrained finite volume model for nonhydrostatic atmospheric dynamics Session: From Dynamical Cores to Operational Models
Two dimensional non‐hydrostatic compressible dynamic cores for atmosphere are developed by using a new nodal type high order conservative method, so called multi‐moment constrained finite volume (MCV) method. Different from conventional finite volume method, the predicted variables (unknowns) in an MCV scheme are the values at the solution points distributed within each mesh cell. The time evolution equations to update the unknown point values are derived from a set of constraint conditions based on multi‐moment concept, where the constraint on the volume integrated average (VIA) for each mesh cell is cast into a flux form and thus guarantees rigorously the numerical conservation. Two important features makes MCV method particularly attractive as an accurate and practical numerical framework for atmospheric and oceanic modelling. 1) Using nodal values at uniformly located solution points as the predicted variables provides great convenience in de aling with complex geometry and source terms, and 2) High order schemes can be built by using constraints in terms of different moments, which makes the numerical schemes more flexible and efficient. We present in this paper the dynamic cores that use the third and the fourth order MCV schemes. We have verified the numerical outputs of both schemes by widely used standard benchmark tests and obtained competitive results. The present numerical cores provide a promising and practical framework for further development of non‐hydrostatic compressible atmospheric models.

INI 1  
12:30 to 13:30  Lunch at Wolfson Court  
13:35 to 14:00 
James Kent (University of Michigan) DCMIP 2012: Tracer Transport Tests in Dynamical Cores Session: From Dynamical Cores to Operational Models
We investigate the accuracy of tracer transport algorithms in dynamical cores of weather and climate models using prescribed velocity test cases. The tests make use of threedimensional flow, with the focus on extreme deformation, horizontalvertical coupling, and flow in the presence of orography. These tests highlight common problems associated with advection in dynamical cores. Standard error measures allow us to assess convergence rates, monotonicity, and the effects of the diffusion mechanisms in tracer transport algorithms. We will present results form the 2012 Dynamical Core Model Intercomparison Project (DCMIP) workshop.

INI 1  
14:00 to 14:25 
Sasa Gabersek (Naval Research Laboratory) Numerical Simulations with a ThreeDimensional Spectral Element Model Session: From Dynamical Cores to Operational Models
Highly accurate numerical methods for solving partial differential equations that have been traditionally used in the computational fluid dynamics have yet to be fully exploited for geophysical fluid dynamics applications, such as weather prediction. We will present results obtained with a fully compressible, nonhydrostatic spectral element (SE) model in three dimensions. All of our results are obtained using an eighth order polynomial (p=8), which provides the best compromise between accuracy and computational cost. The number of elements (h), into which the computational domain is decomposed, is varied among experiments to achieve different effective spatial resolutions.
Introducing physical parameterizations into an SE model is challenging due to a number of aspects including the varying nodal spacing within each element. We highlight some of these challenges and our approaches to address them. The model is applied in a series of idealized experiments that include physical process parameterizations such as: i) flow over complex terrain with subgrid scale mixing, and ii) a sea breeze that includes an evolving planetary boundary layer with surface fluxes.
We discuss the issues that arise when physical parameterizations are introduced into SE models. Also, scalability results over many computational cores will be addressed. In addition, we will compare the efficiency and accuracy of the results with a finite difference model (FD).

INI 1  
14:25 to 14:50 
Celal Konor (Colorado State University) Design of a dynamical core based on the nonhydrostatic unified system of equations Session: From Dynamical Cores to Operational Models
In this talk, we present the design of a dry dynamical core based on the nonhydrostatic “unified system” of equations. The unified system filters vertically propagating acoustic waves. The dynamical core predicts the potential temperature and horizontal momentum. It uses the predicted potential temperature to determine the quasihydrostatic components of the Exner pressure and density. The continuity equation is diagnostic (and used to determine vertical mass flux) because the time derivative of the quasihydrostatic density is obtained from the predicted potential temperature. The nonhydrostatic component of the Exner pressure is obtained from an elliptic equation. The main focus of this paper is on the integration procedure of this unique dynamical core. Height is used as the vertical coordinate, and the equations are vertically discretized on a Lorenztype grid. Cartesian horizontal coordinates are used along with an Arakawa Cgrid in the planar version of the dynamical core. The global version of the dynamical core is based on the vorticity and divergence predicting Zgrid formulation. The performance of the model in simulating a wide range of dynamical scales in the planar and global domains is demonstrated through idealized extratropical cyclogenesis simulations, and warm and cold bubble test cases.

INI 1  
14:50 to 15:15 
Almut Gassmann (LeibnizInstitut für Atmosphärenphysik, Kühlungsborn) ICONIAP: A nonhydrostatic global model designed for energetic consistency Session: From Dynamical Cores to Operational Models
The talk describes a new global nonhydrostatic dynamical core (ICONIAP: Icosahedral Nonhydrostatic model at the Institute for Atmospheric Physics) on a hexagonal Cgrid which is designed to conserve mass and energy. Energy conservation is achieved by discretising the antisymmetric Poisson bracket which mimics correct energy conversions between the different kinds of energy (kinetic, potential, internal). Because of the bracket structure this is even possible in a complicated numerical environment with (i) the occurrence of terrainfollowing coordinates with all the metric terms in it, (ii) the horizontal Cgrid staggering on the Voronoimesh and the complications induced by the need for an acceptable stationary geostrophic mode, and (iii) the necessity for avoiding the Hollingsworthinstability. The model is equipped with a Smagorinskytype nonlinear horizontal diffusion. The associated dissipative heating is accounted for by the application of the discrete product rule for derivatives. The time integration scheme is explicit in the horizontal and implicit in the vertical. In order to ensure energy conservation, the Exner pressure has to be offcentered in the vertical velocity equation and extrapolated in the horizontal velocity equation.
Test simulations are performed for small scale and global scale flows. A test simulation of linear nonhydrostatic flow over a rough mountain range shows the theoretically expected gravity wave propagation. The baroclinic wave test is extended to 40 days in order to check the Lorenz energy cycle. The model exhibits excellent energy conservation properties even in this strongly nonlinear and dissipative case. The HeldSuarez test confirms the reliability of the model over even longer timescales.

INI 1  
15:15 to 15:45  Afternoon Tea  
15:45 to 16:10 
Nils Wedi (European Centre for MediumRange Weather Forecasts (ECMWF)) Global, nonhydrostatic, cloudpermitting, mediumrange forecasts using the spectral transform method: progress and challenges Session: From Dynamical Cores to Operational Models
Numerical weather prediction (NWP) requires an answer in real time with a window of approximately one hour to run a mediumrange global forecast such that it can be delivered in time to its customers worldwide. While computational efficiency remains one of the most pressing needs of NWP, it is an open question how to most efficiently use the computer power available over the next twenty years, while seeking the most accurate and costeffective solution. At the same time, there are significant scientific challenges to increase resolution further, changes to the governing equations, and how subgridscale (SGS) processes are represented. ECMWF plans to implement a global horizontal resolution of approximately 10km by 2015 for its assimilation and deterministic forecast system, and approximately 20km for the ensemble prediction system (EPS). The scales resolved at these resolutions are still hydrostatic and the efficiency of the contemporary hydrostatic, semiLagrangian, semiimp licit solution procedure using the spectral transform method is likely to remain a relevant benchmark. However, due to the relative cost increase of the Legendre transforms compared to the gridpoint computations, very high resolution spectral models may become prohibitively expensive. Moreover, spectraltogridpoint transformations require datarich global communications at every timestep that may become too expensive on massively parallel computers. Recent progress in the development of fast spherical harmonics transforms (Tygert, 2008,2010) based on the butterfly scheme (O’Neil et al, 2010) mitigate the computational expense of the spectral transforms. Results are presented that demonstrate the costeffectiveness of the fast Legendre transforms (FLT) on the ibm_power7 architecture. The FLTs save both memory and computing time enabling the "world's first" successful T7999 (or equivalently ~2.5km horizontal resolution) global weather forecast with a spectral transform model.
Tygert, M., Fast algorithms for spherical harmonic expansions, II, J. of Comput. Physics, Vol. 227 (8), 2008, 42604279.
Tygert, M., Fast algorithms for spherical harmonic expansions, III, J. of Comput. Physics, Vol. 229 (18), 2010, 61816192.
O’Neil, M., F. Woolfe, V. Rohklin, An algorithm for the rapid evaluation of special function transforms, Appl. Comput. Harmon. Anal., Vol. 28(2), 2010, 203226.

INI 1  
16:10 to 16:35 
Michail Diamantakis (European Centre for MediumRange Weather Forecasts (ECMWF)) Evaluation of mass fixing schemes for the ECMWF Integrated Forecasting System (IFS) Session: From Dynamical Cores to Operational Models
The purpose of this talk is to summarize the conservation properties of the ECMWF spectral semiLagrangian, semiimplicit IFS model and present results from some recently implemented algorithms for correcting global mass imbalances due to advection.
Several established algorithms of different complexity, such as the ones described in [1,2,3,4], will be summarized and assessed in the context of NWP forecasts. Results from water species advection and passive tracer advection will be presented. Strengths and limitations of these schemes will be analyzed as well as their ability to preserve monotonicity and the impact they have on forecast skill.
References
Bermejo R., Conde J., 2002: A conservative QuasiMonotone SemiLagrangian Scheme. MWR, 130, 423430.
A. Priestley, 1993: A QuasiConservative Version of the SemiLagrangian Advection Scheme, MWR, 121, 621629.
Zerroukat M., 2010, A simple mass conserving semiLagrangian scheme for transport problems, JCP, 229, 90119019.
Williamson D.L., Rasch P.J., 1994, Water vapor transport in the NCAR CCM2, Tellus, 46A, 3451.

INI 1  
16:35 to 17:00 
Michael Baldauf (Deutscher Wetterdienst (DWD)) An exact analytical solution for gravity wave expansion of the compressible, nonhydrostatic Euler equations on the sphere Session: From Dynamical Cores to Operational Models
For the development and assessment of dynamical cores for atmospheric simulation models, suitable idealized test setups with known solutions are very useful. But only in rare cases exact analytical solutions exist for the underlying equation systems. In this work a slightly modified version of the original idea of Skamarock, Klemp (1994) is proposed: the quasi linear expansion of sound and gravity waves on a sphere induced by a weak warm bubble. For this case an exact analytical solution for the compressible, nonhydrostatic Euler equations was found for a shallow atmosphere and optionally with inclusion of Coriolis effects for a 'spherical fplaneapproximation'.
This solution can be used as reference for convergence studies of global models. 'Small earth' convergence tests with the ICON model of the Deutscher Wetterdienst (DWD) and the MaxPlanck Institut of Meteorology (MPI) are shown.

INI 1  
17:00 to 17:25 
Christiane Jablonowski (University of Michigan) Highlights of the Dynamical Core Model Intercomparison Project (DCMIP) Session: From Dynamical Cores to Operational Models
The talk presents highlights of the Dynamical Core Model Intercomparison Project (DCMIP) and its associated 2week summer school that has been held at NCAR in August 2012. DCMIP paid special attention to nonhydrostatic global models that are an emerging trend, especially in the climate modeling community. These allow highresolution simulations and provide a pathway for embedded variableresolution meshes in regions of interest. The primary goals of DCMIP were to assess the scientific designs of over 16 current and forthcoming dynamical cores, to establish new nonhydrostatic and moist dynamical core test cases in the dynamical core community, and to introduce innovative cyberinfrastructure tools. The talk surveys the DCMIP dynamical core test suite. This test hierarchy with increasing complexity (1) assesses the accuracy of 3D advection schemes via 3D deformational flows, (2) evaluates the evolution of gravity waves in nonrotating model configurations, (3) sheds light on the impact of orography, especially in the presence of orographyfollowing vertical coordinates, (4) assesses dry baroclinic waves with dynamic tracers, and (5) suggests test cases of intermediate complexity that incorporate moisture and very simplified physical parameterizations. The latter include a new moist variant of the baroclinic wave test case, and idealized tropical cyclones. The test suite makes extensive use of smallplanet experiments that are suggested as a computationallyefficient tool for highresolution nonhydrostatic model evaluations. Highlights of the DCMIP model results are shown.

INI 1 
09:00 to 09:25 
Günther Zängl (Deutscher Wetterdienst (DWD)) The Icosahedral Nonhydrostatic (ICON) model: formulation of the dynamical core and physicsdynamics coupling Session: From Dynamical Cores to Operational Models
The nonhydrostatic dynamical core of the ICON model is formulated on an icosahedraltriangular Cgrid and uses the edgenormal wind component, the vertical wind speed, density and virtual potential temperature as prognostic variables. In the vertical, a generalized terrainfollowing coordinate based on height is used that allows for a rapid decay with height of topographic structures, thereby reducing the numerical discretization errors over steep mountains. Moreover, a truly horizontal formulation of the horizontal pressure gradient term is used that greatly improves numerical stability over steep mountains. The spatial discretizations are secondorder, with some slight degradation near the pentagon points of the basic icosahedron, and are primarily optimized for computational efficiency, ensuring strict mass conservation and consistent tracer transport but no exact conservation of energy or vorticityrelated quantities. Time integration is based on a secondorder predictorcorrector scheme that is fully explicit in the horizontal and implicit for the terms entering into vertical sound wave propagation. The dynamics time step is therefore limited by the CFL stability condition for horizontal sound wave propagation, but a longer time step (typically 4x or 5x) is used for tracer transport and fastphysics parameterizations. A important aspect of physicsdynamics coupling is that all terms related to latent heat release have to be converted into temperature changes at constant volume because density is used as a prognostic variable and is kept fixed during the physics call.

INI 1  
09:25 to 09:50 
Mark Taylor (Sandia National Laboratories) CAM highresolution AMIP simulations using a spectral finite element dynamical core Session: From Dynamical Cores to Operational Models
We will present results from highresolution AMIP simulations using CAMSE: The Community Atmosphere Model (CAM), running with the Spectral finite Element dynamical core from NCAR's HighOrder Method Modeling Environment. CAMSE is now supported within the Community Earth System Model. It uses fully unstructured conforming quadrilateral grids, and thus supports global high resolution through quasiuniform grids and localized regions of high resolution with unstructured grids. For global 1/4 and 1/8 degree resolutions, CAMSE runs efficiently on hundreds of thousands of processors on modern Cray and IBM supercomputers and obtains excellent simulation throughput. For variable resolution grids, our simulations benefit from a new mixedformulation tensorbased hyperviscosity operator. The mixedformulation allows for efficient implementation in finite element methods, and a tensor form is used to provide the correct scaleaware dissipation, especially in the distorted elements which appear in resolution transition regions.

INI 1  
09:50 to 10:15 
John McGregor (CSIRO Marine and Atmospheric Research (CMAR)) Comparing the formulations of CCAM and VCAM and their performance as atmospheric GCMs Session: From Dynamical Cores to Operational Models
Two cubebased atmospheric GCMs have been developed at CSIRO, the ConformalCubic Atmospheric Model (CCAM) and, more recently, the Variable Cubic Atmospheric Model (VCAM). The formulations of the dynamical cores of both models will be described and compared. CCAM is formulated on the conformalcubic grid, whereas VCAM is cast on the equiangular gnomoniccubic grid. CCAM is a 2timelevel semiLagrangian semiimplicit Eulerian model, whereas VCAM employs a splitexplicit fluxconserving approach. Both models use reversible staggering for the wind components (McGregor, MWR, 2005) to produce good wave dispersion behavior. CCAM employs several orographic treatments that are not available for use in the VCAM dynamical core, with the interesting consequence that only VCAM requires hybrid vertical coordinates.
Both models include the same efficient messagepassing framework. Although VCAM avoids the messagepassing overheads necessitated by the Helmholtz solver of CCAM, if does have some overheads from more frequent calls to the wind staggering/unstaggering routines.
Aspects of the climatologies of both models will be compared and their overall advantages and disadvantages discussed. VCAM is presently being coupled to the Parallel Cubic Ocean Model (PCOM) from JAMSTEC; both the atmosphere and ocean model employ identical grids, and progress on this activity will be briefly presented.

INI 1  
10:15 to 10:40 
Thomas Melvin & Nigel Wood (Met Office) The dynamical core of the Met Office's Unified Model Session: From Dynamical Cores to Operational Models
The Met Office's Unified Model (MetUM) is used across a wide range of spatial and temporal scales, from very high resolution nowcasting to centennial climate predictions. Key to the accuracy and efficiency of its predictions is the dynamical core. The scalability of the dynamical core is becoming increasingly critical to the efficiency of the model as a whole.
An update on the progress, benefits and plans for the replacement dynamical core of the Met Office's Unified Model (ENDGame) will be given together with an overview of the project to design the next generation dynamical core for the the Unified Model (GungHo).

INI 1  
10:40 to 11:10  Morning Coffee  
11:10 to 11:35 
Abdessamad Qaddouri (Environment Canada) The CMC 15Km operational deterministic global forecast system with YinYang grid Session: From Dynamical Cores to Operational Models
At Canadian meteorological center (CMC), we are currently implementing the future 15km operational deterministic global forecast system. In the horizontal we use spherical coordinates on the overset YinYang grid while in the vertical, we use a loghydrostaticpressure coordinate on the CharneyPhillips grid. The global forecast on the YinYang grid is performed by considering a domain decomposition (a twoway coupling method) between two limitedarea models (LAM), discretized on the two panels of YinYang grid and using the same time step. Each panel of the YinYang grid system is extended by a static halo region that plays the same role as the piloting region in limitedarea modeling. Since the two subgrids of the YinYang grid do not match, the update of the variables in the pilot region is done by interpolation.
Each limitedarea model uses the same fully implicit semiLagrangian method as in the GEM operational model to solve its own dynamic core. Geophysical fields are produced separately for each LAM grid however, The two LAM models use the same parametrization schemes for the physics configuration. For data assimilation, we use the EnsembleVar (EnVar) approach where the 4D ensemble covariances are used to produce 4D analysis without TL/AD models.
In this presentation we will present some preliminary results of this work.

INI 1  
11:35 to 12:00 
William Skamarock (National Center for Atmospheric Research) The global nonhydrostatic atmospheric model MPAS: Preliminary results from variableresolution mesh tests Session: From Dynamical Cores to Operational Models
The Model for Prediction Across Scales (MPAS) is comprised of geophysical fluid flow solvers using horizontallyunstructured spherical centriodal Voronoi meshes and a Cgrid staggering of the prognostic variables. We have constructed a global atmospheric model for these meshes that solves the fully compressible nonhydrostatic equations using a finitevolume formulation coupled with a splitexplicit time integration technique to handle acoustic modes. The Voronoi meshes are unstructured grids that permit variable horizontal resolution, and we plan to make use of variableresolution meshes in our weather and regional climate applications. Towards this end, we have begun testing the robustness of variableresolution global mesh simulations using both idealized test cases (baroclinic waves) and full NWP forecasts. Preliminary results for hydrostaticscale tests (dx > 10 km) indicate that flow features are appropriately resolved relative to the local resolution if the mesh t ransition zones provide a smooth transition from coarser to finer resolution regions and if model filters are appropriately scaled by the local mesh spacing. We will present these results and discuss their implications for variableresolution mesh configurations and for the end applications.

INI 1  
12:00 to 12:25 
Florian Prill (Deutscher Wetterdienst (DWD)) The Icosahedral Nonhydrostatic (ICON) model: Scalability on Massively Parallel Computer Architectures Session: Novel Optimisation Techniques
Simulation in numerical weather prediction and climate forecasting has a fastgrowing demand for memory capacity and processing speed. For the last decade, however, computer technology has shifted towards multicore chip designs while at the same time onchip clock rates have increased only moderately. The parallel implementation of the ICON model's nonhydrostatic dynamical core therefore follows a hybrid distributed/shared memory approach, based on the Message Passing Interface (MPI) and the OpenMP API.
The ICON code couples the different encapsulated components of the earth system model, e.g. dynamics, soil, radiation, and ocean, with highlevel language constructs. Its communication characteristics and programming patterns are designed to meet the main challenges in high performance computing, i.e. load balancing, cache efficiency, and lowlatency networking, and take the unstructured triangular Cgrid into account, which implies indirect addressing. Besides basic optimization strategies such as loop tiling and grid point reordering, the implementation employs special domain decomposition heuristics, parallel rangesearching algorithms with logarithmic complexity, and makes use of asynchronous I/O servers to deal with the potentially prohibitive amount of data generated by earth system models. This facilitates the ICON code to extract an adequate level of performance on a wide range of HPC platforms, targeting large scalar cluster systems with thousands of cores as well as vector computers.

INI 1  
12:30 to 13:30  Lunch at Wolfson Court  
13:30 to 13:55 
Hirofumi Tomita ([RIKEN/AICS]) SCALELES: Strategic development of large eddy simulation suitable to the future HPC Session: Novel Optimisation Techniques
The Large Eddy Simulation is a vital dynamical framework to investigate the cloudaerosolchemistryradiation interaction from the viewpoint of climate problem. So far, the LES using in the meteorological filed are having several problems. One problem was that it is large gridsize used, compromising to the suitability of LES. In addition, the aspect ratio of horizontal and vertical grids was much larger than unity. The gridsize must be reduced to several 10m and it is desirable that the aspect ratio is near unity for the atmospheric LES. The target domain was also narrow for less of computer resources. The largescale computing using the recent powerful supercomputer may enable us to conduct the LES with reasonable gridsize and wide domain. Ultimately, the global LES is one of milestones in near future. Another problem in LES applied on meteorological field is that the heat source owing to water condensation is injected in a grid box. Strictly considering, the gridbox heating collapse the theory of LES that the grid size is in the energy cascade domain. Nevertheless, we have used the dry theory of LES. Beside the above problem that should be resolved in the future, we are now confronting with computational problems for such largescale calculations. The numerical method of fluid dynamical part in the atmospheric model has been shifted from the spectral transform method to the gridpoint method. The former is no longer acceptable on the massively parallel platforms form the limitation of innerconnect communication. On the other hand, the latter also contains a new problem, which is socalled memory bandwidth problem. For example, even on K Computer, the B/F ratio is just 0.5. The key to get high computational performance is the reduction of load/store from and to the main memory and efficient use of cash memory. Similar problem occurs in the communication between computer nodes. The multidisciplinary team (Team SCALE) in RIKE/AICS is now tackling to such problems.

INI 1  
13:55 to 14:20 
Ryuji Yoshida ([RIKEN/AICS]) An InterComparison of Icosahedral Climate Models on the G8 Call: ICOMEX Project Session: Novel Optimisation Techniques
The ICOsahedralgrid Models for EXascale Earth system simulations (ICOMEX) is the consortium for the climate models development toward the exascale computing. It started since October 2011 as one of the G8 call projects. Participated are NICAM (Japan), ICON (Germany), MPAS (UK and US), and DYNAMICO (France) model teams. On the road to the exascale computing, we would find many road blocks, such as file I/O speed, lower byte/flops rate, outer/inner node communication. In order to examine these problems, six working groups are seated together. The Japan team works on the model intercomparison to be synergistic among working groups. Although all the participated models use the icosahedralgrid, the discretization methods are different. For example, the hexagonal shape of a control volume is used in NICAM, while the triangular shape is used in ICON. Through intercomparison of both computational and physical performances between models, we will find which aspects of the model configuration are advantagous toward exascale computing. Two types of experiments were performed until now using NICAM and ICON: the baroclinic wave test (Jablonowski and Williamson, 2006) and the statistical climatology test (Held and Suarez, 1994). The horizontal resolution is from 240km (glevel5) to 14km (glevel9), and higher resolution runs will be tested in future. For the Held and Suarez test case, NICAM and ICON simulated climatology similar to that shown in the original paper. The baroclinic wave test is also compared between the two models. To investigate computational aspects, we examined strong scaling of parallel computing. Tests were performed on the Westmere and Bulldozer machine by using 5 through 40 MPI processes. We found that the two models have a good scaling: the measured scaling efficiency is 0.80.9. We plan to perform the intercomparison experiments using K computer using a larger number of processes O(10^5) to examine detail profiles of very massive parallel cores.

INI 1  
14:20 to 14:45 
Jean Côté (UQAM  Université du Québec à Montréal) Timeparallel algorithms for weather prediction and climate simulation Session: Novel Optimisation Techniques
The forecast of weather relies on computer models that need to be executed in realtime, meaning that a forecast needs to be disseminated to users well before the time period for which it is made. A challenge in the future will be to succeed in using the computing power available in massively parallel highperformance computers and meet the realtime requirement. Until now weather forecast and related climate simulation models have taken advantage of the parallelism of the computers by dividing the task to be performed in the horizontal space dimensions.
The purpose of this work is to develop algorithms that allow also parallelism in the time dimension. This increased parallelism should allow an acceleration of the execution time of weather and climate models. This acceleration in turn permits an increase in the space accuracy of models while still meeting the realtime requirement. The talk will present our preliminary work on the "Parareal" algorithm that has been developed for that purpose (Lions et al., 2001) and whose applications to date have included among others air quality, but ignored weather forecasting.
Weather forecasting presents a challenge for the method because of the presence of waves and advection. An important question is to examine how the traditional way to accelerate models with the semiimplicit semiLagrangian methodology can be advantageously blended with the “Parareal” approach.

INI 1  
14:45 to 15:15  Afternoon Tea 