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 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)
def get_evap_i(lat, elev, wind, srad, tmin, tmax, tavg, month):
if month == 1:
J = 15
else:
J = 15 + (month - 1) * 30
latitude = pyeto.deg2rad(lat)
atmosphericVapourPressure = pyeto.avp_from_tmin(tmin)
saturationVapourPressure = pyeto.svp_from_t(tavg)
ird = pyeto.inv_rel_dist_earth_sun(J)
solarDeclination = pyeto.sol_dec(J)
sha = [pyeto.sunset_hour_angle(l, solarDeclination) for l in latitude]
extraterrestrialRad = [
pyeto.et_rad(x, solarDeclination, y, ird)
for x, y in zip(latitude, sha)
]
clearSkyRad = pyeto.cs_rad(elev, extraterrestrialRad)
netInSolRadnet = pyeto.net_in_sol_rad(srad * 0.001, albedo=0.23)
netOutSolRadnet = pyeto.net_out_lw_rad(tmin, tmax, srad * 0.001,
clearSkyRad,
atmosphericVapourPressure)
netRadiation = pyeto.net_rad(netInSolRadnet, netOutSolRadnet)
tempKelvin = pyeto.celsius2kelvin(tavg)
windSpeed2m = pyeto.wind_speed_2m(wind, 10)
slopeSvp = pyeto.delta_svp(tavg)
atmPressure = pyeto.atm_pressure(elev)
psyConstant = pyeto.psy_const(atmPressure)
return pyeto.fao56_penman_monteith(netRadiation,
tempKelvin,
windSpeed2m,
saturationVapourPressure,
atmosphericVapourPressure,
slopeSvp,
psyConstant,
shf=0.0)
lat=pyeto.deg2rad(lat)
day=day+1
sunshine_hours=(BS["0"][j])
#Radiacion neta
sol_dec=pyeto.sol_dec(day)
sha=pyeto.sunset_hour_angle(lat,sol_dec)
daylight_hours=pyeto.daylight_hours(sha)
ird=pyeto.inv_rel_dist_earth_sun(day)
et_rad=pyeto.et_rad(lat, sol_dec, sha, ird)
sol_rad=pyeto.sol_rad_from_sun_hours(daylight_hours,sunshine_hours,et_rad)
ni_sw_rad=pyeto.net_in_sol_rad(sol_rad, albedo=0.23)
cs_rad=pyeto.cs_rad(altitude, et_rad)
svp_tmin=pyeto.svp_from_t(tmin)
svp_tmax=pyeto.svp_from_t(tmax)
avp=pyeto.avp_from_rhmean(svp_tmin, svp_tmax, rh_mean)
no_lw_rad=pyeto.net_out_lw_rad(tmink, tmaxk, sol_rad, cs_rad, avp)
net_rad=pyeto.net_rad(ni_sw_rad, no_lw_rad)
#Presion de vapor de saturacion
svp=pyeto.svp_from_t(t)
#Delta presion de vapor de saturacion
delta_svp=pyeto.delta_svp(t)
#Constante psicrométrica
atmos_pres=pyeto.atm_pressure(altitude)
psy=pyeto.psy_const(atmos_pres)
#Calculo ETo Fao Penman Monteith
ETo=pyeto.fao56_penman_monteith(net_rad, tk, ws, svp, avp, delta_svp, psy, shf=0.0)
ETO.append(ETo)
if day>=365:
if Ano%4==0:
Ano=Ano+1
else:
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
coastal = True altitude = 147 rh_min = 13 rh_max = 88 ws = 1.3 tmean = pyeto.daily_mean_t(tmin, tmax) atmos_pres = pyeto.atm_pressure(altitude) psy = pyeto.psy_const(atmos_pres) # Humidity svp_tmin = pyeto.svp_from_t(tmin) svp_tmax = pyeto.svp_from_t(tmax) delta_svp = pyeto.delta_svp(tmean) svp = pyeto.mean_svp(tmin, tmax) avp = pyeto.avp_from_rhmin_rhmax(svp_tmin, svp_tmax, rh_min, rh_max) # Radiation sol_dec = pyeto.sol_dec(day_of_year) sha = pyeto.sunset_hour_angle(latitude, sol_dec) ird = pyeto.inv_rel_dist_earth_sun(day_of_year) et_rad = pyeto.et_rad(latitude, sol_dec, sha, ird) cs_rad = pyeto.cs_rad(altitude, et_rad) sol_rad = pyeto.sol_rad_from_t(et_rad, cs_rad, tmin, tmax, coastal) ni_sw_rad = pyeto.net_in_sol_rad(sol_rad) no_lw_rad = pyeto.net_out_lw_rad(pyeto.celsius2kelvin(tmin), pyeto.celsius2kelvin(tmax), sol_rad, cs_rad, avp) net_rad = pyeto.net_rad(ni_sw_rad, no_lw_rad) eto = pyeto.fao56_penman_monteith(net_rad, pyeto.celsius2kelvin(tmean), ws, svp, avp, delta_svp, psy) print eto
# Humidity
svp_tmin = pyeto.svp_from_t(tmin)
svp_tmax = pyeto.svp_from_t(tmax)
delta_svp = pyeto.delta_svp(tmean)
svp = pyeto.mean_svp(tmin, tmax)
avp = pyeto.avp_from_rhmin_rhmax(svp_tmin, svp_tmax, rh_min, rh_max)
# Radiation
sol_dec = pyeto.sol_dec(day_of_year)
sha = pyeto.sunset_hour_angle(latitude, sol_dec)
ird = pyeto.inv_rel_dist_earth_sun(day_of_year)
et_rad = pyeto.et_rad(latitude, sol_dec, sha, ird)
cs_rad = pyeto.cs_rad(altitude, et_rad)
sol_rad = pyeto.sol_rad_from_t(et_rad, cs_rad, tmin, tmax, coastal)
ni_sw_rad = pyeto.net_in_sol_rad(sol_rad)
no_lw_rad = pyeto.net_out_lw_rad(pyeto.celsius2kelvin(tmin), pyeto.celsius2kelvin(tmax), sol_rad, cs_rad, avp)
net_rad = pyeto.net_rad(ni_sw_rad, no_lw_rad)
eto = pyeto.fao56_penman_monteith(net_rad, pyeto.celsius2kelvin(tmean), ws, svp, avp, delta_svp, psy)
print eto
Forever
check the time
if not running
if eto is <0 then
turn on water
var ontime = time.time()
else
time running = ontime - time.time()
if time running > eto then
turn off water