pvlib 0.10.3__tar.gz → 0.10.5__tar.gz

{pvlib-0.10.3/pvlib.egg-info → pvlib-0.10.5}/PKG-INFO

@@ -1,11 +1,10 @@
 Metadata-Version: 2.1
 Name: pvlib
-Version: 0.10.3
+Version: 0.10.5
 Summary: A set of functions and classes for simulating the performance of photovoltaic energy systems.
 Home-page: https://github.com/pvlib/pvlib-python
-Author: pvlib python Developers
-Maintainer-email: pvlib-admin@googlegroups.com
-License: BSD 3-Clause
+Author-email: pvlib python Developers <pvlib-admin@googlegroups.com>
+License: BSD-3-Clause
 Project-URL: Bug Tracker, https://github.com/pvlib/pvlib-python/issues
 Project-URL: Documentation, https://pvlib-python.readthedocs.io/
 Project-URL: Source Code, https://github.com/pvlib/pvlib-python
@@ -16,22 +15,21 @@ Classifier: Intended Audience :: Science/Research
 Classifier: Programming Language :: Python
 Classifier: Programming Language :: Python :: 3
 Classifier: Topic :: Scientific/Engineering
-Requires-Python: >=3.7
+Requires-Python: >=3.8
 Description-Content-Type: text/x-rst
 License-File: LICENSE
 License-File: AUTHORS.md
-Requires-Dist: numpy>=1.16.0
-Requires-Dist: pandas>=0.25.0
+Requires-Dist: numpy>=1.17.3
+Requires-Dist: pandas>=1.3.0
 Requires-Dist: pytz
 Requires-Dist: requests
-Requires-Dist: scipy>=1.5.0
+Requires-Dist: scipy>=1.6.0
 Requires-Dist: h5py
-Requires-Dist: importlib-metadata; python_version < "3.8"
 Provides-Extra: optional
 Requires-Dist: cython; extra == "optional"
 Requires-Dist: ephem; extra == "optional"
 Requires-Dist: nrel-pysam; extra == "optional"
-Requires-Dist: numba; extra == "optional"
+Requires-Dist: numba>=0.17.0; extra == "optional"
 Requires-Dist: solarfactors; extra == "optional"
 Requires-Dist: statsmodels; extra == "optional"
 Provides-Extra: doc
@@ -54,29 +52,7 @@ Requires-Dist: pytest-rerunfailures; extra == "test"
 Requires-Dist: pytest-remotedata; extra == "test"
 Requires-Dist: packaging; extra == "test"
 Provides-Extra: all
-Requires-Dist: cython; extra == "all"
-Requires-Dist: docutils==0.15.2; extra == "all"
-Requires-Dist: ephem; extra == "all"
-Requires-Dist: ipython; extra == "all"
-Requires-Dist: matplotlib; extra == "all"
-Requires-Dist: nrel-pysam; extra == "all"
-Requires-Dist: numba; extra == "all"
-Requires-Dist: packaging; extra == "all"
-Requires-Dist: pillow; extra == "all"
-Requires-Dist: pydata-sphinx-theme==0.8.1; extra == "all"
-Requires-Dist: pytest; extra == "all"
-Requires-Dist: pytest-cov; extra == "all"
-Requires-Dist: pytest-mock; extra == "all"
-Requires-Dist: pytest-remotedata; extra == "all"
-Requires-Dist: pytest-rerunfailures; extra == "all"
-Requires-Dist: pytest-timeout; extra == "all"
-Requires-Dist: requests-mock; extra == "all"
-Requires-Dist: solarfactors; extra == "all"
-Requires-Dist: sphinx==4.5.0; extra == "all"
-Requires-Dist: sphinx-gallery; extra == "all"
-Requires-Dist: sphinx-toggleprompt>=0.0.5; extra == "all"
-Requires-Dist: statsmodels; extra == "all"
-
+Requires-Dist: pvlib[doc,optional,test]; extra == "all"
 
 pvlib python is a community developed toolbox that provides a set of
 functions and classes for simulating the performance of photovoltaic

