A discretization of a hydrostatic primitive equation global atmospheric model on spherical and generalized hybrid vertical coordinates is described in NCEP Office Notes 445. The discretization in the horizontal using a spectral method with spherical transformation is as the same as used in NCEP global model, only the vertical discretization is illustrated (illus. 1). Energy, entropy, and angular momentum conservation are used as constraints to discretize the vertical integration by finite difference scheme. The entire atmosphere is divided into several layers (illus. 2); only pressure and vertical flux are specified at the interfaces, and other variables such as horizontal wind, temperature, specific humidity and specific amount of tracers are specified at each layer. Conservation is a constraint that requires the pressure at each layer to be averaged by the pressures at the immediate neighbor interfaces (the one above and one below a given layer). Since pressures are not in a logarithmic form, the relationship for pressure between layers and interfaces becomes simple, and with pressure equation not in logarithmic form, it provides mass conservation as an extra. Due to the generalized vertical coordinate, vertical flux is solved by applying local changes in the pressure and virtual temperature equations to the definition of the vertical coordinate. It solves vertical fluxes at all interfaces by a simple algebraic equation through tridiagonal inversion. For the sake of time splitting between dynamics and physics processes, the vertical flux obtained in the dynamics is without local change on isentropic surface, but the vertical advection is required in the model physics. The forwardtime uncentered semiimplicit time integration scheme is also given by this finite difference scheme in generalized vertical coordinates under vertical profiles of reference temperature and pressure. A specific definition of a generalized hybrid coordinate, including sigma, pressure and isentropic surfaces, is introduced to define pressure at interface. It is given by surface pressure and virtual temperature. These modifications of the generalized form provide computational saving, instead of solving for the pressure at all interfaces, only the surface pressure equation is needed. Though the elements in the matrixes for semiimplicit computations become more complicated than those in the generalized pressure equation at all levels, the computing time is not increased because the matrixes are of the same degree; also two matrixes reduce to vector computation only because the surface pressure is solved instead of pressure at all the interfaces. The preliminary results indicate the differences among all types of coordinates, sigma, sigmapressure and sigmatheta, are small (illus. 3), but a tendency to have improved performance is shown in hybrid coordinates of sigmapressure and sigmatheta. A bias of zonal mean temperature used to be existed in current operational model. The hybrid coordinates results, especially on sigmatheta, show an indication of a correction over this bias (see the following figure). The noises by divergence or vertical motion over topper layers are smaller or eliminated on hybrid coordinates. More works have to be done to provide proofed impacts before an operational implementation.
