getUWfocs {RFOC}

Get UW focals

Description

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

Usage

getUWfocs(amfile)

Arguments

amfile

character, file name

Details

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

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

Note

Uses UW2 format, so full 4 digit year is required

Author(s)

Jonathan M. Lees<[e-mail]>

References

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

See Also

getCMT

Examples

## Not run: 
#####  uwpickfile is an ascii format file from University of Washington
G1 = getUWfocs(uwpickfile)

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.05)




## End(Not run)