Plot a RMODFLOW riv object

# S3 method for riv
rmf_plot(
  riv,
  dis,
  kper = NULL,
  variable = "id",
  i = NULL,
  j = NULL,
  k = NULL,
  active_only = TRUE,
  fun = ifelse(variable %in% c("stage", "rbot"), mean, sum),
  add = FALSE,
  ...
)

Arguments

riv

an RMODFLOW riv object

dis

a RMODFLOW dis object

kper

integer specifying the stress-period to plot

variable

single character or numeric indicating which column of riv$data to plot. Defaults to 'id', which plots the locations of the cells.

i

row number to plot

j

column number to plot

k

layer number to plot

active_only

logical; indicating if only the active cells should be plotted. Non-active cells are set to NA. Defaults to TRUE.

fun

function to compute values in the case multiple values are defined for the same MODFLOW cell. Typically either mean or sum. Defaults to mean for variables 'stage' & 'rbot' and sum for variable 'cond'.

add

logical; if TRUE, provide ggplot2 layers instead of object, or add 3D plot to existing rgl device; defaults to FALSE

...

additional arguments passed to rmf_plot.rmf_3d_array

Value

ggplot2 object or layer; if plot3D is TRUE, nothing is returned and the plot is made directly