pvlib 0.11.0a1__tar.gz → 0.11.2__tar.gz

{pvlib-0.11.0a1/pvlib.egg-info → pvlib-0.11.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: pvlib
-Version: 0.11.0a1
+Version: 0.11.2
 Summary: A set of functions and classes for simulating the performance of photovoltaic energy systems.
 Home-page: https://github.com/pvlib/pvlib-python
 Author-email: pvlib python Developers <pvlib-admin@googlegroups.com>
@@ -15,11 +15,11 @@ Classifier: Intended Audience :: Science/Research
 Classifier: Programming Language :: Python
 Classifier: Programming Language :: Python :: 3
 Classifier: Topic :: Scientific/Engineering
-Requires-Python: >=3.8
+Requires-Python: >=3.9
 Description-Content-Type: text/x-rst
 License-File: LICENSE
 License-File: AUTHORS.md
-Requires-Dist: numpy>=1.17.3
+Requires-Dist: numpy>=1.19.3
 Requires-Dist: pandas>=1.3.0
 Requires-Dist: pytz
 Requires-Dist: requests
@@ -34,14 +34,17 @@ Requires-Dist: solarfactors; extra == "optional"
 Requires-Dist: statsmodels; extra == "optional"
 Provides-Extra: doc
 Requires-Dist: ipython; extra == "doc"
+Requires-Dist: pickleshare; extra == "doc"
 Requires-Dist: matplotlib; extra == "doc"
-Requires-Dist: sphinx==4.5.0; extra == "doc"
-Requires-Dist: pydata-sphinx-theme==0.8.1; extra == "doc"
+Requires-Dist: sphinx==7.3.7; extra == "doc"
+Requires-Dist: pydata-sphinx-theme==0.15.4; extra == "doc"
 Requires-Dist: sphinx-gallery; extra == "doc"
-Requires-Dist: docutils==0.15.2; extra == "doc"
+Requires-Dist: docutils==0.21; extra == "doc"
 Requires-Dist: pillow; extra == "doc"
-Requires-Dist: sphinx-toggleprompt>=0.0.5; extra == "doc"
+Requires-Dist: sphinx-toggleprompt==0.5.2; extra == "doc"
+Requires-Dist: sphinx-favicon; extra == "doc"
 Requires-Dist: solarfactors; extra == "doc"
+Requires-Dist: sphinx-hoverxref; extra == "doc"
 Provides-Extra: test
 Requires-Dist: pytest; extra == "test"
 Requires-Dist: pytest-cov; extra == "test"

{pvlib-0.11.0a1 → pvlib-0.11.2}/README.md

@@ -81,7 +81,7 @@ Contributing
 ============
 
 We need your help to make pvlib-python a great tool!
-Please see the [Contributing page](http://pvlib-python.readthedocs.io/en/stable/contributing.html) for more on how you can contribute.
+Please see the [Contributing page](https://pvlib-python.readthedocs.io/en/stable/contributing/index.html) for more on how you can contribute.
 The long-term success of pvlib-python requires substantial community support.
 
 

pvlib-0.11.2/docs/examples/bifacial/plot_irradiance_nonuniformity_loss.py

@@ -0,0 +1,129 @@
+"""
+Plot Irradiance Non-uniformity Loss
+===================================
+
+Calculate the DC power lost to irradiance non-uniformity in a bifacial PV
+array.
+"""
+
+# %%
+# The incident irradiance on the backside of a bifacial PV module is
+# not uniform due to neighboring rows, the ground albedo, and site conditions.
+# When each cell works at different irradiance levels, the power produced by
+# the module is less than the sum of the power produced by each cell since the
+# maximum power point (MPP) of each cell is different, but cells connected in
+# series will operate at the same current.
+# This is known as irradiance non-uniformity loss.
+#
+# Calculating the IV curve of each cell and then matching the working point of
+# the whole module is computationally expensive, so a simple model to account
+# for this loss is of interest. Deline et al. [1]_ proposed a model based on
+# the Relative Mean Absolute Difference (RMAD) of the irradiance of each cell.
+# They considered the standard deviation of the cells' irradiances, but they
+# found that the RMAD was a better predictor of the mismatch loss.
+#
+# This example demonstrates how to model the irradiance non-uniformity loss
+# from the irradiance levels of each cell in a PV module.
+#
+# The function
+# :py:func:`pvlib.bifacial.power_mismatch_deline` is
+# used to transform the Relative Mean Absolute Difference (RMAD) of the
+# irradiance into a power loss mismatch. Down below you will find a
+# numpy-based implementation of the RMAD function.
+#
+# References
+# ----------
+# .. [1] C. Deline, S. Ayala Pelaez, S. MacAlpine, and C. Olalla, 'Estimating
+#    and parameterizing mismatch power loss in bifacial photovoltaic
+#    systems', Progress in Photovoltaics: Research and Applications, vol. 28,
+#    no. 7, pp. 691-703, 2020, :doi:`10.1002/pip.3259`.
+#
+# .. sectionauthor:: Echedey Luis <echelual (at) gmail.com>
+
+# %%
+# Problem description
+# -------------------
+# Let's set a fixed irradiance to each cell row of the PV array with the values
+# described in Figure 1 (A), [1]_. We will cover this case for educational
+# purposes, although it can be achieved with the packages
+# `solarfactors <https://github.com/pvlib/solarfactors/>`_ and
+# `bifacial_radiance <https://github.com/NREL/bifacial_radiance>`_.
+#
+# Here we set and plot the global irradiance level of each cell.
+
+import numpy as np
+import matplotlib.pyplot as plt
+from matplotlib.cm import ScalarMappable
+from matplotlib.colors import Normalize
+
+from pvlib.bifacial import power_mismatch_deline
+
+x = np.arange(12, 0, -1)
+y = np.arange(6, 0, -1)
+cells_irrad = np.repeat([1059, 976, 967, 986, 1034, 1128], len(x)).reshape(
+    len(y), len(x)
+)
+
+# plot the irradiance levels of each cell
+color_map = "gray"
+color_norm = Normalize(930, 1150)
+
+fig, ax = plt.subplots()
+fig.suptitle("Global Irradiance Levels of Each Cell")
+fig.colorbar(
+    ScalarMappable(cmap=color_map, norm=color_norm),
+    ax=ax,
+    orientation="vertical",
+    label="$[W/m^2]$",
+)
+ax.set_aspect("equal")
+ax.pcolormesh(
+    x,
+    y,
+    cells_irrad,
+    shading="nearest",
+    edgecolors="black",
+    cmap=color_map,
+    norm=color_norm,
+)
+
+# %%
+# Relative Mean Absolute Difference
+# ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
+# Calculate the Relative Mean Absolute Difference (RMAD) of the cells'
+# irradiances with the following function, Eq. (4) of [1]_:
+#
+# .. math::
+#
+#    \Delta \left[ unitless \right] = \frac{1}{n^2 \bar{G}_{total}}
+#    \sum_{i=1}^{n} \sum_{j=1}^{n} \lvert G_{total,i} - G_{total,j} \rvert
+#
+
+
+def rmad(data):
+    """
+    Relative Mean Absolute Difference. Output is [Unitless]. Eq. (4) of [1]_.
+    """
+    mean = np.mean(data)
+    mad = np.mean(np.absolute(np.subtract.outer(data, data)))
+    return mad / mean
+
+
+rmad_cells = rmad(cells_irrad)
+
+# this is the same as a column's RMAD!
+print(rmad_cells == rmad(cells_irrad[:, 0]))
+
+# %%
+# Mismatch Loss
+# ^^^^^^^^^^^^^
+# Calculate the power loss ratio due to the irradiance non-uniformity
+# with :py:func:`pvlib.bifacial.power_mismatch_deline`.
+
+mismatch_loss = power_mismatch_deline(rmad_cells)
+
+print(f"RMAD of the cells' irradiance: {rmad_cells:.3} [unitless]")
+print(
+    "Power loss due to the irradiance non-uniformity: "
+    + f"{mismatch_loss:.3} [unitless]"
+)

pvlib-0.11.2/docs/examples/floating-pv/README.rst

@@ -0,0 +1,2 @@
+Floating PV Systems Modelling
+-----------------------------

pvlib-0.11.2/docs/examples/floating-pv/plot_floating_pv_cell_temperature.py

@@ -0,0 +1,230 @@
+r"""
+Temperature modeling for floating PV
+================================================
+
+This example uses the PVSyst temperature model to calculate
+cell temperature for floating photovoltaic (FPV) systems.
+
+One of the primary benefits attributed to FPV systems is the potential
+for lower operating temperatures, which are expected to increase the
+operating efficiency. In general, the temperature at which a photovoltaic
+module operates is influenced by various factors including solar radiation,
+ambient temperature, wind speed and direction, and the characteristics of the
+cell and module materials, as well as the mounting structure.
+
+A popular model for calculating PV cell temperature is the
+empirical heat loss factor model suggested by Faiman
+(:py:func:`pvlib.temperature.faiman`). A modified version of this model is
+has been proposed by PVSyst (:py:func:`~pvlib.temperature.pvsyst_cell`).
+The PVSyst model for cell temperature :math:`T_{C}` is given by:
+
+.. math::
+
+    T_{C} = T_{a} + \frac{\alpha \cdot E \cdot (1 - \eta_{m})}{U_{c} + U_{v} \cdot WS},
+
+where :math:`E` is the plane-of-array irradiance, :math:`T_{a}` is the
+ambient air temperature, :math:`WS` is the wind speed, :math:`\alpha` is the
+absorbed fraction of the incident irradiance, :math:`\eta_{m}` is the
+electrical efficiency of the module, :math:`U_{c}` is the wind-independent heat
+loss coefficient, and :math:`U_{v}` is the wind-dependent heat loss coefficient.
+It should be noted that in many cases, similar to land-based PV systems,
+the wind-dependent heat loss coefficient (:math:`U_{v}`) can be set to zero,
+and the denominator is thus reduced to a single combined U-value
+(:math:`U_{c}`).
+
+However, the default heat loss coefficient values of the PVSyst model were
+specified for land-based PV systems and are not necessarily representative
+of FPV systems.
+
+For FPV systems, the module's operating temperature, much like in land-based
+systems, is mainly influenced by the mounting structure (which significantly
+affects both U-value coefficients), wind, and air temperature. Thus, factors
+that help reduce operating temperatures in such systems include lower air
+temperatures and changes in air flow beneath the modules (wind/convection).
+In some designs, where the modules are in direct thermal contact with water,
+cooling effectiveness is largely dictated by the water temperature.
+
+Systems with closely packed modules and restricted airflow behind the modules
+generally exhibit lower heat loss coefficients compared to those with better
+airflow behind the modules.
+
+The table below gives heat loss coefficients derived for different FPV systems
+and locations as found in the literature. It should be noted that, for some
+systems, there are two sets of coefficients, where the second set uses only
+one heat loss coefficient (i.e., only :math:`U_{c}`).
+
+.. table:: Heat transfer coefficients for different PV systems
+   :widths: 40 15 15 15 15
+
+   +--------------------------+-------------+--------------------------------+----------------------------------------+-----------+
+   | System                   | Location    |:math:`U_{c}`                   | :math:`U_{v}`                          | Reference |
+   |                          |             |:math:`[\frac{W}{m^2 \cdot K}]` | :math:`[\frac{W}{m^3 \cdot K \cdot s}]`|           |
+   +==========================+=============+================================+========================================+===========+
+   | - Monofacial module      | Netherlands | 24.4                           | 6.5                                    | [1]_      |
+   | - Open structure         |             |                                |                                        |           |
+   | - Two-axis tracking      |             | 57                             | 0                                      |           |
+   | - Small water footprint  |             |                                |                                        |           |
+   +--------------------------+-------------+--------------------------------+----------------------------------------+-----------+
+   | - Monofacial module      | Netherlands | 25.2                           | 3.7                                    | [1]_      |
+   | - Closed structure       |             |                                |                                        |           |
+   | - Large water footprint  |             | 37                             | 0                                      |           |
+   +--------------------------+-------------+--------------------------------+----------------------------------------+-----------+
+   | - Monofacial module      | Singapore   | 34.8                           | 0.8                                    | [1]_      |
+   | - Closed structure       |             |                                |                                        |           |
+   | - Large water footprint  |             | 36                             | 0                                      |           |
+   +--------------------------+-------------+--------------------------------+----------------------------------------+-----------+
+   | - Monofacial module      | Singapore   | 18.9                           | 8.9                                    | [1]_      |
+   | - Closed structure       |             |                                |                                        |           |
+   | - Medium water footprint |             | 41                             | 0                                      |           |
+   +--------------------------+-------------+--------------------------------+----------------------------------------+-----------+
+   | - Monofacial module      | Singapore   | 35.3                           | 8.9                                    | [1]_      |
+   | - Open structure         |             |                                |                                        |           |
+   | - Free-standing          |             | 55                             | 0                                      |           |
+   +--------------------------+-------------+--------------------------------+----------------------------------------+-----------+
+   | - Monofacial module      | Norway      | 71                             | 0                                      | [2]_      |
+   | - In contact with water  |             |                                |                                        |           |
+   | - Calculated using water |             |                                |                                        |           |
+   |   temperature as         |             |                                |                                        |           |
+   |   :math:`T_{amb}`        |             |                                |                                        |           |
+   +--------------------------+-------------+--------------------------------+----------------------------------------+-----------+
+   | - Monofacial module      | South Italy | 31.9                           | 1.5                                    | [3]_      |
+   | - Open structure         |             |                                |                                        |           |
+   | - Free-standing          |             |                                |                                        |           |
+   +--------------------------+-------------+--------------------------------+----------------------------------------+-----------+
+   | - Bifacial module        | South Italy | 35.2                           | 1.5                                    | [3]_      |
+   | - Open structure         |             |                                |                                        |           |
+   | - Free-standing          |             |                                |                                        |           |
+   +--------------------------+-------------+--------------------------------+----------------------------------------+-----------+
+
+References
+----------
+.. [1] Dörenkämper M., Wahed A., Kumar A., de Jong M., Kroon J., Reindl T.
+    (2021), 'The cooling effect of floating PV in two different climate zones:
+    A comparison of field test data from the Netherlands and Singapore'
+    Solar Energy, vol. 214, pp. 239-247, :doi:`10.1016/j.solener.2020.11.029`.
+
+.. [2] Kjeldstad T., Lindholm D., Marstein E., Selj J. (2021), 'Cooling of
+    floating photovoltaics and the importance of water temperature', Solar
+    Energy, vol. 218, pp. 544-551, :doi:`10.1016/j.solener.2021.03.022`.
+
+.. [3] Tina G.M., Scavo F.B., Merlo L., Bizzarri F. (2021), 'Comparative
+    analysis of monofacial and bifacial photovoltaic modules for floating
+    power plants', Applied Energy, vol 281, pp. 116084,
+    :doi:`10.1016/j.apenergy.2020.116084`.
+"""  # noqa: E501
+
+# %%
+# Read example weather data
+# ^^^^^^^^^^^^^^^^^^^^^^^^^
+# Read weather data from a TMY3 file and calculate the solar position and
+# the plane-of-array irradiance.
+
+import pvlib
+import matplotlib.pyplot as plt
+from pathlib import Path
+
+# Assume a FPV system on a lake with the following specifications
+tilt = 30  # degrees
+azimuth = 180  # south-facing
+
+# Datafile found in the pvlib distribution
+data_file = Path(pvlib.__path__[0]).joinpath('data', '723170TYA.CSV')
+
+tmy, metadata = pvlib.iotools.read_tmy3(
+    data_file, coerce_year=2002, map_variables=True
+)
+tmy = tmy.filter(
+    ['ghi', 'dni', 'dni_extra', 'dhi', 'temp_air', 'wind_speed', 'pressure']
+)  # remaining columns are not needed
+tmy = tmy['2002-06-06 00:00':'2002-06-06 23:59']  # select period
+
+solar_position = pvlib.solarposition.get_solarposition(
+    # TMY timestamp is at end of hour, so shift to center of interval
+    tmy.index.shift(freq='-30T'),
+    latitude=metadata['latitude'],
+    longitude=metadata['longitude'],
+    altitude=metadata['altitude'],
+    pressure=tmy['pressure'] * 100,  # convert from millibar to Pa
+    temperature=tmy['temp_air'],
+)
+solar_position.index = tmy.index  # reset index to end of the hour
+
+# Albedo calculation for inland water bodies
+albedo = pvlib.albedo.inland_water_dvoracek(
+    solar_elevation=solar_position['elevation'],
+    surface_condition='clear_water_no_waves'
+)
+
+# Use transposition model to find plane-of-array irradiance
+irradiance = pvlib.irradiance.get_total_irradiance(
+    surface_tilt=tilt,
+    surface_azimuth=azimuth,
+    solar_zenith=solar_position['apparent_zenith'],
+    solar_azimuth=solar_position['azimuth'],
+    dni=tmy['dni'],
+    dni_extra=tmy['dni_extra'],
+    ghi=tmy['ghi'],
+    dhi=tmy['dhi'],
+    albedo=albedo,
+    model='haydavies'
+)
+
+# %%
+# Calculate and plot cell temperature
+# ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
+# The temperature of the PV cell is calculated for lake-based floating PV
+# systems.
+
+# Make a dictionary containing all the sets of coefficients presented in the
+# above table.
+heat_loss_coeffs = {
+    'open_structure_small_footprint_tracking_NL': [24.4, 6.5, 'C0', 'solid'],
+    'open_structure_small_footprint_tracking_NL_2': [57, 0, 'C0', 'dashed'],
+    'closed_structure_large_footprint_NL': [25.2, 3.7, 'C1', 'solid'],
+    'closed_structure_large_footprint_NL_2': [37, 0, 'C1', 'dashed'],
+    'closed_structure_large_footprint_SG': [34.8, 0.8, 'C2', 'solid'],
+    'closed_structure_large_footprint_SG_2': [36, 0, 'C2', 'dashed'],
+    'closed_structure_medium_footprint_SG': [18.9, 8.9, 'C3', 'solid'],
+    'closed_structure_medium_footprint_SG_2': [41, 0, 'C3', 'dashed'],
+    'open_structure_free_standing_SG': [35.3, 8.9, 'C4', 'solid'],
+    'open_structure_free_standing_SG_2': [55, 0, 'C4', 'dashed'],
+    'in_contact_with_water_NO': [71, 0, 'C5', 'solid'],
+    'open_structure_free_standing_IT': [31.9, 1.5, 'C6', 'solid'],
+    'open_structure_free_standing_bifacial_IT': [35.2, 1.5, 'C7', 'solid'],
+    'default_PVSyst_coeffs_for_land_systems': [29.0, 0, 'C8', 'solid']
+}
+
+# Plot the cell temperature for each set of the above heat loss coefficients
+for coeffs in heat_loss_coeffs:
+    T_cell = pvlib.temperature.pvsyst_cell(
+        poa_global=irradiance['poa_global'],
+        temp_air=tmy['temp_air'],
+        wind_speed=tmy['wind_speed'],
+        u_c=heat_loss_coeffs[coeffs][0],
+        u_v=heat_loss_coeffs[coeffs][1]
+    )
+    # Convert Dataframe Indexes to Hour format to make plotting easier
+    T_cell.index = T_cell.index.strftime("%H")
+    plt.plot(T_cell, label=coeffs, c=heat_loss_coeffs[coeffs][2],
+             ls=heat_loss_coeffs[coeffs][3], alpha=0.8)
+
+plt.xlabel('Hour')
+plt.ylabel('PV cell temperature [°C]')
+plt.ylim(10, 45)
+plt.xlim('06', '20')
+plt.grid()
+plt.legend(loc='upper left', frameon=False, ncols=2, fontsize='x-small',
+           bbox_to_anchor=(0, -0.2))
+plt.tight_layout()
+plt.show()
+
+# %%
+# The figure above illustrates the necessity of choosing appropriate heat loss
+# coefficients when using the PVSyst model for calculating the cell temperature
+# for floating PV systems. A difference of up to 10.3 °C was obtained when
+# using the default PVSyst coefficients versus using coefficients for systems
+# where panels are in contact with water.
+#
+# It should be noted that, using the single combined U-value versus the
+# :math:`U_c` and :math:`U_v` gives significantly different results, even
+# when using the coefficients derived from the same system.

{pvlib-0.11.0a1 → pvlib-0.11.2}/docs/examples/irradiance-transposition/plot_seasonal_tilt.py

@@ -32,6 +32,8 @@ class SeasonalTiltMount(pvsystem.AbstractMount):
     surface_azimuth: float = 180.0
 
     def get_orientation(self, solar_zenith, solar_azimuth):
+        # note: determining tilt based on month may produce different
+        # results depending on the time zone of the timestamps
         tilts = [self.monthly_tilts[m-1] for m in solar_zenith.index.month]
         return pd.DataFrame({
             'surface_tilt': tilts,

{pvlib-0.11.0a1 → pvlib-0.11.2}/docs/examples/irradiance-transposition/plot_transposition_gain.py

@@ -104,6 +104,7 @@ df_monthly = 100 * (df_monthly.divide(ghi_monthly, axis=0) - 1)
 df_monthly.plot()
 plt.xlabel('Month of Year')
 plt.ylabel('Monthly Transposition Gain [%]')
+plt.tight_layout()
 plt.show()
 
 

pvlib-0.11.2/docs/examples/irradiance-transposition/use_perez_modelchain.py

@@ -0,0 +1,126 @@
+"""
+Use different Perez coefficients with the ModelChain
+====================================================
+
+This example demonstrates how to customize the ModelChain
+to use site-specific Perez transposition coefficients.
+"""
+
+# %%
+# The :py:class:`pvlib.modelchain.ModelChain` object provides a useful method
+# for easily constructing a PV system model with a simple, unified interface.
+# However, a user may want to customize the steps
+# in the system model in various ways.
+# One such example is during the irradiance transposition step.
+# The Perez model perform very well on field data, but
+# it requires a set of fitted coefficients from various sites.
+# It has been noted that these coefficients can be specific to
+# various climates, so users may see improved model accuracy
+# when using a site-specific set of coefficients.
+# However, the base  :py:class:`~pvlib.modelchain.ModelChain`
+# only supports the default coefficients.
+# This example shows how the  :py:class:`~pvlib.modelchain.ModelChain` can
+# be adjusted to use a different set of Perez coefficients.
+
+import pandas as pd
+from pvlib.pvsystem import PVSystem
+from pvlib.modelchain import ModelChain
+from pvlib.temperature import TEMPERATURE_MODEL_PARAMETERS
+from pvlib import iotools, location, irradiance
+import pvlib
+import os
+import matplotlib.pyplot as plt
+
+# load in TMY weather data from North Carolina included with pvlib
+PVLIB_DIR = pvlib.__path__[0]
+DATA_FILE = os.path.join(PVLIB_DIR, 'data', '723170TYA.CSV')
+
+tmy, metadata = iotools.read_tmy3(DATA_FILE, coerce_year=1990,
+                                  map_variables=True)
+
+weather_data = tmy[['ghi', 'dhi', 'dni', 'temp_air', 'wind_speed']]
+
+loc = location.Location.from_tmy(metadata)
+
+#%%
+# Now, let's set up a standard PV model using the ``ModelChain``
+
+surface_tilt = metadata['latitude']
+surface_azimuth = 180
+
+# define an example module and inverter
+sandia_modules = pvlib.pvsystem.retrieve_sam('SandiaMod')
+cec_inverters = pvlib.pvsystem.retrieve_sam('cecinverter')
+sandia_module = sandia_modules['Canadian_Solar_CS5P_220M___2009_']
+cec_inverter = cec_inverters['ABB__MICRO_0_25_I_OUTD_US_208__208V_']
+
+temp_params = TEMPERATURE_MODEL_PARAMETERS['sapm']['open_rack_glass_glass']
+
+# define the system and ModelChain
+system = PVSystem(arrays=None,
+                  surface_tilt=surface_tilt,
+                  surface_azimuth=surface_azimuth,
+                  module_parameters=sandia_module,
+                  inverter_parameters=cec_inverter,
+                  temperature_model_parameters=temp_params)
+
+mc = ModelChain(system, location=loc)
+
+# %%
+# Now, let's calculate POA irradiance values outside of the ``ModelChain``.
+# We do this for both the default Perez coefficients and the desired
+# alternative Perez coefficients.  This enables comparison at the end.
+
+# Cape Canaveral seems like the most likely match for climate
+model_perez = 'capecanaveral1988'
+
+solar_position = loc.get_solarposition(times=weather_data.index)
+dni_extra = irradiance.get_extra_radiation(weather_data.index)
+
+POA_irradiance = irradiance.get_total_irradiance(
+    surface_tilt=surface_tilt,
+    surface_azimuth=surface_azimuth,
+    dni=weather_data['dni'],
+    ghi=weather_data['ghi'],
+    dhi=weather_data['dhi'],
+    solar_zenith=solar_position['apparent_zenith'],
+    solar_azimuth=solar_position['azimuth'],
+    model='perez',
+    dni_extra=dni_extra)
+
+POA_irradiance_new_perez = irradiance.get_total_irradiance(
+    surface_tilt=surface_tilt,
+    surface_azimuth=surface_azimuth,
+    dni=weather_data['dni'],
+    ghi=weather_data['ghi'],
+    dhi=weather_data['dhi'],
+    solar_zenith=solar_position['apparent_zenith'],
+    solar_azimuth=solar_position['azimuth'],
+    model='perez',
+    model_perez=model_perez,
+    dni_extra=dni_extra)
+
+# %%
+# Now, run the ``ModelChain`` with both sets of irradiance data and compare
+# (note that to use POA irradiance as input to the ModelChain the method
+# `.run_model_from_poa` is used):
+
+mc.run_model_from_poa(POA_irradiance)
+ac_power_default = mc.results.ac
+
+mc.run_model_from_poa(POA_irradiance_new_perez)
+ac_power_new_perez = mc.results.ac
+
+start, stop = '1990-05-05 06:00:00', '1990-05-05 19:00:00'
+plt.plot(ac_power_default.loc[start:stop],
+         label="Default Composite Perez Model")
+plt.plot(ac_power_new_perez.loc[start:stop],
+         label="Cape Canaveral Perez Model")
+plt.xticks(rotation=90)
+plt.ylabel("AC Power ($W$)")
+plt.legend()
+plt.tight_layout()
+plt.show()
+# %%
+# Note that there is a small, but noticeable difference from the default
+# coefficients that may add up over longer periods of time.

{pvlib-0.11.0a1 → pvlib-0.11.2}/docs/examples/iv-modeling/plot_singlediode.py

@@ -114,6 +114,8 @@ for idx, case in conditions.iterrows():
     # mark the MPP
     plt.plot([v_mp], [i_mp], ls='', marker='o', c='k')
 
+plt.xlim(left=0)
+plt.ylim(bottom=0)
 plt.legend(loc=(1.0, 0))
 plt.xlabel('Module voltage [V]')
 plt.ylabel('Module current [A]')

{pvlib-0.11.0a1 → pvlib-0.11.2}/docs/examples/shading/plot_martinez_shade_loss.py

@@ -47,7 +47,7 @@ import pvlib
 import pandas as pd
 import numpy as np
 import matplotlib.pyplot as plt
-from matplotlib.dates import ConciseDateFormatter
+from matplotlib.dates import DateFormatter
 
 pitch = 4  # meters
 width = 1.5  # meters
@@ -235,10 +235,6 @@ for k, shade_losses in shade_losses_per_module.items():
     ax1.plot(times, shade_losses, label=k, linestyle=linestyle)
 ax1.legend(loc="upper center")
 ax1.grid()
-ax1.set_xlabel("Time")
-ax1.xaxis.set_major_formatter(
-    ConciseDateFormatter("%H:%M", tz="Europe/Madrid")
-)
 ax1.set_ylabel(r"$P_{out}$ losses")
 ax1.set_title("Per module")
 
@@ -248,9 +244,7 @@ for k, shade_losses in shade_losses_per_row.items():
 ax2.legend(loc="upper center")
 ax2.grid()
 ax2.set_xlabel("Time")
-ax2.xaxis.set_major_formatter(
-    ConciseDateFormatter("%H:%M", tz="Europe/Madrid")
-)
+ax2.xaxis.set_major_formatter(DateFormatter("%H:%M", tz="Europe/Madrid"))
 ax2.set_ylabel(r"$P_{out}$ losses")
 ax2.set_title("Per row")
 fig.tight_layout()