getregionals: Extract regional events

In Rquake: Seismic Hypocenter Determination

Description

Extract regional events from a hypocenter list (catalog)

Usage

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<[e-mail]>

See Also

Mine.seis

Examples

## 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 Dec. 16, 2020, 5:06 p.m.