Description Usage Arguments Value Examples
A post-processing tool for DRW
. This function gives a
distributed concentration estimate for a DRW result, using a weighted
kernel density estimate from the particles, using
kde
.
1 2 3 |
DRWmodel |
|
mfdata, wtop |
open NetCDF, character string file path(s) or list of open NetCDFs; MODFLOW data set(s) and saruated water tops |
dkcell |
numeric [1] or [2]; cell spacing for the kernel smooth in x and y directions |
ts, L |
|
smd |
numeric [1]; smoothing distance for kernel smooth: too small and the result will be patchy, too large and the result will be too spread |
H |
matrix [2, 2];
alternative, lower level way of specifying smoothing (see
|
binned |
logical [1];
see |
Kxlim |
|
Kylim |
|
... |
additional arguments for |
ptype |
|
List object of class kDRW with elements:
$conc
numeric [x, y, L, ts]; concentration estimate
$coords
list:
..$x,y
numeric []; x and y grid divide co-ordinates (note not the
cell centres)
..$L,ts
integer []; layers and time steps that are present
$H
numeric [2, 2]; the H matrix used for smoothing (see
kde
)
1 2 3 4 5 6 7 8 9 10 11 12 | setwd(mfdir <- system.file(package = "DRW"))
drmod <- readRDS("DRW_EXAMPLE.rds")
kdr <- ksConc(drmod, "drw_mf_demo.nc", "drw_mf_demo_wtop.nc", 10,
smd = 20, Kxlim = "mf")
# Note here that the result is indexed using strings for the 3rd and 4th
# dimensions. This is much less likely to cause confusion because the
# layers and time steps included are likely not to be a simple sequence
# starting from 1.
with(kdr, image(z = conc[,, "L1", "ts10"], coords$x, coords$y,
col = grDevices::rainbow(51L)))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.