BKpfile2ypx {RSEIS} | R Documentation |
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.
BKpfile2ypx(P)
P |
output of TTteles |
The output of TTteles includes:
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) |
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.
Jonathan M. Lees<jonathan.lees@unc.edu>
TTteles, gwpix2ypx, uwpfile2ypx, PICK.GEN
## 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)