A generalization of the two-dimensional spectral forecast equations is suggested, whereby the atmospheric flow field in horizontal and vertical directions is represented in terms of orthogonal functions, which are eigenfunctions of certain differential operators in the three-dimensional equations. The technique is applied to the quasi-geostrophic potential vorticity equation.