Forest fragmentation leads to changes in scalar concentration and momentum fluxes. These changes create superficial heterogeneity in forests, affecting mass, energy, and meteorological flows across various spatial scales. The abrupt alteration of canopy roughness from clearing-to-forest wind flow triggers instabilities downstream of the forest border, revealing the significant contribution of advective flows. Current methodologies investigate canopy flows with idealized foliage distribution, using the same Leaf Area Index as the experimental data. In these experiments, leaf density is typically higher near the top of the canopy. To simulate clearingto-forest wind flow, we employ Large Eddy Simulation (LES) with momentum drag dependent on foliage density. The maximum value of this drag occurs at the middle of the forest height (0.5h). Our results demonstrate a well-adjusted turbulent profile downstream of the forest. Near the border, increased skewness at the top and in the middle of the canopy indicates a large variability in the wind field. At the top of the canopy, we find