Weather forecasts over mountainous terrain are challenging due to the complex topography that is necessarily smoothed by actual local-area models. As complex mountainous territories represent 20% of the Earth’s surface, accurate forecasts and the numerical resolution of the interaction between the surface and the atmospheric boundary layer are crucial. We present an assessment of the Weather Research and Forecasting model with two different grid spacings (1 km and 0.5 km), using two topography datasets (NASA Shuttle Radar Topography Mission and Global Multi-resolution Terrain Elevation Data 2010, digital elevation models) and four land-cover-description datasets (Corine Land Cover, U.S. Geological Survey land-use, MODIS30 and MODIS15, Moderate Resolution Imaging Spectroradiometer land-use). We investigate the Ortles Cevadale region in the Rhaetian Alps (central Italian Alps), focusing on the upper Forni Glacier proglacial area, where a micrometeorological station operated from 28 August to 11 September 2017. The simulation outputs are compared with observations at this micrometeorological station and four other weather stations distributed around the Forni Glacier with respect to the latent heat, sensible heat and ground heat fluxes, mixing-layer height, soil moisture, 2-m air temperature, and 10-m wind speed. The different model runs make it possible to isolate the contributions of land use, topography, grid spacing, and boundary-layer parametrizations. Among the considered factors, land use proves to have the most significant impact on results.