[175347046.0, 0.0, 0.0],

[3341656.0, 4.6692568, 6283.07585],

[34894.0, 4.6261, 12566.1517],

[3497.0, 2.7441, 5753.3849],

[3418.0, 2.8289, 3.5231],

[3136.0, 3.6277, 77713.7715],

[2676.0, 4.4181, 7860.4194],

[2343.0, 6.1352, 3930.2097],

[1324.0, 0.7425, 11506.7698],

[1273.0, 2.0371, 529.691],

[1199.0, 1.1096, 1577.3435],

[990.0, 5.233, 5884.927],

[902.0, 2.045, 26.298],

[857.0, 3.508, 398.149],

[780.0, 1.179, 5223.694],

[753.0, 2.533, 5507.553],

[505.0, 4.583, 18849.228],

[492.0, 4.205, 775.523],

[357.0, 2.92, 0.067],

[317.0, 5.849, 11790.629],

[284.0, 1.899, 796.298],

[271.0, 0.315, 10977.079],

[243.0, 0.345, 5486.778],

[206.0, 4.806, 2544.314],

[205.0, 1.869, 5573.143],

[202.0, 2.458, 6069.777],

[156.0, 0.833, 213.299],

[132.0, 3.411, 2942.463],

[126.0, 1.083, 20.775],

[115.0, 0.645, 0.98],

[103.0, 0.636, 4694.003],

[102.0, 0.976, 15720.839],

[102.0, 4.267, 7.114],

[99.0, 6.21, 2146.17],

[98.0, 0.68, 155.42],

[86.0, 5.98, 161000.69],

[85.0, 1.3, 6275.96],

[85.0, 3.67, 71430.7],

[80.0, 1.81, 17260.15],

[79.0, 3.04, 12036.46],

[75.0, 1.76, 5088.63],

[74.0, 3.5, 3154.69],

[74.0, 4.68, 801.82],

[70.0, 0.83, 9437.76],

[62.0, 3.98, 8827.39],

[61.0, 1.82, 7084.9],

[57.0, 2.78, 6286.6],

[56.0, 4.39, 14143.5],

[56.0, 3.47, 6279.55],

[52.0, 0.19, 12139.55],

[52.0, 1.33, 1748.02],

[51.0, 0.28, 5856.48],

[49.0, 0.49, 1194.45],

[41.0, 5.37, 8429.24],

[41.0, 2.4, 19651.05],

[39.0, 6.17, 10447.39],

[37.0, 6.04, 10213.29],

[37.0, 2.57, 1059.38],

[36.0, 1.71, 2352.87],

[36.0, 1.78, 6812.77],

[33.0, 0.59, 17789.85],

[30.0, 0.44, 83996.85],

[30.0, 2.74, 1349.87],

[25.0, 3.16, 4690.48]

[628331966747.0, 0.0, 0.0],

[206059.0, 2.678235, 6283.07585],

[4303.0, 2.6351, 12566.1517],

[425.0, 1.59, 3.523],

[119.0, 5.796, 26.298],

[109.0, 2.966, 1577.344],

[93.0, 2.59, 18849.23],

[72.0, 1.14, 529.69],

[68.0, 1.87, 398.15],

[67.0, 4.41, 5507.55],

[59.0, 2.89, 5223.69],

[56.0, 2.17, 155.42],

[45.0, 0.4, 796.3],

[36.0, 0.47, 775.52],

[29.0, 2.65, 7.11],

[21.0, 5.34, 0.98],

[19.0, 1.85, 5486.78],

[19.0, 4.97, 213.3],

[17.0, 2.99, 6275.96],

[16.0, 0.03, 2544.31],

[16.0, 1.43, 2146.17],

[15.0, 1.21, 10977.08],

[12.0, 2.83, 1748.02],

[12.0, 3.26, 5088.63],

[12.0, 5.27, 1194.45],

[12.0, 2.08, 4694.0],

[11.0, 0.77, 553.57],

