docs: def magnetic_lat(inc):
      """
      returns magnetic latitude from inclination
      """
      rad = old_div(np.pi, 180.)
      paleo_lat = old_div(np.arctan(0.5 * np.tan(inc * rad)), rad)
      return paleo_lat