Abstract
A pseudospectral method and its numerical solution algorithm for application to boundary layer problems in the atmosphere are presented. The method introduces the evaluation of a polynomial function when the solution is expressed as the sum of a periodic function and a polynomial function. The periodic function is then treated by Fourier expansion. In the paper, the accuracy of method has been demonstrated. Numerical results for a system of time dependent equations, modeling the atmospheric planetary boundary layer flow and nocturnal flow over terrain are encouraging. The method offers a promising alternative to finite-difference techniques.