faultplane: fault plane projection on focal sphere

Description

given azimuth and dip of fault mechanism, calculate and plot the fault plane.

Usage

faultplane(az, dip, col = par("col"), PLOT = TRUE, UP = FALSE)

Arguments

degrees, strike of the plane (NOT down dip azimuth)

dip

degrees, dip from horizontal

col

color for line

PLOT

option for adding to plot

upper or lower hemisphere

Value

list of points along fault plane
xcoordinates on focal sphere
ycoordinates on focal sphere

Details

Azimuth is the strike in degrees, not the down dip azimuth as described in other routines.

Examples

Run this code

gcol='black'
border='black'
ndiv=36
phi = seq(0,2*pi, by=2*pi/ndiv);
  x = cos(phi);
  y = sin(phi);

plot(x,y, type='n', asp=1)
  lines(x,y, col=border)
  lines(c(-1,1), c(0,0), col=gcol)
  lines(c(0,0), c(-1,1), col=gcol)

faultplane(65, 34)