As the computing power increases, the grid size of atmospheric models moves toward the gray zone of turbulence (the scales on the order of the energy-containing range). Nevertheless, the roughness sublayer, which is a compartment of the inertial sublayer, has not been considered in high-resolution mesoscale models. This study coupled a roughness sublayer parameterization into the Weather Research and Forecasting model and evaluated its performance to simulate climate near the Earth's surface.