RSEIS (version 4.2-4)

getIRIS: 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)

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

Arguments

fn

character, file path name

skip

numeric, number of lines to skip (e.g. for the header)

Author

Jonathan M. Lees<[e-mail]>

Details

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

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)