“Generalized” vertical coordinate systems for nonhydrostatic primitive equations

1995 ◽  
Vol 4 (5) ◽  
pp. 196-202
Author(s):  
Franz Huber-Pock
2010 ◽  
Vol 67 (6) ◽  
pp. 1983-1995 ◽  
Author(s):  
Steven C. Hardiman ◽  
David G. Andrews ◽  
Andy A. White ◽  
Neal Butchart ◽  
Ian Edmond

Abstract Transformed Eulerian mean (TEM) equations and Eliassen–Palm (EP) flux diagnostics are presented for the general nonhydrostatic, fully compressible, deep atmosphere formulation of the primitive equations in spherical geometric coordinates. The TEM equations are applied to a general circulation model (GCM) based on these general primitive equations. It is demonstrated that a naive application in this model of the widely used approximations to the EP diagnostics, valid for the hydrostatic primitive equations using log-pressure as a vertical coordinate and presented, for example, by Andrews et al. in 1987 can lead to misleading features in these diagnostics. These features can be of the same order of magnitude as the diagnostics themselves throughout the winter stratosphere. Similar conclusions are found to hold for “downward control” calculations. The reasons are traced to the change of vertical coordinate from geometric height to log-pressure. Implications for the modeling community, including comparison of model output with that from reanalysis products available only on pressure surfaces, are discussed.


2021 ◽  
Author(s):  
Jerome Chanut ◽  
James Harle ◽  
Tim Graham ◽  
Laurent Debreu

<p>The NEMO platform possesses a versatile block-structured refinement capacity thanks to the AGRIF library. It is however restricted up to versions 4.0x, to the horizontal direction only. In the present work, we explain how we extended the nesting capabilities to the vertical direction, a feature which can appear, in some circumstances, as beneficial as refining the horizontal grid.</p><p>Doing so is not a new concept per se, except that we consider here the general case of child and parent grids with possibly different vertical coordinate systems, hence not logically defined from each other as in previous works. This enables connecting together for instance z (geopotential), s (terrain following) or eventually ALE (Arbitrary Lagrangian Eulerian) coordinate systems. In any cases, two-way exchanges are enabled, which is the other novel aspect tackled here.  </p><p>Considering the vertical nesting procedure itself, we describe the use of high order conservative and monotone polynomial reconstruction operators to remap from parent to child grids and vice versa. Test cases showing the feasibility of the approach are presented, with particular attention on the connection of s and z grids in the context of gravity flow modelling. This work can be considered as a preliminary step towards the application of the vertical nesting concept over major overflow regions in global realistic configurations. The numerical representation of these areas is indeed known to be particularly sensitive to the vertical coordinate formulation. More generally, this work illustrates the typical methodology from the development to the validation of a new feature in the NEMO model.</p>


2007 ◽  
Vol 133 (627) ◽  
pp. 1547-1558 ◽  
Author(s):  
Y. J. Rochon ◽  
L. Garand ◽  
D. S. Turner ◽  
S. Polavarapu

2003 ◽  
Vol 131 (8) ◽  
pp. 1931-1938 ◽  
Author(s):  
Andrew Staniforth ◽  
Nigel Wood

Abstract Previous analysis of the hydrostatic primitive equations using a generalized vertical coordinate is extended to the deep-atmosphere nonhydrostatic Euler equations, and some special vertical coordinates of interest are noted. Energy and axial angular momentum budgets are also derived. This would facilitate the development of conserving finite-difference schemes for deep-atmosphere models. It is found that the implied principles of energy and axial angular momentum conservation depend on the form of the upper boundary. In particular, for a modeled atmosphere of finite extent, global energy conservation is only obtained for a rigid lid, fixed in space and time. To additionally conserve global axial angular momentum, the height of the lid cannot vary with longitude.


2015 ◽  
Vol 8 (2) ◽  
pp. 1749-1800 ◽  
Author(s):  
T. Dubos ◽  
S. Dubey ◽  
M. Tort ◽  
R. Mittal ◽  
Y. Meurdesoif ◽  
...  

