from astroquery.splatalogue import Splatalogue,utils
from astropy import units as u

# find all CO lines between 100 and 1000 GHz
# Note the spaces around CO: chemical_name is a regular 
# expression, so this excludes names like 13CO and C18O
table = Splatalogue.query_lines(100*u.GHz, 1000*u.GHz, 
                                chemical_name=" CO ")
table

print(Splatalogue.query_lines.__doc__)

# find all CO lines *you might expect to see* between 100 and 1000 GHz
# by selecting lines with an upper state energy <= 300 K
table2 = Splatalogue.query_lines(100*u.GHz, 1000*u.GHz,
                                 chemical_name=" CO ", energy_max=300,
                                 energy_type='eu_k')
table2

print(len(table2))

table2.pprint(max_lines=100,max_width=70)

dir(utils)

minitable2 = utils.minimize_table(table2)
minitable2

table3 = Splatalogue.query_lines(0.1*u.GHz, 1000*u.GHz,
                                 chemical_name='Ethylene Glycol')
utils.minimize_table(table3)

table2.write('test.csv',format='ascii.csv')

cat test.csv

pwd

table4 = Splatalogue.query_lines(218.40*u.GHz, 218.50*u.GHz, energy_max=300, energy_type='eu_k', line_lists='JPL')
table4

len(table4)

np.unique(table4['Chemical Name'])

table4[True - (table4['Chemical Name'] == 'Vinyl Cyanide')]

import pylab as pl
pl.plot(table4['Freq-GHz'])

table5 = Splatalogue.query_lines(67*u.GHz,95*u.GHz,chemical_name="D")
len(table5)

mask = ['D' in x for x in table5['Species']]
sum(mask)

table5[np.array(mask)]

from astroquery.simbad import Simbad

Simbad.query_object('L1689N')