code: def lat_from_pole(ref_loc_lon, ref_loc_lat, pole_plon, pole_plat):
      
  
      ref_loc = (ref_loc_lon, ref_loc_lat)
      pole = (pole_plon, pole_plat)
      paleo_lat = 90 - pmag.angle(pole, ref_loc)
      return float(paleo_lat)
docs: Calculate paleolatitude for a reference location based on a paleomagnetic pole
  
      Required Parameters
      ----------
      ref_loc_lon: longitude of reference location in degrees
      ref_loc_lat: latitude of reference location
      pole_plon: paleopole longitude in degrees
      pole_plat: paleopole latitude in degrees