rmf_create_hob creates an RMODFLOW hob object

rmf_create_hob(
  locations,
  time_series,
  dis,
  hydraulic_conductivity = array(1, dim = c(dis$nrow, dis$ncol, dis$nlay)),
  prj = rmf_get_prj(dis),
  iuhobsv = 665,
  hobdry = -888,
  noprint = FALSE,
  tomulth = 1,
  itt = 1,
  unique_obsnam = FALSE
)

Arguments

locations

data.frame or sf POINT object with the name, x and y coordinates (if data.frame), and top and bottom of the monitoring well filter

time_series

data.frame with the name of the monitoring well filter and the observation time and head

dis

RMODFLOW dis object

hydraulic_conductivity

3d array with hydraulic conductivity values, for calculating the weights for a multilayer observations

prj

RMODFLOW prj object; provide if coordinates in locations are real world coordinates

iuhobsv

file unit (> 0) for saving observation data in a file; the file unit should be included as DATA type in the name file; defaults to 665

hobdry

value of the simulated equivalent for a dry cell used in the observation output file

noprint

logical; if TRUE, the input and output are not printed in the listing file

itt

integer flag, either a single value or a value for every observation well as defined in locations. Identifies whether hydraulic heads (1; default) or changes in hydraulic heads (2) are used as observations. See details.

unique_obsnam

logical; should an ID number be added to obsnam for filters with observations at different times prefixed with an underscore? Defaults to FALSE.

tomult

time offset multiplier, to convert the time unit in time_series to the MODFLOW time unit

Value

Object of class hob

Details

Regardless of the value of itt, observed heads should be specified in time_series. MODFLOW will calculate the changes in hydraulic head internally. The locations data.frame must have columns named name, x, y, top, bottom. The time_series data.frame must have columns named name, time, head.

See also