[10.0, 1.3, 6286.6],

[10.0, 4.24, 1349.87],

[9.0, 2.7, 242.73],

[9.0, 5.64, 951.72],

[8.0, 5.3, 2352.87],

[6.0, 2.65, 9437.76],

[6.0, 4.67, 4690.48]

[52919.0, 0.0, 0.0],

[8720.0, 1.0721, 6283.0758],

[309.0, 0.867, 12566.152],

[27.0, 0.05, 3.52],

[16.0, 5.19, 26.3],

[16.0, 3.68, 155.42],

[10.0, 0.76, 18849.23],

[9.0, 2.06, 77713.77],

[7.0, 0.83, 775.52],

[5.0, 4.66, 1577.34],

[4.0, 1.03, 7.11],

[4.0, 3.44, 5573.14],

[3.0, 5.14, 796.3],

[3.0, 6.05, 5507.55],

[3.0, 1.19, 242.73],

[3.0, 6.12, 529.69],

[3.0, 0.31, 398.15],

[3.0, 2.28, 553.57],

[2.0, 4.38, 5223.69],

[2.0, 3.75, 0.98]

[289.0, 5.844, 6283.076],

[35.0, 0.0, 0.0],

[17.0, 5.49, 12566.15],

[3.0, 5.2, 155.42],

[1.0, 4.72, 3.52],

[1.0, 5.3, 18849.23],

[1.0, 5.97, 242.73]

[114.0, 3.142, 0.0],

[8.0, 4.13, 6283.08],

[1.0, 3.84, 12566.15]

[280.0, 3.199, 84334.662],

[102.0, 5.422, 5507.553],

[80.0, 3.88, 5223.69],

[44.0, 3.7, 2352.87],

[32.0, 4.0, 1577.34]

[100013989.0, 0.0, 0.0],

[1670700.0, 3.0984635, 6283.07585],

[13956.0, 3.05525, 12566.1517],

[3084.0, 5.1985, 77713.7715],

[1628.0, 1.1739, 5753.3849],

[1576.0, 2.8469, 7860.4194],

[925.0, 5.453, 11506.77],

[542.0, 4.564, 3930.21],

[472.0, 3.661, 5884.927],

[346.0, 0.964, 5507.553],

[329.0, 5.9, 5223.694],

[307.0, 0.299, 5573.143],

[243.0, 4.273, 11790.629],

[212.0, 5.847, 1577.344],

[186.0, 5.022, 10977.079],

[175.0, 3.012, 18849.228],

[110.0, 5.055, 5486.778],

[98.0, 0.89, 6069.78],

[86.0, 5.69, 15720.84],

[86.0, 1.27, 161000.69],

[65.0, 0.27, 17260.15],

[63.0, 0.92, 529.69],

[57.0, 2.01, 83996.85],

[56.0, 5.24, 71430.7],

[49.0, 3.25, 2544.31],

[47.0, 2.58, 775.52],

[45.0, 5.54, 9437.76],

[43.0, 6.01, 6275.96],

[39.0, 5.36, 4694.0],

[38.0, 2.39, 8827.39],

[37.0, 0.83, 19651.05],

[37.0, 4.9, 12139.55],

[36.0, 1.67, 12036.46],

[35.0, 1.84, 2942.46],

[33.0, 0.24, 7084.9],

[32.0, 0.18, 5088.63],

[32.0, 1.78, 398.15],

[28.0, 1.21, 6286.6],

[28.0, 1.9, 6279.55],

[26.0, 4.59, 10447.39]

[103019.0, 1.10749, 6283.07585],

[1721.0, 1.0644, 12566.1517],

[702.0, 3.142, 0.0],

[32.0, 1.02, 18849.23],

[31.0, 2.84, 5507.55],

[25.0, 1.32, 5223.69],

[18.0, 1.42, 1577.34],

[10.0, 5.91, 10977.08],

[9.0, 1.42, 6275.96],

[9.0, 0.27, 5486.78]

