def test_daily_grid(self):
# We use a 2x1 grid, where point 1, 1 is the same as Example 18, and
# point 1, 2 has some different values.
unit_converters = {
# Eq. 47 p. 56
'wind_speed': lambda x: x * 4.87 / math.log(67.8 * 10 - 5.42),
}
pm = PenmanMonteith(albedo=0.23,
elevation=100,
latitude=50.8,
step_length=timedelta(days=1),
unit_converters=unit_converters,
)
result = pm.calculate(temperature_max=np.array([21.5, 28]),
temperature_min=np.array([12.3, 15]),
humidity_max=np.array([84, 70]),
humidity_min=np.array([63, 60]),
wind_speed=np.array([2.78, 3]),
sunshine_duration=np.array([9.25, 9]),
adatetime=date(2014, 7, 6))
np.testing.assert_almost_equal(result, np.array([3.9, 4.8]),
decimal=1)
# Same thing with solar radiation instead of sunshine duration
result = pm.calculate(temperature_max=np.array([21.5, 28]),
temperature_min=np.array([12.3, 15]),
humidity_max=np.array([84, 70]),
humidity_min=np.array([63, 60]),
wind_speed=np.array([2.78, 3]),
solar_radiation=np.array([22.07, 21.62]),
adatetime=date(2014, 7, 6))
np.testing.assert_almost_equal(result, np.array([3.9, 4.8]),
decimal=1)
def test_daily(self):
# Apply Allen et al. (1998) Example 18 page 72.
unit_converters = {
# Eq. 47 p. 56
'wind_speed': lambda x: x * 4.87 / math.log(67.8 * 10 - 5.42),
}
pm = PenmanMonteith(albedo=0.23,
elevation=100,
latitude=50.8,
step_length=timedelta(days=1),
unit_converters=unit_converters,
)
result = pm.calculate(temperature_max=21.5,
temperature_min=12.3,
humidity_max=84,
humidity_min=63,
wind_speed=2.78,
sunshine_duration=9.25,
adatetime=date(2014, 7, 6))
self.assertAlmostEqual(result, 3.9, places=1)
# We try the same calculation, but instead of sunshine duration we
# provide the solar radiation directly. Should get the same result.
result = pm.calculate(temperature_max=21.5,
temperature_min=12.3,
humidity_max=84,
humidity_min=63,
wind_speed=2.78,
solar_radiation=22.07,
adatetime=date(2014, 7, 6))
self.assertAlmostEqual(result, 3.9, places=1)
def test_hourly(self):
# Apply Allen et al. (1998) Example 19 page 75.
pm = PenmanMonteith(albedo=0.23,
nighttime_solar_radiation_ratio=0.8,
elevation=8,
latitude=16.217,
longitude=-16.25,
step_length=timedelta(hours=1)
)
result = pm.calculate(temperature=38,
humidity=52,
wind_speed=3.3,
pressure=101.3,
solar_radiation=2.450,
adatetime=datetime(2014, 10, 1, 15, 0,
tzinfo=senegal_tzinfo))
self.assertAlmostEqual(result, 0.63, places=2)
result = pm.calculate(temperature=28,
humidity=90,
wind_speed=1.9,
pressure=101.3,
solar_radiation=0,
adatetime=datetime(2014, 10, 1, 2, 30,
tzinfo=senegal_tzinfo))
self.assertAlmostEqual(result, 0.0, places=2)
def test_hourly_grid(self):
# We use a 2x1 grid, where point 1, 1 is the same as Example 19, and
# point 1, 2 has some different values.
pm = PenmanMonteith(albedo=0.23,
nighttime_solar_radiation_ratio=0.8,
elevation=8,
latitude=16.217,
longitude=np.array([-16.25, -15.25]),
step_length=timedelta(hours=1)
)
result = pm.calculate(temperature=np.array([38, 28]),
humidity=np.array([52, 42]),
wind_speed=np.array([3.3, 2.3]),
pressure=101.3,
solar_radiation=np.array([2.450, 1.450]),
adatetime=datetime(2014, 10, 1, 15, 0,
tzinfo=senegal_tzinfo))
np.testing.assert_almost_equal(result, np.array([0.63, 0.36]),
decimal=2)
def test_hourly_with_albedo_grid(self):
# Apply Allen et al. (1998) Example 19 page 75.
pm = PenmanMonteith(albedo=np.array([0.23]),
nighttime_solar_radiation_ratio=0.8,
elevation=8,
latitude=16.217,
longitude=-16.25,
step_length=timedelta(hours=1)
)
result = pm.calculate(temperature=38,
humidity=52,
wind_speed=3.3,
pressure=101.3,
solar_radiation=2.450,
adatetime=datetime(2014, 10, 1, 15, 0,
tzinfo=senegal_tzinfo))
# The following two lines could be written more simply like this:
# self.assertAlmostEqual(result, 0.63, places=2)
# However, it does not work properly on Python 3 because of a numpy
# issue.
self.assertEqual(result.size, 1)
self.assertAlmostEqual(result[0], 0.63, places=2)
