code: def azimuth2rake(strike, dip, azimuth):
      
      plunge, bearing = plane_intersection(strike, dip, azimuth, 90)
      rake = project_onto_plane(strike, dip, plunge, bearing)
      return rake
docs: Projects an azimuth of a linear feature onto a plane as a rake angle.
  
      Parameters
      ----------
      strike, dip : numbers
          The strike and dip of the plane in degrees following the
          right-hand-rule.
      azimuth : numbers
          The azimuth of the linear feature in degrees clockwise from north (i.e.
          a 0-360 azimuth).
  
      Returns
      -------
      rake : number
          A rake angle in degrees measured downwards from horizontal.  Negative
          values correspond to the opposite end of the strike.