faultdip {GEOmap}

Show Fault dip

Description

Show Fault dip

Usage

faultdip(x, y, rot = 0, h = 1, lab = "")

Arguments

x

x-coordinates

y

y-coordinates

rot

cosine and sine of rotation

h

length of mark

lab

labels

Value

Graphical Side effect

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

perpen, PointsAlong, getsplineG

Examples

ff=list()
ff$x=c(-4.850,-4.700,-3.934,-2.528, 0.603, 2.647, 3.861, 2.626)
ff$y=c(-4.045,-2.087,-0.710, 0.172, 1.291, 2.087,-0.753,-4.131)

G =getsplineG(ff$x, ff$y, kdiv=20)
g = PointsAlong(G$x, G$y, N=5)


plot(c(-5,5), c(-5,5), asp=1, type='n' )
lines(G)

angs = 180*atan(g$rot$sn/g$rot$cs)/pi
faultdip(g$x , g$y , rot=angs, h=.5, lab='')