Nothing
`Beachfoc` <-
function(MEC, fcol=gray(0.90), fcolback="white", ALIM = c(-1, -1, +1, +1))
{
if(missing(fcol)) { fcol = gray(0.90) }
if(missing(fcolback)) { fcolback = "white" }
if(missing(ALIM)) {ALIM=c(-1,-1, +1, +1) }
if(is.null(MEC$UP)) MEC$UP=TRUE
if(is.null(MEC$LIM)) MEC$LIM = ALIM
##### net(1)
C = RPMG::circle()
##### print(MEC$LIM)
plot(C$x,C$y, type='n', axes=FALSE, asp=1, xlab="", ylab="", xlim=c(MEC$LIM[1],MEC$LIM[3]), ylim=c(MEC$LIM[2],MEC$LIM[4]) )
lines(C$x,C$y, type='l')
####circ.tics()
pax = focpoint(MEC$P$az, MEC$P$dip, lab="P", UP=MEC$UP, PLOT=FALSE)
PLS = polyfoc(MEC$az1, MEC$dip1, MEC$az2, MEC$dip2, UP=MEC$UP)
### POK = list(x=PLS$Px, y =PLS$Py)
### OLD, does not work: kin = inpoly(pax$x, pax$y,POK)
kin = fields::in.poly(cbind(pax$x, pax$y) ,cbind(PLS$Px, y =PLS$Py) )
if(kin==0)
{
polygon(PLS$Px, PLS$Py, col=fcol )
}
else
{
polygon(C$x,C$y, col=fcol )
polygon(PLS$Px, PLS$Py, col=fcolback )
}
### lines(PLS$Px, PLS$Py, col="black", lwd=2)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.