geophys (version 1.4-1)

draw.brachiopod: Draw Brachiopod

Description

Draw Brachiopod at a particular location

Usage

draw.brachiopod(BK = matrix(), x = 0, y = 0, col = "black", fill = NULL, ...)

Arguments

BK

brachiopod matrix

x

x-coordinate location on plot

y

y-coordinate location on plot

col

col for outline

fill

color for fill

…

other par parameters for plotting

Value

side effects

Details

Used for strain illustrations

See Also

get.brachiopod

Examples

# NOT RUN {
gb = get.brachiopod()

draw.brachiopod(gb)



BB = get.brachiopod()
phi0=0
phi =   rep(phi0, 6)
shears = seq(from=0, to=1.5, length=6)
shearmats = list()
for(i in 1:length(shears))
{
 shr = shears[i]
  shearmats[[i]] = matrix(c(1, shr, 0,  1), ncol=2)
}

phi0=0
phi =   rep(phi0, 6)
phi=runif(6, 0, 90)
shears = rep(0, 6)


 for(i in 1:length(phi))
        {
          shr = shears[i]
          shearmats[[i]] = matrix(c(1, shr, 0,  1), ncol=2)
        }
      

x =  rep(c(25,75), times=3)

y = as.vector( t(matrix( rep(c( 75, 50,25 ), times=2) , ncol=2)) )


siz  = rep(30, 6)


plot(c(0, 100), c(0, 100), asp=1, type='n', ann=FALSE, axes=FALSE)
for(i in 1:length(x))
  {
    cosp=cos(phi[i]*pi/180)
    sinp=sin(phi[i]*pi/180)
#######  make a rotation matrix
#### rmat=matrix(c(cosp,-sinp, x, sinp, cosp, y), ncol=2)

    
 ####
    rmat1=matrix(c(cosp,-sinp, sinp, cosp),  ncol=2)
    

    total =  rmat1 
# }
# NOT RUN {
<!-- %*% shearmats[[i]] -->
# }
# NOT RUN {
    
    rmat= rbind( total, c(x[i],y[i]) )
    
    n = length(BB[,1])
    
    RACK1 = cbind(siz[i]*BB[,1:2], rep(1,n)) 
# }
# NOT RUN {
<!-- %*% rmat -->
# }
# NOT RUN {
   draw.brachiopod(RACK1)


  }


# }