justfocXY: Plot focal mechanism

Description

Add simple focal mechanisms to plot

Usage

justfocXY(MEC, x = x, y = y, focsiz=1 , fcol = gray(0.9),
          fcolback = "white", xpd = TRUE)

Value

Used for graphical side effect

Arguments

MEC: MEC structure
x: x-coordinate of center
y: y-coordinate of center
focsiz: size of focal sphere in inches
fcol: color of shaded region
fcolback: color of background region
xpd: logical, whether to extend the plot beyond, or to clip

Author

Jonathan M. Lees <jonathan.lees@unc.edu>

Details

This routine can be used to add focal mechanisms on geographic map or other plot.

Examples

Run this code


#### read in some data:


Z1 = c(159.33,51.6,206,18,78,
161.89,54.5,257,27,133,
170.03,53.57,-44,13,171,
154.99,50.16,-83,19,-40,
151.09,47.15,123,23,-170,
176.31,51.41,-81,22,122,
153.71,46.63,205,28,59,
178.39,51.21,-77,16,126,
178.27,51.1,-86,15,115,
177.95,51.14,-83,25,126,
178.25,51.18,215,16,27
)

MZ = matrix(Z1, ncol=5, byrow=TRUE)

plot(MZ[,1], MZ[,2], type='n', xlab="LON", ylab="LAT", asp=1)

for(i in 1:length(MZ[,1]))
{
paste(MZ[i,3], MZ[i,4], MZ[i,5])


MEC =  SDRfoc(MZ[i,3], MZ[i,4], MZ[i,5], u=FALSE, ALIM=c(-1,-1, +1, +1), PLOT=FALSE)
fcol =  foc.color(foc.icolor(MEC$rake1), pal=1)
justfocXY(MEC, x=MZ[i,1], y =MZ[i,2] , focsiz=.5, fcol =fcol , fcolback = "white", xpd = TRUE)


}

Run the code above in your browser using DataLab

Last chance! 50% off unlimited learning