[4359.0, 5.7846, 6283.0758],

[124.0, 5.579, 12566.152],

[12.0, 3.14, 0.0],

[9.0, 3.63, 77713.77],

[6.0, 1.87, 5573.14],

[3.0, 5.47, 18849.23]

[-171996, -174.2, 92025, 8.9],

[-13187, -1.6, 5736, -3.1],

[-2274, -0.2, 977, -0.5],

[2062, 0.2, -895, 0.5],

[1426, -3.4, 54, -0.1],

[712, 0.1, -7, 0],

[-517, 1.2, 224, -0.6],

[-386, -0.4, 200, 0],

[-301, 0, 129, -0.1],

[217, -0.5, -95, 0.3],

[-158, 0, 0, 0],

[129, 0.1, -70, 0],

[123, 0, -53, 0],

[63, 0, 0, 0],

[63, 0.1, -33, 0],

[-59, 0, 26, 0],

[-58, -0.1, 32, 0],

[-51, 0, 27, 0],

[48, 0, 0, 0],

[46, 0, -24, 0],

[-38, 0, 16, 0],

[-31, 0, 13, 0],

[29, 0, 0, 0],

[29, 0, -12, 0],

[26, 0, 0, 0],

[-22, 0, 0, 0],

[21, 0, -10, 0],

[17, -0.1, 0, 0],

[16, 0, -8, 0],

[-16, 0.1, 7, 0],

[-15, 0, 9, 0],

[-13, 0, 7, 0],

[-12, 0, 6, 0],

[11, 0, 0, 0],

[-10, 0, 5, 0],

[-8, 0, 3, 0],

[7, 0, -3, 0],

[-7, 0, 0, 0],

[-7, 0, 3, 0],

[-7, 0, 3, 0],

[6, 0, 0, 0],

[6, 0, -3, 0],

[6, 0, -3, 0],

[-6, 0, 3, 0],

[-6, 0, 3, 0],

[5, 0, 0, 0],

[-5, 0, 3, 0],

[-5, 0, 3, 0],

[-5, 0, 3, 0],

[4, 0, 0, 0],

[4, 0, 0, 0],

[4, 0, 0, 0],

[-4, 0, 0, 0],

[-4, 0, 0, 0],

[-4, 0, 0, 0],

[3, 0, 0, 0],

[-3, 0, 0, 0],

[-3, 0, 0, 0],

[-3, 0, 0, 0],

[-3, 0, 0, 0],

[-3, 0, 0, 0],

[-3, 0, 0, 0],

[-3, 0, 0, 0],

[0, 0, 0, 0, 1],

[-2, 0, 0, 2, 2],

[0, 0, 0, 2, 2],

[0, 0, 0, 0, 2],

[0, 1, 0, 0, 0],

[0, 0, 1, 0, 0],

[-2, 1, 0, 2, 2],

[0, 0, 0, 2, 1],

[0, 0, 1, 2, 2],

[-2, -1, 0, 2, 2],

[-2, 0, 1, 0, 0],

[-2, 0, 0, 2, 1],

[0, 0, -1, 2, 2],

[2, 0, 0, 0, 0],

[0, 0, 1, 0, 1],

[2, 0, -1, 2, 2],

[0, 0, -1, 0, 1],

[0, 0, 1, 2, 1],

[-2, 0, 2, 0, 0],

[0, 0, -2, 2, 1],

[2, 0, 0, 2, 2],

[0, 0, 2, 2, 2],

[0, 0, 2, 0, 0],

[-2, 0, 1, 2, 2],

[0, 0, 0, 2, 0],

[-2, 0, 0, 2, 0],

[0, 0, -1, 2, 1],

[0, 2, 0, 0, 0],

[2, 0, -1, 0, 1],

[-2, 2, 0, 2, 2],

[0, 1, 0, 0, 1],

[-2, 0, 1, 0, 1],

[0, -1, 0, 0, 1],

[0, 0, 2, -2, 0],

