getIRIS: get Hypocenters from IRIS web site

In RSEIS: Seismic Time Series Analysis Tools

get Hypocenters from IRIS web site

Description

Convert hypocenters from the IRIS website and prepare for plotting in GEOmap

Usage

getIRIS(fn, skip=0)
getANSS(fn, skip=2)

Arguments

fn	character, file path name
skip	numeric, number of lines to skip (e.g. for the header)

Details

Reads in a file dumped out by the website selection box.

Value

list:

yr	vector year
dom	vector, day of month
mo	vector, mo
hr	vector, hour
mi	vector, minute
sec	vector, sec
lat	vector, latitude
lon	vector, longitude
z	vector, depth
mag	vector, magnitude

Note

Be careful about headers and lines that need to be skipped.

for IRIS: http://www.iris.washington.edu/data/event/eventsearch.htm

For ANSS: http://www.quake.geo.berkeley.edu/anss/catalog-search.html

For NEIC (yet to be added) http://earthquake.usgs.gov/earthquakes/eqarchives/epic/epic_global.php

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

getjul

Examples

fn <- tempfile()
K = c(
'Date       Time             Lat       Lon  Depth   Mag Magt  Nst Gap  Clo  RMS  SRC   Event ID',
'----------------------------------------------------------------------------------------------',
'1994/09/06 09:37:36.48  40.1330  144.6240  33.40  4.60   Mb   28          1.22  NEI 199409064025',
'1994/09/06 10:00:02.97  36.4840  140.5730  66.60  4.90   Mb   39          0.88  NEI 199409064028',
'1994/09/06 10:07:16.53  40.1700  144.5890  33.00  4.70   Mb   49          1.09  NEI 199409064029',
'1994/09/06 17:31:52.27  42.6220  142.7000  33.00  5.00   Mb   13          0.54  NEI 199409064042')

cat(file=fn, K, sep='\n')

### check: z = scan(file=fn, what='', sep='\n')

g <- getANSS(fn, skip=2)
g$jd <- getjul(g$yr, g$mo, g$dom)

RSEIS documentation built on Sept. 13, 2024, 1:09 a.m.