Plot a MODFLOW drawdown file object

# S3 method for ddn
rmf_plot(
  ddn,
  dis,
  i = NULL,
  j = NULL,
  k = NULL,
  l = NULL,
  nstp = NULL,
  kper = NULL,
  kstp = NULL,
  ...
)

Arguments

ddn	`RMODFLOW` ddn object
dis	`RMODFLOW` dis object
i	row number to plot
j	column number to plot
k	layer number to plot
l	time step number to plot; defaults to plotting the final time step.
nstp	integer specifying the flow time step to plot. See details.
kper	integer specifying the stress-period. Use in conjunction with kstp. See details.
kstp	integer specifying the time step of kper. Use in conjunction with kper. See details.
...	additional parameters passed to `rmf_plot.rmf_4d_array`

Value

ggplot2 object or layer

Details

There are three ways to specify which time step to plot. The l argument can be specified which represents the dimension index. This can be smaller than the total number of time steps since output may not be written for all steps. The second option is to specify nstp which is the total time step number. Since output at every time step is not always written, a specified nstp value might not be available in the array. The third option is to specify both kper & kstp which specify the stress-period and corresponding time step in that stress-period. A negative kstp will plot the final time step of the stress-period.

If no output is written for the specified time step, as controlled by the Output Control file in MODFLOW, an error is thrown.