@@ -124,9 +124,10 @@ identifying some spectral lines in the data.
124124
125125 m83table = Alma.query_object(' M83' , public = True )
126126 m83urls = Alma.stage_data(m83table[' Asdm_uid' ])
127- m83files = Alma.download_and_extract_files(m83urls[' URL' ])
128- # Sometimes there can be duplicates
129- m83files = list (set (m83files))
127+ # Sometimes there can be duplicates: avoid them with
128+ # list(set())
129+ m83files = Alma.download_and_extract_files(list (set (m83urls[' URL' ])))
130+ m83files = m83files
130131
131132 Simbad.add_votable_fields(' rvel' )
132133 m83simbad = Simbad.query_object(' M83' )
@@ -148,12 +149,13 @@ identifying some spectral lines in the data.
148149
149150 # Change the cube coordinate system to be in velocity with respect
150151 # to the rest frequency (in the M83 rest frame)
151- rests_frequency = lines[' Freq-GHz' ][0 ]* u.GHz / (1 + rvel/ constants.c)
152+ rest_frequency = lines[' Freq-GHz' ][0 ]* u.GHz / (1 + rvel/ constants.c)
152153 vcube = cube.with_spectral_unit(u.km/ u.s,
153154 rest_value = rest_frequency,
154155 velocity_convention = ' radio' )
155156
156157 # Write the cube with the specified line name
157158 fmt = " {Species} {Resolved QNs}"
159+ row = lines[0 ]
158160 linename = fmt.format(** dict (zip (row.colnames,row.data)))
159161 vcube.write(' M83_ALMA_{linename} .fits' .format(linename = linename))
0 commit comments