R/Comp1Dvel.R

In RSEIS: Seismic Time Series Analysis Tools

`Comp1Dvel` <-
function(v, v2, col=c('blue', 'brown'),   ... )
{
    
      plot(c(v$vp,v$vs, v2$vp, v2$vs), c( -v$zp,-v$zs, -v2$zp,-v2$zs), type='n',
           xlab="Velocity, km/s", ylab="Depth, km")
      
     u = par('usr')
    lines(v$vp, -v$zp, type='s', col=col[1], ...)
    lines(v$vs, -v$zs, type='s', col=col[2], ...)

    lines(v2$vp, -v2$zp, type='s', col=col[1], ...)
    lines(v2$vs, -v2$zs, type='s', col=col[2], ...)

    grid()
    title(v$name)
  }