def estimate_fao56_daily(self, day_of_year, temp_c, temp_c_min, temp_c_max,
tdew, elevation, latitude, rh, wind_m_s,
atmos_pres):
""" Estimate fao56 from weather """
sha = pyeto.sunset_hour_angle(pyeto.deg2rad(latitude),
pyeto.sol_dec(day_of_year))
daylight_hours = pyeto.daylight_hours(sha)
sunshine_hours = 0.8 * daylight_hours
ird = pyeto.inv_rel_dist_earth_sun(day_of_year)
et_rad = pyeto.et_rad(pyeto.deg2rad(latitude),
pyeto.sol_dec(day_of_year), sha, ird)
sol_rad = pyeto.sol_rad_from_sun_hours(daylight_hours, sunshine_hours,
et_rad)
net_in_sol_rad = pyeto.net_in_sol_rad(sol_rad=sol_rad, albedo=0.23)
cs_rad = pyeto.cs_rad(elevation, et_rad)
avp = pyeto.avp_from_tdew(tdew)
net_out_lw_rad = pyeto.net_out_lw_rad(
pyeto.convert.celsius2kelvin(temp_c_min),
pyeto.convert.celsius2kelvin(temp_c_max), sol_rad, cs_rad, avp)
net_rad = pyeto.net_rad(net_in_sol_rad, net_out_lw_rad)
eto = pyeto.fao56_penman_monteith(
net_rad=net_rad,
t=pyeto.convert.celsius2kelvin(temp_c),
ws=wind_m_s,
svp=pyeto.svp_from_t(temp_c),
avp=avp,
delta_svp=pyeto.delta_svp(temp_c),
psy=pyeto.psy_const(atmos_pres))
return eto
def get_data_from_WU():
###array to store the reports
wu_weather_reports = []
## today and last 6 days definition
day1 = now - datetime.timedelta(days=6)
day2 = now - datetime.timedelta(days=5)
day3 = now - datetime.timedelta(days=4)
day4 = now - datetime.timedelta(days=3)
day5 = now - datetime.timedelta(days=2)
day6 = now - datetime.timedelta(days=1)
day7 = now
#### convert dates to WU required format
days = {
'day1': day1.strftime('%Y%m%d'),
'day2': day2.strftime('%Y%m%d'),
'day3': day3.strftime('%Y%m%d'),
'day4': day4.strftime('%Y%m%d'),
'day5': day5.strftime('%Y%m%d'),
'day6': day6.strftime('%Y%m%d'),
'day7': day7.strftime('%Y%m%d')
}
### make API wather hisotry call for each day
for day in days:
url = 'http://api.wunderground.com/api/7c2ab99a0ccee978/history_{0}/q/95316.json'.format(
days[day])
headers = {'content-type': 'application/JSON; charset=utf8'}
response = requests.get(url, headers=headers)
data = json.loads(response.text)
#ETo calculation for the day using FAO-56 Penman-Monteith method
lat = pyeto.deg2rad(37.585652)
altitude = 38
julian_day = datetime.datetime.strptime(days.get(day),
'%Y%m%d').timetuple().tm_yday
sol_dec = pyeto.sol_dec(julian_day)
sha = pyeto.sunset_hour_angle(lat, sol_dec)
ird = pyeto.inv_rel_dist_earth_sun(julian_day)
### net radiation calculator
net_rad = pyeto.et_rad(lat, sol_dec, sha, ird)
temp_c = float(data["history"]["observations"][1]["tempm"])
temp_k = float(data["history"]["observations"][1]["tempi"])
humidity = float(data["history"]["observations"][1]["hum"])
dew_point = float(data["history"]["observations"][1]["dewptm"])
ws = float(data["history"]["observations"][1]["wspdm"])
#actual and saturated vapour pressure in kPH
svp = pyeto.svp_from_t(temp_c)
avp = pyeto.avp_from_tdew(dew_point)
delta_svp = pyeto.delta_svp(temp_c)
atm_pressure = pyeto.atm_pressure(altitude)
psy = pyeto.psy_const(atm_pressure)
#### the ETo plugin retun results in mm, it was converted to inched
ETo_in_mm = pyeto.fao56_penman_monteith(net_rad,
temp_k,
ws,
svp,
avp,
delta_svp,
psy,
shf=0.0)
ETo = ETo_in_mm * 0.039370
## insert eto value to day weather report
data["history"]["observations"][1].update(
{"ETo": "{0:.2f}".format(ETo)})
###add report to report collector array
wu_weather_reports.append(data["history"]["observations"][1])
#return all reports
return wu_weather_reports
def calculate_precipitation(self, d):
if "rain" in d:
self.rain_day = float(d["rain"])
if "snow" in d:
self.snow_day = float(d["snow"])
# def calculate_ev_fao56_factor(self, d):
dt = d['dt']
factor = 0.0
if dt > d['sunrise']:
if dt < d['sunset']:
factor = min(float(dt - d['sunrise'])/3600.0, 1.0)
else:
if dt > d['sunset']:
factor = (dt - d['sunrise'])/3600.0
if factor < 1.0:
factor = 1.0 - factor
return factor
#def estimate_fao56_hourly(self, day_of_year, temp_c, tdew, elevation, latitude, rh, wind_m_s, atmos_pres):
""" Estimate fao56 from weather """
sha = pyeto.sunset_hour_angle(pyeto.deg2rad(latitude),
pyeto.sol_dec(day_of_year))
daylight_hours = pyeto.daylight_hours(sha)
sunshine_hours = 0.8 * daylight_hours
ird = pyeto.inv_rel_dist_earth_sun(day_of_year)
et_rad = pyeto.et_rad(pyeto.deg2rad(latitude),
pyeto.sol_dec(day_of_year), sha, ird)
sol_rad = pyeto.sol_rad_from_sun_hours(daylight_hours, sunshine_hours,
et_rad)
net_in_sol_rad = pyeto.net_in_sol_rad(sol_rad=sol_rad, albedo=0.23)
cs_rad = pyeto.cs_rad(elevation, et_rad)
avp = pyeto.avp_from_tdew(tdew)
#not sure if I trust this net_out_lw_rad calculation here!
net_out_lw_rad = pyeto.net_out_lw_rad(temp_c-1, temp_c, sol_rad,
cs_rad, avp)
net_rad = pyeto.net_rad(net_in_sol_rad, net_out_lw_rad)
eto = pyeto.fao56_penman_monteith(
net_rad=net_rad,
t=pyeto.convert.celsius2kelvin(temp_c),
ws=wind_m_s,
svp=pyeto.svp_from_t(temp_c),
avp=avp,
delta_svp=pyeto.delta_svp(temp_c),
psy=pyeto.psy_const(atmos_pres))
return eto
#def calculate_fao56_hourly(self, d):
day_of_year = datetime.datetime.now().timetuple().tm_yday
T_hr = d['temp']
t_dew = float(d["dew_point"])
pressure = d['pressure']
RH_hr = d['humidity']
u_2 = d['wind_speed']
#print("CALCULATE_FAO56:")
#print("T_hr: {}".format(T_hr))
#print("t_dew: {}".format(t_dew))
#print("RH_hr: {}".format(RH_hr))
#print("u_2: {}".format(u_2))
#print("pressure: {}".format(pressure))
fao56 = self.estimate_fao56_hourly(day_of_year,
T_hr,
t_dew,
self.elevation,
LAT,
RH_hr,
u_2,
pressure)
return fao56
def get_data_from_WU():
###array to store the reports
wu_weather_reports = []
## today and last 6 days definition
day1 = now - datetime.timedelta(days=6)
day2 = now - datetime.timedelta(days=5)
day3 = now - datetime.timedelta(days=4)
day4 = now - datetime.timedelta(days=3)
day5 = now - datetime.timedelta(days=2)
day6 = now - datetime.timedelta(days=1)
day7 = now
#### convert dates to WU required format
days = {
'day1': day1.strftime('%Y%m%d'),
'day2': day2.strftime('%Y%m%d'),
'day3': day3.strftime('%Y%m%d'),
'day4': day4.strftime('%Y%m%d'),
'day5': day5.strftime('%Y%m%d'),
'day6': day6.strftime('%Y%m%d'),
'day7': day7.strftime('%Y%m%d')
}
### make API wather hisotry call for each day
for day in days:
url = 'http://api.wunderground.com/api/7c2ab99a0ccee978/history_{0}/q/95316.json'.format(days[day])
headers = {'content-type': 'application/JSON; charset=utf8'}
response = requests.get(url, headers=headers)
data = json.loads(response.text)
#ETo calculation for the day using FAO-56 Penman-Monteith method
lat = pyeto.deg2rad(37.585652)
altitude = 38
julian_day = datetime.datetime.strptime(days.get(day), '%Y%m%d').timetuple().tm_yday
sol_dec = pyeto.sol_dec(julian_day)
sha = pyeto.sunset_hour_angle(lat, sol_dec)
ird = pyeto.inv_rel_dist_earth_sun(julian_day)
### net radiation calculator
net_rad = pyeto.et_rad(lat, sol_dec, sha, ird)
temp_c = float(data["history"]["observations"][1]["tempm"])
temp_k = float(data["history"]["observations"][1]["tempi"])
humidity = float(data["history"]["observations"][1]["hum"])
dew_point = float(data["history"]["observations"][1]["dewptm"])
ws = float(data["history"]["observations"][1]["wspdm"])
#actual and saturated vapour pressure in kPH
svp = pyeto.svp_from_t(temp_c)
avp = pyeto.avp_from_tdew(dew_point)
delta_svp = pyeto.delta_svp(temp_c)
atm_pressure = pyeto.atm_pressure(altitude)
psy = pyeto.psy_const(atm_pressure)
#### the ETo plugin retun results in mm, it was converted to inched
ETo_in_mm = pyeto.fao56_penman_monteith(net_rad, temp_k, ws, svp, avp, delta_svp, psy, shf=0.0)
ETo = ETo_in_mm * 0.039370
## insert eto value to day weather report
data["history"]["observations"][1].update({"ETo": "{0:.2f}".format(ETo)})
###add report to report collector array
wu_weather_reports.append(data["history"]["observations"][1])
#return all reports
return wu_weather_reports
def test_avp_from_tdew(self):
# Test based on example 20, p.121 of FAO paper
avp = pyeto.avp_from_tdew(14.8)
self.assertAlmostEqual(avp, 1.68, 2)
def test_avp_from_tdew(self):
# Test based on example 20, p.121 of FAO paper
avp = pyeto.avp_from_tdew(14.8)
self.assertAlmostEqual(avp, 1.68, 2)