code: def bc02(data):
      
  
      plate, site_lat, site_lon, age = data[0], data[1], data[2], data[3]
      apwp = get_plate_data(plate)
      recs = apwp.split()
      #
      # put it into  usable form in plate_data
      #
      k, plate_data = 0, []
      while k < len(recs) - 3:
          rec = [float(recs[k]), float(recs[k + 1]), float(recs[k + 2])]
          plate_data.append(rec)
          k = k + 3
  
      #
      # find the right pole for the age
      #
      for i in range(len(plate_data)):
          if age >= plate_data[i][0] and age <= plate_data[i + 1][0]:
              if (age - plate_data[i][0]) < (plate_data[i][0] - age):
                  rec = i
              else:
                  rec = i + 1
              break
      pole_lat = plate_data[rec][1]
      pole_lon = plate_data[rec][2]
      return pole_lat, pole_lon
docs: get APWP from Besse and Courtillot 2002 paper
  
      Parameters
      ----------
      Takes input as [plate, site_lat, site_lon, age]
      plate : string (options: AF, ANT, AU, EU, GL, IN, NA, SA)
      site_lat : float
      site_lon : float
      age : float in Myr
  
      Returns
      ----------