Wall functions have been widely used in computational fluid dynamics (CFD) simulations and can save significant computational costs compared to other near-wall flow treatment strategies. However, most of the existing wall functions were based on the asymptotic characteristics of near-wall flow quantities, which are inapplicable in complex and non-equilibrium flows. A modified wall function is thus derived in this study based on flow over a plate at zero-pressure gradient, instead of on the basis of asymptotic formulations. Turbulent kinetic energy generation ($G_P$), dissipation rate (${varepsilon}$) and shear stress (${ au}_{omega}$) are composed together as the near-wall expressions. Performances of the modified wall function combined with the nonlinear realizable k-${varepsilon}$ turbulence model are investigated in homogeneous equilibrium atmosphere boundary layer (ABL) and flow around a 6 m cube. The computational results and associated comparisons to available full-scale measurements show a clear improvement over the standard wall function, especially in reproducing the boundary layer flow. It is demonstrated through the two case studies that the modified wall function is indeed adaptive and can yield accurate prediction results, in spite of its simplicity.