cosopix: Selection of pickfiles from Coso Geothermal Field

In Rquake: Seismic Hypocenter Determination

Selection of pickfiles from Coso Geothermal Field

Description

Set of selected seismic arrival files with hypocenter locations.

Usage

data("cosopix")

Format

List consisting of:

• PF: original text version of file, as read from disk

• AC: Acard: hypocenter information

• LOC: location

• MC: Fault Mechanizm card

• STAS: Station information

• LIP: Error Ellipse

• E: E-card

• F: F-card

• filename: original file location

• UWFILEID: UW file identification

• comments: Comments on event location

• OSTAS: Station names

• H: High resolution location numbers

• N: Stations Not used in location

Details

Each element of this list is an individual earthquake record.

Examples

data(cosopix)
A = sapply(cosopix, '[[', 'LOC')
###  gather stations

ST.name = vector(mode='character')
ST.lat = vector(mode='numeric')
ST.lon = vector(mode='numeric')
ST.z = vector(mode='numeric')

for(i in 1:length(cosopix))
{
g = cosopix[[i]]
g = data.frame(g$STAS )
w = which(!is.na(g$lat) )
ST.name = c(ST.name, g$name[w])
ST.lat = c(ST.lat, g$lat[w])
ST.lon = c(ST.lon, g$lon[w])
ST.z = c(ST.z, g$z[w])
}

notdup = !duplicated(ST.name)

name = ST.name[notdup ]
lat = ST.lat[notdup ]
lon =ST.lon[notdup ]
z = ST.z[notdup ]

plot(range(c(A[9, ], lon)) , range(c(A[8, ], lat)) , type='n',
xlab='Lon', ylab='Lat')
points(lon, lat, pch=6)

text(lon, lat, labels=name, pos=3)

points(A[9, ], A[8, ])

Rquake documentation built on June 22, 2024, 11:54 a.m.