nipXY: Fault-Slip vector plot

Description Usage Arguments Details Value Author(s) See Also Examples

Description

Plots a fault plane and the slip vector. Used for geographic representation of numerous focal spheres.

Usage

1
nipXY(MEC, x = x, y = y, focsiz=1, fcol = gray(0.9), nipcol = "black", cex = 0.4)

Arguments

MEC

MEC structure

x

coordinate on plot

y

coordinate on plot

focsiz

size in inches

fcol

color for plotting

nipcol

color of slip point

cex

character expansion for slip point

Details

Slip vector is the cross product of the poles to the fault plane and auxilliary planes.

Value

LIST

Q

output of qpoint

N

slip vector

Author(s)

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

See Also

qpoint, CROSSL, lowplane, TOCART

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
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, xlab='km', ylab='km' )
for(i in 1:length(XY$x))
{
  Msdr = CONVERTSDR(MEKS$str1[i], MEKS$dip1[i],MEKS$rake1[i])
     MEC = MRake(Msdr$M)
       MEC$UP = FALSE

         jcol =  foc.color(foc.icolor(MEC$rake1), pal=1)



nipXY(MEC, x = XY$x[i], y = XY$y[i], focsiz=0.5, fcol = jcol, nipcol = 'black' , cex = 1)
}

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