RFOC (version 3.4-14)

getUWfocs: Get UW focals

Description

Get UW focal mechansims from a file. These are often called A and M cards

Usage

getUWfocs(amfile)

Value

List:

lon

numeric, longitude

lat

numeric, latitude

str1

numeric, strike of plane 1

dip1

numeric, dip of plane 1

rake1

numeric, rake of plane 1

str2

numeric, strike of plane 2

dip2

numeric, dip of plane 2

rake2

numeric, rake of plane 2

sc

character, some GMT info for scale

iexp

character, some GMT info for scale

name

character, name

yr

numeric, year

mo

numeric, month

dom

numeric, day of month

jd

numeric, julian day

hr

numeric, hour

mi

numeric, minute

se

numeric, second

z

numeric, depth

mag

numeric, magnitude

Arguments

amfile

character, file name

Author

Jonathan M. Lees<[e-mail]>

Details

UW focal mechanisms are stored as A and M cards. The A card described the hypocenter the M card describes the focal mechanism.

References

http://www.unc.edu/~leesj/XM_DOC/xm_hypo.doc.html

See Also

getCMT

Examples

#####  uwpickfile is an ascii format file from University of Washington

data(GH, package='RSEIS')

pf = GH$pickfile$PF

### save only the A and M cards to temp file:
fn =  tempfile("temp", fileext = ".UW")
cat(file=fn, c(pf[1], pf[34]) , sep='\n')

G1 = getUWfocs(fn)

plot(G1$lon, G1$lat)

 MEKS = list(lon=G1$lon, lat=G1$lat, str1=G1$str1,
dip1=G1$dip1, rake1=G1$rake1, dep=G1$z, name=G1$name)

##   utm projection
 PROJ = GEOmap::setPROJ(type=2, LAT0=mean(G1$lat) , LON0=mean(G1$lon) )   

     XY = GEOmap::GLOB.XY(G1$lat, G1$lon, PROJ)

     plot(range(XY$x), range(XY$y), type='n', asp=1)

     plotmanyfoc(MEKS, PROJ, focsiz=0.5)