origin='almaimf_line_imaging')
continue
else:
logprint("Matched spectral window {0} to line {1}".format(
spw, line_name),
origin='almaimf_line_imaging')
elif line_name in spwnames and line_name.lstrip("spw") != spw:
logprint(
"Skipped spectral window {0} because it's not {1}".format(
spw, line_name),
origin='almaimf_line_imaging')
continue
if exclude_7m:
# don't use variable name 'ms' in python2.7!
vis = [ms_ for ms_ in vis if not is_7m(ms_)]
arrayname = '12M'
else:
arrayname = '7M12M'
lineimagename = os.path.join(
imaging_root,
"{0}_{1}_spw{2}_{3}_{4}".format(field, band, spw, arrayname,
line_name))
logprint("Measurement sets are: " + str(vis),
origin='almaimf_line_imaging')
coosys, racen, deccen = determine_phasecenter(ms=vis, field=field)
phasecenter = "{0} {1}deg {2}deg".format(coosys, racen, deccen)
(dra, ddec, pixscale) = list(
determine_imsizes(
# CASAguides recommend chanchunks=-1, but this resulted in:
# 2018-09-05 23:16:34 SEVERE tclean::task_tclean:: Exception from task_tclean : Invalid Gridding/FTM Parameter set : Must have at least 1 chanchunk
chanchunks = os.getenv('CHANCHUNKS') or 16
# global default: only do robust 0 for lines
robust = 0
for band in band_list:
for field in to_image[band]:
for spw in to_image[band][field]:
vis = list(map(str, to_image[band][field][spw]))
if exclude_7m:
vis = [ms for ms in vis if not(is_7m(ms))]
arrayname = '12M'
else:
arrayname = '7M12M'
lineimagename = os.path.join(imaging_root,
"{0}_{1}_spw{2}_{3}_lines".format(field,
band,
spw,
arrayname))
logprint(str(vis), origin='almaimf_line_imaging')
coosys,racen,deccen = determine_phasecenter(ms=vis, field=field)
phasecenter = "{0} {1}deg {2}deg".format(coosys, racen, deccen)
(dra,ddec,pixscale) = list(determine_imsizes(mses=vis, field=field,