code: def interpolate_delta_t(delta_t_table, tt):
      
      tt_array, delta_t_array = delta_t_table
      delta_t = _to_array(interp(tt, tt_array, delta_t_array, nan, nan))
      missing = isnan(delta_t)
  
      if missing.any():
          # Test if we are dealing with an array and proceed appropriately
          if missing.shape:
              tt = tt[missing]
              delta_t[missing] = delta_t_formula_morrison_and_stephenson_2004(tt)
          else:
              delta_t = delta_t_formula_morrison_and_stephenson_2004(tt)
      return delta_t
docs: Return interpolated Delta T values for the times in `tt`.
  
      The 2xN table should provide TT values as element 0 and
      corresponding Delta T values for element 1.  For times outside the
      range of the table, a long-term formula is used instead.