Mine.Teles: Mine seismic data base for teleseismic/regional events

Description

Mine seismic data base for teleseismic/regional events

Usage

Mine.Teles(DB, ev, sta, ista = 1, usta = "DOM", acomp = "V", velmodel = NULL, IASPpath = "", regdist = 110)

Arguments

Data base of segy/sac files (makeDB)

list of events (lat, lon, date/time, mag)

sta

list of stations (name, lat, lon z)

ista

which stations to predict arrivals for

usta

which stations to show in plot

acomp

which components to show in plot

velmodel

regional/local 1D velocity model

IASPpath

path to IASP91 programs and binary vel models for Taup-P

regdist

radius for regional vs global model calculations

Value

List of predicted arrival times per each event, plotting of events and picks on those events

Details

program must run external codes from Brian Kennets software.

Examples

Run this code

load("/home/Site/lees/SANTIAGUITO/SGRobjs/GUAT09.sta.RData")

sta$name = sta$nam

stinfo  = targetLL(sta, rdist=100)

proj = stinfo$proj



IASPpath = '/home/lees/Progs/TTIMES'
eq1 = scan(file="sg06_eqs_neic.csv", skip=1, sep=",", what=list(yr=0, mo=0, dom=0, t1="", lat=0, lon=0, mag=0, depth=0))

######  parse out the time:
eq1$hr = as.numeric( substr(eq1$t1, 1,2))
eq1$mi = as.numeric( substr(eq1$t1, 3,4))
eq1$sec = as.numeric( substr(eq1$t1, 5,9))

eq1$z = eq1$depth
eq1$jd = getjul(eq1$yr, eq1$mo,  eq1$dom)

w1 = which(eqdists$dist<110)

##########   restrict the global data to only very large events?
w2 = which(eq1$mag>5.5)


##    create a an earthquake structure list

W = c(w1, w2)
##  W = c(w1 )


mylist = list()
for(i in 1:length(W))
{
j = W[i] 
mylist[[i]] = list(yr=eq1$yr[j], jd=eq1$jd[j], mo=eq1$mo[j], dom=eq1$dom[j], hr=eq1$hr[j], 
mi=eq1$mi[j], sec=eq1$sec[j], lat=eq1$lat[j], lon=eq1$lon[j], z=eq1$z[j], mag=eq1$mag[j])


}
NEW.EV = Mine.Teles(DB, mylist, sta, ista=1, usta=sta$name, acomp="V", IASPpath=IASPpath, regdist=110 )