Description Usage Arguments Details Value Author(s) See Also Examples
This routine calculates the value of a given property at the middle of
the grid cells (mid
) and at the interfaces of the grid cells
(int
).
Two possibilities are available: either specifying a mathematical function
(func
) that describes the spatial dependency of the property,
or obtaining the property from interpolation of a data series
(via the input of the data matrix xy
).
For example, in a sediment model, setup.prop.1D
can be used to
specify the porosity, the mixing intensity or some other parameter over
the onedimensional grid. Similarly, in a vertical water column model, setup.prop.1D
can be
used to specify the sinking velocity of particles or other model parameters
changing with water depth.
1 2 3 4 5 
func 
function that describes the spatial dependency. For example, one can use the functions provided in 
value 
constant value given to the property (no spatial dependency) 
xy 
a twocolumn data matrix where the first column ( 
interpolate 
specifies how the interpolation should be done, one
of "spline" or "linear"; only used when 
grid 
list specifying the 1D grid characteristics, see

x 
the object of class 
xyswap 
if 
... 
additional arguments that are passed on to 
There are two options to carry out the data interpolation:
"spline" gives a smooth profile, but sometimes generates strange profiles  always check the result!
"linear" gives a segmented profile
A list of type prop.1D
containing:
mid 
property value in the middle of the grid cells; vector of length N (where N is the number of grid cells) 
int 
property value at the interface of the grid cells; vector of length N+1 
Karline Soetaert <karline.soetaert@nioz.nl>, Filip Meysman <filip.meysman@nioz.nl>
tran.1D
, for a discretisation of the general transport equation in 1D
setup.grid.1D
, the creation of grids in 1D
setup.prop.2D
for defining properties on 2D grids.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22  # Construction of the 1D grid
grid < setup.grid.1D(x.up = 0, L = 10, N = 10)
# Porosity profile via function specification
P.prof < setup.prop.1D(func = p.exp, grid = grid, y.0 = 0.9,
y.inf = 0.5, x.att = 3)
# Porosity profile via data series interpolation
P.data < matrix(ncol = 2, data = c(0,3,6,10,0.9,0.65,0.55,0.5))
P.spline < setup.prop.1D(xy = P.data, grid = grid)
P.linear < setup.prop.1D(xy = P.data, grid = grid, interpolate = "linear")
# Plot different profiles
plot(P.prof, grid = grid, type = "l",
main = "setup.prop, function evaluation")
points(P.data, cex = 1.5, pch = 16)
lines(grid$x.int, P.spline$int, lty = "dashed")
lines(grid$x.int, P.linear$int, lty = "dotdash")

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.