Non-hydrostatic semi-elastic hybrid-coordinate SISL extension of HIRLAM. Part I: numerical scheme

Article English OPEN
Rõõm, Rein ; Männik, Aarne ; Luhamaa, Andres (2007)

Two-time-level, semi-implicit, semi-Lagrangian (SISL) scheme is applied to the non-hydrostatic pressure coordinate equations, constituting a modified Miller–Pearce–White model, in hybrid-coordinate framework. Neutral background is subtracted in the initial continuous dynamics, yielding modified equations for geopotential, temperature and logarithmic surface pressure fluctuation. Implicit Lagrangian marching formulae for single time-step are derived. A disclosure scheme is presented, which results in an uncoupled diagnostic system, consisting of 3-D Poisson equation for omega velocity and 2-D Helmholtz equation for logarithmic pressure fluctuation. The model is discretized to create a non-hydrostatic extension to numerical weather prediction model HIRLAM. The discretization schemes, trajectory computation algorithms and interpolation routines, as well as the physical parametrization package are maintained from parent hydrostatic HIRLAM. For stability investigation, the derived SISL model is linearized with respect to the initial, thermally non-equilibrium resting state. Explicit residuals of the linear model prove to be sensitive to the relative departures of temperature and static stability from the reference state. Relayed on the stability study, the semi-implicit term in the vertical momentum equation is replaced to the implicit term, which results in stability increase of the model.
  • References (40)
    40 references, page 1 of 4

    Bates, J. R. and McDonald, A. 1982. Multiply upstream, semiLagrangian advective schemes: analysis and application to a multilevel primitive equation model. Mon. Wea. Rev. 110, 1831-1842.

    Bates, J. R., Moorthi, S. and Higgins, R. W. 1993. A global multilevel atmospheric model using a vector semi-Lagrangian finite-difference scheme. Part I: adiabatic formulation. Mon. Wea. Rev. 121, 244-263.

    Benard, P., 2003. Stability of semi-implicit and iterative centeredimplicit time discretization for various equation systems used in NWP. Mon. Wea. Rev. 131, 2479-2491.

    Benard, P., 2004. On the use of a wider class of linear systems for the design of constant-coefficient semi-implicit time schemes in NWP. Mon. Wea. Rev. 132, 1319-1324.

    Benoit, R., Desgagne´, M., Pellerin, P., Chartier, Y. and Desjardins, S. 1997. The Canadian MC2: A Semi-Lagrangian, Semi-Implicit Wideband Atmospheric Model Suited for Finescale Process Studies and Simulation. Mon. Wea. Rev. 125, 2382-2415.

    Bubnova, R., Hello, G., Bernard, P. and Geleyn, J.-F. 1995. Integration of the fully elastic equations cast in the hydrostatic pressure terrainfollowing coordinate in the framework of APREGE/Aladin NWP system. Mon. Wea. Rev. 123, 515-535.

    Coˆte´, J. and Staniforth, A. 1988. A two-time-level semi-Lagrangian semi-implicit scheme for spectral models. Mon. Wea. Rev. 116, 2003- 2012.

    Davies, H. C. 1976. A lateral boundary formulation for multilevel prediction models. Q. J. R. Meteorol. Soc. 102, 405-418.

    Girard, C., Benoit, R. and Desgagne´, M. 2005. Finescale Topography and the MC2 Dynamics Kernel. Mon. Wea. Rev. 133, 1463-1477.

    Golding, B. W. 1992. An efficient non-hydrostatic forecast model. Meteorol. Atmos. Phys. 50, 89-103.

  • Similar Research Results (1)
  • Metrics
    No metrics available
Share - Bookmark