from astroquery.alma import Alma

import astroquery
print(astroquery.__version__)

MessierNames = ["M{0:d}".format(ii) for ii in range(1,111)]

MessierALMAData = {name: Alma.query_object(name, cache=True) for name in MessierNames}

MessierALMABands = {name: set(MessierALMAData[name]['Band']) for name in MessierALMAData}

observed_messier_objects = {name: MessierALMABands[name] for name in MessierALMABands if MessierALMABands[name]}
observed_messier_objects

len(observed_messier_objects)

from astropy import coordinates
from astropy import units as u

coord_list = coordinates.SkyCoord([(1.0, 2.0), (-50.0, 1.0), (0.6667,-0.03)]*u.deg, frame='galactic')
coord_list

results = [Alma.query_region(coord, radius=1*u.arcmin) for coord in coord_list]

[len(x) for x in results]

results[-1]

set(results[-1]['Band'])