plotProfiles: Plot dynamic simulation output as vertical profiles

Description Usage Arguments Value Note Author(s)

Description

Plots the output of a dynamic 1D diagenetic model as vertical profiles.

Usage

1
2
3
4
plotProfiles(obj, item, obs, xlab = "", ylab = "depth", xlim = NULL,
  ylim = NULL, pos = "top", timeconv = function(x) {     ISOdatetime(1970,
  1, 1, 0, 0, 0) + x }, width = 4, height = 4, pointsize = 10,
  dir = ".", fmt = "svg", prefix = "profile", replace = TRUE, ...)

Arguments

obj

Object created with sim1D.create.

item

Name of a state variable or process (character string).

obs

Data frame of observations. See notes below.

xlab

Label for the time axis (x).

ylab

Label for the depth axis (y).

xlim

Limits for the value axis (x). Computed if NULL.

ylim

Limits for the depth axis (y). Computed if NULL.

pos

Keyword specifying the position of the legend.

timeconv

A function to convert numeric times to class POSIXct.

width

Plot width in inch.

height

Plot height in inch.

pointsize

Font size.

dir

Output directory for generated graphics.

fmt

Graphics format, either 'svg' or 'pdf'.

prefix

Prefix put before the value of name when constructing file names.

replace

If FALSE existing output files result in an error.

...

Additional arguments to be passed to legend.

Value

A vector holding the names of the generated graphics files.

Note

The data frame of observations passed as obs needs to have 3 columns named 'time', 'depth', and 'value'. If the times of observation do not match exactly with the model output, data for the nearest time points available in dyn will be plotted. Data in the 'time' column must be of the same type as the return value of timeconv. It is OK if the 'depth' and 'value' columns contain NA only and in that case, just the simulated data are plotted.

Author(s)

David Kneis david.kneis@tu-dresden.de


dkneis/diatools documentation built on May 15, 2019, 9:12 a.m.