According to present knowledge, acoustic-gravity waves (AGWs) observed in the upper atmosphere may arise near the Earth surface due to different sources and propagate upwards. Algorithms for two- and three-dimensional numerical simulation of vertical propagation and breaking of nonlinear AGWs from the Earth’s surface to the upper atmosphere were developed recently. The algorithms of the solution of fluid dynamic equations use finite-difference analogues of basic conservation laws. This approach allows us to select physically correct generalized wave solutions of the nonlinear equations. Horizontally moving periodical horizontal sinusoidal structures of vertical velocity on the Earth’s surface serve as AGW sources in the model. Numerical simulation was made in a region of the Earth atmosphere with dimensions up to several thousand kilometers horizontally and 500 km vertically. Vertical profiles of the mean temperature, density, molecular viscosity and thermal conductivity correspond to standard models of the a