Abstract. The design of the icosahedral dynamical core DYNAMICO is presented. DYNAMICO solves the multi-layer rotating shallow-water equations, a compressible variant of the same equivalent to a discretization of the hydrostatic primitive equations in a Lagrangian vertical coordinate, and the primitive equations in a hybrid mass-based vertical coordinate. The common Hamiltonian structure of these sets of equations is exploited to formulate energy-conserving spatial discretizations in a unified way. The horizontal mesh is a quasi-uniform icosahedral C-grid obtained by subdivision of a regular icosahedron. Control volumes for mass, tracers and entropy/potential temperature are the hexagonal cells of the Voronoi mesh to avoid the fast numerical modes of the triangular C-grid. The horizontal discretization is that of Ringler et al. (2010), whose discrete quasi-Hamiltonian structure is identified. The prognostic variables are arranged vertically on a Lorenz grid with all thermodynamical variables collocated with mass. The vertical discretization is obtained from the three-dimensional Hamiltonian formulation. Tracers are transported using a second-order finite volume scheme with slope limiting for positivity. Explicit Runge–Kutta time integration is used for dynamics and forward-in-time integration with horizontal/vertical splitting is used for tracers. Most of the model code is common to the three sets of equations solved, making it easier to develop and validate each piece of the model separately. Representative three-dimensional test cases are run and analyzed, showing correctness of the model. The design permits to consider several extensions in the near future, from higher-order transport to more general dynamics, especially deep-atmosphere and non-hydrostatic equations.


Author(s):  
Vladimir Zeitlin

The chapter gives the foundations of modelling of large-scale atmospheric and oceanic motions and presents the ‘primitive equations’ (PE) model. After a concise reminder on general fluid mechanics, the main hypotheses leading to the PE model are explained, together with the tangent-plane (so-called f and beta plane) approximations, and ‘traditional’ approximation to the hydrodynamical equations on the rotating sphere. PE are derived in parallel for the ocean and for the atmosphere. It is then shown that, with a judicious choice of the vertical coordinate, the ‘pseudo-height’, in the atmosphere, these two sets of equations are practically equivalent. The main properties of PE are derived and the key concepts of wave–vortex dichotomy, and of slow and fast motions, are explained. The essential notion of potential vorticity is introduced and its conservation by fluid masses is demonstrated. Inertia–gravity waves are explained and their properties presented. Limitations of the hydrostatic hypothesis are demonstrated.


2015 ◽  
Vol 8 (10) ◽  
pp. 3131-3150 ◽  
Author(s):  
T. Dubos ◽  
S. Dubey ◽  
M. Tort ◽  
R. Mittal ◽  
Y. Meurdesoif ◽  
...  

Abstract. The design of the icosahedral dynamical core DYNAMICO is presented. DYNAMICO solves the multi-layer rotating shallow-water equations, a compressible variant of the same equivalent to a discretization of the hydrostatic primitive equations in a Lagrangian vertical coordinate, and the primitive equations in a hybrid mass-based vertical coordinate. The common Hamiltonian structure of these sets of equations is exploited to formulate energy-conserving spatial discretizations in a unified way. The horizontal mesh is a quasi-uniform icosahedral C-grid obtained by subdivision of a regular icosahedron. Control volumes for mass, tracers and entropy/potential temperature are the hexagonal cells of the Voronoi mesh to avoid the fast numerical modes of the triangular C-grid. The horizontal discretization is that of Ringler et al. (2010), whose discrete quasi-Hamiltonian structure is identified. The prognostic variables are arranged vertically on a Lorenz grid with all thermodynamical variables collocated with mass. The vertical discretization is obtained from the three-dimensional Hamiltonian formulation. Tracers are transported using a second-order finite-volume scheme with slope limiting for positivity. Explicit Runge–Kutta time integration is used for dynamics, and forward-in-time integration with horizontal/vertical splitting is used for tracers. Most of the model code is common to the three sets of equations solved, making it easier to develop and validate each piece of the model separately. Representative three-dimensional test cases are run and analyzed, showing correctness of the model. The design permits to consider several extensions in the near future, from higher-order transport to more general dynamics, especially deep-atmosphere and non-hydrostatic equations.


