for row in splat:
        closest = np.argmin(np.abs(barton['Freq']-row['Freq']))
        closest_inds.append(closest)

    barton_ = barton[np.array(closest_inds)]
    result = table.hstack([splat, barton_])
    result.add_column(table.Column(data=result['Freq_1'] - result['Freq_2'], name='Splat-Barton'))
    result.add_column(table.Column(data=(result['Freq_1'] - result['Freq_2'])/result['Freq_2']*3e5, name='Splat-Barton_kms'))
    return result

KCls = mt(Splatalogue.query_lines(8*u.GHz, 40000*u.GHz, chemical_name=' KCl', line_lists=['SLAIM']))
KCl = load_barton('39K-35Cl')

KCl_diff = match_splat_barton(KCls, KCl[KCl['vu']<4])
KCl_offset_model = get_offset_model(KCl_diff)


NaCls = mt(Splatalogue.query_lines(8*u.GHz, 40000*u.GHz, chemical_name=' NaCl', line_lists=['SLAIM']))
NaCl = load_barton('23Na-35Cl')

NaCl_diff = match_splat_barton(NaCls, NaCl[NaCl['vu']<4])
NaCl_offset_model = get_offset_model(NaCl_diff)

cabezasNaCl = Vizier(row_limit=1e9).get_catalogs('J/ApJ/825/150')[1]
cabezasNaCl['nuCalc'][cabezasNaCl['x_nuCalc'] == 'cm-1'] = u.Quantity(cabezasNaCl['nuCalc'][cabezasNaCl['x_nuCalc'] == 'cm-1'], u.cm**-1).to(u.MHz, u.spectral())
cabezasNaCl['x_nuCalc'][cabezasNaCl['x_nuCalc'] == 'cm-1'] = 'MHz'
cabezasNaCl.add_column(Column(name='Freq', data=u.Quantity(cabezasNaCl['nuCalc'], u.MHz).to(u.GHz)))

cabezas_match = match_splat_barton(cabezasNaCl, NaCl)
splat_cabezas_match = match_splat_barton(NaCls, cabezasNaCl)
NaClbarton = load_barton('23Na-35Cl')

NaCl_diff = match_splat_barton(NaCls, NaClbarton[NaClbarton['vu'] < 4])
NaCl_offset_model = get_offset_model(NaCl_diff)

NaClbarton['Freq'] = NaClbarton['Freq'] + NaCl_offset_model(
    NaClbarton['vu'], NaClbarton['Ju'])

# use Cabezas instead of Barton
cabezasNaCl = Vizier(row_limit=1e9).get_catalogs('J/ApJ/825/150')[1]
cabezasNaCl['nuCalc'][cabezasNaCl['x_nuCalc'] == 'cm-1'] = u.Quantity(
    cabezasNaCl['nuCalc'][cabezasNaCl['x_nuCalc'] == 'cm-1'],
    u.cm**-1).to(u.MHz, u.spectral())
cabezasNaCl['x_nuCalc'][cabezasNaCl['x_nuCalc'] == 'cm-1'] = 'MHz'
cabezasNaCl.add_column(
    Column(name='Freq',
           data=u.Quantity(cabezasNaCl['nuCalc'], u.MHz).to(u.GHz)))
cabezasNaCl.rename_column('J1', 'Ju')
cabezasNaCl.rename_column('J0', 'Jl')
cabezasNaCl.add_column(Column(name='gu', data=cabezasNaCl['Ju'] * 32 + 16))
cabezasNaCl.add_column(Column(name='gl', data=cabezasNaCl['Jl'] * 32 + 16))
cabezasNaCl.rename_column('V1', 'vu')
cabezasNaCl.rename_column('V0', 'vl')
cabezasNaCl.rename_column('Eup', 'E_U')
cabezasNaCl.add_column(
    Column(name='E_L',
           data=cabezasNaCl['E_U'] - cabezasNaCl['Freq'].quantity.to(
               u.eV, u.spectral()).to(u.K, u.temperature_energy())))
cabezasNaCl.add_column(
    Column(name='QNs',
           data=[