Numerical procedure for planetary wave solution

The newly-developed LLNL two-dimensional chemical-radiative-transport model requires a knowledge of the EP flux divergence as an input momentum forcing. The major contributions for this forcing term come from the synoptic wave in the troposphere, the planetary wave in the stratosphere and mesosphere and the gravity wave in the upper mesosphere. The major source of the zonal momentum forcing in the middle atmosphere is the nonlinear planetary wave breaking. This planetary wave breaking also plays a significant role of mixing the chemical tracers. Garcia suggested a way of parameterizing the planetary wave breaking by using, linear damping of the primary wave. In this note we describe the procedure of obtaining the wave solution for the parameterization