1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
|
import pysolar
from pysolar import radiation, solar
from pysolar import numeric as math
import datetime
import numpy as np
from nose.tools import raises, assert_equal
@raises(TypeError)
def test_fail_with_math():
pysolar.use_math()
lat = np.array([45., 40.])
lon = np.array([3., 4.])
time = datetime.datetime(2018, 5, 8, 12, 0, 0, tzinfo=datetime.timezone.utc)
solar.get_altitude(lat, lon, time)
def test_scalar_with_math():
pysolar.use_math()
lat = 45.
lon = 3.
time = datetime.datetime(2018, 5, 8, 12, 0, 0, tzinfo=datetime.timezone.utc)
print(solar.get_altitude(lat, lon, time))
print(solar.get_azimuth(lat, lon, time))
def test_scalar_with_numpy():
pysolar.use_numpy()
lat = 50.63
lon = 3.05
time = datetime.datetime(2018, 5, 8, 12, 0, 0, tzinfo=datetime.timezone.utc)
print(solar.get_altitude(lat, lon, time))
print(solar.get_azimuth(lat, lon, time))
def test_with_fixed_time():
""" get_altitude and get_azimuth, with scalar date """
pysolar.use_numpy()
lat = np.array([45., 40.])
lon = np.array([3., 4.])
time = datetime.datetime(2018, 5, 8, 12, 0, 0, tzinfo=datetime.timezone.utc)
print(solar.get_altitude(lat, lon, time))
print(solar.get_azimuth(lat, lon, time))
print(solar.get_altitude_fast(lat, lon, time))
print(solar.get_azimuth_fast(lat, lon, time))
def test_with_fixed_position():
""" get_altitude and get_azimuth, with scalar position """
pysolar.use_numpy()
lat = 50.
lon = 3.
time = np.array(['2018-05-08T12:15:00',
'2018-05-08T15:00:00'], dtype='datetime64')
print(solar.get_altitude_fast(lat, lon, time))
print(solar.get_azimuth_fast(lat, lon, time))
def test_datetime_operations():
d0 = datetime.datetime(2018,5,8,12,0,0)
d1 = np.array(d0)
assert_equal(math.tm_yday_math(d0),
math.tm_yday_numpy(d1))
assert_equal(math.tm_hour_math(d0),
math.tm_hour_numpy(d1))
assert_equal(math.tm_min_math(d0),
math.tm_min_numpy(d1))
def test_numpy():
""" get_altitude and get_azimuth, with lat, lon and date arrays """
pysolar.use_numpy()
lat = np.array([45., 40.])
lon = np.array([3., 4.])
time = np.array(['2018-05-08T12:15:00',
'2018-05-08T15:00:00'], dtype='datetime64')
print(solar.get_altitude_fast(lat, lon, time))
print(solar.get_azimuth_fast(lat, lon, time))
def test_numpy_radiation():
"""
get_radiation_direct with lat, lon, and date as arrays
"""
pysolar.use_numpy()
lat = np.array([45., 40., 40.])
lon = np.array([3., 4., 3.])
time = np.array([
'2018-05-08T12:15:00',
'2018-05-08T15:00:00',
'2018-05-08T03:00:00',
], dtype='datetime64')
altitude = solar.get_altitude_fast(lat, lon, time)
rad_results = radiation.get_radiation_direct(time, altitude)
assert rad_results[2] == 0
print(rad_results)
|
