R: Plot Many Focals

plotmanyfoc {RFOC} | R Documentation

## Plot Many Focals

### Description

Plot a long list of focal mechanisms

### Usage

    plotmanyfoc(MEK, PROJ, focsiz = 0.5, foccol = NULL,
    UP=TRUE, focstyle=1, PMAT = NULL, LEG = FALSE, DOBAR = FALSE)

### Arguments

`MEK` | List of Focal Mechanisms, see details
`PROJ` | Projection
`focsiz` | focal size, inches
`foccol` | focal color
`UP` | logical, UP=TRUE means plot upper hemisphere (DEFAULT=TRUE)
`focstyle` | integer, 1=beach ball, 2=nipplot, 3=strike-slip, 4=P-T, 5=P, 6=T
`PMAT` | Projection Matrix from persp
`LEG` | logical, TRUE= add focal legend for color codes
`DOBAR` | add strike dip bar at epicenter

### Details

Input MEK list contains

MEKS = list(lon=0, lat=0, str1=0, dip1=0, rake1=0, dep=0, name="", Elat=0, Elon=0)

### Value

Graphical Side Effects

### Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

### References

Lees, J. M., Geotouch: Software for Three and Four Dimensional GIS in the Earth Sciences, Computers & Geosciences, 26, 7, 751-761, 2000.

### See Also

justfocXY

### Examples

    set.seed(2015)
    lon=runif(20, 268.1563 , 305)
    lat=runif(20, 7.593004,  25.926045)
    str1=runif(20,50,100)
    dip1=runif(20,10, 80)
    rake1=runif(20,5, 180)

    dep=runif(20,1,15)
    name=seq(from=1, to=length(lon), by=1)
    Elat=NULL
    Elon=NULL

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

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

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

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

    plotmanyfoc(MEKS, PROJ, focsiz=0.5)

[Package RFOC version 3.4-3 Index]