[2, 0, -1, 2, 1],

[2, 0, 1, 2, 2],

[0, 1, 0, 2, 2],

[-2, 1, 1, 0, 0],

[0, -1, 0, 2, 2],

[2, 0, 0, 2, 1],

[2, 0, 1, 0, 0],

[-2, 0, 2, 2, 2],

[-2, 0, 1, 2, 1],

[2, 0, -2, 0, 1],

[2, 0, 0, 0, 1],

[0, -1, 1, 0, 0],

[-2, -1, 0, 2, 1],

[-2, 0, 0, 0, 1],

[0, 0, 2, 2, 1],

[-2, 0, 2, 0, 1],

[-2, 1, 0, 2, 1],

[0, 0, 1, -2, 0],

[-1, 0, 1, 0, 0],

[-2, 1, 0, 0, 0],

[1, 0, 0, 0, 0],

[0, 0, 1, 2, 0],

[0, 0, -2, 2, 2],

[-1, -1, 1, 0, 0],

[0, 1, 1, 0, 0],

[0, -1, 1, 2, 2],

[2, -1, -1, 2, 2],

[0, 0, 3, 2, 2],

[2, -1, 0, 2, 2],

"""This is the original way to calculate the julian day from the NREL paper.

if month <= 2:

year = year-1

month = month+12

a = int(year/100)

b = 2 - a + int(a * 0.25)

frac_of_day = (microsecond / 1e6 + (second + minute * 60 + hour * 3600)

) * 1.0 / (3600*24)

d = day + frac_of_day

jd = (int(365.25 * (year + 4716)) + int(30.6001 * (month + 1)) + d +

b - 1524.5)

# shared calculation used for heliocentric longitude, latitude, and radius

s = 0.

for row in range(arr.shape[0]):

s += arr[row, 0] * np.cos(arr[row, 1] + arr[row, 2] * x)

l0 = sum_mult_cos_add_mult(L0, jme)

l1 = sum_mult_cos_add_mult(L1, jme)

l2 = sum_mult_cos_add_mult(L2, jme)

l3 = sum_mult_cos_add_mult(L3, jme)

l4 = sum_mult_cos_add_mult(L4, jme)

l5 = sum_mult_cos_add_mult(L5, jme)

l_rad = (l0 + l1 * jme + l2 * jme**2 + l3 * jme**3 + l4 * jme**4 +

l5 * jme**5)/10**8

l = np.rad2deg(l_rad)

b0 = sum_mult_cos_add_mult(B0, jme)

b1 = sum_mult_cos_add_mult(B1, jme)

b_rad = (b0 + b1 * jme)/10**8

b = np.rad2deg(b_rad)

r0 = sum_mult_cos_add_mult(R0, jme)

r1 = sum_mult_cos_add_mult(R1, jme)

r2 = sum_mult_cos_add_mult(R2, jme)

r3 = sum_mult_cos_add_mult(R3, jme)

r4 = sum_mult_cos_add_mult(R4, jme)

r = (r0 + r1 * jme + r2 * jme**2 + r3 * jme**3 + r4 * jme**4)/10**8

x0 = (297.85036

+ 445267.111480 * julian_ephemeris_century

- 0.0019142 * julian_ephemeris_century**2

+ julian_ephemeris_century**3 / 189474)

x1 = (357.52772

+ 35999.050340 * julian_ephemeris_century

- 0.0001603 * julian_ephemeris_century**2

- julian_ephemeris_century**3 / 300000)

x2 = (134.96298

+ 477198.867398 * julian_ephemeris_century

+ 0.0086972 * julian_ephemeris_century**2

+ julian_ephemeris_century**3 / 56250)

x3 = (93.27191

+ 483202.017538 * julian_ephemeris_century

- 0.0036825 * julian_ephemeris_century**2

+ julian_ephemeris_century**3 / 327270)

x4 = (125.04452

- 1934.136261 * julian_ephemeris_century

+ 0.0020708 * julian_ephemeris_century**2

+ julian_ephemeris_century**3 / 450000)

out):

delta_psi_sum = 0.0

delta_eps_sum = 0.0

for row in range(NUTATION_YTERM_ARRAY.shape[0]):

a = NUTATION_ABCD_ARRAY[row, 0]

b = NUTATION_ABCD_ARRAY[row, 1]

c = NUTATION_ABCD_ARRAY[row, 2]

d = NUTATION_ABCD_ARRAY[row, 3]

arg = np.radians(

NUTATION_YTERM_ARRAY[row, 0]*x0 +

NUTATION_YTERM_ARRAY[row, 1]*x1 +

NUTATION_YTERM_ARRAY[row, 2]*x2 +

NUTATION_YTERM_ARRAY[row, 3]*x3 +

NUTATION_YTERM_ARRAY[row, 4]*x4

)

delta_psi_sum += (a + b * julian_ephemeris_century) * np.sin(arg)

delta_eps_sum += (c + d * julian_ephemeris_century) * np.cos(arg)

delta_psi = delta_psi_sum*1.0/36000000

delta_eps = delta_eps_sum*1.0/36000000

# seems like we ought to be able to return a tuple here instead

# of resorting to `out`, but returning a UniTuple from this

# function caused calculations elsewhere to give the wrong result.

# very difficult to investigate since it did not occur when using

# object mode. issue was observed on numba 0.56.4

out[0] = delta_psi

U = 1.0*julian_ephemeris_millennium/10

e0 = (84381.448 - 4680.93 * U - 1.55 * U**2

+ 1999.25 * U**3 - 51.38 * U**4 - 249.67 * U**5

- 39.05 * U**6 + 7.12 * U**7 + 27.87 * U**8

+ 5.79 * U**9 + 2.45 * U**10)

e0 = mean_ecliptic_obliquity

deleps = obliquity_nutation

e = e0*1.0/3600 + deleps

true_ecliptic_obliquity):

v = mean_sidereal_time + longitude_nutation * np.cos(

np.radians(true_ecliptic_obliquity))

true_ecliptic_obliquity,

geocentric_latitude):

true_ecliptic_obliquity_rad = np.radians(true_ecliptic_obliquity)

apparent_sun_longitude_rad = np.radians(apparent_sun_longitude)

num = (np.sin(apparent_sun_longitude_rad)

* np.cos(true_ecliptic_obliquity_rad)

- np.tan(np.radians(geocentric_latitude))

* np.sin(true_ecliptic_obliquity_rad))

alpha = np.degrees(np.arctan2(num, np.cos(apparent_sun_longitude_rad)))

geocentric_latitude):

geocentric_latitude_rad = np.radians(geocentric_latitude)

true_ecliptic_obliquity_rad = np.radians(true_ecliptic_obliquity)

delta = np.degrees(np.arcsin(np.sin(geocentric_latitude_rad) *

np.cos(true_ecliptic_obliquity_rad) +

np.cos(geocentric_latitude_rad) *

np.sin(true_ecliptic_obliquity_rad) *

np.sin(np.radians(apparent_sun_longitude))))

sun_right_ascension):

"""Measured westward from south"""

H = apparent_sidereal_time + observer_longitude - sun_right_ascension

local_hour_angle, geocentric_sun_declination):

equatorial_horizontal_parallax_rad = \

np.radians(equatorial_horizontal_parallax)

local_hour_angle_rad = np.radians(local_hour_angle)

num = (-xterm * np.sin(equatorial_horizontal_parallax_rad)

* np.sin(local_hour_angle_rad))

denom = (np.cos(np.radians(geocentric_sun_declination))

- xterm * np.sin(equatorial_horizontal_parallax_rad)

* np.cos(local_hour_angle_rad))

delta_alpha = np.degrees(np.arctan2(num, denom))

equatorial_horizontal_parallax,

parallax_sun_right_ascension,

local_hour_angle):

geocentric_sun_declination_rad = np.radians(geocentric_sun_declination)

equatorial_horizontal_parallax_rad = \

np.radians(equatorial_horizontal_parallax)

num = ((np.sin(geocentric_sun_declination_rad) - yterm

* np.sin(equatorial_horizontal_parallax_rad))

* np.cos(np.radians(parallax_sun_right_ascension)))

denom = (np.cos(geocentric_sun_declination_rad) - xterm

* np.sin(equatorial_horizontal_parallax_rad)

* np.cos(np.radians(local_hour_angle)))

delta = np.degrees(np.arctan2(num, denom))

topocentric_sun_declination,

topocentric_local_hour_angle

):

observer_latitude_rad = np.radians(observer_latitude)

topocentric_sun_declination_rad = np.radians(topocentric_sun_declination)

e0 = np.degrees(np.arcsin(

np.sin(observer_latitude_rad)

* np.sin(topocentric_sun_declination_rad)

+ np.cos(observer_latitude_rad)

* np.cos(topocentric_sun_declination_rad)

* np.cos(np.radians(topocentric_local_hour_angle))))

topocentric_elevation_angle_wo_atmosphere,

atmos_refract):

# switch sets delta_e when the sun is below the horizon

switch = topocentric_elevation_angle_wo_atmosphere >= -1.0 * (

0.26667 + atmos_refract)

delta_e = ((local_pressure / 1010.0) * (283.0 / (273 + local_temp))

* 1.02 / (60 * np.tan(np.radians(

topocentric_elevation_angle_wo_atmosphere

+ 10.3 / (topocentric_elevation_angle_wo_atmosphere

+ 5.11))))) * switch

atmospheric_refraction_correction):

e = (topocentric_elevation_angle_without_atmosphere

+ atmospheric_refraction_correction)

topocentric_sun_declination,

observer_latitude):

topocentric_local_hour_angle_rad = np.radians(topocentric_local_hour_angle)

observer_latitude_rad = np.radians(observer_latitude)

num = np.sin(topocentric_local_hour_angle_rad)

denom = (np.cos(topocentric_local_hour_angle_rad)

* np.sin(observer_latitude_rad)

- np.tan(np.radians(topocentric_sun_declination))

* np.cos(observer_latitude_rad))

gamma = np.degrees(np.arctan2(num, denom))

M = (280.4664567 + 360007.6982779 * julian_ephemeris_millennium

+ 0.03032028 * julian_ephemeris_millennium**2

+ julian_ephemeris_millennium**3 / 49931

- julian_ephemeris_millennium**4 / 15300

- julian_ephemeris_millennium**5 / 2000000)

longitude_nutation, true_ecliptic_obliquity):

E = (sun_mean_longitude - 0.0057183 - geocentric_sun_right_ascension +

longitude_nutation * np.cos(np.radians(true_ecliptic_obliquity)))

# limit between 0 and 360

E = E % 360

# convert to minutes

E *= 4

greater = E > 20

less = E < -20

other = (E <= 20) & (E >= -20)

E = greater * (E - 1440) + less * (E + 1440) + other * E

"""Loop through the time array and calculate the solar position"""

lat = loc_args[0]

lon = loc_args[1]

elev = loc_args[2]

pressure = loc_args[3]

temp = loc_args[4]

delta_t = loc_args[5]

atmos_refract = loc_args[6]

sst = loc_args[7]

esd = loc_args[8]

for i in range(unixtime.shape[0]):

utime = unixtime[i]

jd = julian_day(utime)

jde = julian_ephemeris_day(jd, delta_t)

jc = julian_century(jd)

jce = julian_ephemeris_century(jde)

jme = julian_ephemeris_millennium(jce)

R = heliocentric_radius_vector(jme)

if esd:

out[0, i] = R

continue

L = heliocentric_longitude(jme)

B = heliocentric_latitude(jme)

Theta = geocentric_longitude(L)

beta = geocentric_latitude(B)

x0 = mean_elongation(jce)

x1 = mean_anomaly_sun(jce)

x2 = mean_anomaly_moon(jce)

x3 = moon_argument_latitude(jce)

x4 = moon_ascending_longitude(jce)

l_o_nutation = np.empty((2,))

longitude_obliquity_nutation(jce, x0, x1, x2, x3, x4, l_o_nutation)

delta_psi = l_o_nutation[0]

delta_epsilon = l_o_nutation[1]

epsilon0 = mean_ecliptic_obliquity(jme)

epsilon = true_ecliptic_obliquity(epsilon0, delta_epsilon)

delta_tau = aberration_correction(R)

lamd = apparent_sun_longitude(Theta, delta_psi, delta_tau)

v0 = mean_sidereal_time(jd, jc)

v = apparent_sidereal_time(v0, delta_psi, epsilon)

alpha = geocentric_sun_right_ascension(lamd, epsilon, beta)

delta = geocentric_sun_declination(lamd, epsilon, beta)

if sst:

out[0, i] = v

out[1, i] = alpha

out[2, i] = delta

continue

m = sun_mean_longitude(jme)

eot = equation_of_time(m, alpha, delta_psi, epsilon)

H = local_hour_angle(v, lon, alpha)

xi = equatorial_horizontal_parallax(R)

u = uterm(lat)

x = xterm(u, lat, elev)

y = yterm(u, lat, elev)

delta_alpha = parallax_sun_right_ascension(x, xi, H, delta)

delta_prime = topocentric_sun_declination(delta, x, y, xi, delta_alpha,

H)

H_prime = topocentric_local_hour_angle(H, delta_alpha)

e0 = topocentric_elevation_angle_without_atmosphere(lat, delta_prime,

H_prime)

delta_e = atmospheric_refraction_correction(pressure, temp, e0,

atmos_refract)

e = topocentric_elevation_angle(e0, delta_e)

theta = topocentric_zenith_angle(e)

theta0 = topocentric_zenith_angle(e0)

gamma = topocentric_astronomers_azimuth(H_prime, delta_prime, lat)

phi = topocentric_azimuth_angle(gamma)

out[0, i] = theta

out[1, i] = theta0

out[2, i] = e

out[3, i] = e0

out[4, i] = phi

atmos_refract, numthreads, sst=False, esd=False):

"""Calculate the solar position using the numba compiled functions

