R: Plot the daily, monthly and annual aggregations of esimtated...
ETPlot
R Documentation
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.
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 plotting. For average plot the plotting period equals to the calculation period.
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 argument type is 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 argument type is Average, two plots are displayed in the following order
1) Monthly averaged daily estimations of evapotranspiration;
2) Annually averaged daily estiamtions of evapotranspriation.
Author(s)
Danlu Guo
See Also
ETComparison
Examples
# Use processed existing data set and constants from kent Town, Adelaide
data("processeddata")
data("constants")
# Call ET.Penman under the generic function ET
results <- ET.Penman(data, constants, ts="daily", 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")
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(Evapotranspiration)
Loading required package: zoo
Attaching package: 'zoo'
The following objects are masked from 'package:base':
as.Date, as.Date.numeric
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Evapotranspiration/ETPlot.Rd_%03d_medium.png", width=480, height=480)
> ### Name: ETPlot
> ### Title: Plot the daily, monthly and annual aggregations of esimtated
> ### evapotranspiration
> ### Aliases: ETPlot
> ### Keywords: plot evapotranspiration
>
> ### ** Examples
>
> # Use processed existing data set and constants from kent Town, Adelaide
> data("processeddata")
> data("constants")
>
> # Call ET.Penman under the generic function ET
> results <- ET.Penman(data, constants, ts="daily", solar="sunshine hours",
+ wind="yes", windfunction_ver = "1948", alpha = 0.08, z0 = 0.001)
Penman Open-water Evaporation
Evaporative surface: water, albedo = 0.08 ; roughness height = 0.001 m
Sunshine hour data have been used for calculating incoming solar radiation
Wind data have been used for calculating the Penman evaporation. Penman 1948 wind function has been used.
Timestep: daily
Units: mm
Time duration: 2001-03-01 to 2004-08-31
Basic stats
Mean: 4.87
Max: 12.79
Min: 1.04
>
>
> # 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")
>
>
>
>
>
> dev.off()
null device
1
>