{pvlib-0.10.3 → pvlib-0.10.5}/README.md

@@ -7,11 +7,11 @@
     <a href="https://pypi.org/project/pvlib/">
     <img src="https://img.shields.io/pypi/v/pvlib.svg" alt="latest release" />
     </a>
-    <a href="https://anaconda.org/conda-forge/pvlib-python">
-    <img src="https://anaconda.org/conda-forge/pvlib-python/badges/version.svg" />
+    <a href="https://anaconda.org/conda-forge/pvlib">
+    <img src="https://anaconda.org/conda-forge/pvlib/badges/version.svg" />
     </a>
-    <a href="https://anaconda.org/conda-forge/pvlib-python">
-    <img src="https://anaconda.org/conda-forge/pvlib-python/badges/latest_release_date.svg" />
+    <a href="https://anaconda.org/conda-forge/pvlib">
+    <img src="https://anaconda.org/conda-forge/pvlib/badges/latest_release_date.svg" />
     </a>
 </tr>
 <tr>
@@ -50,19 +50,8 @@
     <a href="https://doi.org/10.5281/zenodo.593284">
     <img src="https://zenodo.org/badge/DOI/10.5281/zenodo.593284.svg" alt="zenodo reference">
     </a>
-    <a href="http://joss.theoj.org/papers/41187535cad22dd4b076c89b72f874b1">
-    <img src="http://joss.theoj.org/papers/41187535cad22dd4b076c89b72f874b1/status.svg" alt="JOSS reference" />
-    </a>
-  </td>
-</tr>
-<tr>
-  <td>Downloads</td>
-  <td>
-    <a href="https://pypi.org/project/pvlib/">
-    <img src="https://img.shields.io/pypi/dm/pvlib" alt="PyPI downloads" />
-    </a>
-    <a href="https://anaconda.org/conda-forge/pvlib-python">
-    <img src="https://anaconda.org/conda-forge/pvlib-python/badges/downloads.svg" alt="conda-forge downloads" />
+    <a style="border-width:0" href="https://doi.org/10.21105/joss.05994">
+    <img src="https://joss.theoj.org/papers/10.21105/joss.05994/status.svg" alt="DOI badge" >
     </a>
   </td>
 </tr>
@@ -99,17 +88,32 @@ The long-term success of pvlib-python requires substantial community support.
 Citing
 ======
 
-If you use pvlib-python in a published work, please cite:
+Many of the contributors to pvlib python work in institutions where
+citation metrics are used in performance or career evaluations. If you
+use pvlib python in a published work, please cite:
 
-  William F. Holmgren, Clifford W. Hansen, and Mark A. Mikofski.
+**Recommended citation for the pvlib python project**
+
+  Anderson, K., Hansen, C., Holmgren, W., Jensen, A., Mikofski, M., and Driesse, A.
+  "pvlib python: 2023 project update."
+  Journal of Open Source Software, 8(92), 5994, (2023).
+  https://doi.org/10.21105/joss.05994
+
+
+**Recommended citation for pvlib iotools**
+
+  Jensen, A., Anderson, K., Holmgren, W., Mikofski, M., Hansen, C., Boeman, L., Loonen, R.
+  "pvlib iotools —- Open-source Python functions for seamless access to solar irradiance data."
+  Solar Energy, 266, 112092, (2023).
+  https://doi.org/10.1016/j.solener.2023.112092
+
+**Historical citation for pvlib python**
+
+  Holmgren, W., Hansen, C., and Mikofski, M.
   "pvlib python: a python package for modeling solar energy systems."
   Journal of Open Source Software, 3(29), 884, (2018).
   https://doi.org/10.21105/joss.00884
 
