def test_fao56_penman_monteith(self):
# Test based on example 17, p.110 (monthly calc) and
# example 18, p.113 (daily calc) of FAO paper
# Monthly ETo:
# Note, due to rounding errors in the FAO's example we can only
# test to 1 decimal place here!
eto = pyeto.fao56_penman_monteith(
net_rad=14.33,
t=convert.celsius2kelvin(30.2),
ws=2.0,
svp=4.42,
avp=2.85,
delta_svp=0.246,
psy=0.0674,
shf=0.14,
)
self.assertAlmostEqual(eto, 5.72, 1)
# Daily ETo:
# (Rn, t, ws, es, ea, delta_es, psy, shf=0.0)
eto = pyeto.fao56_penman_monteith(net_rad=13.28,
t=convert.celsius2kelvin(16.9),
ws=2.078,
svp=1.997,
avp=1.409,
delta_svp=0.122,
psy=0.0666)
self.assertAlmostEqual(eto, 3.9, 1)
def test_fao56_penman_monteith(self):
# Test based on example 17, p.110 (monthly calc) and
# example 18, p.113 (daily calc) of FAO paper
# Monthly ETo:
# Note, due to rounding errors in the FAO's example we can only
# test to 1 decimal place here!
eto = pyeto.fao56_penman_monteith(
net_rad=14.33,
t=convert.celsius2kelvin(30.2),
ws=2.0,
svp=4.42,
avp=2.85,
delta_svp=0.246,
psy=0.0674,
shf=0.14,
)
self.assertAlmostEqual(eto, 5.72, 1)
# Daily ETo:
# (Rn, t, ws, es, ea, delta_es, psy, shf=0.0)
eto = pyeto.fao56_penman_monteith(
net_rad=13.28,
t=convert.celsius2kelvin(16.9),
ws=2.078,
svp=1.997,
avp=1.409,
delta_svp=0.122,
psy=0.0666
)
self.assertAlmostEqual(eto, 3.9, 1)
def test_net_out_lw_rad(self):
# Test based on example 11, p.87 of FAO paper
lwrad = pyeto.net_out_lw_rad(tmin=convert.celsius2kelvin(19.1),
tmax=convert.celsius2kelvin(25.1),
sol_rad=14.5,
cs_rad=18.8,
avp=2.1)
self.assertAlmostEqual(lwrad, 3.5, 1)
def test_net_out_lw_rad(self):
# Test based on example 11, p.87 of FAO paper
lwrad = pyeto.net_out_lw_rad(
tmin=convert.celsius2kelvin(19.1),
tmax=convert.celsius2kelvin(25.1),
sol_rad=14.5,
cs_rad=18.8,
avp=2.1
)
self.assertAlmostEqual(lwrad, 3.5, 1)
