R: Get UW focals

getUWfocs {RFOC} | R Documentation

## 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<jonathan.lees@unc.edu>

### References

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

### See Also

getCMT

### Examples

    ## Not run:
    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)

[Package RFOC version 3.4-3 Index]