“…It is based on the wind tunnel experiment by Costa et al presented in [54], in which the authors simulated atmospheric boundary layer conditions with a system of vortex generators, barriers and surface roughness. The same configuration was also used in a numerical study [55]. We use the mean velocity field at the domain inlet based on B flow type configuration defined in [54,55], which is characterized by the following power law profile…”