astrochron (version 1.4)

rankSeries: Create lithofacies rank series from bed thickness data

Description

Create lithofacies rank series from bed thickness data.

Usage

rankSeries(dat,dt,start=0,genplot=T,verbose=T)

Arguments

dat

First column should be bed thickness, and second column should bed lithofacies rank.

dt

Sampling interval for piecewise linear interpolation. By default a grid spacing that is 5 times smaller than the thinnest bed is used. If dt is set to zero, interpolation is skipped.

start

Start at what time/depth/height value?

genplot

Generate summary plots? (T or F)

verbose

Verbose output? (T or F)

Examples

# generate example series with random bed thicknesses
exThick=rnorm(n=20,mean=10,sd=2)
# assign alternating rank of 1 and 2
rank=double(20)
rank[seq(from=1,to=19,by=2)] <- 1
rank[seq(from=2,to=20,by=2)] <- 2

# combine into a dataframe
ex=cb(exThick,rank)

# generate lithofacies rank series
rankSeries(ex)