The great contribution into the formation of the general circulation, thermal regime, and composition of the middle and upper atmosphere comes from Internal gravity waves. One of the most important wave sources is the topography of the earth's surface. Oncoming stream wind interacts with the irregularities of the relief, the result is a stationary gravity waves. Using a mathematical model of the generation of atmospheric waves by the mountains and distribution of orographic waves in the middle and upper atmosphere allowed us to develop computer programs for calculating the wave acceleration and heating in the three-dimensional numerical model of general atmospheric circulation. The main new feature of the parameterization developed is using the theory of stationary waves propagation in an inhomogeneous atmosphere for calculating orographic wave accelerations and heating at different altitudes up to the upper atmosphere. Using the atmospheric general circulation model with the inclusion of the parametrization of the orographic waves allowed us to determine the sensitivity of the atmospheric circulation at all altitudes to the effects of the orographic waves. The parameterization was included in general circulation model MUAM-48 (Middle and upper atmospheric model). The calculated profiles of the wave heating are integrated into the equation of heat and zonal and meridional components of the vector wave acceleration are inserted into the equations of motion. To calculate the topography the database of earth surface ETOPO2 in increments of 2 minutes along the latitudes and longitudes is used. The exponentially increasing vertical profile of turbulent and molecular viscosity is given. The calculations are performed for the winter and summer seasons.