The paper presents a numerical simulation of flow of a turbulent boundary layer, representing a typical wind environment and matching a series of wind tunnel observations. The simulations are carried out at a Reynolds number of 20,000, based on the velocity U at a pseudo-height h, and large enough that the flow be effectively Reynolds number independent. Some wall models are proposed for the LES(Large Eddy Simulation) of the turbulent boundary layer over a rough surface. The Jenson number, $J=h/z_0$, based on the roughness length $z_0$, is 600 to match the wind tunnel data. The computational mesh is uniform with a spacing of h/32, as this aids rapid convergence of the multigrid solver, and the governing equations are discretised using second order finite differences within a parallel multiblock environment. The results presented include the comparison between wind tunnel measurements and LES computations of the turbulent boundary layer over rough surface.