MFimage: Plot MODFLOW variable matrix

In CJBarry/Rflow: Tools for analysing and writing input and output files for USGS MODFLOW

Description

Plot MODFLOW variable matrix

Usage

MFimage(mtx, gccs, grcs, zlim = "auto", col = colorRampPalette(c("blue",
  "white", "red"))(101L), show.range = identical(zlim, "angle"), ...)

Arguments

mtx	matrix; variable to plot
gccs	numeric [NCOL + 1]; column divider x co-ordinates, in ascending order
grcs	numeric [NROW + 1]; row divider y co-ordinates, in ascending order
zlim	numeric [2] or character string; min. and max. values to plot; alternatively use breaks argument as in graphics::image; additionally, three character string options are given: "auto" for automatic range calculation (pretty), "sym" for automatic range that is symmetric about 0 and "angle" for a range between -pi and pi
col	character []; colours
show.range	logical [1]; print computed zlim to screen? (in case that zlim = "angle", shows which directions the colours represent)
...	additional plotting arguments for image

Examples

# plot model domain
fnms <- system.file(c("rflow_mf_demo.bas",
                      "rflow_mf_demo.dis"),
                    package = "Rflow")

dis <- read.DIS(fnms[2L])
bas <- read.BAS(fnms[1L], dis)

# grid divider co-ordinates (the example has a regular grid)
gccs <- with(dis, seq(0, by = DELR, length.out = extent["NCOL"] + 1L))
grcs <- with(dis, seq(0, by = DELC, length.out = extent["NROW"] + 1L))

MFimage(bas$IBOUND[,, 1L], gccs, grcs, c(-1, 1),
        c("blue", "grey", "transparent"))

CJBarry/Rflow documentation built on June 16, 2019, 12:35 p.m.