Hostname: page-component-84b7d79bbc-4hvwz Total loading time: 0 Render date: 2024-07-28T12:30:22.690Z Has data issue: false hasContentIssue false
  • English
  • Français

A New Approach to Implement Sigma Coordinate in a Numerical Model

Published online by Cambridge University Press:  20 August 2015

Yiyuan Li ,
Donghai Wang  and
Bin Wang
Yiyuan Li*
Affiliation:
State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing, China
Donghai Wang*
Affiliation:
State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing, China
Bin Wang*
Affiliation:
State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing, China Ministry of Education Key Laboratory for Earth System Modeling, and Center for Earth System Science, Tsinghua University, Beijing, China
*
Email:liyiyuan@mail.iap.ac.cn
Corresponding author.Email:d.wang@hotmail.com
Email:wab@lasg.iap.ac.cn
Get access

Abstract

This study shows a new way to implement terrain-following σ-coordinate in a numerical model, which does not lead to the well-known “pressure gradient force (PGF)” problem. First, the causes of the PGF problem are analyzed with existing methods that are categorized into two different types based on the causes. Then, the new method that bypasses the PGF problem all together is proposed. By comparing these three methods and analyzing the expression of the scalar gradient in a curvilinear coordinate system, this study finds out that only when using the covariant scalar equations of σ-coordinate will the PGF computational form have one term in each momentum component equation, thereby avoiding the PGF problem completely. A convenient way of implementing the covariant scalar equations of σ-coordinate in a numerical atmospheric model is illustrated, which is to set corresponding parameters in the scalar equations of the Cartesian coordinate. Finally, two idealized experiments manifest that the PGF calculated with the new method is more accurate than using the classic one. This method can be used for oceanic models as well, and needs to be tested in both the atmospheric and oceanic models.

Keywords

