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