getregionals: Extract regional events

Description Usage Arguments Details Value Author(s) See Also Examples

View source: R/getregionals.R

Description

Extract regional events from a hypocenter list (catalog)

Usage

1
getregionals(KAT, Mlat, Mlon, rad = 1000, t1 = 1, t2 = 2)

Arguments

KAT

catalog list

Mlat

central latitude

Mlon

central longitude

rad

radius (km)

t1

start time (julian days)

t2

end time (julian days)

Details

given an earthquake catalog from PDEs, for example, extract the events that are close to a network in a given time frame. The limited data set may be used to help predict arrival times for known hypocenter locations.

Value

Catalog

Author(s)

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

See Also

Mine.seis

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
## Not run: 
fpde = '/Users/lees/Site/CHAC/pde_catalog_NOV_2011.txt'

ipde = prepPDE(fpde)
Mlat = median(stas$lat)
Mlon = median(stas$lon)

Lt1 =  JtimL(upf$LOC)
t1 = Lt1-3/(24*6)
t2 = Lt1+3/(24*6)

localeqs = getregionals(KAT, Mlat, Mlon, rad=100 ,  t1=NULL, t2=NULL)

for(i in 1:length(localeqs))
  {
j = localeqs[i]
 at1 = KAT$jsec[j]
        at2 = at1+20/(24*60)
 GH = Mine.seis(at1, at2, DBnov , gstas, gcomps, kind=2, Iendian=1, BIGLONG=FALSE)

hord = which(GH$COMPS=="V")
        
        gret = swig(GH, sel=hord, PADDLAB=buts)
####  save any picks to a file on disc for later use


  }



## End(Not run)

Rquake documentation built on May 29, 2017, 3:02 p.m.