Modelling incompressible, neutrally stratified, barotropic, horizontally homogeneous and steady-state atmospheric boundary layer (ABL) is an important aspect in computational wind engineering (CWE) applications. The ABL flow can be viewed as a balance of the horizontal pressure gradient force, the Coriolis force and the turbulent stress divergence. While much research has focused on the increase of the wind velocity with height, the Ekman layer effects, entailing veering — the change of the wind velocity direction with height, are far less concerned in wind engineering. In this paper, a modified k-epsilon model is introduced for the ABL simulation considering wind veering. The self-sustainable method is discussed in detail including the precursor simulation, main simulation and near-ground physical quantities adjustment. Comparisons are presented among the simulation results, field measurement values and the wind profiles used in the conventional wind tunnel test. The studies show that the modified k-epsilon model simulation results are consistent with field measurement values. The self-sustainable method is effective to maintain the ABL physical quantities in an empty domain. The wind profiles used in the conventional wind tunnel test have deficiencies in the prediction of upper-level winds. The studies in this paper support future practical super high-rise buildings design in CWE.