When the magnitude of sub-scale ographic forcing is comparable with explicitly ordinary dynamic forcing, the drag effect reduced by ographic gravity wave is to be significant for maintaining dynamic balance of atmo-spheric circulation, as well as the momentum and energy transport. Such sub-scale ographic forcing should be introduced into numerically atmospheric model by means of drag being parameterized. Furthermore, the currently mature ographic gravity wave drag (OGWD) parameterization, i.e., the so-called first-generation(based on lineal single-wave theoretical framework) or the second-generation drag parameterization (includ-ing an important extra forcing by the contribution of critical level absorption), cannot correctly and effectly describe the vertical profile of wave stress under the influence of ambient wind shearing. Based on aforemen-tioned consideration, a new two-wave scheme was proposed to parameterize the ographic gravity wave drag by means of freely propagating gravity waves. It starts with a second order WKB approximation, and treats the wave stress attenuations caused by either the selective critical level absorption or the classical critical level absorption explicitly; while in the regions where critical levels are absent, it transports the wave stress vertically by two sinusoidal waves and deposits them and then damps them according to the wave saturation criteria. This scheme is thus used to conduct some sample computations over the Dabie Mountain region of East China, as an example. The results showed that the new two-wave scheme is able to model the vertical distribution of the wave stress more realistically.