R: Fault-Slip vector plot

nipXY {RFOC} | R Documentation

## Fault-Slip vector plot

### Description

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

### Usage

    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

    set.seed(2015)
    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

     MEKS = list(lon=lon, lat=lat, str1=str1, dip1=dip1,
    rake1=rake1, dep=dep, name=name)

    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)
    }

[Package RFOC version 3.4-3 Index]