Photovoltaic Prediction Software: Evaluation with Real Data from Northern Spain

available and modelling software tools are compared with real data for three different photovoltaic power plants in operation over the past 12 years. Both annual and monthly results are analyzed. Abstract: Prediction of energy production is crucial for the design and installation of PV plants. In this study, ﬁve free and commercial software tools to predict photovoltaic energy production are evaluated: RETScreen, Solar Advisor Model (SAM), PVGIS, PVSyst, and PV*SOL. The evaluation involves a comparison of monthly and annually predicted data on energy supplied to the national grid with real ﬁeld data collected from three real PV plants. All the systems, located in Castile and Leon (Spain), have three different tilting systems: ﬁxed mounting, horizontal-axis tracking, and dual-axis tracking. The last 12 years of operating data, from 2008 to 2020, are used in the evaluation. Although the commercial software tools were easier to use and their installations could be described in detail, their results were not appreciably superior. In annual global terms, the results hid poor estimations throughout the year, where overestimations were compensated by underestimated results. This fact was reﬂected in the monthly results: the software yielded overestimates during the colder months, while the models showed better estimates during the warmer months. In most studies, the deviation was below 10% when the annual results were analyzed. The accuracy of the software was also reduced when the complexity of the dual-axis solar tracking systems replaced the ﬁxed installation.


Introduction
The environmental commitments to which European countries have subscribed to make it essential to drastically reduce greenhouse gas emissions, to reduce energy consumption in both public and residential buildings, and to incorporate renewable energies as the main source of supply. Photovoltaic (PV) energy is a renewable energy that has undergone a faster development than any other and is now a mature and economically competitive technology [1].
Any technical analysis of renewable energy and its implementation requires computer tools. It is time-consuming to create new tools for each analysis when accessible tools exist. Reasonable predictions of PV system energy production levels are vital for continuous supply, from the perspective of both the consumer and the system integrator. These predictions should be attainable based on available or desired system sizing parameters, such as the model of the PV modules, performance rating, orientation, and power loss from wires. Specific commercially available, as well as free, downloadable prediction the study were described: geographical data needed, way to compute solar irradiation, the PV system simulation capacity, etc. The results of the simulation with computer programs are compared with real production data over one complete year of operation of three facilities located in Spain, each one with a different tracking system: fixed installation, installation with one axis sun-tracker and two-axis tracker PV system. The main novelty of the study is that data used to verify simulation came from grid tied facilities in operation: the data correspond with the actual system turnover values and was provided directly by the company that owned the facility. The lengthy data-collection period, from 2008 to 2020, allows us to discount the influence of either warm or cold years. In addition, annual comparisons have been complemented with a monthly analysis where the deviation error throughout the year is discussed.

Methodology
Five different software photovoltaic models: Photovoltaic Geographical Information System (PVGIS) [9], PVsyst [10], RETScreen [11], PV*SOL [12] and System Advisor Model (SAM) [13], are evaluated in this work. The AC energy generation estimates of each software package were compared with the real AC energy production data of three different types of PV plants. Table 1 details a comparison of the main features of each program, whether solar radiation databases can be used, as well as product libraries (PV module and inverter), the number of parameters that are required, difficulty of using the software, number of results calculated, shadow modeling, electrical loses and license costs. One of the greatest benefits of these programs is that preliminary studies of PV project feasibility may be performed even before fully developing the design of the PV plant. The main objective of this type of study is to estimate the energy generated with a low number of known parameters. Although the power plant can be described in detail in some of these programs, the comparison in this study is focused on a preliminary study, where only essential parameters have been defined: location, type of solar tracking technology, installed capacity, number of modules, and efficiency and peak power of the PV module and inverter. The parameters of each station are summarized in Table 2. All other parameters, such as wire and shadow losses, as and when required, can be set as default parameters or with the standard recommended values of the software manual. All the models implement different climate database options, either online or integrated, so the option recommended by each program was selected and compared when different options were available.

