Developing Dynamical Core of a GCM in Sigma-Pressure Vertical Coordinate

Authors

1 M.Sc. Student of Meteorology , University of Yazd

2 Professor of Meteorology , Institute of Geophysics

3 Assistant Professor, Physics Department - Atomic and Molecular

Abstract

Developing atmospheric General Circulation Models (GCMs) and, in particular, their dynamical cores are important steps towards future model improvements. For this reason, most GCMs of the atmosphere are developed in various vertical coordinates and grids. The diabatic contour-advective semi-Lagrangian method (DCASL) is an algorithm based on the use of contour representation for a fundamental dynamic quantity such as a potential vorticity. Previously, DCASL algorithms have been constructed for the shallow-water and multilayer Boussinesq primitive equations and also multilayer non-Boussinesq equations on the sphere using a hybrid terrain-following–isentropic (Sigma-Theta) vertical coordinate on the Charney-Philips grid (CP-grid). In this research, the results of the dynamical core constructed based on the DCASL algorithm are presented for the hybrid terrain-following–pressure (Sigma-P) vertical coordinate and compared with the results in the previously constructed model in coordinate. Also in this study, the performance of the coordinate models constructed using the CP-grid and the Lorenz grid (L-grid) are compared. For assessment and comparison of the dynamical cores in representations of different fields, we have used the baroclinic instability test case introduced by Jablonowski and Williamson in 2006. Qualitative and quantitative comparisons have been made with the results of four dynamical cores of global reference models (three hydrostatic dynamical cores that are parts of NCAR’s Community Atmosphere Model version 3 (CAM3) and also the dynamical core of the operational weather forecast model GME at the German Weather Service (DWD)). The comparisons show that the initialization problem present in the Sigma-Theta vertical coordinate in this test case is removed by the use of the Sigma-P vertical coordinate and thus results become ever closer to the references solutions. Also results of the Sigma-P vertical coordinate in the Lorenz grid for this dynamical core have the least difference with the reference results in the first 15 days of the test case.

Keywords

References

منابع
1- Collins, W. D., P. J. Rasch, B. A. Boville, J. J. Hack, J. R. McCaa, D. L. Williamson, J. T. Kiehl, B. P. Briegleb, C. M. Bitz, S. J. Lin, M. Zhang & Y. Dai, 2004, Description of the NCAR Community Atmosphere Model (CAM3.0), NCAR Technical Note NCAR/TN-464+STR, National Center for Atmospheric Research, Boulder, USA.
 
2- Jablonowski, C. & D. L. Williamson, 2006, A baroclinic instability test case for atmospheric model dynamical cores, Quart. J. Roy. Meteorol. Soc., 132, pp. 2943-2975.
 
3- Konor, C. S. & A. Arakawa, 1997, Design of an atmospheric model based on a generalized vertical coordinate, Mon. We Rev., 125, pp. 1649-1673.
 
4- Lin, S. J., 2004, A ‘vertically Lagrangian’ finite-volume dynamical core for global models. Mon. Weather Rev., 132, 2293-2307.
 
5- Majewski, D., D. Liermann, P. Prohl, B. Ritter, M. Buchhold, T. Hanisch, G. Paul, W. Wergen & J. Baumgardner, 2002, The Operational Global Icosahedral-Hexagonal Gridpoint Model GME: Description and high-resolution tests, Mon. Weather Rev., 130, pp. 319-338.
 
 
6- Mohebalhojeh, A. R., M. Joghataei & D. G. Dritschel, 2016, Toward a PV-based algorithm for the dynamical core of hydrostatic global models, Mon. Wea. Rev., 144, No. 7.
 
7- Simmons, A. J. & D. M. Burridge, 1981: An energy and angular momentum conserving vertical finite-difference scheme and hybrid vertical coordinates, Mon. Wea. Rev., 109, pp. 758-766.
 
 
Nivar
Volume 43, 106-107 - Serial Number 107
September 2019
Pages 108-114

Files

History

  • Receive Date: 02 August 2018
  • Revise Date: 11 November 2018
  • Accept Date: 14 November 2019
  • First Publish Date: 14 November 2019

Share

How to cite

Statistics