Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

IF value: 5.240
IF5.240
IF 5-year value: 5.768
IF 5-year
5.768
CiteScore value: 8.9
CiteScore
8.9
SNIP value: 1.713
SNIP1.713
IPP value: 5.53
IPP5.53
SJR value: 3.18
SJR3.18
Scimago H <br class='widget-line-break'>index value: 71
Scimago H
index
71
h5-index value: 51
h5-index51
Volume 7, issue 4
Geosci. Model Dev., 7, 1767–1778, 2014
https://doi.org/10.5194/gmd-7-1767-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.

Special issue: Isaac Newton Institute programme on multiscale numerics for...

Geosci. Model Dev., 7, 1767–1778, 2014
https://doi.org/10.5194/gmd-7-1767-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.

Development and technical paper 25 Aug 2014

Development and technical paper | 25 Aug 2014

An orthogonal terrain-following coordinate and its preliminary tests using 2-D idealized advection experiments

Y. Li1, B. Wang1,2, D. Wang3, J. Li1, and L. Dong1 Y. Li et al.
  • 1State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing, China
  • 2Ministry of Education Key Laboratory for Earth System Modeling, and Center for Earth System Science, Tsinghua University, Beijing, China
  • 3State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing, China

Abstract. We have designed an orthogonal curvilinear terrain-following coordinate (the orthogonal σ coordinate, or the OS coordinate) to reduce the advection errors in the classic σ coordinate. First, we rotate the basis vectors of the z coordinate in a specific way in order to obtain the orthogonal, terrain-following basis vectors of the OS coordinate, and then add a rotation parameter b to each rotation angle to create the smoother vertical levels of the OS coordinate with increasing height. Second, we solve the corresponding definition of each OS coordinate through its basis vectors; and then solve the 3-D coordinate surfaces of the OS coordinate numerically, therefore the computational grids created by the OS coordinate are not exactly orthogonal and its orthogonality is dependent on the accuracy of a numerical method. Third, through choosing a proper b, we can significantly smooth the vertical levels of the OS coordinate over a steep terrain, and, more importantly, we can create the orthogonal, terrain-following computational grids in the vertical through the orthogonal basis vectors of the OS coordinate, which can reduce the advection errors better than the corresponding hybrid σ coordinate. However, the convergence of the grid lines in the OS coordinate over orography restricts the time step and increases the numerical errors. We demonstrate the advantages and the drawbacks of the OS coordinate relative to the hybrid σ coordinate using two sets of 2-D linear advection experiments.

Publications Copernicus
Download
Citation