86-0886A0586A10Pressure gradient force (PGF)terrain-following sigma coordinatenon-orthogonalbasis vectorsnumerical modelingcomputational errors
Type
Research Article
Information
Communications in Computational Physics , Volume 12 , Issue 4 , October 2012 , pp. 1033 - 1050
Copyright
Copyright © Global Science Press Limited 2012

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Phillips, N. A., A coordinate system having some special advantages for numerical forecasting, J. Meteorol., 14 (1957), 184–185.Google Scholar
[2]Danard, M., Zhang, Q. and Kozlowski, J., On computing the horizontal pressure gradient force in sigma coordinates, Mon. Weather Rev., 121 (1993), 3173–3183.2.0.CO;2>CrossRefGoogle Scholar
[3]Berntsen, J. and Oey, L.-Y., Estimation of the internal pressure gradient in σ-coordinate ocean models: Comparison of second‐, fourth‐ and sixth-order schemes, Ocean Dyn., 60 (2010), 317–330.Google Scholar
[4]Smagorinsky, J., Strickler, R. F., Sangster, W. E., Manabe, S., Halloway, J. L. Jr. and Hem-bree, G. D., Prediction experiments with a general circulation model. Proc. Int. Symp. on Dynamics of Large Scale Atmospheric Processes, Moscow, USSR, 70–134, 1967.Google Scholar
[5]Steppeler, J., Hess, R., Schättler, U. and Bonaventura, L., Review of numerical methods for nonhydrostatic weather prediction models, Meteorol. Atmos. Phys., 82 (2003), 287–301.Google Scholar
[6]Hu, J. and Wang, P., The errors of pressure gradient force in high-resolution mesoscale model with terrain-following coordinate and its revised scheme (in Chinese), Chin. J. Atoms. Sci., 31 (2007), 109–118.Google Scholar
[7]Corby, G. A., Gilchrist, A. and Newson, R. L., A general circulation model of the atmosphere suitable for long period integrations, Q. J. R. Meteorol. Soc., 98 (1972), 809–832.CrossRefGoogle Scholar
[8]Gary, J. M., Estimate of truncation error in transformed coordinate,primitive equation atmospheric models, J. Atmos. Sci., 30 (1973), 223–233.2.0.CO;2>CrossRefGoogle Scholar
[9]Qian, Y. and Zhong, Z., General forms of dynamic equations for atmosphere in numerical models with topography, Adv. Atmos. Sci., 3 (1986), 10–22.Google Scholar
[10]Blumberg, A. F. and Mellor, G. L., A description of a three-dimensional coastal ocean circulation model. Three-dimensional coastal ocean models, N. S. Heaps, Ed., American Geophysical Union, 1987.CrossRefGoogle Scholar
[11]McCalpin, J. D., A comparison of second-order and fourth-order pressure gradient algorithms in a sigma coordinate ocean model, Int. J. Numer. Methods Fluids., 18 (1994), 361–383.Google Scholar
[12]Yang, X. and Qian, Y., The recurrent computational methods of pressure gradient force in the p-σ coordinate numerical models (in Chinese), Chin. J. Atmos. Sci., 27 (2003), 171–181.Google Scholar
[13]Sikiricć, M. D., Janekevicć, I. and Kuzmicć, M., A new approach to bathymetry smoothing in σ-coordinate ocean models, Ocean Modell., 29 (2009), 128–136.Google Scholar
[14]Berntsen, J., A perfectly balanced method for estimating the internal pressure gradients in σ-coordinate ocean models, Ocean Modell., 38 (2011), 85–95.CrossRefGoogle Scholar
[15]Pielke, R. A., Cotton, W. R., Walko, R. L., Tremback, C. J., Lyons, W. A., Grasso, L. D., Nicholls, M. E., Moran, M. D., Wesley, D. A., Lee, T. J. and Copeland, J. H., A comprehensive meteorological modeling system-RAMS, Meteorol. Atmos. Phys., 49 (1992), 69–91.Google Scholar
[16]Hodur, R. M., The naval research laboratory’s coupled ocean/atmosphere mesoscale prediction system (COAMPS), Mon. Weather Rev., 125 (1997), 1414–1430.2.0.CO;2>CrossRefGoogle Scholar
[17]Dudhia, J., cited 2011: WRF modeling system overview. http://www.mmm.ucar.edu/ wrf/users/tutorial/201001/WRF Overview_Dudhia.pdfGoogle Scholar
[18]Haidvogel, D. B., Wilkin, J. L. and Young, R., A semi-spectral primitive equation ocean circulation model using vertical sigma and orthogonal curvilinear horizontal coordinate, J. Comput. Phys., 94 (1991), 151–185.CrossRefGoogle Scholar
[19]Song, Y. and Haidvogel, D., A semi-implicit ocean circulation model using a generalized topography-following coordinate system, J. Comput. Phys., 115 (1994), 228–248.Google Scholar
[20]Haidvogel, D. B., Arango, H. G., Hedstrom, K., Beckmann, A., Malanotte-Rizzoli, P. and Shchepetkin, A. F., Model evaluation experiments in the north atlantic basin: Simulations in nonlinear terrain-following coordinates, Dyn. Atmos. Oceans., 32 (2000), 239–281.Google Scholar
[21]Gal-Chen, T. and Somerville, R. C. J., On the use of a coordinate transformation for the solution of the Navier-stokes equations, J. Comput. Phys., 17 (1975), 209–228.Google Scholar
[22]Xue, M., Droegemeier, K. K. and Wong, V., The advanced regional prediction system (ARPS) - a multi-scale nonhydrostatic atmospheric simulation and prediction model. Part I: Model dynamics and verification, Meteorol. Atmos. Phys., 75 (2000), 161–193.Google Scholar
[23]Wang, X. and Xiong, A., Advanced Fluid Dynamics (in Chinese), Huazhong University of Science and Technology Press, 2003.Google Scholar
[24] W Zdunkowski and Bott, A., Dynamics of the Atmosphere: A Course in Theoretical Meteorology, Cambridge University Press, 2003.Google Scholar
[25]Mellor, G. L., -Y Oey, L. and Ezer, T., Sigma coordinate pressure gradient errors and the seamount problem, J. Atmos. Oceanic Technol., 15 (1998), 1122–1131.Google Scholar