# Use processed existing data set from kent Town, Adelaide
data("processeddata")
data("constants")
# Set data class to be "Penman" to call function
funname <- "Penman"
class(data) <- funname
# Call generic function Evapotranspiration(data, constants, ...) with class "Penman"
results <- ET(data, constants, solar="sunshine hours", wind=
"yes", windfunction_ver = "1948", alpha = 0.08, z0 = 0.001)
# Read evaporation data
data("E_OBS")
OBS <- ReadOBSEvaporation(E_OBS, data)
# Plot the aggregation of estimated Penman open-water evaporation with observed evaporation
ETPlot(results, type = "Aggregation", OBS, OBSplot = TRUE, Sdate = "2001-05-01",
Edate = "2004-05-01")
Run the code above in your browser using DataLab