# draw.brachiopod: Draw Brachiopod In geophys: Geophysics, Continuum Mechanics, Gravity Modeling

## Description

Draw Brachiopod at a particular location

## Usage

 `1` ```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

## Details

Used for strain illustrations

side effects

## Author(s)

Jonathan M. Lees<[email protected]>

 ``` 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 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63``` ```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 rmat= rbind( total, c(x[i],y[i]) ) n = length(BB[,1]) RACK1 = cbind(siz[i]*BB[,1:2], rep(1,n)) draw.brachiopod(RACK1) } ```