RFOC (version 3.4-10)

plotmanyfoc: 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)

Value

Graphical Side Effects

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

Author

Jonathan M. Lees<[e-mail]>

Details

Input MEK list contains

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

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)
N = 20
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
yr = rep(2017, times=N)
jd = runif(N, min=1, max=365)

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

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)