Fault-Slip vector plot
Plots a fault plane and the slip vector. Used for geographic representation of numerous focal spheres.
nipXY(MEC, x = x, y = y, focsiz=1, fcol = gray(0.9), nipcol = "black", cex = 0.4)
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 |
Slip vector is the cross product of the poles to the fault plane and auxilliary planes.
LIST
Q |
output of qpoint |
N |
slip vector |
Jonathan M. Lees<jonathan.lees@unc.edu>
qpoint, CROSSL, lowplane, TOCART
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)
}Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.