Learn R Programming

Evapotranspiration (version 1.7)

ETPlot: Plot the daily, monthly and annual aggregations of esimtated evapotranspiration

Description

Produce plot of aggregated estimations of evapotranspiration in daily, monthly and annual steps, or averaged daily estimations in monthly or annual steps.

Usage

ETPlot(results, type = "Aggregation", OBS, OBSplot, Sdate = time(results$ET.Daily)[1], 
Edate = time(results$ET.Daily)[length(results$ET.Daily)])

Arguments

results
A list named "results" which has been derived from function ET..().
type
A character string of either "Aggregation" or "Average" to indicate the type of plot required. The default is "aggregation". For aggregation plot the user can define the start and end date of plotting or by default using the calculation period for plottin
OBS
A list named "OBS" which has been derived from function ReadOBSEvaporation().
OBSplot
Must be eith TRUE or FALSE. TRUE indicates that the observed evaporation will be plotted together with the estimations and FALSE indicates that the observations will not be shown on the plots.
Sdate
Only used when type = "Aggregation" to define the start date for the plotting windows, the default is the first day for the estimate evapotranspiration, but can be defined by user in the format "YYYY-MM-DD".
Edate
Only used when type = "Aggregation" to define the end date for the plotting windows, the default is the last day for the estimate evapotranspiration, but can be defined by user in the format "YYYY-MM-DD".

Value

  • If type = "Aggregation", three plots are displayed in the following order (the next one appears after pressing enter): 1) Daily evapotranspiration estimates; 2) Monthly evapotranspiration estimates aggregated from daily estimates; 3) Annual evapotranspiration estimates aggregated from daily estimates. If type = "Average", two plots are displayed in the following order 1) Monthly averaged daily estimations of evapotranspiration; 2) Annually averaged daily estiamtions of evapotranspriation.

Examples

Run this code
# 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