pvlib 0.11.2__py3-none-any.whl → 0.12.0__py3-none-any.whl

pvlib/tests/test_spa.py

@@ -1,454 +0,0 @@
-import os
-import datetime as dt
-import warnings
-import pytest
-import pvlib
-
-try:
-    from importlib import reload
-except ImportError:
-    try:
-        from imp import reload
-    except ImportError:
-        pass
-
-import numpy as np
-from numpy.testing import assert_almost_equal
-import pandas as pd
-
-import unittest
-from .conftest import requires_numba
-
-
-times = (pd.date_range('2003-10-17 12:30:30', periods=1, freq='D')
-           .tz_localize('MST'))
-unixtimes = np.array(times.tz_convert('UTC').view(np.int64)*1.0/10**9)
-
-lat = 39.742476
-lon = -105.1786
-elev = 1830.14
-pressure = 820
-temp = 11
-delta_t = 67.0
-atmos_refract = 0.5667
-
-JD = 2452930.312847
-JC = 0.0379277986858
-JDE = 2452930.313623
-JCE = 0.037927819916852
-JME = 0.003792781991685
-L = 24.0182616917
-B = -0.0001011219
-R = 0.9965422974
-Theta = 204.0182616917
-beta = 0.0001011219
-X0 = 17185.861179
-X1 = 1722.893218
-X2 = 18234.075703
-X3 = 18420.071012
-X4 = 51.686951
-dPsi = -0.00399840
-dEpsilon = 0.00166657
-epsilon0 = 84379.672625
-epsilon = 23.440465
-dTau = -0.005711
-lamd = 204.0085519281
-v0 = 318.515579
-v = 318.511910
-alpha = 202.227408
-delta = -9.31434
-H = 11.10590
-xi = 0.002451
-dAlpha = -0.000369
-alpha_prime = 202.22704
-delta_prime = -9.316179
-H_prime = 11.10627
-e0 = 39.872046
-de = 0.016332
-e = 39.888378
-theta = 50.11162
-theta0 = 90 - e0
-Gamma = 14.340241
-Phi = 194.340241
-year = 1985
-month = 2
-year_array = np.array([-499, 500, 1000, 1500, 1800, 1860, 1900, 1950,
-                       1970, 1985, 1990, 2000, 2005, 2050, 2150])
-# `month_array` is used with `year_array` in `test_calculate_deltat`.
-# Both arrays need to have the same length for the test, hence the duplicates.
-month_array = np.array([1, 2, 3, 4, 5, 6, 6, 7, 8, 9, 10, 11, 12, 12, 12])
-dt_actual = 54.413442486
-dt_actual_array = np.array([1.7184831e+04, 5.7088051e+03, 1.5730419e+03,
-                            1.9801820e+02, 1.3596506e+01, 7.8316585e+00,
-                            -2.1171894e+00, 2.9289261e+01, 4.0824887e+01,
-                            5.4724581e+01, 5.7426651e+01, 6.4108015e+01,
-                            6.5038015e+01, 9.4952955e+01, 3.3050693e+02])
-mix_year_array = np.full((10), year)
-mix_month_array = np.full((10), month)
-mix_year_actual = np.full((10), dt_actual)
-mix_month_actual = mix_year_actual
-
-
-class SpaBase:
-    """Test functions common to numpy and numba spa"""
-    def test_julian_day_dt(self):
-        # add 1us manually to the test timestamp (GH #940)
-        dt = times.tz_convert('UTC')[0] + pd.Timedelta(1, unit='us')
-        year = dt.year
-        month = dt.month
-        day = dt.day
-        hour = dt.hour
-        minute = dt.minute
-        second = dt.second
-        microsecond = dt.microsecond
-        assert_almost_equal(JD + 1e-6 / (3600*24),  # modify expected JD by 1us
-                            self.spa.julian_day_dt(
-                                year, month, day, hour,
-                                minute, second, microsecond), 6)
-
-    def test_julian_ephemeris_day(self):
-        assert_almost_equal(JDE, self.spa.julian_ephemeris_day(JD, delta_t), 5)
-
-    def test_julian_century(self):
-        assert_almost_equal(JC, self.spa.julian_century(JD), 6)
-
-    def test_julian_ephemeris_century(self):
-        assert_almost_equal(JCE, self.spa.julian_ephemeris_century(JDE), 10)
-
-    def test_julian_ephemeris_millenium(self):
-        assert_almost_equal(JME, self.spa.julian_ephemeris_millennium(JCE), 10)
-
-    def test_heliocentric_longitude(self):
-        assert_almost_equal(L, self.spa.heliocentric_longitude(JME), 6)
-
-    def test_heliocentric_latitude(self):
-        assert_almost_equal(B, self.spa.heliocentric_latitude(JME), 6)
-
-    def test_heliocentric_radius_vector(self):
-        assert_almost_equal(R, self.spa.heliocentric_radius_vector(JME), 6)
-
-    def test_geocentric_longitude(self):
-        assert_almost_equal(Theta, self.spa.geocentric_longitude(L), 6)
-
-    def test_geocentric_latitude(self):
-        assert_almost_equal(beta, self.spa.geocentric_latitude(B), 6)
-
-    def test_mean_elongation(self):
-        assert_almost_equal(X0, self.spa.mean_elongation(JCE), 5)
-
-    def test_mean_anomaly_sun(self):
-        assert_almost_equal(X1, self.spa.mean_anomaly_sun(JCE), 5)
-
-    def test_mean_anomaly_moon(self):
-        assert_almost_equal(X2, self.spa.mean_anomaly_moon(JCE), 5)
-
-    def test_moon_argument_latitude(self):
-        assert_almost_equal(X3, self.spa.moon_argument_latitude(JCE), 5)
-
-    def test_moon_ascending_longitude(self):
-        assert_almost_equal(X4, self.spa.moon_ascending_longitude(JCE), 6)
-
-    def test_longitude_obliquity_nutation(self):
-        out = np.empty((2,))
-        self.spa.longitude_obliquity_nutation(JCE, X0, X1, X2, X3, X4, out)
-        _dPsi, _dEpsilon = out[0], out[1]
-        assert_almost_equal(dPsi, _dPsi, 6)
-        assert_almost_equal(dEpsilon, _dEpsilon, 6)
-
-    def test_mean_ecliptic_obliquity(self):
-        assert_almost_equal(epsilon0, self.spa.mean_ecliptic_obliquity(JME), 6)
-
-    def test_true_ecliptic_obliquity(self):
-        assert_almost_equal(epsilon, self.spa.true_ecliptic_obliquity(
-            epsilon0, dEpsilon), 6)
-
-    def test_aberration_correction(self):
-        assert_almost_equal(dTau, self.spa.aberration_correction(R), 6)
-
-    def test_apparent_sun_longitude(self):
-        assert_almost_equal(lamd, self.spa.apparent_sun_longitude(
-            Theta, dPsi, dTau), 6)
-
-    def test_mean_sidereal_time(self):
-        assert_almost_equal(v0, self.spa.mean_sidereal_time(JD, JC), 3)
-
-    def test_apparent_sidereal_time(self):
-        assert_almost_equal(v, self.spa.apparent_sidereal_time(
-            v0, dPsi, epsilon), 5)
-
-    def test_geocentric_sun_right_ascension(self):
-        assert_almost_equal(alpha, self.spa.geocentric_sun_right_ascension(
-            lamd, epsilon, beta), 6)
-
-    def test_geocentric_sun_declination(self):
-        assert_almost_equal(delta, self.spa.geocentric_sun_declination(
-            lamd, epsilon, beta), 6)
-
-    def test_local_hour_angle(self):
-        assert_almost_equal(H, self.spa.local_hour_angle(v, lon, alpha), 4)
-
-    def test_equatorial_horizontal_parallax(self):
-        assert_almost_equal(xi, self.spa.equatorial_horizontal_parallax(R), 6)
-
-    def test_parallax_sun_right_ascension(self):
-        u = self.spa.uterm(lat)
-        x = self.spa.xterm(u, lat, elev)
-        assert_almost_equal(dAlpha, self.spa.parallax_sun_right_ascension(
-            x, xi, H, delta), 4)
-
-    def test_topocentric_sun_right_ascension(self):
-        assert_almost_equal(alpha_prime,
-                            self.spa.topocentric_sun_right_ascension(
-                                alpha, dAlpha), 5)
-
-    def test_topocentric_sun_declination(self):
-        u = self.spa.uterm(lat)
-        x = self.spa.xterm(u, lat, elev)
-        y = self.spa.yterm(u, lat, elev)
-        assert_almost_equal(delta_prime, self.spa.topocentric_sun_declination(
-            delta, x, y, xi, dAlpha, H), 5)
-
-    def test_topocentric_local_hour_angle(self):
-        assert_almost_equal(H_prime, self.spa.topocentric_local_hour_angle(
-            H, dAlpha), 5)
-
-    def test_topocentric_elevation_angle_without_atmosphere(self):
-        assert_almost_equal(
-            e0, self.spa.topocentric_elevation_angle_without_atmosphere(
-                lat, delta_prime, H_prime), 6)
-
-    def test_atmospheric_refraction_correction(self):
-        assert_almost_equal(de, self.spa.atmospheric_refraction_correction(
-            pressure, temp, e0, atmos_refract), 6)
-
-    def test_topocentric_elevation_angle(self):
-        assert_almost_equal(e, self.spa.topocentric_elevation_angle(e0, de), 6)
-
-    def test_topocentric_zenith_angle(self):
-        assert_almost_equal(theta, self.spa.topocentric_zenith_angle(e), 5)
-
-    def test_topocentric_astronomers_azimuth(self):
-        assert_almost_equal(Gamma, self.spa.topocentric_astronomers_azimuth(
-            H_prime, delta_prime, lat), 5)
-
-    def test_topocentric_azimuth_angle(self):
-        assert_almost_equal(Phi, self.spa.topocentric_azimuth_angle(Gamma), 5)
-
-    def test_solar_position(self):
-        with warnings.catch_warnings():
-            # don't warn on method reload
-            warnings.simplefilter("ignore")
-            spa_out_0 = self.spa.solar_position(
-                unixtimes, lat, lon, elev, pressure, temp, delta_t,
-                atmos_refract)[:-1]
-            spa_out_1 = self.spa.solar_position(
-                unixtimes, lat, lon, elev, pressure, temp, delta_t,
-                atmos_refract, sst=True)[:3]
-        assert_almost_equal(np.array([[theta, theta0, e, e0, Phi]]).T,
-                            spa_out_0, 5)
-        assert_almost_equal(np.array([[v, alpha, delta]]).T, spa_out_1, 5)
-
-    def test_equation_of_time(self):
-        eot = 14.64
-        M = self.spa.sun_mean_longitude(JME)
-        assert_almost_equal(eot, self.spa.equation_of_time(
-            M, alpha, dPsi, epsilon), 2)
-
-    def test_transit_sunrise_sunset(self):
-        # tests at greenwich
-        times = pd.DatetimeIndex([dt.datetime(1996, 7, 5, 0),
-                                  dt.datetime(2004, 12, 4, 0)]
-                                 ).tz_localize(
-                                     'UTC').view(np.int64)*1.0/10**9
-        sunrise = pd.DatetimeIndex([dt.datetime(1996, 7, 5, 7, 8, 15),
-                                    dt.datetime(2004, 12, 4, 4, 38, 57)]
-                                   ).tz_localize(
-                                       'UTC').view(np.int64)*1.0/10**9
-        sunset = pd.DatetimeIndex([dt.datetime(1996, 7, 5, 17, 1, 4),
-                                   dt.datetime(2004, 12, 4, 19, 2, 2)]
-                                  ).tz_localize(
-                                      'UTC').view(np.int64)*1.0/10**9
-        times = np.array(times)
-        sunrise = np.array(sunrise)
-        sunset = np.array(sunset)
-        result = self.spa.transit_sunrise_sunset(times, -35.0, 0.0, 64.0, 1)
-        assert_almost_equal(sunrise/1e3, result[1]/1e3, 3)
-        assert_almost_equal(sunset/1e3, result[2]/1e3, 3)
-
-        times = pd.DatetimeIndex([dt.datetime(1994, 1, 2), ]
-                                 ).tz_localize(
-                                     'UTC').view(np.int64)*1.0/10**9
-        sunset = pd.DatetimeIndex([dt.datetime(1994, 1, 2, 16, 59, 55), ]
-                                  ).tz_localize(
-                                      'UTC').view(np.int64)*1.0/10**9
-        sunrise = pd.DatetimeIndex([dt.datetime(1994, 1, 2, 7, 8, 12), ]
-                                   ).tz_localize(
-                                       'UTC').view(np.int64)*1.0/10**9
-        times = np.array(times)
-        sunrise = np.array(sunrise)
-        sunset = np.array(sunset)
-        result = self.spa.transit_sunrise_sunset(times, 35.0, 0.0, 64.0, 1)
-        assert_almost_equal(sunrise/1e3, result[1]/1e3, 3)
-        assert_almost_equal(sunset/1e3, result[2]/1e3, 3)
-
-        # tests from USNO
-        # Golden
-        times = pd.DatetimeIndex([dt.datetime(2015, 1, 2),
-                                  dt.datetime(2015, 4, 2),
-                                  dt.datetime(2015, 8, 2),
-                                  dt.datetime(2015, 12, 2)],
-                                 ).tz_localize(
-                                     'UTC').view(np.int64)*1.0/10**9
-        sunrise = pd.DatetimeIndex([dt.datetime(2015, 1, 2, 7, 19),
-                                    dt.datetime(2015, 4, 2, 5, 43),
-                                    dt.datetime(2015, 8, 2, 5, 1),
-                                    dt.datetime(2015, 12, 2, 7, 1)],
-                                   ).tz_localize(
-                                       'MST').view(np.int64)*1.0/10**9
-        sunset = pd.DatetimeIndex([dt.datetime(2015, 1, 2, 16, 49),
-                                   dt.datetime(2015, 4, 2, 18, 24),
-                                   dt.datetime(2015, 8, 2, 19, 10),
-                                   dt.datetime(2015, 12, 2, 16, 38)],
-                                  ).tz_localize(
-                                      'MST').view(np.int64)*1.0/10**9
-        times = np.array(times)
-        sunrise = np.array(sunrise)
-        sunset = np.array(sunset)
-        result = self.spa.transit_sunrise_sunset(times, 39.0, -105.0, 64.0, 1)
-        assert_almost_equal(sunrise/1e3, result[1]/1e3, 1)
-        assert_almost_equal(sunset/1e3, result[2]/1e3, 1)
-
-        # Beijing
-        times = pd.DatetimeIndex([dt.datetime(2015, 1, 2),
-                                  dt.datetime(2015, 4, 2),
-                                  dt.datetime(2015, 8, 2),
-                                  dt.datetime(2015, 12, 2)],
-                                 ).tz_localize(
-                                     'UTC').view(np.int64)*1.0/10**9
-        sunrise = pd.DatetimeIndex([dt.datetime(2015, 1, 2, 7, 36),
-                                    dt.datetime(2015, 4, 2, 5, 58),
-                                    dt.datetime(2015, 8, 2, 5, 13),
-                                    dt.datetime(2015, 12, 2, 7, 17)],
-                                   ).tz_localize('Asia/Shanghai').view(
-                                       np.int64)*1.0/10**9
-        sunset = pd.DatetimeIndex([dt.datetime(2015, 1, 2, 17, 0),
-                                   dt.datetime(2015, 4, 2, 18, 39),
-                                   dt.datetime(2015, 8, 2, 19, 28),
-                                   dt.datetime(2015, 12, 2, 16, 50)],
-                                  ).tz_localize('Asia/Shanghai').view(
-                                      np.int64)*1.0/10**9
-        times = np.array(times)
-        sunrise = np.array(sunrise)
-        sunset = np.array(sunset)
-        result = self.spa.transit_sunrise_sunset(
-            times, 39.917, 116.383, 64.0, 1)
-        assert_almost_equal(sunrise/1e3, result[1]/1e3, 1)
-        assert_almost_equal(sunset/1e3, result[2]/1e3, 1)
-
-    def test_earthsun_distance(self):
-        times = (pd.date_range('2003-10-17 12:30:30', periods=1, freq='D')
-                   .tz_localize('MST'))
-        unixtimes = times.tz_convert('UTC').view(np.int64)*1.0/10**9
-        unixtimes = np.array(unixtimes)
-        result = self.spa.earthsun_distance(unixtimes, 64.0, 1)
-        assert_almost_equal(R, result, 6)
-
-    def test_calculate_deltat(self):
-        result_mix_year = self.spa.calculate_deltat(mix_year_array, month)
-        assert_almost_equal(mix_year_actual, result_mix_year)
-
-        result_mix_month = self.spa.calculate_deltat(year, mix_month_array)
-        assert_almost_equal(mix_month_actual, result_mix_month)
-
-        result_array = self.spa.calculate_deltat(year_array, month_array)
-        assert_almost_equal(dt_actual_array, result_array, 3)
-
-        result_scalar = self.spa.calculate_deltat(year, month)
-        assert_almost_equal(dt_actual, result_scalar)
-
-
-class NumpySpaTest(unittest.TestCase, SpaBase):
-    """Import spa without compiling to numba then run tests"""
-    @classmethod
-    def setUpClass(self):
-        os.environ['PVLIB_USE_NUMBA'] = '0'
-        import pvlib.spa as spa
-        spa = reload(spa)
-        self.spa = spa
-
-    @classmethod
-    def tearDownClass(self):
-        del os.environ['PVLIB_USE_NUMBA']
-
-    def test_julian_day(self):
-        assert_almost_equal(JD, self.spa.julian_day(unixtimes)[0], 6)
-
-
-@requires_numba
-class NumbaSpaTest(unittest.TestCase, SpaBase):
-    """Import spa, compiling to numba, and run tests"""
-    @classmethod
-    def setUpClass(self):
-        os.environ['PVLIB_USE_NUMBA'] = '1'
-        import pvlib.spa as spa
-        spa = reload(spa)
-        self.spa = spa
-
-    @classmethod
-    def tearDownClass(self):
-        del os.environ['PVLIB_USE_NUMBA']
-
-    def test_julian_day(self):
-        assert_almost_equal(JD, self.spa.julian_day(unixtimes[0]), 6)
-
-    def test_solar_position_singlethreaded(self):
-        assert_almost_equal(
-            np.array([[theta, theta0, e, e0, Phi]]).T, self.spa.solar_position(
-                unixtimes, lat, lon, elev, pressure, temp, delta_t,
-                atmos_refract, numthreads=1)[:-1], 5)
-        assert_almost_equal(
-            np.array([[v, alpha, delta]]).T, self.spa.solar_position(
-                unixtimes, lat, lon, elev, pressure, temp, delta_t,
-                atmos_refract, numthreads=1, sst=True)[:3], 5)
-
-    def test_solar_position_multithreaded(self):
-        result = np.array([theta, theta0, e, e0, Phi])
-        nresult = np.array([result, result, result]).T
-        times = np.array([unixtimes[0], unixtimes[0], unixtimes[0]])
-        assert_almost_equal(
-            nresult, self.spa.solar_position(
-                times, lat, lon, elev, pressure, temp, delta_t,
-                atmos_refract, numthreads=3)[:-1], 5)
-        result = np.array([v, alpha, delta])
-        nresult = np.array([result, result, result]).T
-        assert_almost_equal(
-            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.0),
-]
-
-
-@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}"