Plot a RMODFLOW hfb object

# S3 method for hfb
rmf_plot(
  hfb,
  dis,
  variable = "id",
  all_hfb = FALSE,
  i = NULL,
  j = NULL,
  k = NULL,
  type = "fill",
  active_only = TRUE,
  fun = ifelse(variable == "hydchr", sum, mean),
  prj = rmf_get_prj(dis),
  crs = NULL,
  size = 1,
  colour = "orange",
  crop = FALSE,
  add = FALSE,
  ...
)

Arguments

hfb

an RMODFLOW hfb object

dis

a RMODFLOW dis object

variable

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

all_hfb

logical; should all horizontal-flow barriers be plotted (TRUE) or only the ones that are active in the simulation (FALSE; default)

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 sum for variable 'hydchr'

...

additional arguments passed to rmf_plot.rmf_3d_array

Value

ggplot2 object or layer