According to present views, acoustic-gravity waves (AGWs) existing in the upper
atmosphere may be generated near the Earth’s surface due to different sources and propagate upwards. An algorithm for three-dimensional numerical simulation of vertical propagation of nonlinear AGWs from the Earth’s surface to the upper atmosphere was developed recently. The algorithm of the solution of fluid dynamic equations uses finite-difference analogues of basic conservation laws. This approach allows us to select physically correct generalized wave solutions of the nonlinear equations.