# these args are the same for each thread

loc_args = np.array([lat, lon, elev, pressure, temp, delta_t,

atmos_refract, sst, esd])

# construct dims x ulength array to put the results in

ulength = unixtime.shape[0]

if sst:

dims = 3

elif esd:

dims = 1

else:

dims = 6

result = np.empty((dims, ulength), dtype=np.float64)

if unixtime.dtype != np.float64:

unixtime = unixtime.astype(np.float64)

if ulength < numthreads:

warnings.warn('The number of threads is more than the length of '

'the time array. Only using %s threads.'.format(ulength))

numthreads = ulength

if numthreads <= 1:

solar_position_loop(unixtime, loc_args, result)

return result

# split the input and output arrays into numthreads chunks

split0 = np.array_split(unixtime, numthreads)

split2 = np.array_split(result, numthreads, axis=1)

chunks = [[a0, loc_args, split2[i]] for i, a0 in enumerate(split0)]

# Spawn one thread per chunk

threads = [threading.Thread(target=solar_position_loop, args=chunk)

for chunk in chunks]

for thread in threads:

thread.start()

for thread in threads:

thread.join()

atmos_refract, numthreads, sst=False, esd=False):

"""Calculate the solar position assuming unixtime is a numpy array. Note

jd = julian_day(unixtime)

jde = julian_ephemeris_day(jd, delta_t)

jc = julian_century(jd)

jce = julian_ephemeris_century(jde)

jme = julian_ephemeris_millennium(jce)

R = heliocentric_radius_vector(jme)

if esd:

return (R, )

L = heliocentric_longitude(jme)