RFOC (version 2.0-00)

getUWfocs: 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

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

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

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)


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

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

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

     plotmanyfoc(MEKS, PROJ, focsiz=0.05)