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)