DDFacet command:
DDF.py --Data-MS T16R02C02-corr-scan3-chan512-t16-XXYY.ms,T16R02C04-corr-scan4-chan512-t16-XXYY.ms,T16R02C06-corr-scan5-chan512-t16-XXYY.ms --Data-ColName DATA --Data-ChunkHours 3.0 --Data-Sort 1 --Predict-ColName POINTING_MODEL --Output-Mode Clean --Output-Name mgpls_im_01/T16R02-avg-mosaic-di1/im-di1 --Output-Cubes all --Output-Images all --Image-NPix 18000 --Image-Cell 1.0 --Facets-NFacets 9 --Weight-Robust -1.0 --Weight-SigmoidTaperInnerRolloffStrength 0 --Weight-SigmoidTaperOuterRolloffStrength 0 --RIME-DecorrMode FT --Parallel-NCPU 60 --Cache-Reset 1 --Beam-Model FITS --Beam-NBand 3 --Beam-CenterNorm 1 --Beam-Smooth 1 --Beam-FITSFile mktat_beams/new_mkt_pb_avg_$(xy)_$(reim).fits --Beam-FITSFeed xy --Beam-DtBeamMin 4.0 --Beam-FITSParAngleIncDeg 0.5 --Freq-NBand 3 --Freq-NDegridBand 12 --Deconv-Mode SSD2 --Deconv-MaxMajorIter 4 --Deconv-RMSFactor 3.0 --Deconv-PeakFactor 0.01 --Mask-Auto 1 --SSD2-PolyFreqOrder 3 --Log-Memory 1 --Log-Boring 0
Name= im-dil.app.restored.fits
HDU = 0
Shape = [18225, 18225, 1, 1]
Number of channels = 1
Number of stokes = 1
Coordinate type = Right Ascension, Declination
Projection = SIN
Image reference pixels = [9113, 9113]
Image reference coords [14:17:01.3800, -061.10.37.9000]
Image ref coords (deg) = [214.256 deg, -61.1772 deg]
Pixel increment = -1′′, 1′′
Celestial frame = ICRS
Spectral frame = TOPO
Velocity definition = RADIO
Pixel unit = Jy/beam
Restoring beam = 21.8848′′ X 10.2914", 50.8834 deg
Dirty Images for each pointing