new_continuum_ms = continuum_ms.replace(".cal.ms", "_12M.cal.ms")
split(vis=continuum_ms, outputvis=new_continuum_ms, antenna=antennae,
field=field, datacolumn='data')
continuum_ms = new_continuum_ms
else:
antennae = ""
arrayname = '7M12M'
if os.getenv("USE_SELFCAL_MS"):
selfcal_ms = basename+"_"+arrayname+"_selfcal.ms"
continuum_ms = selfcal_ms
logprint("Imaging MS {0} with array {1}".format(continuum_ms, arrayname),
origin='almaimf_cont_imaging')
coosys,racen,deccen = determine_phasecenter(ms=continuum_ms, field=field)
phasecenter = "{0} {1}deg {2}deg".format(coosys, racen, deccen)
(dra,ddec,pixscale) = list(determine_imsize(ms=continuum_ms, field=field,
phasecenter=(racen,deccen),
exclude_7m=exclude_7m,
spw=0, pixfraction_of_fwhm=1/4.))
imsize = [dra, ddec]
cellsize = ['{0:0.2f}arcsec'.format(pixscale)] * 2
contimagename = os.path.join(imaging_root, basename) + "_" + arrayname
# SKIP the plot
#if not os.path.exists(contimagename+".uvwave_vs_amp.png"):
# # make a diagnostic plot to show the UV distribution
# plotms(vis=continuum_ms,
# xaxis='uvwave',
split(
vis=continuum_ms,
outputvis=selfcal_ms,
datacolumn=datacolumn,
antenna=antennae,
spw=spwstr,
width=width,
field=field,
)
logprint("Selfcal MS is: "
"{0}".format(selfcal_ms),
origin='contim_selfcal')
coosys, racen, deccen = determine_phasecenter(ms=selfcal_ms, field=field)
phasecenter = "{0} {1}deg {2}deg".format(coosys, racen, deccen)
(dra, ddec, pixscale) = list(
determine_imsize(ms=selfcal_ms,
field=field,
phasecenter=(racen, deccen),
exclude_7m=exclude_7m,
only_7m=only_7m,
spw='all',
pixfraction_of_fwhm=1 / 4.))
imsize = [dra, ddec]
cellsize = ['{0:0.2f}arcsec'.format(pixscale)] * 2
contimagename = os.path.join(imaging_root, basename) + "_" + arrayname
if not os.path.exists(contimagename + ".uvwave_vs_amp.png"):
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(
mses=vis,
field=field,
phasecenter=(racen, deccen),
spw=0,
pixfraction_of_fwhm=1 / 3.,
exclude_7m=exclude_7m,
min_pixscale=0.1, # arcsec
))
imsize = [int(dra), int(ddec)]
cellsize = ['{0:0.2f}arcsec'.format(pixscale)] * 2
dirty_tclean_made_residual = False
