code: def pinc(lat):
      
      tanl = np.tan(np.radians(lat))
      inc = np.arctan(2. * tanl)
      return np.degrees(inc)
docs: calculate paleoinclination from latitude using dipole formula: tan(I) = 2tan(lat)
      Parameters
      ________________
  
      lat : either a single value or an array of latitudes
  
      Returns
      -------
  
      array of inclinations