# 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 ET(data, constants, ...) with class "Penman"
results_Penman <- ET(data, constants, solar="sunshine hours", wind=
"yes", windfunction_ver = "1948", alpha = 0.26, z0 = 0.02)
# Set data class to be "PenmanMonteith" to call function
funname <- "PenmanMonteith"
class(data) <- funname
# Call generic function ET(data, constants, ...) with class "PenmanMonteith"
results_PenmanMonteith <- ET(data, constants, solar="sunshine hours", wind="yes", crop="short")
# Plot the estimated Penman open-water evaporation against average temperature
ETComparison(results_Penman, results_PenmanMonteith, type = "Monthly", ymin = NULL, ymax = NULL)Run the code above in your browser using DataLab