MFncimage: Plot MODFLOW variable matrix from NetCDF

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

Description

Plot MODFLOW variable matrix from NetCDF

Usage

MFncimage(nc, variable, start = c(NA, NA, 1L, 1L), count = c(NA, NA, 1L,
  1L), zlim = "auto", col = colorRampPalette(c("blue", "white",
  "red"))(101L), show.range = identical(zlim, "angle"), ...)

Arguments

nc	NetCDF object; an open connection to a MODFLOW/ MT3DMS data set
variable	character string; name of the variable to be plotted
start	integer [4]; index in each dimension at which to start reading (see var.get.nc, pkg = RNetCDF); NA implies start
count	integer [4]; numbers of indices to read in each dimension (see var.get.nc, pkg = RNetCDF); NA implies to end
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)
...
...	Arguments passed on to MFimage 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)

Examples

library("RNetCDF")
fnm <- system.file("rflow_mf_demo.nc", package = "Rflow")
mfdata <- open.nc(fnm)

# plot the flux from river (blue is negative)
MFncimage(mfdata, "RiverLeakage",
          c(NA, NA, 1L, 10L),
          c(NA, NA, 1L, 1L),
          "sym", show.range = TRUE)

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