synthStrat {astrochron} R Documentation

## 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)


[Package astrochron version 1.1 Index]