Maeswrap (version 1.7)

plotuspar: Plot the understorey PAR points

Description

Reads the 'uspar.dat' file in the current working directory, and plots the incident or absorbed (or diffuse, or direct) PAR at the understorey points. Either produces a plot, or makes a pdf with a plot for each hour of the selected day.

Reads the point-wise output file when the understorey was simulated.

Usage

plotuspar(what = c("PARbeam", "PARtotal", "PARdiffuse", "APAR"), dataset = NULL, day = 1, hour = NA, xlim = NULL, ylim = NULL, makepdf = FALSE, outputfile = "aparunderstorey.pdf", scaleeach = TRUE, addNarrow = TRUE)
readuspar(filename = "uspar.dat")

Arguments

what

Either 'diff', 'apar', 'ipar', or 'beam' (the default).

dataset

If left alone, reads the uspar dataset.

day

Which day to use, if left alone uses the first day only.

hour

Which hour to plot. If left alone, makes a plot for each hour.

xlim,ylim

X- and Y-axis limits.

makepdf

Logical. If TRUE, produces a pdf in the working directory.

outputfile

Name of the pdf file.

scaleeach

Logical. Rescale grey scale for each plot, or same for all hours?

addNarrow

Logical. Add an arrow pointing North.

filename

The understorey file

Value

A lattice device, or a pdf.

Details

If addNarrow is TRUE, attempts to read the trees.dat file in the current working directory as well. Prints a warning when this file cannot be opened.

Examples

## Not run: 
# 
# # Plot one hour of the first day, showing incident PAR on understorey:
# plotuspar("ipar", day=1,hour=12,makepdf=FALSE)
# 
# # Make pdf of the whole day, plotting beam radiation:
# plotuspar("beam", day=1, outputfile="beam uspar")
# 
# ## End(Not run)