getseis24: Get 24 Hours of Seismic Data
In RSEIS: Seismic Time Series Analysis Tools

getseis24

R Documentation

Get 24 Hours of Seismic Data

Description

Get 24 Hours of Seismic Data

Usage

getseis24(DB, iyear = 2009, iday = 1, usta = "",
 acomp = "", kind = 1,  Iendian=1, BIGLONG=FALSE)

Arguments

`DB`	Data base of meta-data about the seismic trace files
`iyear`	Year for extraction
`iday`	Julian day for extraction
`usta`	station to show
`acomp`	component to show
`kind`	kind of data, default=1, 0="RDATA" , -1="RDS", 0="RDATA", 1 = "segy", 2 = "sac"
`Iendian`	Endian-ness of the data: 1,2,3: "little", "big", "swap". Default = 1 (little)
`BIGLONG`	logical, TRUE=long=8 bytes

Details

The DB file consists of a list of information on where to find the data and what times are covered. DB is

fn: full path to file
yr: year
jd: julian day
hr: hour
mi: minute
sec: second
dur: duration, seconds
origyr: origin time for epoch calculations

Value

`yr`	start year
`jd`	start julian day
`t1`	start t1 (with epoch)
`t2`	start t2 (with epoch day)
`ed`	epoch day
`hr`	start hour
`mi`	start minute
`sec`	start seconds
`gamp`	Amplitude of each trace
`gdt`	delta-t, sample interval, in seconds
`gnam`	station name
`gfile`	file information
`sigs`	List of time series
`zna`	List of NA values in each time series

Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

Examples


## Not run: 
JJ <- getseis24(DB, 2008, 366, usta, acomp, kind = 1)
 
pjj <- plotseis24(JJ, dy=1/18, FIX=24, SCALE=1,
FILT=list(ON<-FALSE, fl=0.05 , fh=20.0, type="BP", proto="BU"),
RCOLS=c(rgb(0.2, .2, 1), rgb(.2, .2, .2))  )


## End(Not run)

RSEIS documentation built on Aug. 19, 2023, 5:07 p.m.