R/radiateSH.R

In RFOC: Graphics for Spherical Distributions and Earthquake Focal Mechanisms

`radiateSH` <-
function(MEC, SCALE=FALSE, col=col, TIT=FALSE)
{
  if(missing(SCALE)) { SCALE=FALSE }
   if(missing(col)) { col=heat.colors(20) }
  if(missing(TIT)) {  TIT=FALSE }
 
  updown = "Lower"

  if(MEC$UP==TRUE){updown = " Upper"} else {updown = " Lower"}

  imageSH(MEC$az1, MEC$dip1, MEC$rake1, SCALE=SCALE, UP=MEC$UP, col=col )
  net(add=TRUE)
  # PlotPlanes(MEC)
  if(TIT==TRUE)
    {
      
      title(main=paste(" ", "SH-wave",MEC$name, updown),
            sub=paste(sep="", "dip=", MEC$dip1, " strike=",
              MEC$az1, " rake=", format.default(MEC$rake1, digits=3)))
    }
}