getGAP: Get Seismic Gap

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

View source: R/getGAP.R

Description

Given an earthquake and a set of stations, return the maximum angle subtended between adjacent stations relative to the epicenter.

Usage

1
getGAP(EQ, Ldat, PLOT = FALSE)

Arguments

EQ

List, Earthequake location, elements (lat, lon) must be present

Ldat

List, station information, (lat, lon) must be present

PLOT

logical, plot the stations and show the gap

Details

Theangles are calculated in cartesian coordinates with the epicenter at the origin using a UTM projection.

Value

numeric, gap in degrees

Author(s)

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

See Also

eqwrapup

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
set.seed(0)

N = 10
snames = paste(sep="", "A", as.character(1:N))
stas = list(name=snames, lat=runif(N, 35.9823, 36.1414), lon=runif(N, -118.0031, -117.6213))

NEQ = 3
WEQ = list(lat=runif(NEQ, 35.9823, 36.1414), lon=runif(NEQ, -118.0031, -117.6213))

 MLAT = median(stas$lat)
  MLON = median(stas$lon)
  proj = GEOmap::setPROJ(type=2, LAT0=MLAT, LON0=MLON)

  XYSTAS = GEOmap::GLOB.XY(stas$lat,  stas$lon , proj)
  eqxy = GEOmap::GLOB.XY(WEQ$lat, WEQ$lon, proj)


plot(range(c(XYSTAS$x, eqxy$x)), range(c(XYSTAS$y, eqxy$y)), type='n', asp=1, xlab="km", ylab="km" )
points(XYSTAS$x, XYSTAS$y, pch=6)

for(i in 1:NEQ)
{
EQ = list(lat=WEQ$lat[i], lon=WEQ$lon[i])


g = getGAP(EQ, stas, PLOT=FALSE)

points(eqxy$x[i], eqxy$y[i], pch=8, col='red')
text(eqxy$x[i], eqxy$y[i], labels=paste("gap=", format(g)), pos=3)

}

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

Search within the Rquake package
Search all R packages, documentation and source code