2023-11-04 07:15:18 INFO ::casa Using configuration file ~/.casa/config.py 2023-11-04 07:15:18 INFO ::casa 2023-11-04 07:15:18 INFO ::casa Checking Measures tables in data repository sub-directory /wheeldirectory/casa-CAS-14200-7-py3.8.el8/lib/py/lib/python3.8/site-packages/casadata/__data__/geodetic 2023-11-04 07:15:18 INFO ::casa IERSeop2000 (version date, last date in table (UTC)): 2023/10/06/15:00, 2023/09/06/00:00:00 2023-11-04 07:15:18 INFO ::casa IERSeop97 (version date, last date in table (UTC)): 2023/10/06/15:00, 2023/09/06/00:00:00 2023-11-04 07:15:18 INFO ::casa IERSpredict (version date, last date in table (UTC)): 2023/10/29/15:00, 2024/01/27/00:00:00 2023-11-04 07:15:18 INFO ::casa TAI_UTC (version date, last date in table (UTC)): 2023/10/13/15:00, 2017/01/01/00:00:00 2023-11-04 07:15:19 INFO fixvis::::casa ########################################## 2023-11-04 07:15:19 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:19 INFO fixvis::::casa fixvis( vis='0420+417.ms', outputvis='output.ms', field=[], refcode='J2000', reuse=True, phasecenter='', distances=[], datacolumn='all' ) 2023-11-04 07:15:19 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:19 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:19 INFO fixvis::::casa Will leave phase centers unchanged. 2023-11-04 07:15:19 INFO fixvis::::casa Recalculating the UVW coordinates ... 2023-11-04 07:15:19 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/output.ms 2023-11-04 07:15:19 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-04 07:15:19 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:19 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:18.593717 End time: 2023-11-04 07:15:18.820978 2023-11-04 07:15:19 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:19 INFO fixvis::::casa ########################################## .2023-11-04 07:15:19 INFO split::::casa ########################################## 2023-11-04 07:15:19 INFO split::::casa ##### Begin Task: split ##### 2023-11-04 07:15:19 INFO split::::casa split( vis='testx.ms', outputvis='test9tmp.ms', keepmms=True, field='', spw='', scan='', antenna='', correlation='', timerange='', intent='', array='', uvrange='', observation='', feed='', datacolumn='corrected', keepflags=True, width=1, timebin='0s', combine='' ) 2023-11-04 07:15:19 INFO MSTransformManager::parseMsSpecParams Input file name is testx.ms 2023-11-04 07:15:19 INFO MSTransformManager::parseMsSpecParams Data column is CORRECTED 2023-11-04 07:15:19 INFO MSTransformManager::parseMsSpecParams Output file name is test9tmp.ms 2023-11-04 07:15:19 INFO MSTransformManager::colCheckInfo Adding DATA column to output MS from input CORRECTED_DATA column 2023-11-04 07:15:19 INFO MSTransformManager::open Select data 2023-11-04 07:15:19 INFO MSTransformManager::createOutputMSStructure Create output MS structure 2023-11-04 07:15:20 INFO split::::casa Task split complete. Start time: 2023-11-04 07:15:18.978713 End time: 2023-11-04 07:15:19.666595 2023-11-04 07:15:20 INFO split::::casa ##### End Task: split ##### 2023-11-04 07:15:20 INFO split::::casa ########################################## 2023-11-04 07:15:20 INFO fixvis::::casa ########################################## 2023-11-04 07:15:20 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:20 INFO fixvis::::casa fixvis( vis='test9tmp.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='-0.2333266794430661s 0deg', distances=[], datacolumn='CORRECTED' ) 2023-11-04 07:15:20 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:20 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:20 SEVERE fixvis::::casa Exception Reported: Error in fixvis: Datacolumn CORRECTED_DATA not present 2023-11-04 07:15:20 INFO fixvis::::casa Traceback (most recent call last): 2023-11-04 07:15:20 INFO fixvis::::casa+ File "/wheeldirectory/casa-CAS-14200-7-py3.8.el8/lib/py/lib/python3.8/site-packages/casashell/private/fixvis.py", line 472, in __call__ 2023-11-04 07:15:20 INFO fixvis::::casa+ task_result = _fixvis_t( _pc.document['vis'],_pc.document['outputvis'],_pc.document['field'],_pc.document['refcode'],_pc.document['reuse'],_pc.document['phasecenter'],_pc.document['distances'],_pc.document['datacolumn'] ) 2023-11-04 07:15:20 INFO fixvis::::casa+ File "/wheeldirectory/casa-CAS-14200-7-py3.8.el8/lib/py/lib/python3.8/site-packages/casatasks/private/task_fixvis.py", line 219, in fixvis 2023-11-04 07:15:20 INFO fixvis::::casa+ raise RuntimeError(msg) 2023-11-04 07:15:20 INFO fixvis::::casa+ RuntimeError: Datacolumn CORRECTED_DATA not present 2023-11-04 07:15:20 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:19.673430 End time: 2023-11-04 07:15:19.707595 2023-11-04 07:15:20 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:20 INFO fixvis::::casa ########################################## .2023-11-04 07:15:20 INFO split::::casa ########################################## 2023-11-04 07:15:20 INFO split::::casa ##### Begin Task: split ##### 2023-11-04 07:15:20 INFO split::::casa split( vis='testx.ms', outputvis='test9tmp.ms', keepmms=True, field='', spw='', scan='', antenna='', correlation='', timerange='', intent='', array='', uvrange='', observation='', feed='', datacolumn='corrected', keepflags=True, width=1, timebin='0s', combine='' ) 2023-11-04 07:15:20 INFO MSTransformManager::parseMsSpecParams Input file name is testx.ms 2023-11-04 07:15:20 INFO MSTransformManager::parseMsSpecParams Data column is CORRECTED 2023-11-04 07:15:20 INFO MSTransformManager::parseMsSpecParams Output file name is test9tmp.ms 2023-11-04 07:15:20 INFO MSTransformManager::colCheckInfo Adding DATA column to output MS from input CORRECTED_DATA column 2023-11-04 07:15:20 INFO MSTransformManager::open Select data 2023-11-04 07:15:20 INFO MSTransformManager::createOutputMSStructure Create output MS structure 2023-11-04 07:15:20 INFO split::::casa Task split complete. Start time: 2023-11-04 07:15:19.825076 End time: 2023-11-04 07:15:20.484475 2023-11-04 07:15:20 INFO split::::casa ##### End Task: split ##### 2023-11-04 07:15:20 INFO split::::casa ########################################## 2023-11-04 07:15:20 INFO fixvis::::casa ########################################## 2023-11-04 07:15:20 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:20 INFO fixvis::::casa fixvis( vis='test9tmp.ms', outputvis='testx.ms', field='1', refcode='J2000', reuse=True, phasecenter='-0.2333266794430661s 0deg', distances=[], datacolumn='DATA' ) 2023-11-04 07:15:20 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:20 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:21 INFO fixvis::::casa Will only modify the visibilities in the columns DATA 2023-11-04 07:15:21 INFO fixvis::::casa field: non-blanktest2_0 2023-11-04 07:15:21 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3331 -030.59.29.9987 2023-11-04 07:15:21 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.57062666] [-0.54090662] 2023-11-04 07:15:21 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:21 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2023-11-04 07:15:21 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-04 07:15:21 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0997 -030.59.29.9987 2023-11-04 07:15:21 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706436299021347 -0.5409066175416467 2023-11-04 07:15:21 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 1. 2023-11-04 07:15:21 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:21 INFO SubMS::parseColumnNames() Using DATA column. 2023-11-04 07:15:21 INFO FixVis::setFields Selecting fields 1 2023-11-04 07:15:21 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:21 INFO FixVis::makeSelection() Processing field 1 2023-11-04 07:15:21 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:20.493980 End time: 2023-11-04 07:15:20.695899 2023-11-04 07:15:21 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:21 INFO fixvis::::casa ########################################## 2023-11-04 07:15:21 INFO tclean::::casa ########################################## 2023-11-04 07:15:21 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:21 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:21 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:21 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:21 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:21 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:21 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:21 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:21 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:21 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.97 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:21 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:21 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:21 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:21 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:21 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.014245 2023-11-04 07:15:21 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:21 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:21 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:22 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:22 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:22 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:22 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:22 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:20.774854 End time: 2023-11-04 07:15:21.729735 2023-11-04 07:15:22 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:22 INFO tclean::::casa ########################################## 2023-11-04 07:15:22 INFO imstat::::casa ########################################## 2023-11-04 07:15:22 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:22 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:22 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:22 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:22 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:22 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:22 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:22 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:22 INFO imstat:::: Values --- 2023-11-04 07:15:22 INFO imstat:::: -- flux density [flux]: 0.0112617 Jy 2023-11-04 07:15:22 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:22 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-04 07:15:22 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-04 07:15:22 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-04 07:15:22 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-04 07:15:22 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:22 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-04 07:15:22 INFO imstat:::: -- Sum of pixel values [sum]: 5.84977 Jy/beam 2023-11-04 07:15:22 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-04 07:15:22 INFO imstat:::: Statistics --- 2023-11-04 07:15:22 INFO imstat:::: -- Mean of the pixel values [mean]: 0.000357042 Jy/beam 2023-11-04 07:15:22 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-04 07:15:22 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-04 07:15:22 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-04 07:15:22 INFO imstat:::: -- Median of the pixel values [median]: -0.0988178 Jy/beam 2023-11-04 07:15:22 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649163 Jy/beam 2023-11-04 07:15:22 INFO imstat:::: -- IQR [quartile]: 0.146872 Jy/beam 2023-11-04 07:15:22 INFO imstat:::: -- First quartile [q1]: -0.155428 Jy/beam 2023-11-04 07:15:22 INFO imstat:::: -- Third quartile [q3]: -0.0085559 Jy/beam 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:22 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 5.849774e+00 3.570419e-04 3.646062e-01 3.646172e-01 -1.857701e-01 2.768178e+00 2023-11-04 07:15:22 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:21.753182 End time: 2023-11-04 07:15:21.782164 2023-11-04 07:15:22 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:22 INFO imstat::::casa ########################################## 2023-11-04 07:15:22 INFO tclean::::casa ########################################## 2023-11-04 07:15:22 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:22 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='1', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testz', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:22 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:22 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 1 | [Opened in readonly mode] 2023-11-04 07:15:22 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:22 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testz] : 2023-11-04 07:15:22 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:22 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56667e+07] 2023-11-04 07:15:22 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testz] with ftmachine : gridft 2023-11-04 07:15:22 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:22 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:22 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testz] : clark 2023-11-04 07:15:22 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:22 INFO SIImageStore::calcSensitivity [testz] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:22 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010299 2023-11-04 07:15:22 INFO SIImageStore::printBeamSet Beam : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2023-11-04 07:15:22 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:22 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:23 INFO task_tclean::SDAlgorithmBase::restore [testz] : Restoring model image. 2023-11-04 07:15:23 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:23 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2023-11-04 07:15:23 INFO tclean::::casa Searching for images with prefix 'testz'... Found these, writing history into them: ['testz.model', 'testz.image', 'testz.pb', 'testz.sumwt', 'testz.residual', 'testz.psf'] 2023-11-04 07:15:23 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:21.830226 End time: 2023-11-04 07:15:22.667542 2023-11-04 07:15:23 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:23 INFO tclean::::casa ########################################## 2023-11-04 07:15:23 INFO imstat::::casa ########################################## 2023-11-04 07:15:23 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:23 INFO imstat::::casa imstat( imagename='testz.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:23 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:23 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:23 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:23 INFO imstat::image::statistics Determining stats for image testz.image 2023-11-04 07:15:23 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:23 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.597, -30.59.36.399, I, 2.25982e+11Hz to 18:00:01.610, -30.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.597, -30.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.610, -30.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:23 INFO imstat:::: Values --- 2023-11-04 07:15:23 INFO imstat:::: -- flux density [flux]: 0.187451 Jy 2023-11-04 07:15:23 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:23 INFO imstat:::: -- maximum value [max]: 2.76828 Jy/beam 2023-11-04 07:15:23 INFO imstat:::: -- minimum value [min]: -0.185749 Jy/beam 2023-11-04 07:15:23 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2023-11-04 07:15:23 INFO imstat:::: -- position of min value (pixel) [minpos]: [31, 109, 0, 0] 2023-11-04 07:15:23 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.333, -30.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:23 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.356, -30.59.25.499, I, 2.25982e+11Hz 2023-11-04 07:15:23 INFO imstat:::: -- Sum of pixel values [sum]: 97.5737 Jy/beam 2023-11-04 07:15:23 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2189.16 Jy/beam.Jy/beam 2023-11-04 07:15:23 INFO imstat:::: Statistics --- 2023-11-04 07:15:23 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00595543 Jy/beam 2023-11-04 07:15:23 INFO imstat:::: -- Variance of the pixel values : 0.133588 Jy/beam 2023-11-04 07:15:23 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365497 Jy/beam 2023-11-04 07:15:23 INFO imstat:::: -- Root mean square [rms]: 0.365535 Jy/beam 2023-11-04 07:15:23 INFO imstat:::: -- Median of the pixel values [median]: -0.0939644 Jy/beam 2023-11-04 07:15:23 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0676308 Jy/beam 2023-11-04 07:15:23 INFO imstat:::: -- IQR [quartile]: 0.155104 Jy/beam 2023-11-04 07:15:23 INFO imstat:::: -- First quartile [q1]: -0.153538 Jy/beam 2023-11-04 07:15:23 INFO imstat:::: -- Third quartile [q3]: 0.00156572 Jy/beam 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:23 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 9.757369e+01 5.955426e-03 3.655346e-01 3.654972e-01 -1.857485e-01 2.768277e+00 2023-11-04 07:15:23 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:22.669455 End time: 2023-11-04 07:15:22.686104 2023-11-04 07:15:23 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:23 INFO imstat::::casa ########################################## .2023-11-04 07:15:23 INFO split::::casa ########################################## 2023-11-04 07:15:23 INFO split::::casa ##### Begin Task: split ##### 2023-11-04 07:15:23 INFO split::::casa split( vis='testx.ms', outputvis='test9tmp.ms', keepmms=True, field='', spw='', scan='', antenna='', correlation='', timerange='', intent='', array='', uvrange='', observation='', feed='', datacolumn='corrected', keepflags=True, width=1, timebin='0s', combine='' ) 2023-11-04 07:15:23 INFO MSTransformManager::parseMsSpecParams Input file name is testx.ms 2023-11-04 07:15:23 INFO MSTransformManager::parseMsSpecParams Data column is CORRECTED 2023-11-04 07:15:23 INFO MSTransformManager::parseMsSpecParams Output file name is test9tmp.ms 2023-11-04 07:15:23 INFO MSTransformManager::colCheckInfo Adding DATA column to output MS from input CORRECTED_DATA column 2023-11-04 07:15:23 INFO MSTransformManager::open Select data 2023-11-04 07:15:23 INFO MSTransformManager::createOutputMSStructure Create output MS structure 2023-11-04 07:15:23 INFO split::::casa Task split complete. Start time: 2023-11-04 07:15:22.812383 End time: 2023-11-04 07:15:23.472072 2023-11-04 07:15:23 INFO split::::casa ##### End Task: split ##### 2023-11-04 07:15:23 INFO split::::casa ########################################## 2023-11-04 07:15:23 INFO fixvis::::casa ########################################## 2023-11-04 07:15:23 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:23 INFO fixvis::::casa fixvis( vis='test9tmp.ms', outputvis='testx.ms', field='1', refcode='J2000', reuse=True, phasecenter='-0.2333266794430661s 0deg', distances=[], datacolumn='all' ) 2023-11-04 07:15:23 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:23 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:24 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-04 07:15:24 INFO fixvis::::casa field: non-blanktest2_0 2023-11-04 07:15:24 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3331 -030.59.29.9987 2023-11-04 07:15:24 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.57062666] [-0.54090662] 2023-11-04 07:15:24 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:24 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2023-11-04 07:15:24 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-04 07:15:24 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0997 -030.59.29.9987 2023-11-04 07:15:24 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706436299021347 -0.5409066175416467 2023-11-04 07:15:24 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 1. 2023-11-04 07:15:24 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:24 INFO SubMS::parseColumnNames() Using DATA column. 2023-11-04 07:15:24 INFO FixVis::setFields Selecting fields 1 2023-11-04 07:15:24 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:24 INFO FixVis::makeSelection() Processing field 1 2023-11-04 07:15:24 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:23.479704 End time: 2023-11-04 07:15:23.673885 2023-11-04 07:15:24 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:24 INFO fixvis::::casa ########################################## 2023-11-04 07:15:24 INFO tclean::::casa ########################################## 2023-11-04 07:15:24 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:24 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:24 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:24 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:24 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:24 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:24 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:24 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:24 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:24 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:24 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:24 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:24 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:24 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:24 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010348 2023-11-04 07:15:24 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:24 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:24 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:24 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:24 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:25 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:25 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:25 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:23.726774 End time: 2023-11-04 07:15:24.619085 2023-11-04 07:15:25 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:25 INFO tclean::::casa ########################################## 2023-11-04 07:15:25 INFO imstat::::casa ########################################## 2023-11-04 07:15:25 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:25 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:25 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:25 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:25 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:25 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:25 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:25 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:25 INFO imstat:::: Values --- 2023-11-04 07:15:25 INFO imstat:::: -- flux density [flux]: 0.0112617 Jy 2023-11-04 07:15:25 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:25 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-04 07:15:25 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-04 07:15:25 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-04 07:15:25 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-04 07:15:25 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:25 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-04 07:15:25 INFO imstat:::: -- Sum of pixel values [sum]: 5.84977 Jy/beam 2023-11-04 07:15:25 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-04 07:15:25 INFO imstat:::: Statistics --- 2023-11-04 07:15:25 INFO imstat:::: -- Mean of the pixel values [mean]: 0.000357042 Jy/beam 2023-11-04 07:15:25 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-04 07:15:25 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-04 07:15:25 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-04 07:15:25 INFO imstat:::: -- Median of the pixel values [median]: -0.0988178 Jy/beam 2023-11-04 07:15:25 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649163 Jy/beam 2023-11-04 07:15:25 INFO imstat:::: -- IQR [quartile]: 0.146872 Jy/beam 2023-11-04 07:15:25 INFO imstat:::: -- First quartile [q1]: -0.155428 Jy/beam 2023-11-04 07:15:25 INFO imstat:::: -- Third quartile [q3]: -0.0085559 Jy/beam 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:25 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 5.849774e+00 3.570419e-04 3.646062e-01 3.646172e-01 -1.857701e-01 2.768178e+00 2023-11-04 07:15:25 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:24.621176 End time: 2023-11-04 07:15:24.637658 2023-11-04 07:15:25 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:25 INFO imstat::::casa ########################################## 2023-11-04 07:15:25 INFO tclean::::casa ########################################## 2023-11-04 07:15:25 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:25 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='1', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testz', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:25 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:25 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 1 | [Opened in readonly mode] 2023-11-04 07:15:25 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:25 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testz] : 2023-11-04 07:15:25 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:25 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56667e+07] 2023-11-04 07:15:25 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testz] with ftmachine : gridft 2023-11-04 07:15:25 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:25 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:25 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testz] : clark 2023-11-04 07:15:25 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:25 INFO SIImageStore::calcSensitivity [testz] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:25 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010299 2023-11-04 07:15:25 INFO SIImageStore::printBeamSet Beam : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2023-11-04 07:15:25 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:25 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:25 INFO task_tclean::SDAlgorithmBase::restore [testz] : Restoring model image. 2023-11-04 07:15:25 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:25 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2023-11-04 07:15:25 INFO tclean::::casa Searching for images with prefix 'testz'... Found these, writing history into them: ['testz.model', 'testz.image', 'testz.pb', 'testz.sumwt', 'testz.residual', 'testz.psf'] 2023-11-04 07:15:25 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:24.678141 End time: 2023-11-04 07:15:25.495147 2023-11-04 07:15:25 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:25 INFO tclean::::casa ########################################## 2023-11-04 07:15:25 INFO imstat::::casa ########################################## 2023-11-04 07:15:25 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:25 INFO imstat::::casa imstat( imagename='testz.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:26 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:26 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:26 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:26 INFO imstat::image::statistics Determining stats for image testz.image 2023-11-04 07:15:26 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:26 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.597, -30.59.36.399, I, 2.25982e+11Hz to 18:00:01.610, -30.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.597, -30.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.610, -30.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:26 INFO imstat:::: Values --- 2023-11-04 07:15:26 INFO imstat:::: -- flux density [flux]: 0.187451 Jy 2023-11-04 07:15:26 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:26 INFO imstat:::: -- maximum value [max]: 2.76828 Jy/beam 2023-11-04 07:15:26 INFO imstat:::: -- minimum value [min]: -0.185749 Jy/beam 2023-11-04 07:15:26 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2023-11-04 07:15:26 INFO imstat:::: -- position of min value (pixel) [minpos]: [31, 109, 0, 0] 2023-11-04 07:15:26 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.333, -30.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:26 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.356, -30.59.25.499, I, 2.25982e+11Hz 2023-11-04 07:15:26 INFO imstat:::: -- Sum of pixel values [sum]: 97.5737 Jy/beam 2023-11-04 07:15:26 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2189.16 Jy/beam.Jy/beam 2023-11-04 07:15:26 INFO imstat:::: Statistics --- 2023-11-04 07:15:26 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00595543 Jy/beam 2023-11-04 07:15:26 INFO imstat:::: -- Variance of the pixel values : 0.133588 Jy/beam 2023-11-04 07:15:26 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365497 Jy/beam 2023-11-04 07:15:26 INFO imstat:::: -- Root mean square [rms]: 0.365535 Jy/beam 2023-11-04 07:15:26 INFO imstat:::: -- Median of the pixel values [median]: -0.0939644 Jy/beam 2023-11-04 07:15:26 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0676308 Jy/beam 2023-11-04 07:15:26 INFO imstat:::: -- IQR [quartile]: 0.155104 Jy/beam 2023-11-04 07:15:26 INFO imstat:::: -- First quartile [q1]: -0.153538 Jy/beam 2023-11-04 07:15:26 INFO imstat:::: -- Third quartile [q3]: 0.00156572 Jy/beam 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:26 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 9.757369e+01 5.955426e-03 3.655346e-01 3.654972e-01 -1.857485e-01 2.768277e+00 2023-11-04 07:15:26 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:25.496991 End time: 2023-11-04 07:15:25.523917 2023-11-04 07:15:26 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:26 INFO imstat::::casa ########################################## .2023-11-04 07:15:26 INFO fixvis::::casa ########################################## 2023-11-04 07:15:26 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:26 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='outms3.ms', field='0', refcode='J2000', reuse=True, phasecenter='J2000 18h00m02.3092s -29d49m29.9987s', distances=[], datacolumn='all' ) 2023-11-04 07:15:26 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:26 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:26 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-04 07:15:26 INFO fixvis::::casa field: non-blanktest_0 2023-11-04 07:15:26 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:26 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-04 07:15:26 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:26 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.49.29.9987 2023-11-04 07:15:26 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5205444431047309 2023-11-04 07:15:26 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-04 07:15:26 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/outms3.ms 2023-11-04 07:15:26 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-04 07:15:26 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:26 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:26 INFO FixVis::makeSelection() Processing field 0 2023-11-04 07:15:26 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:25.652071 End time: 2023-11-04 07:15:25.960564 2023-11-04 07:15:26 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:26 INFO fixvis::::casa ########################################## 2023-11-04 07:15:26 INFO fixvis::::casa ########################################## 2023-11-04 07:15:26 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:26 INFO fixvis::::casa fixvis( vis='outms3.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='J2000 18h00m02.3092s -29d59m29.9987s', distances=[], datacolumn='all' ) 2023-11-04 07:15:26 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:26 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:26 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-04 07:15:26 INFO fixvis::::casa field: non-blanktest_0 2023-11-04 07:15:26 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.49.29.9987 2023-11-04 07:15:26 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52054444] 2023-11-04 07:15:26 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:26 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:26 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2023-11-04 07:15:26 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-04 07:15:26 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:26 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-04 07:15:26 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:26 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:26 INFO FixVis::makeSelection() Processing field 0 2023-11-04 07:15:26 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:25.961804 End time: 2023-11-04 07:15:26.177522 2023-11-04 07:15:26 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:26 INFO fixvis::::casa ########################################## 2023-11-04 07:15:26 INFO tclean::::casa ########################################## 2023-11-04 07:15:26 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:26 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:26 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:26 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:26 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:26 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:26 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:26 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:26 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:26 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:26 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:26 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:26 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:27 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:27 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.01063 2023-11-04 07:15:27 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:27 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:27 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:27 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:27 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:27 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:27 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:27 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:26.216369 End time: 2023-11-04 07:15:27.200597 2023-11-04 07:15:27 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:27 INFO tclean::::casa ########################################## 2023-11-04 07:15:27 INFO imstat::::casa ########################################## 2023-11-04 07:15:27 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:27 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:27 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:27 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:27 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:27 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:27 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:27 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:27 INFO imstat:::: Values --- 2023-11-04 07:15:27 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2023-11-04 07:15:27 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:27 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-04 07:15:27 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-04 07:15:27 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-04 07:15:27 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-04 07:15:27 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:27 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-04 07:15:27 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2023-11-04 07:15:27 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-04 07:15:27 INFO imstat:::: Statistics --- 2023-11-04 07:15:27 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2023-11-04 07:15:27 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-04 07:15:27 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-04 07:15:27 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-04 07:15:27 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2023-11-04 07:15:27 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2023-11-04 07:15:27 INFO imstat:::: -- IQR [quartile]: 0.14686 Jy/beam 2023-11-04 07:15:27 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2023-11-04 07:15:27 INFO imstat:::: -- Third quartile [q3]: -0.00856261 Jy/beam 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:27 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 5.849415e+00 3.570200e-04 3.646062e-01 3.646172e-01 -1.857702e-01 2.768178e+00 2023-11-04 07:15:27 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:27.204449 End time: 2023-11-04 07:15:27.236511 2023-11-04 07:15:27 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:27 INFO imstat::::casa ########################################## 2023-11-04 07:15:27 INFO fixvis::::casa ########################################## 2023-11-04 07:15:27 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:27 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='J2000 18h00m02.3092s -29d59m29.9987s', distances='', datacolumn='all' ) 2023-11-04 07:15:27 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:27 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:27 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-04 07:15:27 INFO fixvis::::casa field: non-blanktest_0 2023-11-04 07:15:27 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:27 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-04 07:15:27 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:27 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:27 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2023-11-04 07:15:27 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-04 07:15:27 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:27 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-04 07:15:27 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:27 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:27 INFO FixVis::makeSelection() Processing field 0 2023-11-04 07:15:28 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:27.254471 End time: 2023-11-04 07:15:27.519545 2023-11-04 07:15:28 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:28 INFO fixvis::::casa ########################################## 2023-11-04 07:15:28 INFO tclean::::casa ########################################## 2023-11-04 07:15:28 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:28 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:28 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:28 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:28 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:28 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:28 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:28 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:28 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:28 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:28 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:28 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:28 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:28 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:28 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010361 2023-11-04 07:15:28 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:28 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:28 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:28 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:28 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:28 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:28 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:28 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:27.560669 End time: 2023-11-04 07:15:28.491477 2023-11-04 07:15:28 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:28 INFO tclean::::casa ########################################## 2023-11-04 07:15:28 INFO imstat::::casa ########################################## 2023-11-04 07:15:28 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:28 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:28 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:28 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:28 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:29 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:29 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:29 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:29 INFO imstat:::: Values --- 2023-11-04 07:15:29 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2023-11-04 07:15:29 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:29 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-04 07:15:29 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-04 07:15:29 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-04 07:15:29 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-04 07:15:29 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:29 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-04 07:15:29 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2023-11-04 07:15:29 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-04 07:15:29 INFO imstat:::: Statistics --- 2023-11-04 07:15:29 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2023-11-04 07:15:29 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-04 07:15:29 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-04 07:15:29 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-04 07:15:29 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2023-11-04 07:15:29 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2023-11-04 07:15:29 INFO imstat:::: -- IQR [quartile]: 0.146861 Jy/beam 2023-11-04 07:15:29 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2023-11-04 07:15:29 INFO imstat:::: -- Third quartile [q3]: -0.0085626 Jy/beam 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:29 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 5.849419e+00 3.570202e-04 3.646062e-01 3.646172e-01 -1.857702e-01 2.768178e+00 2023-11-04 07:15:29 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:28.493438 End time: 2023-11-04 07:15:28.510490 2023-11-04 07:15:29 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:29 INFO imstat::::casa ########################################## .2023-11-04 07:15:29 INFO fixvis::::casa ########################################## 2023-11-04 07:15:29 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:29 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='J2000 18h00m02.3092s -29d59m29.9987s', distances=['100AU'], datacolumn='all' ) 2023-11-04 07:15:29 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:29 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:29 INFO fixvis::::casa Will refocus to the given distances: ['100AU'] 2023-11-04 07:15:29 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-04 07:15:29 INFO fixvis::::casa field: non-blanktest_0 2023-11-04 07:15:29 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:29 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-04 07:15:29 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:29 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:29 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2023-11-04 07:15:29 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-04 07:15:29 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:29 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-04 07:15:29 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:29 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:29 INFO FixVis::makeSelection() Processing field 0 2023-11-04 07:15:29 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:28.588010 End time: 2023-11-04 07:15:28.803993 2023-11-04 07:15:29 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:29 INFO fixvis::::casa ########################################## 2023-11-04 07:15:29 INFO tclean::::casa ########################################## 2023-11-04 07:15:29 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:29 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:29 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:29 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:29 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:29 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:29 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:29 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:29 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:29 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:29 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:29 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:29 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:29 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:29 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.025019 2023-11-04 07:15:29 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:29 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:30 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:30 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:30 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:30 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:30 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:30 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:28.837325 End time: 2023-11-04 07:15:29.918465 2023-11-04 07:15:30 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:30 INFO tclean::::casa ########################################## 2023-11-04 07:15:30 INFO imstat::::casa ########################################## 2023-11-04 07:15:30 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:30 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:30 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:30 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:30 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:30 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:30 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:30 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:30 INFO imstat:::: Values --- 2023-11-04 07:15:30 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2023-11-04 07:15:30 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:30 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-04 07:15:30 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-04 07:15:30 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-04 07:15:30 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-04 07:15:30 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:30 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-04 07:15:30 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2023-11-04 07:15:30 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-04 07:15:30 INFO imstat:::: Statistics --- 2023-11-04 07:15:30 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2023-11-04 07:15:30 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-04 07:15:30 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-04 07:15:30 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-04 07:15:30 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2023-11-04 07:15:30 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2023-11-04 07:15:30 INFO imstat:::: -- IQR [quartile]: 0.146861 Jy/beam 2023-11-04 07:15:30 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2023-11-04 07:15:30 INFO imstat:::: -- Third quartile [q3]: -0.0085626 Jy/beam 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:30 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 5.849419e+00 3.570202e-04 3.646062e-01 3.646172e-01 -1.857702e-01 2.768178e+00 2023-11-04 07:15:30 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:29.920473 End time: 2023-11-04 07:15:29.937316 2023-11-04 07:15:30 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:30 INFO imstat::::casa ########################################## .2023-11-04 07:15:30 INFO fixvis::::casa ########################################## 2023-11-04 07:15:30 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:30 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='J2000 18h00m02.3092s -29d59m29.9987s', distances='100AU', datacolumn='all' ) 2023-11-04 07:15:30 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:30 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:30 INFO fixvis::::casa Will refocus to the given distances: 100AU 2023-11-04 07:15:30 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-04 07:15:30 INFO fixvis::::casa field: non-blanktest_0 2023-11-04 07:15:30 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:30 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-04 07:15:30 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:30 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:30 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2023-11-04 07:15:30 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-04 07:15:30 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:30 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-04 07:15:30 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:30 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:30 INFO FixVis::makeSelection() Processing field 0 2023-11-04 07:15:30 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:30.023206 End time: 2023-11-04 07:15:30.245669 2023-11-04 07:15:30 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:30 INFO fixvis::::casa ########################################## 2023-11-04 07:15:30 INFO tclean::::casa ########################################## 2023-11-04 07:15:30 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:30 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:30 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:30 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:30 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:30 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:30 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:30 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:30 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:30 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:30 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:30 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:30 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:31 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:31 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010271 2023-11-04 07:15:31 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:31 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:31 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:31 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:31 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:31 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:31 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:31 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:30.284790 End time: 2023-11-04 07:15:31.202147 2023-11-04 07:15:31 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:31 INFO tclean::::casa ########################################## 2023-11-04 07:15:31 INFO imstat::::casa ########################################## 2023-11-04 07:15:31 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:31 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:31 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:31 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:31 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:31 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:31 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:31 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:31 INFO imstat:::: Values --- 2023-11-04 07:15:31 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2023-11-04 07:15:31 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:31 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-04 07:15:31 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-04 07:15:31 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-04 07:15:31 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-04 07:15:31 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:31 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-04 07:15:31 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2023-11-04 07:15:31 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-04 07:15:31 INFO imstat:::: Statistics --- 2023-11-04 07:15:31 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2023-11-04 07:15:31 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-04 07:15:31 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-04 07:15:31 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-04 07:15:31 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2023-11-04 07:15:31 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2023-11-04 07:15:31 INFO imstat:::: -- IQR [quartile]: 0.146861 Jy/beam 2023-11-04 07:15:31 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2023-11-04 07:15:31 INFO imstat:::: -- Third quartile [q3]: -0.0085626 Jy/beam 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:31 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 5.849419e+00 3.570202e-04 3.646062e-01 3.646172e-01 -1.857702e-01 2.768178e+00 2023-11-04 07:15:31 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:31.205742 End time: 2023-11-04 07:15:31.236496 2023-11-04 07:15:31 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:31 INFO imstat::::casa ########################################## .2023-11-04 07:15:31 INFO fixvis::::casa ########################################## 2023-11-04 07:15:31 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:31 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='J2000 18h00m02.3092s -29d59m29.9987s', distances='', datacolumn='all' ) 2023-11-04 07:15:31 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:31 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:31 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-04 07:15:31 INFO fixvis::::casa field: non-blanktest_0 2023-11-04 07:15:31 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:31 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-04 07:15:31 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:31 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:31 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2023-11-04 07:15:31 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-04 07:15:31 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:31 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-04 07:15:31 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:31 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:31 INFO FixVis::makeSelection() Processing field 0 2023-11-04 07:15:32 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:31.346837 End time: 2023-11-04 07:15:31.573342 2023-11-04 07:15:32 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:32 INFO fixvis::::casa ########################################## 2023-11-04 07:15:32 INFO tclean::::casa ########################################## 2023-11-04 07:15:32 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:32 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:32 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:32 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:32 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:32 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:32 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:32 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:32 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:32 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:32 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:32 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:32 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:32 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:32 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.022831 2023-11-04 07:15:32 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:32 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:32 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:33 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:33 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:33 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:33 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:33 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:31.612801 End time: 2023-11-04 07:15:32.652592 2023-11-04 07:15:33 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:33 INFO tclean::::casa ########################################## 2023-11-04 07:15:33 INFO imstat::::casa ########################################## 2023-11-04 07:15:33 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:33 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:33 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:33 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:33 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:33 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:33 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:33 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:33 INFO imstat:::: Values --- 2023-11-04 07:15:33 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2023-11-04 07:15:33 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:33 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-04 07:15:33 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-04 07:15:33 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-04 07:15:33 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-04 07:15:33 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:33 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-04 07:15:33 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2023-11-04 07:15:33 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-04 07:15:33 INFO imstat:::: Statistics --- 2023-11-04 07:15:33 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2023-11-04 07:15:33 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-04 07:15:33 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-04 07:15:33 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-04 07:15:33 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2023-11-04 07:15:33 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2023-11-04 07:15:33 INFO imstat:::: -- IQR [quartile]: 0.146861 Jy/beam 2023-11-04 07:15:33 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2023-11-04 07:15:33 INFO imstat:::: -- Third quartile [q3]: -0.0085626 Jy/beam 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:33 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 5.849419e+00 3.570202e-04 3.646062e-01 3.646172e-01 -1.857702e-01 2.768178e+00 2023-11-04 07:15:33 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:32.654474 End time: 2023-11-04 07:15:32.671507 2023-11-04 07:15:33 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:33 INFO imstat::::casa ########################################## .2023-11-04 07:15:33 INFO fixvis::::casa ########################################## 2023-11-04 07:15:33 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:33 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='J2000 18h00m02.3092s -29d59m26.9987s', distances='', datacolumn='all' ) 2023-11-04 07:15:33 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:33 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:33 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-04 07:15:33 INFO fixvis::::casa field: non-blanktest_0 2023-11-04 07:15:33 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:33 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-04 07:15:33 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:33 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.26.9987 2023-11-04 07:15:33 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234387807809548 2023-11-04 07:15:33 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-04 07:15:33 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:33 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-04 07:15:33 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:33 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:33 INFO FixVis::makeSelection() Processing field 0 2023-11-04 07:15:33 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:32.762190 End time: 2023-11-04 07:15:32.981980 2023-11-04 07:15:33 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:33 INFO fixvis::::casa ########################################## 2023-11-04 07:15:33 INFO tclean::::casa ########################################## 2023-11-04 07:15:33 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:33 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:33 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:33 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:33 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:33 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:33 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:33 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:33 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:33 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:33 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:33 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:33 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:34 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:34 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.01035 2023-11-04 07:15:34 INFO SIImageStore::printBeamSet Beam : 2.2794 arcsec, 2.01116 arcsec, 30.5337 deg 2023-11-04 07:15:34 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:34 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:34 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:34 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:34 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.2794 arcsec, 2.01116 arcsec, 30.5337 deg 2023-11-04 07:15:34 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:34 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:33.020169 End time: 2023-11-04 07:15:34.083487 2023-11-04 07:15:34 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:34 INFO tclean::::casa ########################################## 2023-11-04 07:15:34 INFO imstat::::casa ########################################## 2023-11-04 07:15:34 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:34 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:34 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:34 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:34 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:34 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:34 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:34 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.33.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.20.699, I, 2.25982e+11Hz) 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.33.399, I, 2.25982e+11Hz 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.20.699, I, 2.25982e+11Hz 2023-11-04 07:15:34 INFO imstat:::: Values --- 2023-11-04 07:15:34 INFO imstat:::: -- flux density [flux]: 0.249791 Jy 2023-11-04 07:15:34 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:34 INFO imstat:::: -- maximum value [max]: 2.76821 Jy/beam 2023-11-04 07:15:34 INFO imstat:::: -- minimum value [min]: -0.185766 Jy/beam 2023-11-04 07:15:34 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 34, 0, 0] 2023-11-04 07:15:34 INFO imstat:::: -- position of min value (pixel) [minpos]: [61, 79, 0, 0] 2023-11-04 07:15:34 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:34 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.332, -29.59.25.499, I, 2.25982e+11Hz 2023-11-04 07:15:34 INFO imstat:::: -- Sum of pixel values [sum]: 129.75 Jy/beam 2023-11-04 07:15:34 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2161.57 Jy/beam.Jy/beam 2023-11-04 07:15:34 INFO imstat:::: Statistics --- 2023-11-04 07:15:34 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00791933 Jy/beam 2023-11-04 07:15:34 INFO imstat:::: -- Variance of the pixel values : 0.131877 Jy/beam 2023-11-04 07:15:34 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.363149 Jy/beam 2023-11-04 07:15:34 INFO imstat:::: -- Root mean square [rms]: 0.363224 Jy/beam 2023-11-04 07:15:34 INFO imstat:::: -- Median of the pixel values [median]: -0.0875305 Jy/beam 2023-11-04 07:15:34 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0698687 Jy/beam 2023-11-04 07:15:34 INFO imstat:::: -- IQR [quartile]: 0.153301 Jy/beam 2023-11-04 07:15:34 INFO imstat:::: -- First quartile [q1]: -0.149306 Jy/beam 2023-11-04 07:15:34 INFO imstat:::: -- Third quartile [q3]: 0.00399534 Jy/beam 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:34 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 1.297503e+02 7.919332e-03 3.632243e-01 3.631491e-01 -1.857659e-01 2.768206e+00 2023-11-04 07:15:34 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:34.085391 End time: 2023-11-04 07:15:34.102291 2023-11-04 07:15:34 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:34 INFO imstat::::casa ########################################## .2023-11-04 07:15:34 INFO fixvis::::casa ########################################## 2023-11-04 07:15:34 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:34 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='0h -0d0m3s', distances='', datacolumn='all' ) 2023-11-04 07:15:34 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:34 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:34 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-04 07:15:34 INFO fixvis::::casa field: non-blanktest_0 2023-11-04 07:15:34 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:34 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-04 07:15:34 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:34 INFO fixvis::::casa No equinox given in parameter 'phasecenter': 0h -0d0m3s 2023-11-04 07:15:34 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-04 07:15:34 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.32.9987 2023-11-04 07:15:34 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283963048133 -0.5234678696793916 2023-11-04 07:15:34 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-04 07:15:34 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:34 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-04 07:15:34 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:34 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:34 INFO FixVis::makeSelection() Processing field 0 2023-11-04 07:15:34 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:34.186472 End time: 2023-11-04 07:15:34.472503 2023-11-04 07:15:34 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:34 INFO fixvis::::casa ########################################## 2023-11-04 07:15:35 INFO tclean::::casa ########################################## 2023-11-04 07:15:35 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:35 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:35 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:35 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:35 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:35 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:35 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:35 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:35 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:35 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:35 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:35 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:35 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:35 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:35 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.021684 2023-11-04 07:15:35 INFO SIImageStore::printBeamSet Beam : 2.27942 arcsec, 2.01116 arcsec, 30.5324 deg 2023-11-04 07:15:35 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:35 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:36 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:36 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:36 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27942 arcsec, 2.01116 arcsec, 30.5324 deg 2023-11-04 07:15:36 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:36 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:34.515236 End time: 2023-11-04 07:15:35.729215 2023-11-04 07:15:36 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:36 INFO tclean::::casa ########################################## 2023-11-04 07:15:36 INFO imstat::::casa ########################################## 2023-11-04 07:15:36 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:36 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:36 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:36 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:36 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:36 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:36 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:36 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.39.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.26.699, I, 2.25982e+11Hz) 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.39.399, I, 2.25982e+11Hz 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.26.699, I, 2.25982e+11Hz 2023-11-04 07:15:36 INFO imstat:::: Values --- 2023-11-04 07:15:36 INFO imstat:::: -- flux density [flux]: 0.274462 Jy 2023-11-04 07:15:36 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:36 INFO imstat:::: -- maximum value [max]: 2.76821 Jy/beam 2023-11-04 07:15:36 INFO imstat:::: -- minimum value [min]: -0.185768 Jy/beam 2023-11-04 07:15:36 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 94, 0, 0] 2023-11-04 07:15:36 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 49, 0, 0] 2023-11-04 07:15:36 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:36 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-04 07:15:36 INFO imstat:::: -- Sum of pixel values [sum]: 142.567 Jy/beam 2023-11-04 07:15:36 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2160.11 Jy/beam.Jy/beam 2023-11-04 07:15:36 INFO imstat:::: Statistics --- 2023-11-04 07:15:36 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00870157 Jy/beam 2023-11-04 07:15:36 INFO imstat:::: -- Variance of the pixel values : 0.131775 Jy/beam 2023-11-04 07:15:36 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.363008 Jy/beam 2023-11-04 07:15:36 INFO imstat:::: -- Root mean square [rms]: 0.363102 Jy/beam 2023-11-04 07:15:36 INFO imstat:::: -- Median of the pixel values [median]: -0.0857827 Jy/beam 2023-11-04 07:15:36 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0706394 Jy/beam 2023-11-04 07:15:36 INFO imstat:::: -- IQR [quartile]: 0.154155 Jy/beam 2023-11-04 07:15:36 INFO imstat:::: -- First quartile [q1]: -0.148851 Jy/beam 2023-11-04 07:15:36 INFO imstat:::: -- Third quartile [q3]: 0.00530347 Jy/beam 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:36 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 1.425666e+02 8.701574e-03 3.631017e-01 3.630085e-01 -1.857684e-01 2.768206e+00 2023-11-04 07:15:36 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:35.733320 End time: 2023-11-04 07:15:35.767374 2023-11-04 07:15:36 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:36 INFO imstat::::casa ########################################## .2023-11-04 07:15:36 INFO fixvis::::casa ########################################## 2023-11-04 07:15:36 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:36 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='J2000 18h00m02.5401s -29d59m29.9987s', distances='', datacolumn='all' ) 2023-11-04 07:15:36 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:36 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:36 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-04 07:15:36 INFO fixvis::::casa field: non-blanktest_0 2023-11-04 07:15:36 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:36 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-04 07:15:36 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:36 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.5401 -029.59.29.9987 2023-11-04 07:15:36 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706116055101893 -0.5234533251913881 2023-11-04 07:15:36 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-04 07:15:36 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:36 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-04 07:15:36 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:36 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:36 INFO FixVis::makeSelection() Processing field 0 2023-11-04 07:15:36 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:35.887491 End time: 2023-11-04 07:15:36.149908 2023-11-04 07:15:36 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:36 INFO fixvis::::casa ########################################## 2023-11-04 07:15:36 INFO tclean::::casa ########################################## 2023-11-04 07:15:36 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:36 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:36 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:36 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:36 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:36 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:36 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:36 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:36 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:36 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:36 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:36 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:36 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:37 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:37 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010334 2023-11-04 07:15:37 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01117 arcsec, 30.5334 deg 2023-11-04 07:15:37 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:37 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:37 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:37 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:37 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01117 arcsec, 30.5334 deg 2023-11-04 07:15:37 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:37 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:36.197148 End time: 2023-11-04 07:15:37.174650 2023-11-04 07:15:37 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:37 INFO tclean::::casa ########################################## 2023-11-04 07:15:37 INFO imstat::::casa ########################################## 2023-11-04 07:15:37 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:37 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:37 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:37 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:37 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:37 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:37 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:37 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:03.033, -29.59.36.399, I, 2.25982e+11Hz to 18:00:02.055, -29.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:03.033, -29.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:02.055, -29.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:37 INFO imstat:::: Values --- 2023-11-04 07:15:37 INFO imstat:::: -- flux density [flux]: 0.216271 Jy 2023-11-04 07:15:37 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:37 INFO imstat:::: -- maximum value [max]: 2.76821 Jy/beam 2023-11-04 07:15:37 INFO imstat:::: -- minimum value [min]: -0.185787 Jy/beam 2023-11-04 07:15:37 INFO imstat:::: -- position of max value (pixel) [maxpos]: [94, 64, 0, 0] 2023-11-04 07:15:37 INFO imstat:::: -- position of min value (pixel) [minpos]: [97, 19, 0, 0] 2023-11-04 07:15:37 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:37 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-04 07:15:37 INFO imstat:::: -- Sum of pixel values [sum]: 112.34 Jy/beam 2023-11-04 07:15:37 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2183.54 Jy/beam.Jy/beam 2023-11-04 07:15:37 INFO imstat:::: Statistics --- 2023-11-04 07:15:37 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00685667 Jy/beam 2023-11-04 07:15:37 INFO imstat:::: -- Variance of the pixel values : 0.133234 Jy/beam 2023-11-04 07:15:37 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365012 Jy/beam 2023-11-04 07:15:37 INFO imstat:::: -- Root mean square [rms]: 0.365066 Jy/beam 2023-11-04 07:15:37 INFO imstat:::: -- Median of the pixel values [median]: -0.092606 Jy/beam 2023-11-04 07:15:37 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0686032 Jy/beam 2023-11-04 07:15:37 INFO imstat:::: -- IQR [quartile]: 0.156016 Jy/beam 2023-11-04 07:15:37 INFO imstat:::: -- First quartile [q1]: -0.152775 Jy/beam 2023-11-04 07:15:37 INFO imstat:::: -- Third quartile [q3]: 0.00324101 Jy/beam 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:37 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 1.123397e+02 6.856670e-03 3.650655e-01 3.650123e-01 -1.857871e-01 2.768210e+00 2023-11-04 07:15:37 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:37.178784 End time: 2023-11-04 07:15:37.211411 2023-11-04 07:15:37 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:37 INFO imstat::::casa ########################################## .2023-11-04 07:15:37 INFO fixvis::::casa ########################################## 2023-11-04 07:15:37 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:37 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='-0d0m3s 0deg', distances='', datacolumn='all' ) 2023-11-04 07:15:37 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:37 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-04 07:15:37 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-04 07:15:37 INFO fixvis::::casa field: non-blanktest_0 2023-11-04 07:15:37 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:37 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-04 07:15:37 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:37 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0d0m3s 0deg 2023-11-04 07:15:37 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-04 07:15:37 INFO fixvis::::casa RA offset is an angle (not a time). Will divide by cos(DEC) to compute time offset. 2023-11-04 07:15:37 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0783 -029.59.29.9987 2023-11-04 07:15:37 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706451893538216 -0.5234533252689584 2023-11-04 07:15:37 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-04 07:15:37 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:37 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-04 07:15:37 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:37 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:37 INFO FixVis::makeSelection() Processing field 0 2023-11-04 07:15:38 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:37.319696 End time: 2023-11-04 07:15:37.545787 2023-11-04 07:15:38 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:38 INFO fixvis::::casa ########################################## 2023-11-04 07:15:38 INFO tclean::::casa ########################################## 2023-11-04 07:15:38 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:38 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:38 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:38 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:38 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:38 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:38 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:38 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:38 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:38 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:38 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:38 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:38 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:38 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:38 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.020443 2023-11-04 07:15:38 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2023-11-04 07:15:38 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:38 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:39 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:39 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:39 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2023-11-04 07:15:39 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:39 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:37.585053 End time: 2023-11-04 07:15:38.721744 2023-11-04 07:15:39 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:39 INFO tclean::::casa ########################################## 2023-11-04 07:15:39 INFO imstat::::casa ########################################## 2023-11-04 07:15:39 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:39 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:39 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:39 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:39 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:39 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:39 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:39 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.571, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.593, -29.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.571, -29.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.593, -29.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:39 INFO imstat:::: Values --- 2023-11-04 07:15:39 INFO imstat:::: -- flux density [flux]: 0.179787 Jy 2023-11-04 07:15:39 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:39 INFO imstat:::: -- maximum value [max]: 2.76821 Jy/beam 2023-11-04 07:15:39 INFO imstat:::: -- minimum value [min]: -0.185809 Jy/beam 2023-11-04 07:15:39 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2023-11-04 07:15:39 INFO imstat:::: -- position of min value (pixel) [minpos]: [31, 109, 0, 0] 2023-11-04 07:15:39 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:39 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.332, -29.59.25.499, I, 2.25982e+11Hz 2023-11-04 07:15:39 INFO imstat:::: -- Sum of pixel values [sum]: 93.3882 Jy/beam 2023-11-04 07:15:39 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2184.71 Jy/beam.Jy/beam 2023-11-04 07:15:39 INFO imstat:::: Statistics --- 2023-11-04 07:15:39 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00569996 Jy/beam 2023-11-04 07:15:39 INFO imstat:::: -- Variance of the pixel values : 0.13332 Jy/beam 2023-11-04 07:15:39 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.36513 Jy/beam 2023-11-04 07:15:39 INFO imstat:::: -- Root mean square [rms]: 0.365163 Jy/beam 2023-11-04 07:15:39 INFO imstat:::: -- Median of the pixel values [median]: -0.0940252 Jy/beam 2023-11-04 07:15:39 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.067585 Jy/beam 2023-11-04 07:15:39 INFO imstat:::: -- IQR [quartile]: 0.154968 Jy/beam 2023-11-04 07:15:39 INFO imstat:::: -- First quartile [q1]: -0.153589 Jy/beam 2023-11-04 07:15:39 INFO imstat:::: -- Third quartile [q3]: 0.00137894 Jy/beam 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:39 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 9.338820e+01 5.699964e-03 3.651633e-01 3.651300e-01 -1.858088e-01 2.768210e+00 2023-11-04 07:15:39 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:38.725166 End time: 2023-11-04 07:15:38.752951 2023-11-04 07:15:39 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:39 INFO imstat::::casa ########################################## .2023-11-04 07:15:39 INFO fixvis::::casa ########################################## 2023-11-04 07:15:39 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:39 INFO fixvis::::casa fixvis( vis='testx.ms', outputvis='testx.ms', field='1', refcode='J2000', reuse=True, phasecenter='-0.2333266794430661s 0deg', distances=[], datacolumn='all' ) 2023-11-04 07:15:39 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:39 INFO fixvis::::casa Will overwrite original MS ... 2023-11-04 07:15:39 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-04 07:15:39 INFO fixvis::::casa field: non-blanktest2_0 2023-11-04 07:15:39 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3331 -030.59.29.9987 2023-11-04 07:15:39 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.57062666] [-0.54090662] 2023-11-04 07:15:39 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:39 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2023-11-04 07:15:39 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-04 07:15:39 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0997 -030.59.29.9987 2023-11-04 07:15:39 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706436299021347 -0.5409066175416467 2023-11-04 07:15:39 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 1. 2023-11-04 07:15:39 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:39 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-04 07:15:39 INFO FixVis::setFields Selecting fields 1 2023-11-04 07:15:39 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:39 INFO FixVis::makeSelection() Processing field 1 2023-11-04 07:15:39 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:38.907315 End time: 2023-11-04 07:15:39.177551 2023-11-04 07:15:39 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:39 INFO fixvis::::casa ########################################## 2023-11-04 07:15:39 INFO tclean::::casa ########################################## 2023-11-04 07:15:39 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:39 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:39 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:39 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:39 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:39 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:39 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:39 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:39 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:39 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:39 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:39 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:39 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:40 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:40 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.021557 2023-11-04 07:15:40 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:40 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:40 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:40 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:40 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:40 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-04 07:15:40 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:40 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:39.228275 End time: 2023-11-04 07:15:40.319421 2023-11-04 07:15:40 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:40 INFO tclean::::casa ########################################## 2023-11-04 07:15:40 INFO imstat::::casa ########################################## 2023-11-04 07:15:40 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:40 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:40 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:40 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:40 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:40 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:40 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:40 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:40 INFO imstat:::: Values --- 2023-11-04 07:15:40 INFO imstat:::: -- flux density [flux]: 0.0112617 Jy 2023-11-04 07:15:40 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:40 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-04 07:15:40 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-04 07:15:40 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-04 07:15:40 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-04 07:15:40 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:40 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-04 07:15:40 INFO imstat:::: -- Sum of pixel values [sum]: 5.84977 Jy/beam 2023-11-04 07:15:40 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-04 07:15:40 INFO imstat:::: Statistics --- 2023-11-04 07:15:40 INFO imstat:::: -- Mean of the pixel values [mean]: 0.000357042 Jy/beam 2023-11-04 07:15:40 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-04 07:15:40 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-04 07:15:40 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-04 07:15:40 INFO imstat:::: -- Median of the pixel values [median]: -0.0988178 Jy/beam 2023-11-04 07:15:40 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649163 Jy/beam 2023-11-04 07:15:40 INFO imstat:::: -- IQR [quartile]: 0.146872 Jy/beam 2023-11-04 07:15:40 INFO imstat:::: -- First quartile [q1]: -0.155428 Jy/beam 2023-11-04 07:15:40 INFO imstat:::: -- Third quartile [q3]: -0.0085559 Jy/beam 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:40 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 5.849774e+00 3.570419e-04 3.646062e-01 3.646172e-01 -1.857701e-01 2.768178e+00 2023-11-04 07:15:40 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:40.323544 End time: 2023-11-04 07:15:40.355995 2023-11-04 07:15:40 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:40 INFO imstat::::casa ########################################## 2023-11-04 07:15:40 INFO tclean::::casa ########################################## 2023-11-04 07:15:40 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:40 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='1', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testz', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:40 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:40 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 1 | [Opened in readonly mode] 2023-11-04 07:15:40 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:40 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testz] : 2023-11-04 07:15:41 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:41 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56667e+07] 2023-11-04 07:15:41 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testz] with ftmachine : gridft 2023-11-04 07:15:41 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:41 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:41 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testz] : clark 2023-11-04 07:15:41 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:41 INFO SIImageStore::calcSensitivity [testz] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:41 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.01036 2023-11-04 07:15:41 INFO SIImageStore::printBeamSet Beam : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2023-11-04 07:15:41 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:41 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:41 INFO task_tclean::SDAlgorithmBase::restore [testz] : Restoring model image. 2023-11-04 07:15:41 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:41 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2023-11-04 07:15:41 INFO tclean::::casa Searching for images with prefix 'testz'... Found these, writing history into them: ['testz.model', 'testz.image', 'testz.pb', 'testz.sumwt', 'testz.residual', 'testz.psf'] 2023-11-04 07:15:41 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:40.408100 End time: 2023-11-04 07:15:41.460624 2023-11-04 07:15:41 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:41 INFO tclean::::casa ########################################## 2023-11-04 07:15:41 INFO imstat::::casa ########################################## 2023-11-04 07:15:41 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:41 INFO imstat::::casa imstat( imagename='testz.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:41 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:41 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:41 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:41 INFO imstat::image::statistics Determining stats for image testz.image 2023-11-04 07:15:41 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:41 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.597, -30.59.36.399, I, 2.25982e+11Hz to 18:00:01.610, -30.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:41 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:41 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:41 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:41 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:41 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.597, -30.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:41 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.610, -30.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:41 INFO imstat:::: Values --- 2023-11-04 07:15:41 INFO imstat:::: -- flux density [flux]: 0.187451 Jy 2023-11-04 07:15:41 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:41 INFO imstat:::: -- maximum value [max]: 2.76828 Jy/beam 2023-11-04 07:15:41 INFO imstat:::: -- minimum value [min]: -0.185749 Jy/beam 2023-11-04 07:15:41 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2023-11-04 07:15:41 INFO imstat:::: -- position of min value (pixel) [minpos]: [31, 109, 0, 0] 2023-11-04 07:15:41 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.333, -30.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:41 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.356, -30.59.25.499, I, 2.25982e+11Hz 2023-11-04 07:15:41 INFO imstat:::: -- Sum of pixel values [sum]: 97.5737 Jy/beam 2023-11-04 07:15:41 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2189.16 Jy/beam.Jy/beam 2023-11-04 07:15:41 INFO imstat:::: Statistics --- 2023-11-04 07:15:41 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00595543 Jy/beam 2023-11-04 07:15:41 INFO imstat:::: -- Variance of the pixel values : 0.133588 Jy/beam 2023-11-04 07:15:41 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365497 Jy/beam 2023-11-04 07:15:41 INFO imstat:::: -- Root mean square [rms]: 0.365535 Jy/beam 2023-11-04 07:15:41 INFO imstat:::: -- Median of the pixel values [median]: -0.0939644 Jy/beam 2023-11-04 07:15:41 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0676308 Jy/beam 2023-11-04 07:15:41 INFO imstat:::: -- IQR [quartile]: 0.155104 Jy/beam 2023-11-04 07:15:41 INFO imstat:::: -- First quartile [q1]: -0.153538 Jy/beam 2023-11-04 07:15:41 INFO imstat:::: -- Third quartile [q3]: 0.00156572 Jy/beam 2023-11-04 07:15:41 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:41 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:41 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:41 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:41 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:41 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:42 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 9.757369e+01 5.955426e-03 3.655346e-01 3.654972e-01 -1.857485e-01 2.768277e+00 2023-11-04 07:15:42 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:41.463736 End time: 2023-11-04 07:15:41.501900 2023-11-04 07:15:42 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:42 INFO imstat::::casa ########################################## .2023-11-04 07:15:42 INFO fixvis::::casa ########################################## 2023-11-04 07:15:42 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:42 INFO fixvis::::casa fixvis( vis='testx.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='-0.2333266794430661s 0deg', distances=[], datacolumn='DATA,CORRECTED' ) 2023-11-04 07:15:42 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:42 INFO fixvis::::casa Will overwrite original MS ... 2023-11-04 07:15:42 INFO fixvis::::casa Will only modify the visibilities in the columns DATA,CORRECTED 2023-11-04 07:15:42 INFO fixvis::::casa field: non-blanktest_0 2023-11-04 07:15:42 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:42 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-04 07:15:42 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:42 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2023-11-04 07:15:42 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-04 07:15:42 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0759 -029.59.29.9987 2023-11-04 07:15:42 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706453643230793 -0.5234533252689584 2023-11-04 07:15:42 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-04 07:15:42 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:42 INFO SubMS::parseColumnNames() Using DATA CORRECTED_DATA columns. 2023-11-04 07:15:42 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:42 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:42 INFO FixVis::makeSelection() Processing field 0 2023-11-04 07:15:42 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:41.664517 End time: 2023-11-04 07:15:41.909812 2023-11-04 07:15:42 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:42 INFO fixvis::::casa ########################################## 2023-11-04 07:15:42 INFO tclean::::casa ########################################## 2023-11-04 07:15:42 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:42 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:42 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:42 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:42 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:42 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:42 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:42 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:42 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:42 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:42 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:42 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:42 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:42 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:42 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010461 2023-11-04 07:15:42 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2023-11-04 07:15:42 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:43 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:43 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:43 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:43 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2023-11-04 07:15:43 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:43 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:41.945539 End time: 2023-11-04 07:15:42.981400 2023-11-04 07:15:43 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:43 INFO tclean::::casa ########################################## 2023-11-04 07:15:43 INFO imstat::::casa ########################################## 2023-11-04 07:15:43 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:43 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:43 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:43 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:43 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:43 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:43 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:43 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.569, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.591, -29.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.569, -29.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.591, -29.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:43 INFO imstat:::: Values --- 2023-11-04 07:15:43 INFO imstat:::: -- flux density [flux]: 0.186776 Jy 2023-11-04 07:15:43 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:43 INFO imstat:::: -- maximum value [max]: 2.76655 Jy/beam 2023-11-04 07:15:43 INFO imstat:::: -- minimum value [min]: -0.185788 Jy/beam 2023-11-04 07:15:43 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2023-11-04 07:15:43 INFO imstat:::: -- position of min value (pixel) [minpos]: [30, 110, 0, 0] 2023-11-04 07:15:43 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.307, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:43 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.338, -29.59.25.399, I, 2.25982e+11Hz 2023-11-04 07:15:43 INFO imstat:::: -- Sum of pixel values [sum]: 97.0186 Jy/beam 2023-11-04 07:15:43 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2184.63 Jy/beam.Jy/beam 2023-11-04 07:15:43 INFO imstat:::: Statistics --- 2023-11-04 07:15:43 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00592155 Jy/beam 2023-11-04 07:15:43 INFO imstat:::: -- Variance of the pixel values : 0.133312 Jy/beam 2023-11-04 07:15:43 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365119 Jy/beam 2023-11-04 07:15:43 INFO imstat:::: -- Root mean square [rms]: 0.365156 Jy/beam 2023-11-04 07:15:43 INFO imstat:::: -- Median of the pixel values [median]: -0.0935363 Jy/beam 2023-11-04 07:15:43 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0680586 Jy/beam 2023-11-04 07:15:43 INFO imstat:::: -- IQR [quartile]: 0.155433 Jy/beam 2023-11-04 07:15:43 INFO imstat:::: -- First quartile [q1]: -0.153311 Jy/beam 2023-11-04 07:15:43 INFO imstat:::: -- Third quartile [q3]: 0.00212144 Jy/beam 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:43 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 9.701861e+01 5.921546e-03 3.651561e-01 3.651192e-01 -1.857883e-01 2.766552e+00 2023-11-04 07:15:43 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:42.984950 End time: 2023-11-04 07:15:43.015323 2023-11-04 07:15:43 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:43 INFO imstat::::casa ########################################## .2023-11-04 07:15:43 INFO fixvis::::casa ########################################## 2023-11-04 07:15:43 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-04 07:15:43 INFO fixvis::::casa fixvis( vis='testx.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='-0.2333266794430661s 0deg', distances=[], datacolumn='DATA' ) 2023-11-04 07:15:43 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-04 07:15:43 INFO fixvis::::casa Will overwrite original MS ... 2023-11-04 07:15:43 INFO fixvis::::casa Will only modify the visibilities in the columns DATA 2023-11-04 07:15:43 INFO fixvis::::casa field: non-blanktest_0 2023-11-04 07:15:43 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-04 07:15:43 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-04 07:15:43 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-04 07:15:43 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2023-11-04 07:15:43 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-04 07:15:43 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0759 -029.59.29.9987 2023-11-04 07:15:43 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706453643230793 -0.5234533252689584 2023-11-04 07:15:43 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-04 07:15:43 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-04 07:15:43 INFO SubMS::parseColumnNames() Using DATA column. 2023-11-04 07:15:43 INFO FixVis::setFields Selecting fields 0 2023-11-04 07:15:43 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-04 07:15:43 INFO FixVis::makeSelection() Processing field 0 2023-11-04 07:15:43 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-04 07:15:43.165393 End time: 2023-11-04 07:15:43.460664 2023-11-04 07:15:43 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-04 07:15:43 INFO fixvis::::casa ########################################## 2023-11-04 07:15:44 INFO tclean::::casa ########################################## 2023-11-04 07:15:44 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-04 07:15:44 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, fullsummary=False, nmajor=-1, usemask='user', mask='user', pbmask=0.0, sidelobethreshold=3.0, noisethreshold=5.0, lownoisethreshold=1.5, negativethreshold=0.0, smoothfactor=1.0, minbeamfrac=0.3, cutthreshold=0.01, growiterations=75, dogrowprune=True, minpercentchange=-1.0, verbose=False, fastnoise=True, restart=True, savemodel='none', calcres=True, calcpsf=True, psfcutoff=0.35, parallel=False ) 2023-11-04 07:15:44 INFO tclean::::casa Verifying Input Parameters 2023-11-04 07:15:44 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-04 07:15:44 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-04 07:15:44 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-04 07:15:44 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-04 07:15:44 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-04 07:15:44 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-04 07:15:44 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.96 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-04 07:15:44 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-04 07:15:44 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-04 07:15:44 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-04 07:15:44 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-04 07:15:44 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.021985 2023-11-04 07:15:44 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2023-11-04 07:15:44 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-04 07:15:44 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-04 07:15:45 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-04 07:15:45 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2265) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-04 07:15:45 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2023-11-04 07:15:45 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.sumwt', 'testy.psf', 'testy.pb', 'testy.model', 'testy.image', 'testy.residual'] 2023-11-04 07:15:45 INFO tclean::::casa Task tclean complete. Start time: 2023-11-04 07:15:43.504528 End time: 2023-11-04 07:15:44.660149 2023-11-04 07:15:45 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-04 07:15:45 INFO tclean::::casa ########################################## 2023-11-04 07:15:45 INFO imstat::::casa ########################################## 2023-11-04 07:15:45 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-04 07:15:45 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2023-11-04 07:15:45 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-04 07:15:45 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-04 07:15:45 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-04 07:15:45 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-04 07:15:45 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-04 07:15:45 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.569, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.591, -29.59.23.699, I, 2.25982e+11Hz) 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.569, -29.59.36.399, I, 2.25982e+11Hz 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.591, -29.59.23.699, I, 2.25982e+11Hz 2023-11-04 07:15:45 INFO imstat:::: Values --- 2023-11-04 07:15:45 INFO imstat:::: -- flux density [flux]: 0.0113138 Jy 2023-11-04 07:15:45 INFO imstat:::: -- number of points [npts]: 16384 2023-11-04 07:15:45 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-04 07:15:45 INFO imstat:::: -- minimum value [min]: -0.185778 Jy/beam 2023-11-04 07:15:45 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-04 07:15:45 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-04 07:15:45 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.076, -29.59.29.999, I, 2.25982e+11Hz 2023-11-04 07:15:45 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.053, -29.59.34.499, I, 2.25982e+11Hz 2023-11-04 07:15:45 INFO imstat:::: -- Sum of pixel values [sum]: 5.87683 Jy/beam 2023-11-04 07:15:45 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.04 Jy/beam.Jy/beam 2023-11-04 07:15:45 INFO imstat:::: Statistics --- 2023-11-04 07:15:45 INFO imstat:::: -- Mean of the pixel values [mean]: 0.000358693 Jy/beam 2023-11-04 07:15:45 INFO imstat:::: -- Variance of the pixel values : 0.132945 Jy/beam 2023-11-04 07:15:45 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364616 Jy/beam 2023-11-04 07:15:45 INFO imstat:::: -- Root mean square [rms]: 0.364605 Jy/beam 2023-11-04 07:15:45 INFO imstat:::: -- Median of the pixel values [median]: -0.0988264 Jy/beam 2023-11-04 07:15:45 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.064913 Jy/beam 2023-11-04 07:15:45 INFO imstat:::: -- IQR [quartile]: 0.146844 Jy/beam 2023-11-04 07:15:45 INFO imstat:::: -- First quartile [q1]: -0.155414 Jy/beam 2023-11-04 07:15:45 INFO imstat:::: -- Third quartile [q3]: -0.00856965 Jy/beam 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-04 07:15:45 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 5.876831e+00 3.586933e-04 3.646055e-01 3.646164e-01 -1.857778e-01 2.768178e+00 2023-11-04 07:15:45 INFO imstat::::casa Task imstat complete. Start time: 2023-11-04 07:15:44.663384 End time: 2023-11-04 07:15:44.690946 2023-11-04 07:15:45 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-04 07:15:45 INFO imstat::::casa ########################################## . Using configuration file ~/.casa/config.py Using matplotlib backend: agg CASA 6.6.0.11a14200.dev7 -- Common Astronomy Software Applications [6.6.0.11a14200.dev7] Testing non-existing datacolumn - expected error *** Expected error *** ---------------------------------------------------------------------- Ran 15 tests in 27.090s OK