We present a numerical simulation of flow around a surface mounted cube placed in a turbulent boundary layer which, although representing a typical wind environment, has been specifically tailored to match a series of wind tunnel observations. The simulations were carried out at a Reynolds number, based on the velocity U at the cube height h, of 20, 000 -large enough that many aspects of the flow are effectively Reynolds number independent. The turbulence intensity was about 18% at the cube height, and the integral length scale L u x was about 0.8 times the cube height h. The Jenson number Je = h/z0, based on the approach flow roughness length z0, was 600, to match the wind tunnel situation. The computational mesh was uniform with a spacing of h/32, aiding rapid convergence of the multigrid solver, and the governing equations were discretised using second order finite differences within a parallel multiblock environment. The results presented include detailed comparison between measurements and LES computations of both the inflow boundary layer and the flow field around the cube including mean and fluctuating surface pressures. It is concluded that provided properly formulated inflow and surface boundary conditions are used, LES is now a viable tool for use in wind engineering problems concerning flow over isolated bodies. In particular, both mean and fluctuating surface pressures can be obtained with a similar degree of uncertainty as usually associated with wind tunnel modelling.