synthStrat: Synthesize stratigraphy from forcing function

In astrochron: A Computational Tool for Astrochronology

Synthesize stratigraphy from forcing function

Description

Synthesize stratigraphy from forcing function.

Usage

synthStrat(signal=NULL,nfacies=4,clip=T,flip=F,fmax=0.1,output=F,genplot=2,verbose=T)

Arguments

signal	Forcing signal. First column should be time (in ka), second column should be forcing.
nfacies	Number of sedimentary facies to model.
clip	Clip forcing signal at mean value? (T or F)
flip	Reverse the sign of the forcing? (T or F)
fmax	Maximum frequency for spectra (if genplot=2).
output	Output facies series? (T or F)
genplot	Generate summary plots? (1) plot stratigraphy, (2) plot statigraphy and spectra.
verbose	Verbose output? (T or F)

Value

modeled stratigraphic series.

Examples

 ## Not run: 
# EX.1: precession, unclipped
signal=etp(tmin=8400,tmax=8900,pWt=1,oWt=0,eWt=0)
synthStrat(signal,nfacies=4,clip=FALSE,genplot=2)

# EX.2: more finely resolved facies
#synthStrat(signal,nfacies=15,clip=FALSE,genplot=2)

# EX.3: couplets
#synthStrat(signal,nfacies=2,clip=FALSE,genplot=2)

# EX.4: precession, clipped
#synthStrat(signal,nfacies=4,genplot=2)

# EX.5: noise
noise=ar1(npts=501,rho=0.8)
#synthStrat(noise,nfacies=4,genplot=2)

# EX.6: precession + noise 
#signal2=signal
#signal2[2]=signal2[2]+0.75*noise[2]
#synthStrat(signal2,nfacies=4,genplot=2)

# EX.7: p-0.5t, clipped (demonstrates interference pattern; compare with EX.4
#signal3=etp(tmin=8400,tmax=8900,pWt=1,oWt=-0.5,eWt=0)
#synthStrat(signal3,nfacies=4,genplot=2)

# EX.8: ice sheet model, using p-0.5t
#ice=imbrie()
#synthStrat(ice,nfacies=5,clip=FALSE,genplot=2)

# EX.9: precession, clipped, ramping sedimentation rate
#synthStrat(linterp(sedRamp(signal,genplot=FALSE),genplot=FALSE),nfacies=6,
# clip=TRUE,genplot=2,fmax=10)

 
## End(Not run)

astrochron documentation built on Sept. 30, 2024, 9:14 a.m.