plotDiurnal: Plot average diurnal patterns
In FREddyPro: Post-Processing EddyPro Full Output File

Description Usage Arguments Details Author(s) Examples

Calculating and plotting the diurnal pattern of any variable with the option to use mean or both mean and median. There is also an option to plot quantiles.

plotDiurnal(data, xVar, hourX, dataY = NULL, yVar = NULL, hourY = NULL,
median = FALSE, quantiles = FALSE, probs = c(0.05, 0.95),
stdev = FALSE, sterr = FALSE, legend = FALSE, legendSide = NULL,
legendText = NULL, type = "l", lty = c(1, 4), col = c(1, 2),
cex.legend = 1, horiz.legend = TRUE, ylab = "Mean", xlab = "Hour", ...)

`data`	The data frame to plot
`xVar`	The variable to plot
`hourX`	The hour variable associated with the X variable
`dataY`	A second data frame to plot (optional)
`yVar`	A second optional variable to plot.
`hourY`	The hour variable associated with the Y variable, also optional.
`median`	Logical. If TRUE the median is also plotted with the mean.
`quantiles`	Logical. If TRUE the quantiles will be plotted as a colour band
`probs`	An array of the quantiles to use. The default is 0.05 and 0.95 (5th and 95th).
`stdev`	Logical. If True the standard deviation is plotted as a colour band
`sterr`	Logical. If True the standard error is plotted as a colour band
`legend`	Logical. If TRUE a legend will be plotted.
`legendSide`	The side of the plot to place the legend.
`legendText`	The text to use for legend.
`type`	The type of plot. The default is a line plot.
`lty`	The type of line when a line plot is used.
`col`	The colour to use for the lines or points. It can be a single values or an array.
`cex.legend`	The size of the legend labels.
`horiz.legend`	A logical parameter. If TRUE the legend is will be placed horizontally.
`ylab`	The y-axis label.
`xlab`	The x-axis label.
`...`	Further graphical options.

The function first aggregates the given variable to half hourly and plots by default the mean. The median can also be plotted with the mean. There is an optional graphical output for quantiles. By default is the 5th and 95th but other quantiles can be used. There are two separate hourly inputs to the functions to allow plotting variables from different data frames.

Georgios Xenakis

## Load the data
data(fluxes)

## Clean and despike the fluxes
fluxes=cleanFluxes(fluxes,sdCor=TRUE,sdTimes=3,distCor=TRUE,timesList=3,
thresholdList=list(H=c(-100,1000),LE=c(-100,1000)))	

## Create timestamp
fluxes=createTimestamp(fluxes)

## Plot the diurnal patterns
plotDiurnal(fluxes, xVar="H", hourX="hour",yVar="LE",hourY="hour",
median=TRUE,stdev=TRUE,ylim=c(-100,400),lty=c(1,2,3,4),col=c(1,2,3,4),
legend=TRUE,legendSide='topleft',legendText=c("sensible heat","latent heat"),
lwd=2,ylab="Heat (W/m"^2~")",main="Harwood forest",cex.legend=0.9,
horiz.legend=FALSE)

abline(h=0,lty=2)