Plot a RMODFLOW ghb object

# S3 method for ghb
rmf_plot(
  ghb,
  dis,
  kper = NULL,
  variable = "id",
  i = NULL,
  j = NULL,
  k = NULL,
  active_only = TRUE,
  fun = ifelse(variable == "bhead", mean, sum),
  add = FALSE,
  ...
)

Arguments

ghb

an RMODFLOW ghb object

dis

a RMODFLOW dis object

kper

integer specifying the stress-period to plot

variable

single character or numeric indicating which column of ghb$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 variable 'bhead' 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