add-test_julian_day_issue_2207

Author: dgapitts

pvlib/tests/test_spa.py

@@ -1,6 +1,7 @@
 import os
 import datetime as dt
 import warnings
+import pytest
 
 try:
     from importlib import reload
@@ -423,3 +424,28 @@ def test_solar_position_multithreaded(self):
             nresult, self.spa.solar_position(
                 times, lat, lon, elev, pressure, temp, delta_t,
                 atmos_refract, numthreads=3, sst=True)[:3], 5)
+
+    # Define extra test cases for issue #2077
+    test_cases_issue_2207 = [
+        ((2000, 1, 1, 12, 0, 0), 2451545.0),
+        ((1999, 1, 1, 0, 0, 0), 2451179.5),
+        ((1987, 1, 27, 0, 0, 0), 2446822.5),
+        ((1987, 6, 19, 12, 0, 0), 2446966.0),
+        ((1988, 1, 27, 0, 0, 0), 2447187.5),
+        ((1988, 6, 19, 12, 0, 0), 2447332.0),
+        ((1900, 1, 1, 0, 0, 0), 2415020.5),
+        ((1600, 1, 1, 0, 0, 0), 2305447.5),
+        ((1600, 12, 31, 0, 0, 0), 2305812.5),
+        ((837, 4, 10, 7, 12, 0), 2026871.8),
+        ((-123, 12, 31, 0, 0, 0), 1676496.5),
+        ((-122, 1, 1, 0, 0, 0), 1676497.5),
+        ((-1000, 7, 12, 12, 0, 0), 1356001.0),
+        ((-1000, 2, 29, 0, 0, 0), 1355866.5),
+        ((-1001, 8, 17, 21, 36, 0), 1355671.4),
+        ((-4712, 1, 1, 12, 0, 0), 0.1),
+    ]
+
+    @pytest.mark.parametrize("inputs, expected", test_cases_issue_2207)
+    def test_julian_day_issue_2207(inputs, expected):
+        result = pvlib.spa.julian_day_dt(*inputs, microsecond=0)
+        assert result == expected, f"Failed for inputs {inputs}: expected {expected}, got {result}"