Plotting functions {FESDIA}R Documentation

Utility functions for plotting FESDIA dynamic or steady-state model output.

Description

image2D generates a 2-D image (x=time, y=depth) of the dynamic output.

matplot1D plots how profiles change with time. These plots will have the depth on the y-axis, and will plot all profiles.

plot plots steady-state profiles, with the depth on the y-axis.

Usage


## S3 method for class 'FESDIAdyn'
image2D(z, which, ylim = c(20, 0), 
   colkey = list(cex.clab = 0.8, line.clab = 0.5, cex.axis = 0.8), ...) 

## Default S3 method:
matplot1D(z, ...) 
## S3 method for class 'FESDIAdyn'
matplot1D(z, which, ylim = c(20, 0), 
  type = "l", col = "grey",  lty = 1, ...) 

matplot1D (z, ...) 

## S3 method for class 'FESDIAstd'
plot(x, ..., which, ylim = c(20, 0)) 
 

Arguments

z

object of class FESDIAdyn generated by FESDIAdyna.

x

object of class FESDIAstd generated by FESDIAsolve.

which

The name(s) of the 1-dimensional output to be plotted versus time.

ylim

The ranges of the y-axis. default it to have the y-axis extend downward, from 0 to 20 cm

colkey

A list with specifications of the color key. The default is to have a smaller title (cex.clab) and axis labels (cex.axis), the key title positioned close to the color key (line.clab), See colkey

type, col, lty

The type of plot, the color and type of the lines used for plotting

...

Any argument passed to the functions image2D and matplot.1D .

Author(s)

Karline Soetaert

Examples


# =============================
# Plotting dynamic output
# =============================

  out <- FESDIAdyna(CfluxForc = list(amp = 1))
  image2D(out, which = 1:9)


  image2D(out, which = "O2", ylim = c(5, 0), mfrow = c(1,2))
  matplot1D(out, which = 3, col = "grey", ylim = c(5,0), mfrow = NULL) 
  O2mean <- FESDIA1D(out, "O2")
  head(O2mean)
  lines(O2mean$O2, O2mean$x) 
  
  
# The same but more elaborate
  O2 <- subset(out, which = "O2")
  image2D(O2, x = out[,"time"], y = FESDIAdepth(out), ylim = c(10,0))

# =============================
# Plotting steady-state output
# =============================

  std <- FESDIAsolve()
  plot(std, which = c("TOC", "O2"))

[Package FESDIA version 1.0 Index]