code: def plat(inc):
      
      tani = np.tan(np.radians(inc))
      lat = np.arctan(tani/2.)
      return np.degrees(lat)
docs: calculate paleolatitude from inclination using dipole formula: tan(I) = 2tan(lat)
      Parameters
      ________________
  
      inc : either a single value or an array of inclinations
  
      Returns
      -------
  
      array of latitudes