-Please also cite the DOI corresponding to the specific version of
-pvlib-python that you used. pvlib-python DOIs are listed at
-[Zenodo.org](https://zenodo.org/search?page=1&size=20&q=conceptrecid:593284&all_versions&sort=-version)
-
 If you use pvlib-python in a commercial or publicly-available application, please
 consider displaying one of the "powered by pvlib" logos:
 
@@ -141,11 +145,10 @@ BSD 3-clause.
 History and acknowledgement
 ===========================
 
-pvlib python started out as a Python translation of the PVLIB MATLAB
-toolbox (developed by the `PVPMC <https://pvpmc.sandia.gov/>`_ at
-Sandia National Laboratories) in 2013 and has grown substantially since then.
+pvlib python began in 2013 as a Python translation of the [PVLIB for Matlab](https://github.com/sandialabs/MATLAB_PV_LIB)
+toolbox developed by Sandia National Laboratories. pvlib python has grown substantially since then.
 Today it contains code contributions from over a hundred individuals worldwide
-and is maintained by a core group of PV modelers from a variety institutions.
+and is maintained by a core group of PV modelers from a variety of institutions.
 
 pvlib has been supported directly and indirectly by DOE, NumFOCUS, and
 Google Summer of Code funding, university research projects,

{pvlib-0.10.3 → pvlib-0.10.5}/docs/examples/shading/plot_partial_module_shading_simple.py

@@ -317,7 +317,7 @@ for diffuse_fraction in np.linspace(0, 1, 11):
 
 results = pd.DataFrame(data)
 results['pmp'] /= results['pmp'].max()  # normalize power to 0-1
-results_pivot = results.pivot('fd', 'fs', 'pmp')
+results_pivot = results.pivot(index='fd', columns='fs', values='pmp')
 plt.figure()
 plt.imshow(results_pivot, origin='lower', aspect='auto')
 plt.xlabel('shaded fraction')

{pvlib-0.10.3 → pvlib-0.10.5}/docs/examples/solar-position/plot_sunpath_diagrams.py

@@ -103,8 +103,7 @@ import matplotlib.pyplot as plt
 
 tz = 'Asia/Calcutta'
 lat, lon = 28.6, 77.2
-times = pd.date_range('2019-01-01 00:00:00', '2020-01-01', closed='left',
-                      freq='H', tz=tz)
+times = pd.date_range('2019-01-01 00:00:00', '2020-01-01', freq='H', tz=tz)
 
 solpos = solarposition.get_solarposition(times, lat, lon)
 # remove nighttime

pvlib-0.10.5/docs/examples/system-models/README.rst

@@ -0,0 +1,2 @@
+System Models
+-------------

pvlib-0.10.5/docs/examples/system-models/plot_oedi_9068.py

@@ -0,0 +1,283 @@
+"""
+4.7 MW CdTe single-axis tracking (OEDI System 9068)
+===================================================
+
+A basic model of a 4.7 MW single-axis tracking CdTe system located in
+Colorado, United States.
+"""
+# %%
+# This example model uses satellite-based solar resource data from the
+# NSRDB PSM3. This approach is useful for pre-construction energy modeling
+# and in retrospective analyses where the system’s own irradiance
+# measurements are not present or unreliable.
+#
+# The system has public monitoring data available at the Open Energy Data
+# Initiative (OEDI) under `System ID
+# 9068 <https://data.openei.org/s3_viewer?bucket=oedi-data-lake&prefix=pvdaq%2F2023-solar-data-prize%2F9068_OEDI%2F>`__.
+# For more information about the system, see its `OEDI
+# page <https://openei.org/wiki/PVDAQ/Sites/SR_CO>`__.
+
+# sphinx_gallery_thumbnail_path = "_images/OEDI_9068_daily_timeseries.png"
+import pvlib
+import pandas as pd
+import matplotlib.pyplot as plt
+
+# %%
+# System parameters
+# -----------------
+#
+# The system description on the OEDI provides some high-level system
+# information, but unfortunately we have to make some guesses about other
+# aspects of the system’s configuration.
+#
+# The cells below define the system parameter values required in the
+# simulation.
+
+# information provided by system description on OEDI
+
+latitude = 40.3864
+longitude = -104.5512
+
+# the inverters have identical PV array topologies:
+modules_per_string = 15
+strings_per_inverter = 1344
+
+# %%
+
+# "unofficial" information
+
+# We know the system uses 117.5 W CdTe modules. Based on the system vintage
+# (data begins in 2017), it seems likely that the array uses First Solar
+# Series 4 modules (FS-4117).
+cec_module_db = pvlib.pvsystem.retrieve_sam('cecmod')
+module_parameters = cec_module_db['First_Solar__Inc__FS_4117_3']
+# ensure that correct spectral correction is applied
+module_parameters['Technology'] = 'CdTe'
+
+# default Faiman model parameters:
+temperature_model_parameters = dict(u0=25.0, u1=6.84)
+module_unit_mass = 12 / 0.72  # kg/m^2, taken from datasheet values
+
+# The OEDI metadata says the inverters have AC capacities of 1910 kW,
+# but the clipping level in the measured inverter output is more like 1840 kW.
+# It's not clear what specific model is installed, so let's just assume
+# this inverter, which the CEC database lists as having a nominal AC
+# capacity of 1833 kW:
+cec_inverter_db = pvlib.pvsystem.retrieve_sam('cecinverter')
+inverter_parameters = cec_inverter_db['TMEIC__PVL_L1833GRM']
+
+# We'll use the PVWatts v5 losses model. Set shading to zero as it is
+# accounted for elsewhere in the model, and disable availability loss since
+# we want a "clean" simulation.
+# Leaving the other pvwatts loss types (mismatch, wiring, etc) unspecified
+# causes them to take their default values.
+losses_parameters = dict(shading=0, availability=0)
+
+# Google Street View images show that each row is four modules high, in
+# landscape orientation. Assuming the modules are First Solar Series 4,
+# each of them is 600 mm wide.
+# Assume ~1 centimeter gap between modules (three gaps total).
+# And from Google Earth, the array's pitch is estimated to be about 7.0 meters.
+# From these we calculate the ground coverage ratio (GCR):
+pitch = 7  # meters
+gcr = (4 * 0.6 + 3 * 0.01) / pitch
+
+# The tracker rotation measurements reveal that the tracker rotation limits
+# are +/- 60 degrees, and backtracking is not enabled:
+max_angle = 60  # degrees
+backtrack = False
+
+# Google Earth shows that the tracker axes are very close to north-south:
+axis_azimuth = 180  # degrees
+
+# Estimated from Google Street View images
+axis_height = 1.5  # meters
+
+# %%
+# Create system objects
+# ---------------------
+#
+# The system has two inverters which seem to have identical specifications
+# and arrays. To save some code and computation repetition, we will just
+# model one inverter.
+
+location = pvlib.location.Location(latitude, longitude)
+mount = pvlib.pvsystem.SingleAxisTrackerMount(
+    gcr=gcr,
+    backtrack=backtrack,
+    max_angle=max_angle,
+    axis_azimuth=axis_azimuth
+)
+array = pvlib.pvsystem.Array(
+    mount,
+    module_parameters=module_parameters,
+    modules_per_string=modules_per_string,
+    temperature_model_parameters=temperature_model_parameters,
+    strings=strings_per_inverter
+)
+system = pvlib.pvsystem.PVSystem(
+    array,
+    inverter_parameters=inverter_parameters,
+    losses_parameters=losses_parameters
+)
+
+model = pvlib.modelchain.ModelChain(
+    system,
+    location,
+    spectral_model='first_solar',
+    aoi_model='physical',
+    losses_model='pvwatts'
+)
+
+# %%
+# Fetch weather data
+# ------------------
+#
+# The system does have measured plane-of-array irradiance data, but the
+# measurements suffer from row-to-row shading and tracker stalls. In this
+# example, we will use weather data taken from the NSRDB PSM3 for the year
+# 2019.
+
+api_key = 'DEMO_KEY'
+email = 'your_email@domain.com'
+
+keys = ['ghi', 'dni', 'dhi', 'temp_air', 'wind_speed',
+        'albedo', 'precipitable_water']
+psm3, psm3_metadata = pvlib.iotools.get_psm3(latitude, longitude, api_key,
+                                             email, interval=5, names=2019,
+                                             map_variables=True, leap_day=True,
+                                             attributes=keys)
+
+# %%
+# Pre-generate some model inputs
+# ------------------------------
+#
+# This system’s trackers are configured to not backtrack, meaning the
+# array shades itself when the sun is low in the sky. pvlib’s
+# ``ModelChain`` currently has no shade modeling ability, so we will model
+# it separately.
+#
+# Since this system uses thin-film modules, oriented in such a way that
+# row-to-row shadows affect each cell in the module equally, we can assume
+# that the effect of shading is linear with the reduction in incident beam
+# irradiance. That means we can use pvlib’s infinite sheds model, which
+# penalizes incident beam irradiance according to the calculated shaded
+# module fraction and returns the average irradiance over the total module
+# surface.
+
+solar_position = location.get_solarposition(psm3.index, latitude, longitude)
+tracker_angles = mount.get_orientation(
+    solar_position['apparent_zenith'],
+    solar_position['azimuth']
+)
+dni_extra = pvlib.irradiance.get_extra_radiation(psm3.index)
+
+# note: this system is monofacial, so only calculate irradiance for the
+# front side:
+averaged_irradiance = pvlib.bifacial.infinite_sheds.get_irradiance_poa(
+    tracker_angles['surface_tilt'], tracker_angles['surface_azimuth'],
+    solar_position['apparent_zenith'], solar_position['azimuth'],
+    gcr, axis_height, pitch,
+    psm3['ghi'], psm3['dhi'], psm3['dni'], psm3['albedo'],
+    model='haydavies', dni_extra=dni_extra,
+)
+
+# %%
+# ``ModelChain`` does not consider thermal transience either, so since we
+# are using 5-minute weather data, we will precalculate the cell
+# temperature as well:
+
+cell_temperature_steady_state = pvlib.temperature.faiman(
+    poa_global=averaged_irradiance['poa_global'],
+    temp_air=psm3['temp_air'],
+    wind_speed=psm3['wind_speed'],
+    **temperature_model_parameters,
+)
+
+cell_temperature = pvlib.temperature.prilliman(
+    cell_temperature_steady_state,
+    psm3['wind_speed'],
+    unit_mass=module_unit_mass
+)
+
+# %%
+# Run the model
+# -------------
+#
+# Finally, we are ready to run the rest of the system model. Since we want
+# to use pre-calculated plane-of-array irradiance, we will use
+# :py:meth:`~pvlib.modelchain.ModelChain.run_model_from_poa`:
+
+weather_inputs = pd.DataFrame({
+    'poa_global': averaged_irradiance['poa_global'],
+    'poa_direct': averaged_irradiance['poa_direct'],
+    'poa_diffuse': averaged_irradiance['poa_diffuse'],
+    'cell_temperature': cell_temperature,
+    'precipitable_water': psm3['precipitable_water'],  # for the spectral model
+})
+model.run_model_from_poa(weather_inputs)
+
+
+# %%
+# Compare with measured production
+# --------------------------------
+#
+# Now, let’s compare our modeled AC power with the system’s actual
+# inverter-level AC power measurements:
+
+fn = r"path/to/9068_ac_power_data.csv"
+df_inverter_measured = pd.read_csv(fn, index_col=0, parse_dates=True)
+df_inverter_measured = df_inverter_measured.tz_localize('US/Mountain',
+                                                        ambiguous='NaT',
+                                                        nonexistent='NaT')
+# convert to standard time to match the NSRDB-based simulation
+df_inverter_measured = df_inverter_measured.tz_convert('Etc/GMT+7')
+
+# %%
+
+inverter_ac_powers = [
+    'inverter_1_ac_power_(kw)_inv_150143',
+    'inverter_2_ac_power_(kw)_inv_150144'
+]
+df = df_inverter_measured.loc['2019', inverter_ac_powers]
+df['model'] = model.results.ac / 1000  # convert W to kW
+
+# %%
+
+for column_name in inverter_ac_powers:
+    fig, axes = plt.subplots(1, 3, figsize=(12, 4))
+    df.plot.scatter('model', column_name, ax=axes[0], s=1, alpha=0.1)
+    axes[0].axline((0, 0), slope=1, c='k')
+    axes[0].set_ylabel('Measured 5-min power [kW]')
+    axes[0].set_xlabel('Modeled 5-min power [kW]')
+
+    hourly_average = df.resample('h').mean()
+    hourly_average.plot.scatter('model', column_name, ax=axes[1], s=2)
+    axes[1].axline((0, 0), slope=1, c='k')
+    axes[1].set_ylabel('Measured hourly energy [kWh]')
+    axes[1].set_xlabel('Modeled hourly energy [kWh]')
+
+    daily_total = hourly_average.resample('d').sum()
+    daily_total.plot.scatter('model', column_name, ax=axes[2], s=5)
+    axes[2].axline((0, 0), slope=1, c='k')
+    axes[2].set_ylabel('Measured daily energy [kWh]')
+    axes[2].set_xlabel('Modeled daily energy [kWh]')
+
+    fig.suptitle(column_name)
+    fig.tight_layout()
+
+
+# %%
+# .. image:: ../../_images/OEDI_9068_inverter1_comparison.png
+#
+# .. image:: ../../_images/OEDI_9068_inverter2_comparison.png
+
+fig, ax = plt.subplots(figsize=(12, 4))
+daily_energy = df.clip(lower=0).resample('h').mean().resample('d').sum()
+daily_energy.plot(ax=ax)
+plt.ylim(bottom=0)
+plt.ylabel('Daily Production [kWh]')
+plt.tight_layout()
+
+# %%
+# .. image:: ../../_images/OEDI_9068_daily_timeseries.png

{pvlib-0.10.3 → pvlib-0.10.5}/docs/sphinx/source/conf.py

@@ -357,15 +357,18 @@ suppress_warnings = ['ref.footnote']
 sphinx_gallery_conf = {
     'examples_dirs': ['../../examples'],  # location of gallery scripts
     'gallery_dirs': ['gallery'],  # location of generated output
-    # sphinx-gallery only shows plots from plot_*.py files by default:
-    # 'filename_pattern': '*.py',
+    # execute all scripts except for ones in the "system-models" directory:
+    'filename_pattern': '^((?!system-models).)*$',
 
     # directory where function/class granular galleries are stored
     'backreferences_dir': 'reference/generated/gallery_backreferences',
 
     # Modules for which function/class level galleries are created. In
-    # this case only pvlib, could include others though.  must be tuple of str
+    # this case only pvlib, could include others though. Must be tuple of str
     'doc_module': ('pvlib',),
+
+    # https://sphinx-gallery.github.io/dev/configuration.html#removing-config-comments  # noqa: E501
+    'remove_config_comments': True,
 }
 # supress warnings in gallery output
 # https://sphinx-gallery.github.io/stable/configuration.html

{pvlib-0.10.3 → pvlib-0.10.5}/docs/sphinx/source/contributing.rst

@@ -29,12 +29,8 @@ pvlib-python, git, or Python:
   `easy <https://github.com/pvlib/pvlib-python/labels/easy>`_,
   or `help wanted <https://github.com/pvlib/pvlib-python/labels/help%20wanted>`_.
 * Improve the documentation and the unit tests.
-* Improve the IPython/Jupyter Notebook tutorials or write new ones that
+* Improve the Example Gallery or add new examples that
   demonstrate how to use pvlib-python in your area of expertise.
-* If you have MATLAB experience, you can help us keep pvlib-python
-  up to date with PVLIB_MATLAB or help us develop common unit tests.
-  For more, see `Issue #2 <https://github.com/pvlib/pvlib-python/issues/2>`_
-  and `Issue #3 <https://github.com/pvlib/pvlib-python/issues/3>`_.
 * Tell your friends and colleagues about pvlib-python.
 * Add your project to our
   `Projects and publications that use pvlib-python wiki