2019 ◽  
Vol 49 (7) ◽  
pp. 1809-1826 ◽  
Author(s):  
Knut Klingbeil ◽  
Johannes Becherer ◽  
Elisabeth Schulz ◽  
Huib E. de Swart ◽  
Henk M. Schuttelaars ◽  
...  

AbstractThis paper presents thickness-weighted averaging (TWA) in generalized vertical coordinates as a unified framework for a variety of existing tidal-averaging concepts in seas and estuaries. Vertical profiles of resulting residual quantities depend on the specific vertical coordinate, which is held fixed during the averaging process. This dependence is demonstrated through the application to one-dimensional analytical tidal flow with sediment transport, to field observations from a tidal channel, and to model results from a two-dimensional estuary. The use of different coordinate systems provides complementary views on the residual dynamics and stresses the importance of a correct interpretation of residual quantities obtained by tidal averaging.


1975 ◽  
Vol 26 ◽  
pp. 87-92
Author(s):  
P. L. Bender

AbstractFive important geodynamical quantities which are closely linked are: 1) motions of points on the Earth’s surface; 2)polar motion; 3) changes in UT1-UTC; 4) nutation; and 5) motion of the geocenter. For each of these we expect to achieve measurements in the near future which have an accuracy of 1 to 3 cm or 0.3 to 1 milliarcsec.From a metrological point of view, one can say simply: “Measure each quantity against whichever coordinate system you can make the most accurate measurements with respect to”. I believe that this statement should serve as a guiding principle for the recommendations of the colloquium. However, it also is important that the coordinate systems help to provide a clear separation between the different phenomena of interest, and correspond closely to the conceptual definitions in terms of which geophysicists think about the phenomena.In any discussion of angular motion in space, both a “body-fixed” system and a “space-fixed” system are used. Some relevant types of coordinate systems, reference directions, or reference points which have been considered are: 1) celestial systems based on optical star catalogs, distant galaxies, radio source catalogs, or the Moon and inner planets; 2) the Earth’s axis of rotation, which defines a line through the Earth as well as a celestial reference direction; 3) the geocenter; and 4) “quasi-Earth-fixed” coordinate systems.When a geophysicists discusses UT1 and polar motion, he usually is thinking of the angular motion of the main part of the mantle with respect to an inertial frame and to the direction of the spin axis. Since the velocities of relative motion in most of the mantle are expectd to be extremely small, even if “substantial” deep convection is occurring, the conceptual “quasi-Earth-fixed” reference frame seems well defined. Methods for realizing a close approximation to this frame fortunately exist. Hopefully, this colloquium will recommend procedures for establishing and maintaining such a system for use in geodynamics. Motion of points on the Earth’s surface and of the geocenter can be measured against such a system with the full accuracy of the new techniques.The situation with respect to celestial reference frames is different. The various measurement techniques give changes in the orientation of the Earth, relative to different systems, so that we would like to know the relative motions of the systems in order to compare the results. However, there does not appear to be a need for defining any new system. Subjective figures of merit for the various system dependon both the accuracy with which measurements can be made against them and the degree to which they can be related to inertial systems.The main coordinate system requirement related to the 5 geodynamic quantities discussed in this talk is thus for the establishment and maintenance of a “quasi-Earth-fixed” coordinate system which closely approximates the motion of the main part of the mantle. Changes in the orientation of this system with respect to the various celestial systems can be determined by both the new and the conventional techniques, provided that some knowledge of changes in the local vertical is available. Changes in the axis of rotation and in the geocenter with respect to this system also can be obtained, as well as measurements of nutation.


Sign in / Sign up

Export Citation Format

Share Document