code: def geographic2pole(lon, lat):
      
      plunge, bearing = geographic2plunge_bearing(lon, lat)
      strike = bearing + 90
      strike[strike >= 360] -= 360
      dip = 90 - plunge
      return strike, dip
docs: Converts a longitude and latitude (from a stereonet) into the strike and dip
      of the plane whose pole lies at the given longitude(s) and latitude(s).
  
      Parameters
      ----------
      lon : array-like
          A sequence of longitudes (or a single longitude) in radians
      lat : array-like
          A sequence of latitudes (or a single latitude) in radians
  
      Returns
      -------
      strike : array
          A sequence of strikes in degrees
      dip : array
          A sequence of dips in degrees