# PlotTernfoc: Ternary Distribution of focal mechanisms In RFOC: Graphics for Spherical Distributions and Earthquake Focal Mechanisms

## Description

Create and plot a ternary diagram using rake angle to distribute focal mechanisms on a ternary diagram.

## Usage

 `1` ```PlotTernfoc(h, v, x = 0, y = 0, siz = 1, fcols = "black", LABS = FALSE, add = FALSE) ```

## Arguments

 `h` x-coordinate on ternary plot `v` y-coordinate of ternary plot `x` x Location of center of Ternary plot `y` y Location of center of Ternary plot `siz` size of plot in user coordinates `fcols` vector of colors associated with each focal mechanism `LABS` logical, TRUE=add labels at vertices of Ternary plot `add` logical, add to plot=TRUE

## Value

Used for graphical side effect.

## Author(s)

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

## References

J. M. Lees. Geotouch: Software for three and four dimensional gis in the earth sciences. Computers & Geosciences, 26(7):751–761, 2000

## See Also

ternfoc.point, Bfocvec

## Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50``` ```Z1 = c(159.33,51.6,206,18,78, 161.89,54.5,257,27,133, 170.03,53.57,-44,13,171, 154.99,50.16,-83,19,-40, 151.09,47.15,123,23,-170, 176.31,51.41,-81,22,122, 153.71,46.63,205,28,59, 178.39,51.21,-77,16,126, 178.27,51.1,-86,15,115, 177.95,51.14,-83,25,126, 178.25,51.18,215,16,27 ) MZ = matrix(Z1, ncol=5, byrow=TRUE) h = vector() v = vector() Fcol = vector() for(i in 1:length(MZ[,3])) { Msdr = CONVERTSDR(MZ[i,3], MZ[i,4], MZ[i,5]) MEC = MRake(Msdr\$M) MEC\$UP = FALSE az1 = Msdr\$M\$az1 dip1 = Msdr\$M\$d1 az2 = Msdr\$M\$az2 dip2 = Msdr\$M\$d2 BBB = Bfocvec(az1, dip1, az2, dip2) V = ternfoc.point(BBB\$Bdip, Msdr\$M\$pd, Msdr\$M\$td ) h[i] = V\$h v[i] = V\$v Fcol[i] = foc.color(foc.icolor(MEC\$rake1), pal=1) } PlotTernfoc(h,v,x=0, y=0, siz=1, fcols=Fcol, add=FALSE, LAB=TRUE) MFOC1 = SDRfoc(65,90,1, u=FALSE, ALIM=c(-1,-1, +1, +1), PLOT=FALSE) Fcol1 = foc.color(foc.icolor(MFOC1\$rake1), pal=1) MFOC2 = SDRfoc(135,45,-90, u=FALSE, ALIM=c(-1,-1, +1, +1), PLOT=FALSE) Fcol2 = foc.color(foc.icolor(MFOC2\$rake1), pal=1) MFOC3 = SDRfoc(135,45,90, u=FALSE, ALIM=c(-1,-1, +1, +1), PLOT=FALSE) Fcol3 = foc.color(foc.icolor(MFOC3\$rake1), pal=1) justfocXY( MFOC3, fcol = Fcol3, 1.2, -0.9, focsiz = 0.4 ) justfocXY( MFOC2, fcol = Fcol2, -1.2, -0.9, focsiz = 0.4 ) justfocXY( MFOC1, fcol = Fcol1, 0, 1.414443+.2, focsiz = 0.4 ) ```

RFOC documentation built on May 2, 2019, 1:38 p.m.