Numerical modeling of troposphere-induced gravity wave propagation

Sources of internal gravity waves (IGW) in the upper atmosphere are assumed to be meteorological processes in the troposphere. These sources are vertically and horizontally inhomogeneous and time dependent. In order to describe the IGW propagation from such sources, a numerical solution of a system of hydrodynamical equations is required. In addition, it is necessary to take into account the influence of the altitude latitude inhomogeneity of the temperature and wind fields on the IGW propagation as well as the processes of dissipation. An algorithm is proposed for numerical modelling of the IGW propagation over a limited area from tropospheric local sources to the upper atmosphere. The algorithm takes into account all the above features. A spectral grid method is used with the expansion of wave fields into the Fourier series over longitude. The upper limit conditions were obtained from the requirement of a limited energy dissipation rate in an atmospheric column. The no slip (zero velocity) condition was used at the Earth's surface.