This study took place in the Pyrenees Range, in the northeastern Iberian Peninsula. The Pyrenees extend longitudinally, separating the Iberian Peninsula from the rest of Europe, and high peaks around 3000 m arise from deep valleys. As a mountain range it creates a barrier to advection, in this case from the north and south, and typical meteorological phenomena of mountainous areas occur within it (inversions, Foehn effect, extreme wind-chill, snow storms). Thus, two specific valleys in Catalonia were considered, Val d'Aran and Cerdanya. In both valleys automatic weather stations (AWSs) are available at similar heights. Although these valleys are only 100 km apart, they have different climates. However, the main reason for developing the study was that Numerical Weather Prediction (NWP) has problems when forecasting temperatures in complex terrain areas, mainly in the valley floor in winter season.Firstly, different equations based on a multilinear regression were obtained for each weather station. Multilinear regression was considered in this case as the most suitable downscaling method and data used were provided by the AWSs and MM5 (PSU/NCAR mesoscale model) numerical weather prediction model outputs.These equations were obtained to set up a Geographically Weighted Regression (GWR) method, although this one was modified and changed to a Vertically Weighted Regression (VWR) in order to create vertical temperature profiles.