BKpfile2ypx {RSEIS}R Documentation

Read output of ttimes

Description

Convert a teleseismic pickfile to structure usefule for PICK.GEN. Takes as input the output file of program ttimes, by Brian Kennett, which uses the TauP method.

Usage

BKpfile2ypx(P)

Arguments

P output of TTteles

Details

The output of TTteles includes:

num
phase arrival index
phase
phase name
sec
seconds
tmin
minutes
tsec
fractional seconds after minutes
dtdd
derivative of time w.r.t distance dt/dd
dtdh
derivative of time w.r.t depth dt/dh
dddp
derivative of distance w.r.t ray parameter dd/dp

Value

List of PICKS:

tag tag for unique identification of trace
name station name
comp component name
c3 three element id
phase phase name
err uncertainity of arrival time
pol polarity
flg flag for location program
res residual fom predicted model
yr year
mo month
dom day of month
jd Julian day
hr hour
mi minute
sec second
col color for plotting
onoff logical, flag whether to plot or not
cat category (usually equal to phase)

Note

This routine assumes you have a program called ttimes, or some equivalent. The output of that routine is dumped to a file by ttBK a routine I wrote, and R then reads this in.

Author(s)

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

See Also

TTteles, gwpix2ypx, uwpfile2ypx, PICK.GEN

Examples


## Not run: 

 cmdd = paste(sep=' ', "ttBK 0 ", i , " > ttBK.tempoutput")
system(cmdd )

 A = scan(file="ttBK.tempoutput", list(num=0, phase="", sec=0, tmin=0, tsec=0, dtdd=0, dtdh=0, dddp=0), skip=1)
      attr(A, 'name') = "NPL"
       attr(A, 'input') = list(delta=i , depth=0)
     
## End(Not run)


[Package RSEIS version 2.1-4 Index]