R/TOSPHERE.R

In RFOC: Graphics for Spherical Distributions and Earthquake Focal Mechanisms

`TOSPHERE` <-
function(x, y, z)
{
#######  /* Convert from Cartesian to spherical (az, dip in degrees) coordinates */
  RAD2DEG = 180/pi;
  vlength = sqrt( x*x + y*y + z*z);
  if(any(vlength==0.0)) return(list(az=NaN,dip=NaN, x=NaN, y=NaN, z=NaN));
  diprad = acos(z/vlength);
  azrad = atan2(y, x);
  az = azrad * RAD2DEG;
  dip = diprad * RAD2DEG;
  return(list(az=az,dip=dip, x=x, y=y, z=z))
  
}