A finite‐volume integration method for computing pressure gradient force in general vertical coordinates

A finite-volume integration method is proposed for computing the pressure gradient force in general vertical coordinates. It is based on fundamental physical principles in the discrete physical space, rather than on the common approach of transforming analytically the pressure gradient terms in differential form from the vertical physical (i.e., height or pressure) coordinate to one following the bottom topography. The finite-volume discretization is compact, involving only the four vertices of the finite volume. The accuracy of the method is evaluated statically in a two-dimensional environment and dynamically in three-dimensional dynamical cores for general circulation models. The errors generated by the proposed method are demonstrated to be very low in these tests.

[1]  A. Arakawa,et al.  A Potential Enstrophy and Energy Conserving Scheme for the Shallow Water Equations , 1981 .

[2]  Shian‐Jiann Lin,et al.  Multidimensional Flux-Form Semi-Lagrangian Transport Schemes , 1996 .

[3]  A. Simmons,et al.  An Energy and Angular-Momentum Conserving Vertical Finite-Difference Scheme and Hybrid Vertical Coordinates , 1981 .

[4]  G. Stelling,et al.  On the approximation of horizontal gradients in sigma co‐ordinates for bathymetry with steep bottom slopes , 1994 .

[5]  Richard C. J. Somerville,et al.  On the use of a coordinate transformation for the solution of the Navier-Stokes equations , 1975 .

[6]  Shian-Jiann Lin,et al.  A Class of the van Leer-type Transport Schemes and Its Application to the Moisture Transport in a General Circulation Model , 1994 .

[7]  Shian-Jiann Lin,et al.  An explicit flux‐form semi‐lagrangian shallow‐water model on the sphere , 1997 .

[8]  A. Arakawa,et al.  Vertical Differencing of the Primitive Equations in Sigma Coordinates , 1983 .

[9]  A. Kasahara.  Various Vertical Coordinate Systems Used for Numerical Weather Prediction , 1974 .

[10]  A. Arakawa,et al.  Vertical Differencing of the Primitive Equations Based on the Charney–Phillips Grid in Hybrid &sigma–p Vertical Coordinates , 1996 .

[11]  A. Gilchrist,et al.  A general circulation model of the atmosphere suitable for long period integrations , 1972 .

[12]  Z. Janjic,et al.  Comparison of methods for reducing the error of the pressure gradient force in sigma coordinate models , 1986 .

[13]  N. A. Phillips,et al.  A COORDINATE SYSTEM HAVING SOME SPECIAL ADVANTAGES FOR NUMERICAL FORECASTING , 1957 .