PVGIS
The Photovoltaic Geographical Information System (PVGIS) provides a map-based inventory of solar energy resource and an assessment of electricity generation from photovoltaic systems in Europe, Asia, Africa, and South America. The database includes monthly and yearly average values of the global irradiation on horizontal and inclined surfaces, as well as climatic parameters needed for an assessment of the potential PV electricity generation [14,15]. With the following applications, it is possible to choose between databases collected by satellite that includes data from 2005 to 2016: SARAH, NREL, CMSAF, ECMWF, and COSMO-REA6. It uses an interpolation of data models for geographic resolution.
Using PV-GIS, the PV facility can be described by selecting the type of solar tracking system, its integration or otherwise within the buildings, installed power, system performance, and PV panel technology (Si Crystalline CdTe, CIS o CIGS). PVGIS provides daily and monthly average global irradiation data and daily and monthly average PV power production. It can also be used to optimize the angles of orientation of the PV panels and their tilt.
The simplicity and free use of PVGIS mean that it is frequently used for initial energy estimations of PV facilities. The databases are widely implemented in Europe, Africa, and Asia with regular updates, and their reliability has been demonstrated in a large number of scientific papers [16][17][18][19][20].

PVSYST
PVSYST software is a European-based PV system predictor developed by the University of Geneva for the European Energy Center. It is suitable for studying, simulating and analyzing stand-alone and grid-connected PV systems. Additionally, it is also useful for solar-thermal systems and solar applications. PVSYST is a well-designed tool with an extensive number of features for detailing the PV system and includes an extensive database of different PV related products. PVSYST includes algorithms for thermal and wind effects calculation, wire losses and incidence angle modifier losses, albedo values, and a horizon adjustment capability. This software has received significantly high scores and recommendations in various comparative studies [15,[21][22][23].
PVSYST imports meteorological and global irradiation data from Meteonorm and NASA-SSE. In addition, other external databases can be introduced in the study. PVsyst implements theoretical models for the calculation of both global irradiation on an inclined plane and diffuse irradiation [24,25]. The database of available PV panels, inverters, and batteries is very broad, is directly updated by the manufacturers, but can also be imported from the PHOTON database and even permits the definition of "ad hoc" panels of different technologies: mono-Si, poli-Si, and thin film. It is compulsory to introduce the design of the system: distribution of PV panels, and connection with the inverter, tilted angle and sun-tracking system specification. The incidence angle can be optimized and shadows and other elements can be introduced with PVSYST that may have an impact on the performance of the installation. Although the PV plant can be precisely described, when looking for a preliminary study, some parameters are not pre-defined, which complicates any definition of the basic study.
The main data outputs are annual PV production (MW) and the specific PV production (kWh/kW p year), as well as the performance factor. A monthly diagram of standard production per installed kW p is generated, including the following information: power loss, total energy produced, the effective power to the output of the generator, the energy injected into the network, the electrical efficiency of the facility, and a table with the solar irradiation data (horizontal global irradiation, incident irradiation at the receptor). In addition, the tool generates a diagram of losses during the year, indicating the estimated losses due to temperature, wiring, inverter system, and shadows. It is a payment program, but a complete version is free for one month. Subsequently, the program functions in DEMO mode, with very limited options.

System Advisor Model (SAM)
System Advisor Model (SAM) [8] is provided by the National Renewable Energy Laboratory (NREL) of the U.S. Department of Energy. SAM provides a consistent framework both for analyzing and for comparing power system costs and performance across the range of solar technologies and markets, from PV systems for residential and commercial markets to concentrating solar power and large PV systems for utility markets. SAM combines an hourly simulation model with performance, cost, and finance models to calculate energy output, energy costs, and cash flows. The software can also account for the effect of incentives on cash flows. Most of SAM's inputs can be used as parametric variables for sensitivity studies [15,[26][27][28].
SAM input consists of hourly weather data and the software can read files in either TMY2, TMY3 (typical meteorological year [29]), and EPW (energy plus weather data file) formats from over 20 sources [30] and data that the user may introduce. The program provides horizontal infrared radiation from sky, global horizontal radiation, direct normal radiation, and diffuse horizontal radiation, directly from the database, or calculated through various user-selected models, such as Hay [25], Perez [24,31], and Reindl [32].
SAM includes a wide database of PV panels and inverters and it is possible to define new elements. The complete design of the facility cannot be entered, but the estimated losses in each section can be indicated: DC and AC wiring, transformer efficiency, thermal Appl. Sci. 2021, 11, 5025 6 of 14 and power losses from wires and the sun-tracking system. It also implements a position of the panel model to predict shadows in the solar plant.
As a result, the program provides a summary table with major energy and economic data of the facility and its temporal evolution. It is a very complete free program with great possibilities for the analysis of solar photovoltaic, thermal, concentration, wind, and geothermal energy projects of all sizes and complexity.

RETSCREEN
RETScreen is managed by the research center of Natural Resources Canada (NRCan). This software can be used to evaluate the energy production and savings, costs, emission reductions, financial viability, and risk for various types of renewable-energy and energyefficient technologies (RETs) [32][33][34][35][36][37]. In addition, the RETScreen software integrates a series of databases that help to overcome the costs and difficulties associated with gathering meteorological data, product performance data, etc. Hence, global meteorological data are incorporated in the RETScreen software. This meteorological database includes both the ground-based meteorological data and NASA's satellite-derived meteorological data sets.
RETScreen has a basic user-expandable PV panel. The database covers no other element of the PV facility. The following information must be inputted: type of panel (poly-Si, mono-Si, etc.), total PV installed power of facility, efficiency, estimated electrical losses, losses caused by shadows. The inverter data, defined by the user, are input power, performance and estimated losses.
The economic analysis of the system requires the definition of estimated operational life, initial investment, loans, and grants, and the program can be used to analyze various economic scenarios.
As output, it generates different tables and graphs: daily solar horizontal irradiation by month, the estimation of annual electricity production (MWh), environmental benefits (CO 2 emissions savings, equivalent in barrels of petroleum), return on investment period, evolution of the benefits over the operating life of the installation, internal rate of return (IRR), net present value (NPV), and unit cost of energy. In the most recent updated version, RETScreen Expert, has become a paid program.

PV*SOL
PV*SOL is software for planning, design, and simulation of photovoltaic systems developed by Valentin Software. In addition to PV*SOL, this company has other software tools for solar thermal systems (T*SOL) and geothermal solar systems (GeoT*SOL). PV*SOL offers the possibility of configuring different types of PV systems, i.e., stand-alone systems with or without backup generator, grid-connected systems, and grid-connected systems with electrical appliances, which can include battery systems and electric vehicles [38].
Site-dependent climate data supplied by PV*SOL is based on Meteonorm database. It is also possible to load user-defined climate data measured at a specific meteorological station. There is a premium version of the program, PV*SOL premium, for 3D design of the system including environmental obstacles that can influence the irradiation received on the PV module.
When simulating the electrical production of the system and its performance, there is a choice between 8 models for the calculation of the diffuse component of the solar irradiation and 5 models (Liu and Jordan, Hay and Davies, Klucher, Perez et al., and Reindl) [24,31,[39][40][41][42] for the estimation of irradiance received on the tilted plane. A single annual albedo value may be included or modified for each month of the year. Both annual and monthly soiling losses can also be included. The software can be used to select the temporal resolution of the simulation: hourly (faster) or each minute (greater accurate).
The software also has an extensive database of PV modules, inverters, storage systems, and other devices. New devices can be defined for inclusion in the project. Regarding the PV modules, PV*SOL allows the configuration of different types of installation, i.e., fixed, with a single-axis tracking system, and with a dual-axis tracking system. The inverter configuration includes the definition of the interconnection of the PV modules. Although not mandatory, it is also possible to consider wiring losses.
Once all the elements of the installation are defined, the program checks that the system is correct, and the simulation can proceed. The output data includes the predicted electrical production of the system, its performance ratio, the irradiance received on the module's plane and its temperature, a general energy balance of the system, and a financial analysis of the installation. Finally, it is possible to generate a report detailing the configuration of the system, its electrical diagrams, and the results. PV*SOL is subscription software and has a 30-day trial version.

The Facilities
This case study involves three facilities (System 1, System 2, and System 3) located slightly to the north from the center of the Spanish autonomous region of Castile and Leon, at the city of Palencia (System 1), at Torquemada (Palencia) (System 2), and at Los Balbases (Burgos) (System 3). The map locations are shown in Figure 1. The facilities are subject to very similar environmental conditions, in terms of temperature, solar irradiation, humidity, and wind speed. Solar irradiation is estimated at approximately 1450 kWh/m 2 year [9]. The ambient temperature range is between 4 • C and 20 • C and the number of cloudy days is very low [43]. Figures 2 and 3 and Figure 4 show photographs of the three facilities.
As can be seen in Figures 2 and 3 and Figure 4, all the facilities in this study have different types of sun-tracking technology. System 1, located at Palencia City, is a 38.88 kW p peak power fixed roof solar facility. It is composed of 218 PV panels (model FOTONA 180D) tilted at 20 • and has an electrical performance of 12%. System 1 uses six 7 kW p model SMC 7000TL inverters.
When simulating the electrical production of the system and its performance, there is a choice between 8 models for the calculation of the diffuse component of the solar irradiation and 5 models (Liu and Jordan, Hay and Davies, Klucher, Perez et al., and Reindl) [24,31,[39][40][41][42] for the estimation of irradiance received on the tilted plane. A single annual albedo value may be included or modified for each month of the year. Both annual and monthly soiling losses can also be included. The software can be used to select the temporal resolution of the simulation: hourly (faster) or each minute (greater accurate).
The software also has an extensive database of PV modules, inverters, storage systems, and other devices. New devices can be defined for inclusion in the project. Regarding the PV modules, PV*SOL allows the configuration of different types of installation, i.e., fixed, with a single-axis tracking system, and with a dual-axis tracking system. The inverter configuration includes the definition of the interconnection of the PV modules.
Although not mandatory, it is also possible to consider wiring losses.
Once all the elements of the installation are defined, the program checks that the system is correct, and the simulation can proceed. The output data includes the predicted electrical production of the system, its performance ratio, the irradiance received on the module's plane and its temperature, a general energy balance of the system, and a financial analysis of the installation. Finally, it is possible to generate a report detailing the configuration of the system, its electrical diagrams, and the results. PV*SOL is subscription software and has a 30-day trial version.

The Facilities
This case study involves three facilities (System 1, System 2, and System 3) located slightly to the north from the center of the Spanish autonomous region of Castile and Leon, at the city of Palencia (System 1), at Torquemada (Palencia) (System 2), and at Los Balbases (Burgos) (System 3). The map locations are shown in Figure 1. The facilities are subject to very similar environmental conditions, in terms of temperature, solar irradiation, humidity, and wind speed. Solar irradiation is estimated at approximately 1450 kWh/m 2 year [9]. The ambient temperature range is between 4 °C and 20 °C and the number of cloudy days is very low [43]. Figures 2-4 show photographs of the three facilities.      As can be seen in Figures 2-4, all the facilities in this study have different types of sun-tracking technology. System 1, located at Palencia City, is a 38.88 kW peak power fixed roof solar facility. It is composed of 218 PV panels (model FOTONA 180D) tilted at 20° and has an electrical performance of 12%. System 1 uses six 7 kW model SMC 7000TL inverters. System 2, located at Torquemada (Palencia), is a 101.01 kW power plant on the ground. It has a support structure, with which the panels can be manually adjusted in steps of 5° to inclinations between 15º and 55º thereby optimizing the angle of incidence throughout the year, approximately once per month with the optimal tile angle for each month. System 2 is composed of 546 PV panels (model BP-7185S), 15% electrical performance. Only one 100 kW inverter was selected, Ingecom Sun 100 Model. System 3, located at Los Balbases (Burgos), is a 2-axis sun-tracker facility with 51 modules and an electrical performance of 14% that can generate 14.28 kW (model Sunport 72M). It has three 4.4 kW Sunny Boy 4200TL HC inverters. The technical specifications of each of the three systems are summarized in Table 2.

Real Exported Energy Data
As previously mentioned, the real values of the energy exported to the grid were used as a reference. All the systems were set between 2006 and 2008, therefore the data period used for the study was from 2008 to 2020. Due to the length of the analysis, the    As can be seen in Figures 2-4, all the facilities in this study have different types of sun-tracking technology. System 1, located at Palencia City, is a 38.88 kW peak power fixed roof solar facility. It is composed of 218 PV panels (model FOTONA 180D) tilted at 20° and has an electrical performance of 12%. System 1 uses six 7 kW model SMC 7000TL inverters. System 2, located at Torquemada (Palencia), is a 101.01 kW power plant on the ground. It has a support structure, with which the panels can be manually adjusted in steps of 5° to inclinations between 15º and 55º thereby optimizing the angle of incidence throughout the year, approximately once per month with the optimal tile angle for each month. System 2 is composed of 546 PV panels (model BP-7185S), 15% electrical performance. Only one 100 kW inverter was selected, Ingecom Sun 100 Model. System 3, located at Los Balbases (Burgos), is a 2-axis sun-tracker facility with 51 modules and an electrical performance of 14% that can generate 14.28 kW (model Sunport 72M). It has three 4.4 kW Sunny Boy 4200TL HC inverters. The technical specifications of each of the three systems are summarized in Table 2.

Real Exported Energy Data
As previously mentioned, the real values of the energy exported to the grid were used as a reference. All the systems were set between 2006 and 2008, therefore the data period used for the study was from 2008 to 2020. Due to the length of the analysis, the System 2, located at Torquemada (Palencia), is a 101.01 kW p power plant on the ground. It has a support structure, with which the panels can be manually adjusted in steps of 5 • to inclinations between 15 • and 55 • thereby optimizing the angle of incidence throughout the year, approximately once per month with the optimal tile angle for each month. System 2 is composed of 546 PV panels (model BP-7185S), 15% electrical performance. Only one 100 kW inverter was selected, Ingecom Sun 100 Model. System 3, located at Los Balbases (Burgos), is a 2-axis sun-tracker facility with 51 modules and an electrical performance of 14% that can generate 14.28 kW p (model Sunport 72M). It has three 4.4 kW p Sunny Boy 4200TL HC inverters. The technical specifications of each of the three systems are summarized in Table 2.

Real Exported Energy Data
As previously mentioned, the real values of the energy exported to the grid were used as a reference. All the systems were set between 2006 and 2008, therefore the data period used for the study was from 2008 to 2020. Due to the length of the analysis, the influence of warm or cold years is not relevant. For that reason, average values can be considered as a typical average year of each power plant.
For each system, daily AC energy exported to the grid is available. However, as only monthly production may be modelled with most programs, the monthly average represented in Figure 5 was calculated for each installation. As mentioned above, the facilities have different nominal power capacities, therefore, to homogenize the data, all the values were divided by the installed power capacity of each system.
The energy generated is higher in systems that have some type of solar tracking. However, in the case of the two-axis tracker (System 3), it is penalized during cold months due to windy days. These elements go into safe mode when the wind is high, reducing the amount of energy that is generated. However, during summer months, this fact is not so frequent achieving considerably higher generation levels. It should also be noted that tracking facilities require more maintenance and have higher failure rates. Therefore, longterm production differences are never as high as may be expected. In Table 3, total energy generated in a year and the total production for each installed kW p are both included. considered as a typical average year of each power plant.
For each system, daily AC energy exported to the grid is available. However, as only monthly production may be modelled with most programs, the monthly average represented in Figure 5 was calculated for each installation. As mentioned above, the facilities have different nominal power capacities, therefore, to homogenize the data, all the values were divided by the installed power capacity of each system. The energy generated is higher in systems that have some type of solar tracking. However, in the case of the two-axis tracker (System 3), it is penalized during cold months due to windy days. These elements go into safe mode when the wind is high, reducing the amount of energy that is generated. However, during summer months, this fact is not so frequent achieving considerably higher generation levels. It should also be noted that tracking facilities require more maintenance and have higher failure rates. Therefore, long-term production differences are never as high as may be expected. In Table 3, total energy generated in a year and the total production for each installed kW are both included.

Results and Discussion
The results of the different models under analysis were compared with the real data from the facilities. The comparison includes results obtained from all climate models available in SAM (isotropic, Perez, HDKR) and all climate databases included in PVGIS and PVSYST. PV*SOL was set with the Perez radiation model, since it is optimal in this study. In total, 11 studies were performed for each PV installation.
One of the most widely used values when assessing the suitability of a PV project, is the amount of electrical energy that can be generated throughout a year. Table 4

Results and Discussion
The results of the different models under analysis were compared with the real data from the facilities. The comparison includes results obtained from all climate models available in SAM (isotropic, Perez, HDKR) and all climate databases included in PVGIS and PVSYST. PV*SOL was set with the Perez radiation model, since it is optimal in this study. In total, 11 studies were performed for each PV installation.
One of the most widely used values when assessing the suitability of a PV project, is the amount of electrical energy that can be generated throughout a year. Table 4 shows the total annual results obtained from the different models, in addition to the percentage deviation from the actual value for each system.
The results showed large differences depending on the type of sun-tracking technology. In general, the deviation obtained for System 1 (fixed) was low. The simplicity of this installation, in addition to the low operating uncertainty, meant that the simulations could be adjusted with deviations of below 5%. Except for PV*SOL and RETScreen, the results obtained were slightly conservative, estimating lower energy productions.
In System 2 (1 axis), the results showed greater deviation. This installation incorporates a manual tilt adjustment that increases uncertainty due to system operation errors. However, despite this factor, the deviation was still lower than 10%.
Finally, the most complex system had the 2-axis solar tracker (System 3). This type of facility is highly conditioned by the wind velocity, which, together with the greater probability of failures in these devices, reduces the availability of the facility. Therefore, the results obtained by the models considerably overestimated the energy generated, even above 20%.
When the program results were individually analyzed, the results of the free software were generally better. In general, these models are simpler, benefiting from the type of preliminary study carried out (with a low number of predefined parameters). Although PVGIS is the simplest software, its results were excellent when using the COSMO database. The results provided by SAM showed the best global fit of the three systems under analysis. Specifically, regardless of the irradiation model in use, the deviation was less than 6%. On the contrary, PV*SOL showed very high deviation values, overestimating by around 10% for all installations, but this deviation was consistent for the three systems under analysis.

Monthly Results
In addition, a comparison with monthly results was also performed. In this way, the results could be analyzed in more detail, so as to establish whether the deviation was balanced throughout the year. The deviation percentage from the actual monthly values are shown in Table 5 for the three systems. So that the interpretation of the large volume of data is easier, the cells in this table are color-coded. Red cells show results that were overestimated, blue cells show the underestimated ones.
An analysis of the monthly results showed that most models overestimated the energy generated during the winter months (when daylight hours and solar irradiation are lower), while they underestimated the results during summer, compensating the deviation throughout the year.
The model for fixed installations (System 1), as previously mentioned, showed very good results. The models that showed better results were SAM isotropic and PVSYST new. However, the results were consistent for most months. From November to January the results showed high deviations.
Despite the greater complexity of System 2 and therefore its higher operational uncertainty, due to possible orientation mismatch, in general, all the programs yielded accurate predictions in the warmer months. On the other hand, errors, even higher than 60%, were obtained during the winter months. In this case, the model with the best results was SAM. It can therefore be deduced that none of the climate databases were accurate in wintertime.
Lastly, the results for System 3 were shown to be inconsistent. In this type of installation, the models provided very different values throughout the year. In May, the deviation increased in all the models. Spring is usually windy; therefore, the solar tracker automatically activated save mode and was not correctly oriented. RETScreen and PVGIS COSMO provided the best results during the different months of the year.

Conclusions
The values of the energy supplied to the grid from three photovoltaic power plants with different sun-tracking systems have been compared, over a period of 11 years, using the results from five computer software programs with which the energy production of photovoltaic systems can be modeled. In addition, the use of the different climate databases available by default in the programs have been included in the study. Thus, eleven different simulations were available. The computer tools under analysis included subscription and open access (free) versions.
The results of the simulations carried out using the subscription tools were not appreciably better than those obtained through free-to-use programs. However, these commercial programs present a more intuitive and simpler interface, in addition to having more pos-sibilities of incorporating different formats of climatic bases, and extensive and updated elements libraries (PV panels, inverters, etc.). One of the main advantages of these systems is the possibility of modeling the installation in detail, considering many aspects, and with very good results.
The most influential parameter in the simulation results was the climatic database for the area. PV*SOL presented high but very consistent errors, above 10%. The model may therefore be calibrated to adjust to a set resolution.
The higher failure rate and safe position of the solar trackers reduce the amount of energy that is generated when there is a lot of wind. This aspect, which is difficult to predict, means that the results of all the simulations showed large differences with respect to the real production data. A model should be included with a coefficient of losses due to wind and system failures.
In addition, System 2 has a manual tilt change mechanism, which increases the losses due to errors with the adjustment of the tilt angle. Prediction of this aspect is very difficult in a global model and reduces the accuracy of these PV models.
In annual global terms, the results can occasionally hide bad predictions throughout the year, where overestimations are compensated by underestimated results. This fact was reflected in the monthly results of the different models that were not consistent. In general, during the colder months they were overestimated, while during the warmer months the models showed a better fit.
The best global results were from SAM and RETScreen. The deviations of both programs with respect to real data were lower than 10%.