A planetary boundary-layer numerical model containing third-order derivative terms

In this paper, the third-order derivative of velocity with respect to height is included in the traditional motion equations of the neutral PBL. The nonlinear equations are solved numerically to obtain the vertical distribution of wind in the PBL and some PBL characteristic parameters. Reasonable simulations of the Leipzig wind profile using these parameters show the success of this kind of nonlocal closure in a real PBL simulation.