The understanding of atmospheric flows is crucial in the analysis of dispersion of a contaminant or pollutant, wind energy and air-quality assessment to name a few. Additionally, the effects of complex terrain and associated orographic forcing are crucial in wind energy production. Furthermore, the use of the Reynolds-averaged Navier-Stokes (RANS) equations in the analysis of complex terrain is still considered the "workhorse" since millions of mesh points are required to accurately capture the details of the surface. On the other hand, solving the same problem by means of the instantaneous governing equations of the flow (i.e., in a suite of DNS or LES) would imply almost prohibitive computational resources. In this study, numerical predictions of atmospheric boundary layers are performed over a complex topography located in Nygårdsfjell, Norway. The Nygårdsfjell wind farm is located in a valley at approximately 420 meters above sea level surrounded by mountains in the north and south near the Swedish border. Majority of the winds are believed to be originated from Torneträsk lake in the east which is covered with ice during the winter time. The air closest to the surface on surrounding mountains gets colder and denser. The air then slides down the hill and accumulates over the lake. Later, the air spills out westward towards Ofotfjord through the broader channel that directs and transforms it into highly accelerated winds. Consequently, one of the objectives of the present article is to study the influence of local terrain on shaping these winds over the wind farm. It is worth mentioning that we are not considering any wind turbine model in the present investigation, being the main purpose to understand the influence of the local surface topography and roughness on the wind flow. Nevertheless, future research will include modeling the presence of a wind turbine and will be published elsewhere. The governing equations of the flow are solved by using a RANS approach and by considering three different two-equation turbulence models: k-omega (k-ω), k-epsilon (k-ε) and shear stress transport (SST). Furthermore, the real topographical characteristics of the terrain have been modeled by extracting the required area from the larger digital elevation model (DEM) spanning over 100 km square. The geometry is then extruded using Rhino and meshed in ANSYS Fluent. The terrain dimensions are approximately 2000x1000 meter square.