makeGrid: Depth grid for 1D models

In dkneis/diatools: Utilities for Diagenetic Modeling

Description

Creates a grid with high resolution near surface and lower resolution at greater depth.

Usage

makeGrid(dz0, dzMax, zMax, beta = 1)

Arguments

dz0	Resolution at surface.
dzMax	Resolution at suggested maximum depth.
zMax	Suggested maximum depth.
beta	Shape parameter. The default of 1 creates a grid whose resolution declines linearly with depth.

Value

A data frame with 1 row per layer and the following columns

• dz Thickness of layer.

• zUp Depth of layer's upper boundary.

• zLw Depth of layer's lower boundary.

Note

The maximum depth of the generated grid (i.e. the lower boundary of the last layer) is guaranteed to be >= zMax. For a grid with high resolution near the surface but low resolution at greater depth one typically wants to set beta to a value > 1.

Author(s)

David Kneis [e-mail]

Examples

# 0.3 m total depth, resolution of 2 mm near surface and 2 cm at depth
gr= makeGrid(dz0=0.002, dzMax=0.02, zMax=0.3, beta=2)
plot(gr$zUp, gr$dz, type="S", xlab="Depth", ylab="Thickness of layer")

dkneis/diatools documentation built on May 15, 2019, 9:12 a.m.