series input for ghi in test_infinite_sheds.py

Author: jason-rpkt

docs/sphinx/source/whatsnew/v0.11.3.rst

@@ -10,7 +10,8 @@ Deprecations
 
 Enhancements
 ~~~~~~~~~~~~
-Fix syntax for unshaded ground fraction (:issue:`2245`, :pull:`2359`)
+Fix bug in :py:func:`pvlib.bifacial.get_irradiance_poa` which may had yielded non-zero ground
+irradiance with sun below horizon (:issue:`2245`, :pull:`2359`)
 
 Documentation
 ~~~~~~~~~~~~~

pvlib/bifacial/infinite_sheds.py

@@ -260,9 +260,8 @@ def get_irradiance_poa(surface_tilt, surface_azimuth, solar_zenith,
     Returns
     -------
     output : dict or DataFrame
-        Output is a "pandas.DataFrame" when "ghi" is a Series.
-        Otherwise it is a dict of "numpy.ndarray"
-        See Issue #2245
+        Output is a ``pandas.DataFrame`` when ``ghi`` is a Series.
+        Otherwise it is a dict of ``numpy.ndarray``
         See Notes for descriptions of content.
 
     Notes

pvlib/bifacial/utils.py

@@ -87,8 +87,8 @@ def _unshaded_ground_fraction(surface_tilt, surface_azimuth, solar_zenith,
                                         surface_azimuth)
     f_gnd_beam = 1.0 - np.minimum(
         1.0, gcr * np.abs(cosd(surface_tilt) + sind(surface_tilt) * tan_phi))
-    f_gnd_beam = np.where(solar_zenith > max_zenith, 0., f_gnd_beam)
     # [1], Eq. 4
+    f_gnd_beam = np.where(solar_zenith > max_zenith, 0., f_gnd_beam)
     return f_gnd_beam  # 1 - min(1, abs()) < 1 always
 
 

pvlib/tests/bifacial/test_infinite_sheds.py

@@ -130,6 +130,7 @@ def test_get_irradiance_poa():
     # series inputs
     surface_tilt = pd.Series(surface_tilt)
     surface_azimuth = pd.Series(data=surface_azimuth, index=surface_tilt.index)
+    ghi = pd.Series(data=ghi, index=surface_tilt.index)
     solar_zenith = pd.Series(solar_zenith, index=surface_tilt.index)
     solar_azimuth = pd.Series(data=solar_azimuth, index=surface_tilt.index)
     expected_diffuse = pd.Series(