Reproduction of atmospheric boundary layer wind tunnel experiments by numerical simulation is achieved in this work by directly modeling with immersed boundary method the geometrical elements placed in the wind tunnel's floor to induce the desired characteristics to the boundary layer.The wind tunnel has a cross section of 2.2 m x 2.25 m, with an inlet region 14 m long and a working region 2 m long. Boundary layer development is shaped up with a series of cubical elements, 3 cm in side, placed in a regular staggered arrangement with a 15 cm spacement. Vortex induction, Standen spires type elements, of 13,4 cm height, and a wall of 31.5 cm height are placed at the inlet. This arrangement is used to reproduce a representative urban site boundary layer (figure 1).The numerical model is implemented on the basis of the open source modelcaffa3d.MBRi [Usera et al 2008], which uses a finite volume method over block structured grids, coupled with various LES approaches for turbulence modeling and parallelization through domain decomposition with MPI [Mendina et al 2013]. Simulations were setup with approximately 2 million cells per block, with a 26 block arrangement. The computational grid is horizontally uniform with a resolution of 1.04 cm x 1.04 cm and nonuniform in vertical direction with the grid points concentrated near the floor . The grid spacing is geometrically stretched away from the floor with a minimum value of 1mm. The time step is 0.1 second and the computation is distributed in 26 cores on the Cluster-FING infraestructure [www.fing.edu.uy/cluster]. The Immersed boundary method approach followed the work of [Liao et al 2009]. Numerical simulation results are compared to wind tunnel measurements for the mean velocity profiles (figure 2), rms profiles and spectrums, providing good overall agreement. We conclude that the Immersed Boundary Condition method is a promising approach to numerically reproduce ABL Boundary Layer development methods used in physical modeling.