plotRF {SoPhy}R Documentation

Plotting random fields

Description

the function plots random fields including stones and roots using the definition of the horizons in xswms2d.

Usage

  plotRF(h, col.txt='black', col.stones='white',
         col.rf=if (is.null(h$col.rf)) rainbow(100) else h$col.rf, 
         titl=TRUE, line = 0.7, lim=1, quadratic=TRUE,
         cex=1, cex.leg=0.8, pch=5,
         cex.pch = max(0.05, 3/min(length(h$grid.x), length(h$grid.y))),
         root.col=grey(1 - (i-1) / (1.6 * length(h$plants))),
         what=c('Root.RF', 'Stone.RF', 'RF'),
         transf=c("K", "H", "theta", "none"), legend = TRUE, ylim, zlim)
 

Arguments

h a list of the same format as the output of xswms2d
col.txt colour of the title
col.stones colour of the stones and the air
col.rf colour spectrum of the random field; plants are always grey
titl logical. If TRUE a title is plotted using the colour col.txt
line place where the titel is plottes, see title.
lim value in [0,1] that is used only if zlim is missing; quantile for the upper bound zlim[2]; the value of lim should be less than 1 if the random field contains some extreme values, since the colour scale is linear
quadratic logical. If TRUE the figure matrix is enlarged and filled with NA symmetrically in x-direction or on the bottom in y-direction such the matrix (and the figure) become quadratic
cex A numerical value giving the amount by which the title text, the axes and the labels of the axes should be scaled relative to the default
cex.leg A numerical value giving the amount by which the legend text should be scaled relative to the default
pch the symbol for a root segment
cex.pch size of the symbol
root.col the colours for a root segments
what chooses between the random field that includes modifications by roots and stones or by stones allone, or the random field without modifications
transf Kind of random field: saturated hydraulic conductivity "K", pressure head "H", the water content "theta", or the original Gaussian random field ("none")
legend logical. If (TRUE) a legend is added.
ylim missing or vector of two components. If missing the range of the simulation is used.
zlim missing or vector of two components. Values below zlim[1] are plotted in white, values above zlim[2] are plotted also in white. If missing then zlim[1] is the minimal value and zlim[2] is calculated by means of lim.

Details

If transf="K" then the function plots the random field of hydraulic conductivity, i.e., the Gaussian random field after h$miller.link and h$millerK have been applied.

Value

the matrix of the random field values, always excluding the plants, if no error has occured and an error message otherwise.

Author(s)

Martin Schlather, martin.schlather@math.uni-goettingen.de http://www.stochastik.math.uni-goettingen.de/institute

See Also

plotWater, xswms2d

Examples

   ## see simulateHorizons
  

[Package SoPhy version 1.0.34 Index]