2023-11-08 07:11:57 INFO ::casa Using configuration file ~/.casa/config.py 2023-11-08 07:11:57 INFO ::casa 2023-11-08 07:11:57 INFO ::casa Checking Measures tables in data repository sub-directory /wheeldirectory/casa-CAS-14239-1-py3.8.el8/lib/py/lib/python3.8/site-packages/casadata/__data__/geodetic 2023-11-08 07:11:57 INFO ::casa IERSeop2000 (version date, last date in table (UTC)): 2023/10/06/15:00, 2023/09/06/00:00:00 2023-11-08 07:11:57 INFO ::casa IERSeop97 (version date, last date in table (UTC)): 2023/10/06/15:00, 2023/09/06/00:00:00 2023-11-08 07:11:57 INFO ::casa IERSpredict (version date, last date in table (UTC)): 2023/11/05/15:00, 2024/02/03/00:00:00 2023-11-08 07:11:57 INFO ::casa TAI_UTC (version date, last date in table (UTC)): 2023/10/13/15:00, 2017/01/01/00:00:00 2023-11-08 07:11:58 INFO fixvis::::casa ########################################## 2023-11-08 07:11:58 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:11:58 INFO fixvis::::casa fixvis( vis='0420+417.ms', outputvis='output.ms', field=[], refcode='J2000', reuse=True, phasecenter='', distances=[], datacolumn='all' ) 2023-11-08 07:11:58 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:11:58 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:11:58 INFO fixvis::::casa Will leave phase centers unchanged. 2023-11-08 07:11:58 INFO fixvis::::casa Recalculating the UVW coordinates ... 2023-11-08 07:11:58 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/output.ms 2023-11-08 07:11:58 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-08 07:11:58 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:11:58 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:11:57.898527 End time: 2023-11-08 07:11:58.087873 2023-11-08 07:11:58 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:11:58 INFO fixvis::::casa ########################################## .2023-11-08 07:11:58 INFO split::::casa ########################################## 2023-11-08 07:11:58 INFO split::::casa ##### Begin Task: split ##### 2023-11-08 07:11:58 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-08 07:11:58 INFO MSTransformManager::parseMsSpecParams Input file name is testx.ms 2023-11-08 07:11:58 INFO MSTransformManager::parseMsSpecParams Data column is CORRECTED 2023-11-08 07:11:58 INFO MSTransformManager::parseMsSpecParams Output file name is test9tmp.ms 2023-11-08 07:11:58 INFO MSTransformManager::colCheckInfo Adding DATA column to output MS from input CORRECTED_DATA column 2023-11-08 07:11:58 INFO MSTransformManager::open Select data 2023-11-08 07:11:58 INFO MSTransformManager::createOutputMSStructure Create output MS structure 2023-11-08 07:11:59 INFO split::::casa Task split complete. Start time: 2023-11-08 07:11:58.216032 End time: 2023-11-08 07:11:58.885787 2023-11-08 07:11:59 INFO split::::casa ##### End Task: split ##### 2023-11-08 07:11:59 INFO split::::casa ########################################## 2023-11-08 07:11:59 INFO fixvis::::casa ########################################## 2023-11-08 07:11:59 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:11:59 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-08 07:11:59 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:11:59 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:11:59 SEVERE fixvis::::casa Exception Reported: Error in fixvis: Datacolumn CORRECTED_DATA not present 2023-11-08 07:11:59 INFO fixvis::::casa Traceback (most recent call last): 2023-11-08 07:11:59 INFO fixvis::::casa+ File "/wheeldirectory/casa-CAS-14239-1-py3.8.el8/lib/py/lib/python3.8/site-packages/casashell/private/fixvis.py", line 472, in __call__ 2023-11-08 07:11:59 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-08 07:11:59 INFO fixvis::::casa+ File "/wheeldirectory/casa-CAS-14239-1-py3.8.el8/lib/py/lib/python3.8/site-packages/casatasks/private/task_fixvis.py", line 219, in fixvis 2023-11-08 07:11:59 INFO fixvis::::casa+ raise RuntimeError(msg) 2023-11-08 07:11:59 INFO fixvis::::casa+ RuntimeError: Datacolumn CORRECTED_DATA not present 2023-11-08 07:11:59 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:11:58.894529 End time: 2023-11-08 07:11:58.949195 2023-11-08 07:11:59 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:11:59 INFO fixvis::::casa ########################################## .2023-11-08 07:11:59 INFO split::::casa ########################################## 2023-11-08 07:11:59 INFO split::::casa ##### Begin Task: split ##### 2023-11-08 07:11:59 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-08 07:11:59 INFO MSTransformManager::parseMsSpecParams Input file name is testx.ms 2023-11-08 07:11:59 INFO MSTransformManager::parseMsSpecParams Data column is CORRECTED 2023-11-08 07:11:59 INFO MSTransformManager::parseMsSpecParams Output file name is test9tmp.ms 2023-11-08 07:11:59 INFO MSTransformManager::colCheckInfo Adding DATA column to output MS from input CORRECTED_DATA column 2023-11-08 07:11:59 INFO MSTransformManager::open Select data 2023-11-08 07:11:59 INFO MSTransformManager::createOutputMSStructure Create output MS structure 2023-11-08 07:12:00 INFO split::::casa Task split complete. Start time: 2023-11-08 07:11:59.088585 End time: 2023-11-08 07:11:59.764362 2023-11-08 07:12:00 INFO split::::casa ##### End Task: split ##### 2023-11-08 07:12:00 INFO split::::casa ########################################## 2023-11-08 07:12:00 INFO fixvis::::casa ########################################## 2023-11-08 07:12:00 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:00 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-08 07:12:00 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:00 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:12:00 INFO fixvis::::casa Will only modify the visibilities in the columns DATA 2023-11-08 07:12:00 INFO fixvis::::casa field: non-blanktest2_0 2023-11-08 07:12:00 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3331 -030.59.29.9987 2023-11-08 07:12:00 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.57062666] [-0.54090662] 2023-11-08 07:12:00 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:00 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2023-11-08 07:12:00 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-08 07:12:00 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0997 -030.59.29.9987 2023-11-08 07:12:00 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706436299021347 -0.5409066175416467 2023-11-08 07:12:00 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 1. 2023-11-08 07:12:00 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:00 INFO SubMS::parseColumnNames() Using DATA column. 2023-11-08 07:12:00 INFO FixVis::setFields Selecting fields 1 2023-11-08 07:12:00 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:00 INFO FixVis::makeSelection() Processing field 1 2023-11-08 07:12:00 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:11:59.772180 End time: 2023-11-08 07:11:59.947478 2023-11-08 07:12:00 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:00 INFO fixvis::::casa ########################################## 2023-11-08 07:12:00 INFO tclean::::casa ########################################## 2023-11-08 07:12:00 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:00 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, nmajor=-1, fullsummary=False, 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-08 07:12:00 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:00 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:00 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:00 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:00 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:00 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:00 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:00 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:00 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:00 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:00 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:01 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:01 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.015357 2023-11-08 07:12:01 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:01 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:01 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:01 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:01 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:01 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:01 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-08 07:12:01 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:00.054975 End time: 2023-11-08 07:12:01.150928 2023-11-08 07:12:01 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:01 INFO tclean::::casa ########################################## 2023-11-08 07:12:01 INFO imstat::::casa ########################################## 2023-11-08 07:12:01 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:01 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-08 07:12:01 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:01 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:01 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:01 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:01 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:01 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-08 07:12:01 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:01 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:01 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:01 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:01 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:01 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:01 INFO imstat:::: Values --- 2023-11-08 07:12:01 INFO imstat:::: -- flux density [flux]: 0.0112617 Jy 2023-11-08 07:12:01 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:01 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-08 07:12:01 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-08 07:12:01 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-08 07:12:01 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-08 07:12:01 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:01 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-08 07:12:01 INFO imstat:::: -- Sum of pixel values [sum]: 5.84977 Jy/beam 2023-11-08 07:12:01 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-08 07:12:01 INFO imstat:::: Statistics --- 2023-11-08 07:12:01 INFO imstat:::: -- Mean of the pixel values [mean]: 0.000357042 Jy/beam 2023-11-08 07:12:01 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-08 07:12:01 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-08 07:12:01 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-08 07:12:01 INFO imstat:::: -- Median of the pixel values [median]: -0.0988178 Jy/beam 2023-11-08 07:12:01 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649163 Jy/beam 2023-11-08 07:12:01 INFO imstat:::: -- IQR [quartile]: 0.146872 Jy/beam 2023-11-08 07:12:01 INFO imstat:::: -- First quartile [q1]: -0.155428 Jy/beam 2023-11-08 07:12:01 INFO imstat:::: -- Third quartile [q3]: -0.0085559 Jy/beam 2023-11-08 07:12:01 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:01 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:01 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:01 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:01 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:01 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:01 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-08 07:12:01 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:01.171248 End time: 2023-11-08 07:12:01.196678 2023-11-08 07:12:01 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:01 INFO imstat::::casa ########################################## 2023-11-08 07:12:01 INFO tclean::::casa ########################################## 2023-11-08 07:12:01 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:01 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, nmajor=-1, fullsummary=False, 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-08 07:12:01 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:01 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 1 | [Opened in readonly mode] 2023-11-08 07:12:01 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:01 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testz] : 2023-11-08 07:12:01 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:01 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56667e+07] 2023-11-08 07:12:01 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testz] with ftmachine : gridft 2023-11-08 07:12:01 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:01 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:01 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testz] : clark 2023-11-08 07:12:01 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:02 INFO SIImageStore::calcSensitivity [testz] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:02 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010838 2023-11-08 07:12:02 INFO SIImageStore::printBeamSet Beam : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2023-11-08 07:12:02 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:02 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:02 INFO task_tclean::SDAlgorithmBase::restore [testz] : Restoring model image. 2023-11-08 07:12:02 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:02 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2023-11-08 07:12:02 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-08 07:12:02 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:01.255789 End time: 2023-11-08 07:12:02.332007 2023-11-08 07:12:02 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:02 INFO tclean::::casa ########################################## 2023-11-08 07:12:02 INFO imstat::::casa ########################################## 2023-11-08 07:12:02 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:02 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-08 07:12:02 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:02 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:02 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:02 INFO imstat::image::statistics Determining stats for image testz.image 2023-11-08 07:12:02 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:02 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-08 07:12:02 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:02 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:02 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:02 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:02 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.597, -30.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:02 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.610, -30.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:02 INFO imstat:::: Values --- 2023-11-08 07:12:02 INFO imstat:::: -- flux density [flux]: 0.187451 Jy 2023-11-08 07:12:02 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:02 INFO imstat:::: -- maximum value [max]: 2.76828 Jy/beam 2023-11-08 07:12:02 INFO imstat:::: -- minimum value [min]: -0.185749 Jy/beam 2023-11-08 07:12:02 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2023-11-08 07:12:02 INFO imstat:::: -- position of min value (pixel) [minpos]: [31, 109, 0, 0] 2023-11-08 07:12:02 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.333, -30.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:02 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.356, -30.59.25.499, I, 2.25982e+11Hz 2023-11-08 07:12:02 INFO imstat:::: -- Sum of pixel values [sum]: 97.5737 Jy/beam 2023-11-08 07:12:02 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2189.16 Jy/beam.Jy/beam 2023-11-08 07:12:02 INFO imstat:::: Statistics --- 2023-11-08 07:12:02 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00595543 Jy/beam 2023-11-08 07:12:02 INFO imstat:::: -- Variance of the pixel values : 0.133588 Jy/beam 2023-11-08 07:12:02 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365497 Jy/beam 2023-11-08 07:12:02 INFO imstat:::: -- Root mean square [rms]: 0.365535 Jy/beam 2023-11-08 07:12:02 INFO imstat:::: -- Median of the pixel values [median]: -0.0939644 Jy/beam 2023-11-08 07:12:02 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0676308 Jy/beam 2023-11-08 07:12:02 INFO imstat:::: -- IQR [quartile]: 0.155104 Jy/beam 2023-11-08 07:12:02 INFO imstat:::: -- First quartile [q1]: -0.153538 Jy/beam 2023-11-08 07:12:02 INFO imstat:::: -- Third quartile [q3]: 0.00156572 Jy/beam 2023-11-08 07:12:02 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:02 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:02 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:02 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:02 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:02 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:02 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-08 07:12:02 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:02.335424 End time: 2023-11-08 07:12:02.360717 2023-11-08 07:12:02 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:02 INFO imstat::::casa ########################################## .2023-11-08 07:12:03 INFO split::::casa ########################################## 2023-11-08 07:12:03 INFO split::::casa ##### Begin Task: split ##### 2023-11-08 07:12:03 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-08 07:12:03 INFO MSTransformManager::parseMsSpecParams Input file name is testx.ms 2023-11-08 07:12:03 INFO MSTransformManager::parseMsSpecParams Data column is CORRECTED 2023-11-08 07:12:03 INFO MSTransformManager::parseMsSpecParams Output file name is test9tmp.ms 2023-11-08 07:12:03 INFO MSTransformManager::colCheckInfo Adding DATA column to output MS from input CORRECTED_DATA column 2023-11-08 07:12:03 INFO MSTransformManager::open Select data 2023-11-08 07:12:03 INFO MSTransformManager::createOutputMSStructure Create output MS structure 2023-11-08 07:12:03 INFO split::::casa Task split complete. Start time: 2023-11-08 07:12:02.504771 End time: 2023-11-08 07:12:03.160120 2023-11-08 07:12:03 INFO split::::casa ##### End Task: split ##### 2023-11-08 07:12:03 INFO split::::casa ########################################## 2023-11-08 07:12:03 INFO fixvis::::casa ########################################## 2023-11-08 07:12:03 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:03 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-08 07:12:03 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:03 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:12:03 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-08 07:12:03 INFO fixvis::::casa field: non-blanktest2_0 2023-11-08 07:12:03 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3331 -030.59.29.9987 2023-11-08 07:12:03 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.57062666] [-0.54090662] 2023-11-08 07:12:03 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:03 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2023-11-08 07:12:03 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-08 07:12:03 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0997 -030.59.29.9987 2023-11-08 07:12:03 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706436299021347 -0.5409066175416467 2023-11-08 07:12:03 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 1. 2023-11-08 07:12:03 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:03 INFO SubMS::parseColumnNames() Using DATA column. 2023-11-08 07:12:03 INFO FixVis::setFields Selecting fields 1 2023-11-08 07:12:03 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:03 INFO FixVis::makeSelection() Processing field 1 2023-11-08 07:12:03 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:03.168069 End time: 2023-11-08 07:12:03.366182 2023-11-08 07:12:03 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:03 INFO fixvis::::casa ########################################## 2023-11-08 07:12:03 INFO tclean::::casa ########################################## 2023-11-08 07:12:03 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:03 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, nmajor=-1, fullsummary=False, 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-08 07:12:03 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:03 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:03 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:03 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:04 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:04 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:04 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:04 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:04 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:04 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:04 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:04 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:04 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.020351 2023-11-08 07:12:04 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:04 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:04 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:04 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:04 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:04 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:04 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-08 07:12:04 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:03.406048 End time: 2023-11-08 07:12:04.362662 2023-11-08 07:12:04 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:04 INFO tclean::::casa ########################################## 2023-11-08 07:12:04 INFO imstat::::casa ########################################## 2023-11-08 07:12:04 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:04 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-08 07:12:04 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:04 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:04 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:04 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:04 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:04 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-08 07:12:04 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:04 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:04 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:04 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:04 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:04 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:04 INFO imstat:::: Values --- 2023-11-08 07:12:04 INFO imstat:::: -- flux density [flux]: 0.0112617 Jy 2023-11-08 07:12:04 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:04 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-08 07:12:04 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-08 07:12:04 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-08 07:12:04 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-08 07:12:04 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:04 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-08 07:12:04 INFO imstat:::: -- Sum of pixel values [sum]: 5.84977 Jy/beam 2023-11-08 07:12:04 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-08 07:12:04 INFO imstat:::: Statistics --- 2023-11-08 07:12:04 INFO imstat:::: -- Mean of the pixel values [mean]: 0.000357042 Jy/beam 2023-11-08 07:12:04 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-08 07:12:04 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-08 07:12:04 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-08 07:12:04 INFO imstat:::: -- Median of the pixel values [median]: -0.0988178 Jy/beam 2023-11-08 07:12:04 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649163 Jy/beam 2023-11-08 07:12:04 INFO imstat:::: -- IQR [quartile]: 0.146872 Jy/beam 2023-11-08 07:12:04 INFO imstat:::: -- First quartile [q1]: -0.155428 Jy/beam 2023-11-08 07:12:04 INFO imstat:::: -- Third quartile [q3]: -0.0085559 Jy/beam 2023-11-08 07:12:04 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:04 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:04 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:04 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:04 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:04 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:04 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-08 07:12:04 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:04.367913 End time: 2023-11-08 07:12:04.404785 2023-11-08 07:12:04 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:04 INFO imstat::::casa ########################################## 2023-11-08 07:12:04 INFO tclean::::casa ########################################## 2023-11-08 07:12:04 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:04 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, nmajor=-1, fullsummary=False, 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-08 07:12:04 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:04 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 1 | [Opened in readonly mode] 2023-11-08 07:12:04 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:05 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testz] : 2023-11-08 07:12:05 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:05 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56667e+07] 2023-11-08 07:12:05 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testz] with ftmachine : gridft 2023-11-08 07:12:05 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:05 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:05 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testz] : clark 2023-11-08 07:12:05 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:05 INFO SIImageStore::calcSensitivity [testz] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:05 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010003 2023-11-08 07:12:05 INFO SIImageStore::printBeamSet Beam : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2023-11-08 07:12:05 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:05 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:05 INFO task_tclean::SDAlgorithmBase::restore [testz] : Restoring model image. 2023-11-08 07:12:05 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:05 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2023-11-08 07:12:05 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-08 07:12:05 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:04.461295 End time: 2023-11-08 07:12:05.341937 2023-11-08 07:12:05 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:05 INFO tclean::::casa ########################################## 2023-11-08 07:12:05 INFO imstat::::casa ########################################## 2023-11-08 07:12:05 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:05 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-08 07:12:05 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:05 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:05 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:05 INFO imstat::image::statistics Determining stats for image testz.image 2023-11-08 07:12:05 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:05 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-08 07:12:05 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:05 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:05 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:05 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:05 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.597, -30.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:05 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.610, -30.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:05 INFO imstat:::: Values --- 2023-11-08 07:12:05 INFO imstat:::: -- flux density [flux]: 0.187451 Jy 2023-11-08 07:12:05 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:05 INFO imstat:::: -- maximum value [max]: 2.76828 Jy/beam 2023-11-08 07:12:05 INFO imstat:::: -- minimum value [min]: -0.185749 Jy/beam 2023-11-08 07:12:05 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2023-11-08 07:12:05 INFO imstat:::: -- position of min value (pixel) [minpos]: [31, 109, 0, 0] 2023-11-08 07:12:05 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.333, -30.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:05 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.356, -30.59.25.499, I, 2.25982e+11Hz 2023-11-08 07:12:05 INFO imstat:::: -- Sum of pixel values [sum]: 97.5737 Jy/beam 2023-11-08 07:12:05 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2189.16 Jy/beam.Jy/beam 2023-11-08 07:12:05 INFO imstat:::: Statistics --- 2023-11-08 07:12:05 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00595543 Jy/beam 2023-11-08 07:12:05 INFO imstat:::: -- Variance of the pixel values : 0.133588 Jy/beam 2023-11-08 07:12:05 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365497 Jy/beam 2023-11-08 07:12:05 INFO imstat:::: -- Root mean square [rms]: 0.365535 Jy/beam 2023-11-08 07:12:05 INFO imstat:::: -- Median of the pixel values [median]: -0.0939644 Jy/beam 2023-11-08 07:12:05 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0676308 Jy/beam 2023-11-08 07:12:05 INFO imstat:::: -- IQR [quartile]: 0.155104 Jy/beam 2023-11-08 07:12:05 INFO imstat:::: -- First quartile [q1]: -0.153538 Jy/beam 2023-11-08 07:12:05 INFO imstat:::: -- Third quartile [q3]: 0.00156572 Jy/beam 2023-11-08 07:12:05 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:05 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:05 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:05 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:05 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:05 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:05 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-08 07:12:05 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:05.344285 End time: 2023-11-08 07:12:05.360884 2023-11-08 07:12:05 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:05 INFO imstat::::casa ########################################## .2023-11-08 07:12:05 INFO fixvis::::casa ########################################## 2023-11-08 07:12:05 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:05 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-08 07:12:05 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:05 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:12:06 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-08 07:12:06 INFO fixvis::::casa field: non-blanktest_0 2023-11-08 07:12:06 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:06 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-08 07:12:06 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:06 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.49.29.9987 2023-11-08 07:12:06 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5205444431047309 2023-11-08 07:12:06 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-08 07:12:06 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/outms3.ms 2023-11-08 07:12:06 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-08 07:12:06 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:12:06 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:06 INFO FixVis::makeSelection() Processing field 0 2023-11-08 07:12:06 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:05.450819 End time: 2023-11-08 07:12:05.724903 2023-11-08 07:12:06 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:06 INFO fixvis::::casa ########################################## 2023-11-08 07:12:06 INFO fixvis::::casa ########################################## 2023-11-08 07:12:06 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:06 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-08 07:12:06 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:06 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:12:06 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-08 07:12:06 INFO fixvis::::casa field: non-blanktest_0 2023-11-08 07:12:06 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.49.29.9987 2023-11-08 07:12:06 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52054444] 2023-11-08 07:12:06 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:06 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:06 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2023-11-08 07:12:06 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-08 07:12:06 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:06 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-08 07:12:06 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:12:06 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:06 INFO FixVis::makeSelection() Processing field 0 2023-11-08 07:12:06 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:05.726038 End time: 2023-11-08 07:12:05.931904 2023-11-08 07:12:06 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:06 INFO fixvis::::casa ########################################## 2023-11-08 07:12:06 INFO tclean::::casa ########################################## 2023-11-08 07:12:06 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:06 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, nmajor=-1, fullsummary=False, 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-08 07:12:06 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:06 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:06 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:06 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:06 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:06 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:06 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:06 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:06 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:06 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:06 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:06 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:06 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.01 2023-11-08 07:12:06 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:06 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:06 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:07 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:07 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:07 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:07 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-08 07:12:07 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:05.967769 End time: 2023-11-08 07:12:06.884492 2023-11-08 07:12:07 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:07 INFO tclean::::casa ########################################## 2023-11-08 07:12:07 INFO imstat::::casa ########################################## 2023-11-08 07:12:07 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:07 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-08 07:12:07 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:07 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:07 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:07 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:07 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:07 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-08 07:12:07 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:07 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:07 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:07 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:07 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:07 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:07 INFO imstat:::: Values --- 2023-11-08 07:12:07 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2023-11-08 07:12:07 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:07 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-08 07:12:07 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-08 07:12:07 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-08 07:12:07 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-08 07:12:07 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:07 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-08 07:12:07 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2023-11-08 07:12:07 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-08 07:12:07 INFO imstat:::: Statistics --- 2023-11-08 07:12:07 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2023-11-08 07:12:07 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-08 07:12:07 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-08 07:12:07 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-08 07:12:07 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2023-11-08 07:12:07 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2023-11-08 07:12:07 INFO imstat:::: -- IQR [quartile]: 0.14686 Jy/beam 2023-11-08 07:12:07 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2023-11-08 07:12:07 INFO imstat:::: -- Third quartile [q3]: -0.00856261 Jy/beam 2023-11-08 07:12:07 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:07 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:07 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:07 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:07 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:07 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:07 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-08 07:12:07 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:06.888240 End time: 2023-11-08 07:12:06.916314 2023-11-08 07:12:07 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:07 INFO imstat::::casa ########################################## 2023-11-08 07:12:07 INFO fixvis::::casa ########################################## 2023-11-08 07:12:07 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:07 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-08 07:12:07 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:07 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:12:07 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-08 07:12:07 INFO fixvis::::casa field: non-blanktest_0 2023-11-08 07:12:07 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:07 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-08 07:12:07 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:07 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:07 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2023-11-08 07:12:07 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-08 07:12:07 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:07 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-08 07:12:07 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:12:07 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:07 INFO FixVis::makeSelection() Processing field 0 2023-11-08 07:12:07 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:06.933027 End time: 2023-11-08 07:12:07.170406 2023-11-08 07:12:07 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:07 INFO fixvis::::casa ########################################## 2023-11-08 07:12:07 INFO tclean::::casa ########################################## 2023-11-08 07:12:07 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:07 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, nmajor=-1, fullsummary=False, 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-08 07:12:07 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:07 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:07 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:07 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:07 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:07 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:07 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:07 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:07 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:07 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:07 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:08 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:08 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010184 2023-11-08 07:12:08 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:08 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:08 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:08 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:08 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:08 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:08 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-08 07:12:08 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:07.228712 End time: 2023-11-08 07:12:08.226145 2023-11-08 07:12:08 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:08 INFO tclean::::casa ########################################## 2023-11-08 07:12:08 INFO imstat::::casa ########################################## 2023-11-08 07:12:08 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:08 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-08 07:12:08 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:08 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:08 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:08 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:08 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:08 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-08 07:12:08 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:08 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:08 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:08 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:08 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:08 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:08 INFO imstat:::: Values --- 2023-11-08 07:12:08 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2023-11-08 07:12:08 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:08 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-08 07:12:08 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-08 07:12:08 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-08 07:12:08 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-08 07:12:08 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:08 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-08 07:12:08 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2023-11-08 07:12:08 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-08 07:12:08 INFO imstat:::: Statistics --- 2023-11-08 07:12:08 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2023-11-08 07:12:08 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-08 07:12:08 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-08 07:12:08 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-08 07:12:08 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2023-11-08 07:12:08 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2023-11-08 07:12:08 INFO imstat:::: -- IQR [quartile]: 0.146861 Jy/beam 2023-11-08 07:12:08 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2023-11-08 07:12:08 INFO imstat:::: -- Third quartile [q3]: -0.0085626 Jy/beam 2023-11-08 07:12:08 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:08 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:08 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:08 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:08 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:08 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:08 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-08 07:12:08 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:08.229551 End time: 2023-11-08 07:12:08.255185 2023-11-08 07:12:08 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:08 INFO imstat::::casa ########################################## .2023-11-08 07:12:08 INFO fixvis::::casa ########################################## 2023-11-08 07:12:08 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:08 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-08 07:12:08 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:08 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:12:08 INFO fixvis::::casa Will refocus to the given distances: ['100AU'] 2023-11-08 07:12:08 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-08 07:12:08 INFO fixvis::::casa field: non-blanktest_0 2023-11-08 07:12:08 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:08 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-08 07:12:08 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:08 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:08 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2023-11-08 07:12:08 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-08 07:12:08 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:08 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-08 07:12:08 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:12:08 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:08 INFO FixVis::makeSelection() Processing field 0 2023-11-08 07:12:09 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:08.372165 End time: 2023-11-08 07:12:08.665854 2023-11-08 07:12:09 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:09 INFO fixvis::::casa ########################################## 2023-11-08 07:12:09 INFO tclean::::casa ########################################## 2023-11-08 07:12:09 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:09 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, nmajor=-1, fullsummary=False, 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-08 07:12:09 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:09 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:09 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:09 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:09 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:09 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:09 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:09 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:09 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:09 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:09 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:09 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:09 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010319 2023-11-08 07:12:09 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:09 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:09 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:10 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:10 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:10 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:10 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-08 07:12:10 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:08.699356 End time: 2023-11-08 07:12:09.667797 2023-11-08 07:12:10 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:10 INFO tclean::::casa ########################################## 2023-11-08 07:12:10 INFO imstat::::casa ########################################## 2023-11-08 07:12:10 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:10 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-08 07:12:10 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:10 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:10 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:10 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:10 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:10 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-08 07:12:10 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:10 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:10 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:10 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:10 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:10 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:10 INFO imstat:::: Values --- 2023-11-08 07:12:10 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2023-11-08 07:12:10 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:10 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-08 07:12:10 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-08 07:12:10 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-08 07:12:10 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-08 07:12:10 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:10 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-08 07:12:10 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2023-11-08 07:12:10 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-08 07:12:10 INFO imstat:::: Statistics --- 2023-11-08 07:12:10 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2023-11-08 07:12:10 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-08 07:12:10 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-08 07:12:10 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-08 07:12:10 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2023-11-08 07:12:10 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2023-11-08 07:12:10 INFO imstat:::: -- IQR [quartile]: 0.146861 Jy/beam 2023-11-08 07:12:10 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2023-11-08 07:12:10 INFO imstat:::: -- Third quartile [q3]: -0.0085626 Jy/beam 2023-11-08 07:12:10 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:10 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:10 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:10 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:10 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:10 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:10 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-08 07:12:10 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:09.670151 End time: 2023-11-08 07:12:09.697237 2023-11-08 07:12:10 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:10 INFO imstat::::casa ########################################## .2023-11-08 07:12:10 INFO fixvis::::casa ########################################## 2023-11-08 07:12:10 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:10 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-08 07:12:10 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:10 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:12:10 INFO fixvis::::casa Will refocus to the given distances: 100AU 2023-11-08 07:12:10 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-08 07:12:10 INFO fixvis::::casa field: non-blanktest_0 2023-11-08 07:12:10 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:10 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-08 07:12:10 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:10 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:10 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2023-11-08 07:12:10 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-08 07:12:10 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:10 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-08 07:12:10 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:12:10 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:10 INFO FixVis::makeSelection() Processing field 0 2023-11-08 07:12:10 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:09.789946 End time: 2023-11-08 07:12:10.011867 2023-11-08 07:12:10 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:10 INFO fixvis::::casa ########################################## 2023-11-08 07:12:10 INFO tclean::::casa ########################################## 2023-11-08 07:12:10 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:10 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, nmajor=-1, fullsummary=False, 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-08 07:12:10 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:10 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:10 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:10 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:10 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:10 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:10 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:10 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:10 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:10 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:10 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:11 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:11 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010217 2023-11-08 07:12:11 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:11 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:11 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:11 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:11 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:11 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:11 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-08 07:12:11 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:10.053783 End time: 2023-11-08 07:12:11.054335 2023-11-08 07:12:11 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:11 INFO tclean::::casa ########################################## 2023-11-08 07:12:11 INFO imstat::::casa ########################################## 2023-11-08 07:12:11 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:11 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-08 07:12:11 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:11 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:11 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:11 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:11 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:11 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-08 07:12:11 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:11 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:11 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:11 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:11 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:11 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:11 INFO imstat:::: Values --- 2023-11-08 07:12:11 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2023-11-08 07:12:11 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:11 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-08 07:12:11 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-08 07:12:11 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-08 07:12:11 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-08 07:12:11 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:11 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-08 07:12:11 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2023-11-08 07:12:11 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-08 07:12:11 INFO imstat:::: Statistics --- 2023-11-08 07:12:11 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2023-11-08 07:12:11 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-08 07:12:11 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-08 07:12:11 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-08 07:12:11 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2023-11-08 07:12:11 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2023-11-08 07:12:11 INFO imstat:::: -- IQR [quartile]: 0.146861 Jy/beam 2023-11-08 07:12:11 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2023-11-08 07:12:11 INFO imstat:::: -- Third quartile [q3]: -0.0085626 Jy/beam 2023-11-08 07:12:11 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:11 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:11 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:11 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:11 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:11 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:11 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-08 07:12:11 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:11.056564 End time: 2023-11-08 07:12:11.073024 2023-11-08 07:12:11 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:11 INFO imstat::::casa ########################################## .2023-11-08 07:12:11 INFO fixvis::::casa ########################################## 2023-11-08 07:12:11 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:11 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-08 07:12:11 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:11 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:12:11 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-08 07:12:11 INFO fixvis::::casa field: non-blanktest_0 2023-11-08 07:12:11 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:11 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-08 07:12:11 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:11 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:11 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2023-11-08 07:12:11 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-08 07:12:11 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:11 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-08 07:12:11 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:12:11 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:11 INFO FixVis::makeSelection() Processing field 0 2023-11-08 07:12:11 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:11.146871 End time: 2023-11-08 07:12:11.362809 2023-11-08 07:12:11 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:11 INFO fixvis::::casa ########################################## 2023-11-08 07:12:11 INFO tclean::::casa ########################################## 2023-11-08 07:12:11 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:11 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, nmajor=-1, fullsummary=False, 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-08 07:12:11 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:11 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:11 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:11 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:12 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:12 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:12 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:12 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:12 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:12 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:12 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:12 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:12 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.009949 2023-11-08 07:12:12 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:12 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:12 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:12 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:12 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:12 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:12 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-08 07:12:12 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:11.403315 End time: 2023-11-08 07:12:12.339062 2023-11-08 07:12:12 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:12 INFO tclean::::casa ########################################## 2023-11-08 07:12:12 INFO imstat::::casa ########################################## 2023-11-08 07:12:12 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:12 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-08 07:12:12 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:12 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:12 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:12 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:12 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:12 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-08 07:12:12 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:12 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:12 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:12 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:12 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:12 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:12 INFO imstat:::: Values --- 2023-11-08 07:12:12 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2023-11-08 07:12:12 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:12 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-08 07:12:12 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-08 07:12:12 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-08 07:12:12 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-08 07:12:12 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:12 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-08 07:12:12 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2023-11-08 07:12:12 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-08 07:12:12 INFO imstat:::: Statistics --- 2023-11-08 07:12:12 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2023-11-08 07:12:12 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-08 07:12:12 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-08 07:12:12 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-08 07:12:12 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2023-11-08 07:12:12 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2023-11-08 07:12:12 INFO imstat:::: -- IQR [quartile]: 0.146861 Jy/beam 2023-11-08 07:12:12 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2023-11-08 07:12:12 INFO imstat:::: -- Third quartile [q3]: -0.0085626 Jy/beam 2023-11-08 07:12:12 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:12 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:12 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:12 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:12 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:12 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:12 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-08 07:12:12 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:12.342875 End time: 2023-11-08 07:12:12.370865 2023-11-08 07:12:12 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:12 INFO imstat::::casa ########################################## .2023-11-08 07:12:12 INFO fixvis::::casa ########################################## 2023-11-08 07:12:12 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:12 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-08 07:12:12 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:12 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:12:12 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-08 07:12:12 INFO fixvis::::casa field: non-blanktest_0 2023-11-08 07:12:12 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:12 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-08 07:12:12 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:13 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.26.9987 2023-11-08 07:12:13 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234387807809548 2023-11-08 07:12:13 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-08 07:12:13 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:13 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-08 07:12:13 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:12:13 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:13 INFO FixVis::makeSelection() Processing field 0 2023-11-08 07:12:13 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:12.466651 End time: 2023-11-08 07:12:12.698356 2023-11-08 07:12:13 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:13 INFO fixvis::::casa ########################################## 2023-11-08 07:12:13 INFO tclean::::casa ########################################## 2023-11-08 07:12:13 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:13 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, nmajor=-1, fullsummary=False, 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-08 07:12:13 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:13 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:13 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:13 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:13 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:13 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:13 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:13 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:13 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:13 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:13 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:13 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:13 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.013129 2023-11-08 07:12:13 INFO SIImageStore::printBeamSet Beam : 2.2794 arcsec, 2.01116 arcsec, 30.5337 deg 2023-11-08 07:12:13 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:13 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:14 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:14 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:14 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.2794 arcsec, 2.01116 arcsec, 30.5337 deg 2023-11-08 07:12:14 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-08 07:12:14 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:12.741070 End time: 2023-11-08 07:12:13.668007 2023-11-08 07:12:14 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:14 INFO tclean::::casa ########################################## 2023-11-08 07:12:14 INFO imstat::::casa ########################################## 2023-11-08 07:12:14 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:14 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-08 07:12:14 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:14 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:14 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:14 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:14 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:14 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-08 07:12:14 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:14 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:14 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:14 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:14 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.33.399, I, 2.25982e+11Hz 2023-11-08 07:12:14 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.20.699, I, 2.25982e+11Hz 2023-11-08 07:12:14 INFO imstat:::: Values --- 2023-11-08 07:12:14 INFO imstat:::: -- flux density [flux]: 0.249791 Jy 2023-11-08 07:12:14 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:14 INFO imstat:::: -- maximum value [max]: 2.76821 Jy/beam 2023-11-08 07:12:14 INFO imstat:::: -- minimum value [min]: -0.185766 Jy/beam 2023-11-08 07:12:14 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 34, 0, 0] 2023-11-08 07:12:14 INFO imstat:::: -- position of min value (pixel) [minpos]: [61, 79, 0, 0] 2023-11-08 07:12:14 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:14 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.332, -29.59.25.499, I, 2.25982e+11Hz 2023-11-08 07:12:14 INFO imstat:::: -- Sum of pixel values [sum]: 129.75 Jy/beam 2023-11-08 07:12:14 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2161.57 Jy/beam.Jy/beam 2023-11-08 07:12:14 INFO imstat:::: Statistics --- 2023-11-08 07:12:14 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00791933 Jy/beam 2023-11-08 07:12:14 INFO imstat:::: -- Variance of the pixel values : 0.131877 Jy/beam 2023-11-08 07:12:14 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.363149 Jy/beam 2023-11-08 07:12:14 INFO imstat:::: -- Root mean square [rms]: 0.363224 Jy/beam 2023-11-08 07:12:14 INFO imstat:::: -- Median of the pixel values [median]: -0.0875305 Jy/beam 2023-11-08 07:12:14 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0698687 Jy/beam 2023-11-08 07:12:14 INFO imstat:::: -- IQR [quartile]: 0.153301 Jy/beam 2023-11-08 07:12:14 INFO imstat:::: -- First quartile [q1]: -0.149306 Jy/beam 2023-11-08 07:12:14 INFO imstat:::: -- Third quartile [q3]: 0.00399534 Jy/beam 2023-11-08 07:12:14 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:14 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:14 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:14 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:14 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:14 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:14 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-08 07:12:14 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:13.670194 End time: 2023-11-08 07:12:13.686585 2023-11-08 07:12:14 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:14 INFO imstat::::casa ########################################## .2023-11-08 07:12:14 INFO fixvis::::casa ########################################## 2023-11-08 07:12:14 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:14 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='0h -0d0m3s', distances='', datacolumn='all' ) 2023-11-08 07:12:14 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:14 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:12:14 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-08 07:12:14 INFO fixvis::::casa field: non-blanktest_0 2023-11-08 07:12:14 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:14 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-08 07:12:14 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:14 INFO fixvis::::casa No equinox given in parameter 'phasecenter': 0h -0d0m3s 2023-11-08 07:12:14 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-08 07:12:14 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.32.9987 2023-11-08 07:12:14 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283963048133 -0.5234678696793916 2023-11-08 07:12:14 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-08 07:12:14 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:14 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-08 07:12:14 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:12:14 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:14 INFO FixVis::makeSelection() Processing field 0 2023-11-08 07:12:14 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:13.768478 End time: 2023-11-08 07:12:13.984590 2023-11-08 07:12:14 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:14 INFO fixvis::::casa ########################################## 2023-11-08 07:12:14 INFO tclean::::casa ########################################## 2023-11-08 07:12:14 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:14 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, nmajor=-1, fullsummary=False, 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-08 07:12:14 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:14 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:14 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:14 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:14 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:14 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:14 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:14 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:14 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:14 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:14 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:15 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:15 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.019855 2023-11-08 07:12:15 INFO SIImageStore::printBeamSet Beam : 2.27942 arcsec, 2.01116 arcsec, 30.5324 deg 2023-11-08 07:12:15 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:15 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:15 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:15 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:15 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27942 arcsec, 2.01116 arcsec, 30.5324 deg 2023-11-08 07:12:15 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-08 07:12:15 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:14.020923 End time: 2023-11-08 07:12:15.121881 2023-11-08 07:12:15 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:15 INFO tclean::::casa ########################################## 2023-11-08 07:12:15 INFO imstat::::casa ########################################## 2023-11-08 07:12:15 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:15 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-08 07:12:15 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:15 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:15 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:15 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:15 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:15 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-08 07:12:15 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:15 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:15 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:15 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:15 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.39.399, I, 2.25982e+11Hz 2023-11-08 07:12:15 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.26.699, I, 2.25982e+11Hz 2023-11-08 07:12:15 INFO imstat:::: Values --- 2023-11-08 07:12:15 INFO imstat:::: -- flux density [flux]: 0.274462 Jy 2023-11-08 07:12:15 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:15 INFO imstat:::: -- maximum value [max]: 2.76821 Jy/beam 2023-11-08 07:12:15 INFO imstat:::: -- minimum value [min]: -0.185768 Jy/beam 2023-11-08 07:12:15 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 94, 0, 0] 2023-11-08 07:12:15 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 49, 0, 0] 2023-11-08 07:12:15 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:15 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-08 07:12:15 INFO imstat:::: -- Sum of pixel values [sum]: 142.567 Jy/beam 2023-11-08 07:12:15 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2160.11 Jy/beam.Jy/beam 2023-11-08 07:12:15 INFO imstat:::: Statistics --- 2023-11-08 07:12:15 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00870157 Jy/beam 2023-11-08 07:12:15 INFO imstat:::: -- Variance of the pixel values : 0.131775 Jy/beam 2023-11-08 07:12:15 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.363008 Jy/beam 2023-11-08 07:12:15 INFO imstat:::: -- Root mean square [rms]: 0.363102 Jy/beam 2023-11-08 07:12:15 INFO imstat:::: -- Median of the pixel values [median]: -0.0857827 Jy/beam 2023-11-08 07:12:15 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0706394 Jy/beam 2023-11-08 07:12:15 INFO imstat:::: -- IQR [quartile]: 0.154155 Jy/beam 2023-11-08 07:12:15 INFO imstat:::: -- First quartile [q1]: -0.148851 Jy/beam 2023-11-08 07:12:15 INFO imstat:::: -- Third quartile [q3]: 0.00530347 Jy/beam 2023-11-08 07:12:15 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:15 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:15 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:15 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:15 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:15 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:15 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-08 07:12:15 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:15.124108 End time: 2023-11-08 07:12:15.141457 2023-11-08 07:12:15 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:15 INFO imstat::::casa ########################################## .2023-11-08 07:12:15 INFO fixvis::::casa ########################################## 2023-11-08 07:12:15 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:15 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-08 07:12:15 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:15 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:12:15 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-08 07:12:15 INFO fixvis::::casa field: non-blanktest_0 2023-11-08 07:12:15 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:15 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-08 07:12:15 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:15 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.5401 -029.59.29.9987 2023-11-08 07:12:15 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706116055101893 -0.5234533251913881 2023-11-08 07:12:15 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-08 07:12:15 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:15 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-08 07:12:15 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:12:15 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:15 INFO FixVis::makeSelection() Processing field 0 2023-11-08 07:12:15 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:15.225854 End time: 2023-11-08 07:12:15.454999 2023-11-08 07:12:15 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:15 INFO fixvis::::casa ########################################## 2023-11-08 07:12:15 INFO tclean::::casa ########################################## 2023-11-08 07:12:15 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:15 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, nmajor=-1, fullsummary=False, 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-08 07:12:15 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:15 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:16 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:16 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:16 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:16 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:16 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:16 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:16 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:16 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:16 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:16 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:16 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010186 2023-11-08 07:12:16 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01117 arcsec, 30.5334 deg 2023-11-08 07:12:16 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:16 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:16 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:16 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:16 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01117 arcsec, 30.5334 deg 2023-11-08 07:12:16 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-08 07:12:16 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:15.490047 End time: 2023-11-08 07:12:16.490423 2023-11-08 07:12:16 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:16 INFO tclean::::casa ########################################## 2023-11-08 07:12:16 INFO imstat::::casa ########################################## 2023-11-08 07:12:16 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:16 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-08 07:12:16 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:16 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:16 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:16 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:17 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:17 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-08 07:12:17 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:17 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:17 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:17 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:17 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:03.033, -29.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:17 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:02.055, -29.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:17 INFO imstat:::: Values --- 2023-11-08 07:12:17 INFO imstat:::: -- flux density [flux]: 0.216271 Jy 2023-11-08 07:12:17 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:17 INFO imstat:::: -- maximum value [max]: 2.76821 Jy/beam 2023-11-08 07:12:17 INFO imstat:::: -- minimum value [min]: -0.185787 Jy/beam 2023-11-08 07:12:17 INFO imstat:::: -- position of max value (pixel) [maxpos]: [94, 64, 0, 0] 2023-11-08 07:12:17 INFO imstat:::: -- position of min value (pixel) [minpos]: [97, 19, 0, 0] 2023-11-08 07:12:17 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:17 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-08 07:12:17 INFO imstat:::: -- Sum of pixel values [sum]: 112.34 Jy/beam 2023-11-08 07:12:17 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2183.54 Jy/beam.Jy/beam 2023-11-08 07:12:17 INFO imstat:::: Statistics --- 2023-11-08 07:12:17 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00685667 Jy/beam 2023-11-08 07:12:17 INFO imstat:::: -- Variance of the pixel values : 0.133234 Jy/beam 2023-11-08 07:12:17 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365012 Jy/beam 2023-11-08 07:12:17 INFO imstat:::: -- Root mean square [rms]: 0.365066 Jy/beam 2023-11-08 07:12:17 INFO imstat:::: -- Median of the pixel values [median]: -0.092606 Jy/beam 2023-11-08 07:12:17 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0686032 Jy/beam 2023-11-08 07:12:17 INFO imstat:::: -- IQR [quartile]: 0.156016 Jy/beam 2023-11-08 07:12:17 INFO imstat:::: -- First quartile [q1]: -0.152775 Jy/beam 2023-11-08 07:12:17 INFO imstat:::: -- Third quartile [q3]: 0.00324101 Jy/beam 2023-11-08 07:12:17 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:17 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:17 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:17 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:17 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:17 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:17 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-08 07:12:17 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:16.492630 End time: 2023-11-08 07:12:16.509053 2023-11-08 07:12:17 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:17 INFO imstat::::casa ########################################## .2023-11-08 07:12:17 INFO fixvis::::casa ########################################## 2023-11-08 07:12:17 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:17 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='-0d0m3s 0deg', distances='', datacolumn='all' ) 2023-11-08 07:12:17 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:17 INFO fixvis::::casa Copying original MS to outputvis ... 2023-11-08 07:12:17 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-08 07:12:17 INFO fixvis::::casa field: non-blanktest_0 2023-11-08 07:12:17 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:17 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-08 07:12:17 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:17 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0d0m3s 0deg 2023-11-08 07:12:17 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-08 07:12:17 INFO fixvis::::casa RA offset is an angle (not a time). Will divide by cos(DEC) to compute time offset. 2023-11-08 07:12:17 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0783 -029.59.29.9987 2023-11-08 07:12:17 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706451893538216 -0.5234533252689584 2023-11-08 07:12:17 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-08 07:12:17 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:17 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-08 07:12:17 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:12:17 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:17 INFO FixVis::makeSelection() Processing field 0 2023-11-08 07:12:17 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:16.597091 End time: 2023-11-08 07:12:16.816499 2023-11-08 07:12:17 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:17 INFO fixvis::::casa ########################################## 2023-11-08 07:12:17 INFO tclean::::casa ########################################## 2023-11-08 07:12:17 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:17 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, nmajor=-1, fullsummary=False, 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-08 07:12:17 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:17 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:17 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:17 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:17 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:17 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:17 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:17 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:17 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:17 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:17 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:17 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:17 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.01012 2023-11-08 07:12:17 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2023-11-08 07:12:17 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:17 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:18 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:18 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:18 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2023-11-08 07:12:18 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-08 07:12:18 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:16.855882 End time: 2023-11-08 07:12:17.895820 2023-11-08 07:12:18 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:18 INFO tclean::::casa ########################################## 2023-11-08 07:12:18 INFO imstat::::casa ########################################## 2023-11-08 07:12:18 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:18 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-08 07:12:18 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:18 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:18 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:18 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:18 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:18 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-08 07:12:18 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:18 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:18 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:18 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:18 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.571, -29.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:18 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.593, -29.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:18 INFO imstat:::: Values --- 2023-11-08 07:12:18 INFO imstat:::: -- flux density [flux]: 0.179787 Jy 2023-11-08 07:12:18 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:18 INFO imstat:::: -- maximum value [max]: 2.76821 Jy/beam 2023-11-08 07:12:18 INFO imstat:::: -- minimum value [min]: -0.185809 Jy/beam 2023-11-08 07:12:18 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2023-11-08 07:12:18 INFO imstat:::: -- position of min value (pixel) [minpos]: [31, 109, 0, 0] 2023-11-08 07:12:18 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:18 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.332, -29.59.25.499, I, 2.25982e+11Hz 2023-11-08 07:12:18 INFO imstat:::: -- Sum of pixel values [sum]: 93.3882 Jy/beam 2023-11-08 07:12:18 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2184.71 Jy/beam.Jy/beam 2023-11-08 07:12:18 INFO imstat:::: Statistics --- 2023-11-08 07:12:18 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00569996 Jy/beam 2023-11-08 07:12:18 INFO imstat:::: -- Variance of the pixel values : 0.13332 Jy/beam 2023-11-08 07:12:18 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.36513 Jy/beam 2023-11-08 07:12:18 INFO imstat:::: -- Root mean square [rms]: 0.365163 Jy/beam 2023-11-08 07:12:18 INFO imstat:::: -- Median of the pixel values [median]: -0.0940252 Jy/beam 2023-11-08 07:12:18 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.067585 Jy/beam 2023-11-08 07:12:18 INFO imstat:::: -- IQR [quartile]: 0.154968 Jy/beam 2023-11-08 07:12:18 INFO imstat:::: -- First quartile [q1]: -0.153589 Jy/beam 2023-11-08 07:12:18 INFO imstat:::: -- Third quartile [q3]: 0.00137894 Jy/beam 2023-11-08 07:12:18 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:18 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:18 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:18 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:18 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:18 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:18 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-08 07:12:18 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:17.899888 End time: 2023-11-08 07:12:17.930470 2023-11-08 07:12:18 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:18 INFO imstat::::casa ########################################## .2023-11-08 07:12:18 INFO fixvis::::casa ########################################## 2023-11-08 07:12:18 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:18 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-08 07:12:18 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:18 INFO fixvis::::casa Will overwrite original MS ... 2023-11-08 07:12:18 INFO fixvis::::casa Will modify the visibilities in all columns 2023-11-08 07:12:18 INFO fixvis::::casa field: non-blanktest2_0 2023-11-08 07:12:18 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3331 -030.59.29.9987 2023-11-08 07:12:18 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.57062666] [-0.54090662] 2023-11-08 07:12:18 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:18 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2023-11-08 07:12:18 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-08 07:12:18 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0997 -030.59.29.9987 2023-11-08 07:12:18 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706436299021347 -0.5409066175416467 2023-11-08 07:12:18 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 1. 2023-11-08 07:12:18 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:18 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2023-11-08 07:12:18 INFO FixVis::setFields Selecting fields 1 2023-11-08 07:12:18 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:18 INFO FixVis::makeSelection() Processing field 1 2023-11-08 07:12:18 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:18.088116 End time: 2023-11-08 07:12:18.413925 2023-11-08 07:12:18 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:18 INFO fixvis::::casa ########################################## 2023-11-08 07:12:18 INFO tclean::::casa ########################################## 2023-11-08 07:12:18 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:18 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, nmajor=-1, fullsummary=False, 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-08 07:12:18 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:18 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:18 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:19 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:19 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:19 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:19 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:19 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:19 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:19 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:19 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:19 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:19 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010179 2023-11-08 07:12:19 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:19 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:19 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:19 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:19 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:19 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2023-11-08 07:12:19 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-08 07:12:20 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:18.456614 End time: 2023-11-08 07:12:19.552456 2023-11-08 07:12:20 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:20 INFO tclean::::casa ########################################## 2023-11-08 07:12:20 INFO imstat::::casa ########################################## 2023-11-08 07:12:20 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:20 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-08 07:12:20 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:20 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:20 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:20 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:20 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:20 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-08 07:12:20 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:20 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:20 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:20 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:20 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:20 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:20 INFO imstat:::: Values --- 2023-11-08 07:12:20 INFO imstat:::: -- flux density [flux]: 0.0112617 Jy 2023-11-08 07:12:20 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:20 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-08 07:12:20 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2023-11-08 07:12:20 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-08 07:12:20 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-08 07:12:20 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:20 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2023-11-08 07:12:20 INFO imstat:::: -- Sum of pixel values [sum]: 5.84977 Jy/beam 2023-11-08 07:12:20 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2023-11-08 07:12:20 INFO imstat:::: Statistics --- 2023-11-08 07:12:20 INFO imstat:::: -- Mean of the pixel values [mean]: 0.000357042 Jy/beam 2023-11-08 07:12:20 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2023-11-08 07:12:20 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2023-11-08 07:12:20 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2023-11-08 07:12:20 INFO imstat:::: -- Median of the pixel values [median]: -0.0988178 Jy/beam 2023-11-08 07:12:20 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649163 Jy/beam 2023-11-08 07:12:20 INFO imstat:::: -- IQR [quartile]: 0.146872 Jy/beam 2023-11-08 07:12:20 INFO imstat:::: -- First quartile [q1]: -0.155428 Jy/beam 2023-11-08 07:12:20 INFO imstat:::: -- Third quartile [q3]: -0.0085559 Jy/beam 2023-11-08 07:12:20 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:20 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:20 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:20 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:20 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:20 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:20 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-08 07:12:20 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:19.557283 End time: 2023-11-08 07:12:19.592836 2023-11-08 07:12:20 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:20 INFO imstat::::casa ########################################## 2023-11-08 07:12:20 INFO tclean::::casa ########################################## 2023-11-08 07:12:20 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:20 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, nmajor=-1, fullsummary=False, 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-08 07:12:20 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:20 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 1 | [Opened in readonly mode] 2023-11-08 07:12:20 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:20 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testz] : 2023-11-08 07:12:20 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:20 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56667e+07] 2023-11-08 07:12:20 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testz] with ftmachine : gridft 2023-11-08 07:12:20 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:20 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:20 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testz] : clark 2023-11-08 07:12:20 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:20 INFO SIImageStore::calcSensitivity [testz] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:20 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010223 2023-11-08 07:12:20 INFO SIImageStore::printBeamSet Beam : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2023-11-08 07:12:20 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:20 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:20 INFO task_tclean::SDAlgorithmBase::restore [testz] : Restoring model image. 2023-11-08 07:12:20 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:21 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2023-11-08 07:12:21 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-08 07:12:21 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:19.642915 End time: 2023-11-08 07:12:20.610842 2023-11-08 07:12:21 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:21 INFO tclean::::casa ########################################## 2023-11-08 07:12:21 INFO imstat::::casa ########################################## 2023-11-08 07:12:21 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:21 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-08 07:12:21 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:21 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:21 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:21 INFO imstat::image::statistics Determining stats for image testz.image 2023-11-08 07:12:21 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:21 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-08 07:12:21 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:21 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:21 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:21 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:21 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.597, -30.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:21 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.610, -30.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:21 INFO imstat:::: Values --- 2023-11-08 07:12:21 INFO imstat:::: -- flux density [flux]: 0.187451 Jy 2023-11-08 07:12:21 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:21 INFO imstat:::: -- maximum value [max]: 2.76828 Jy/beam 2023-11-08 07:12:21 INFO imstat:::: -- minimum value [min]: -0.185749 Jy/beam 2023-11-08 07:12:21 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2023-11-08 07:12:21 INFO imstat:::: -- position of min value (pixel) [minpos]: [31, 109, 0, 0] 2023-11-08 07:12:21 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.333, -30.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:21 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.356, -30.59.25.499, I, 2.25982e+11Hz 2023-11-08 07:12:21 INFO imstat:::: -- Sum of pixel values [sum]: 97.5737 Jy/beam 2023-11-08 07:12:21 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2189.16 Jy/beam.Jy/beam 2023-11-08 07:12:21 INFO imstat:::: Statistics --- 2023-11-08 07:12:21 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00595543 Jy/beam 2023-11-08 07:12:21 INFO imstat:::: -- Variance of the pixel values : 0.133588 Jy/beam 2023-11-08 07:12:21 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365497 Jy/beam 2023-11-08 07:12:21 INFO imstat:::: -- Root mean square [rms]: 0.365535 Jy/beam 2023-11-08 07:12:21 INFO imstat:::: -- Median of the pixel values [median]: -0.0939644 Jy/beam 2023-11-08 07:12:21 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0676308 Jy/beam 2023-11-08 07:12:21 INFO imstat:::: -- IQR [quartile]: 0.155104 Jy/beam 2023-11-08 07:12:21 INFO imstat:::: -- First quartile [q1]: -0.153538 Jy/beam 2023-11-08 07:12:21 INFO imstat:::: -- Third quartile [q3]: 0.00156572 Jy/beam 2023-11-08 07:12:21 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:21 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:21 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:21 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:21 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:21 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:21 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-08 07:12:21 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:20.613204 End time: 2023-11-08 07:12:20.630480 2023-11-08 07:12:21 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:21 INFO imstat::::casa ########################################## .2023-11-08 07:12:21 INFO fixvis::::casa ########################################## 2023-11-08 07:12:21 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:21 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-08 07:12:21 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:21 INFO fixvis::::casa Will overwrite original MS ... 2023-11-08 07:12:21 INFO fixvis::::casa Will only modify the visibilities in the columns DATA,CORRECTED 2023-11-08 07:12:21 INFO fixvis::::casa field: non-blanktest_0 2023-11-08 07:12:21 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:21 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-08 07:12:21 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:21 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2023-11-08 07:12:21 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-08 07:12:21 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0759 -029.59.29.9987 2023-11-08 07:12:21 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706453643230793 -0.5234533252689584 2023-11-08 07:12:21 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-08 07:12:21 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:21 INFO SubMS::parseColumnNames() Using DATA CORRECTED_DATA columns. 2023-11-08 07:12:21 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:12:21 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:21 INFO FixVis::makeSelection() Processing field 0 2023-11-08 07:12:21 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:20.759042 End time: 2023-11-08 07:12:20.960470 2023-11-08 07:12:21 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:21 INFO fixvis::::casa ########################################## 2023-11-08 07:12:21 INFO tclean::::casa ########################################## 2023-11-08 07:12:21 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12: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, nmajor=-1, fullsummary=False, 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-08 07:12:21 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:21 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:21 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:21 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:21 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:21 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:21 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:21 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:21 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:21 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:21 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:21 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:21 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010364 2023-11-08 07:12:21 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2023-11-08 07:12:21 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:22 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:22 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:22 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:22 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2023-11-08 07:12: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-08 07:12:22 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:20.997817 End time: 2023-11-08 07:12:21.927409 2023-11-08 07:12:22 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:22 INFO tclean::::casa ########################################## 2023-11-08 07:12:22 INFO imstat::::casa ########################################## 2023-11-08 07:12:22 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12: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-08 07:12:22 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:22 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:22 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:22 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:22 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:22 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-08 07:12:22 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:22 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:22 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:22 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:22 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.569, -29.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:22 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.591, -29.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:22 INFO imstat:::: Values --- 2023-11-08 07:12:22 INFO imstat:::: -- flux density [flux]: 0.186776 Jy 2023-11-08 07:12:22 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:22 INFO imstat:::: -- maximum value [max]: 2.76655 Jy/beam 2023-11-08 07:12:22 INFO imstat:::: -- minimum value [min]: -0.185788 Jy/beam 2023-11-08 07:12:22 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2023-11-08 07:12:22 INFO imstat:::: -- position of min value (pixel) [minpos]: [30, 110, 0, 0] 2023-11-08 07:12:22 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.307, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:22 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.338, -29.59.25.399, I, 2.25982e+11Hz 2023-11-08 07:12:22 INFO imstat:::: -- Sum of pixel values [sum]: 97.0186 Jy/beam 2023-11-08 07:12:22 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2184.63 Jy/beam.Jy/beam 2023-11-08 07:12:22 INFO imstat:::: Statistics --- 2023-11-08 07:12:22 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00592155 Jy/beam 2023-11-08 07:12:22 INFO imstat:::: -- Variance of the pixel values : 0.133312 Jy/beam 2023-11-08 07:12:22 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365119 Jy/beam 2023-11-08 07:12:22 INFO imstat:::: -- Root mean square [rms]: 0.365156 Jy/beam 2023-11-08 07:12:22 INFO imstat:::: -- Median of the pixel values [median]: -0.0935363 Jy/beam 2023-11-08 07:12:22 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0680586 Jy/beam 2023-11-08 07:12:22 INFO imstat:::: -- IQR [quartile]: 0.155433 Jy/beam 2023-11-08 07:12:22 INFO imstat:::: -- First quartile [q1]: -0.153311 Jy/beam 2023-11-08 07:12:22 INFO imstat:::: -- Third quartile [q3]: 0.00212144 Jy/beam 2023-11-08 07:12:22 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:22 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:22 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:22 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:22 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:22 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:22 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-08 07:12:22 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:21.929596 End time: 2023-11-08 07:12:21.946108 2023-11-08 07:12:22 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:22 INFO imstat::::casa ########################################## .2023-11-08 07:12:22 INFO fixvis::::casa ########################################## 2023-11-08 07:12:22 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2023-11-08 07:12:22 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-08 07:12:22 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2023-11-08 07:12:22 INFO fixvis::::casa Will overwrite original MS ... 2023-11-08 07:12:22 INFO fixvis::::casa Will only modify the visibilities in the columns DATA 2023-11-08 07:12:22 INFO fixvis::::casa field: non-blanktest_0 2023-11-08 07:12:22 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2023-11-08 07:12:22 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2023-11-08 07:12:22 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2023-11-08 07:12:22 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2023-11-08 07:12:22 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2023-11-08 07:12:22 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0759 -029.59.29.9987 2023-11-08 07:12:22 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706453643230793 -0.5234533252689584 2023-11-08 07:12:22 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2023-11-08 07:12:22 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2023-11-08 07:12:22 INFO SubMS::parseColumnNames() Using DATA column. 2023-11-08 07:12:22 INFO FixVis::setFields Selecting fields 0 2023-11-08 07:12:22 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2023-11-08 07:12:22 INFO FixVis::makeSelection() Processing field 0 2023-11-08 07:12:22 INFO fixvis::::casa Task fixvis complete. Start time: 2023-11-08 07:12:22.076484 End time: 2023-11-08 07:12:22.294768 2023-11-08 07:12:22 INFO fixvis::::casa ##### End Task: fixvis ##### 2023-11-08 07:12:22 INFO fixvis::::casa ########################################## 2023-11-08 07:12:22 INFO tclean::::casa ########################################## 2023-11-08 07:12:22 INFO tclean::::casa ##### Begin Task: tclean ##### 2023-11-08 07:12:22 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, nmajor=-1, fullsummary=False, 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-08 07:12:22 INFO tclean::::casa Verifying Input Parameters 2023-11-08 07:12:22 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2023-11-08 07:12:22 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2023-11-08 07:12:22 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2023-11-08 07:12:22 INFO SynthesisImagerVi2::defineImage Impars : start 2023-11-08 07:12:22 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2023-11-08 07:12:22 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2023-11-08 07:12:22 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.94 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2023-11-08 07:12:22 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2023-11-08 07:12:22 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2023-11-08 07:12:22 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2023-11-08 07:12:23 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2023-11-08 07:12:23 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.0102 2023-11-08 07:12:23 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2023-11-08 07:12:23 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2023-11-08 07:12:23 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2023-11-08 07:12:23 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2023-11-08 07:12:23 WARN task_tclean::SIImageStore::restore (file /source/casa6/casatools/src/code/synthesis/ImagerObjects/SIImageStore.cc, line 2284) Restoring with an empty model image. Only residuals will be processed to form the output restored image. 2023-11-08 07:12:23 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2023-11-08 07:12:23 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-08 07:12:23 INFO tclean::::casa Task tclean complete. Start time: 2023-11-08 07:12:22.333585 End time: 2023-11-08 07:12:23.332925 2023-11-08 07:12:23 INFO tclean::::casa ##### End Task: tclean ##### 2023-11-08 07:12:23 INFO tclean::::casa ########################################## 2023-11-08 07:12:23 INFO imstat::::casa ########################################## 2023-11-08 07:12:23 INFO imstat::::casa ##### Begin Task: imstat ##### 2023-11-08 07:12:23 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-08 07:12:23 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2023-11-08 07:12:23 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2023-11-08 07:12:23 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2023-11-08 07:12:23 INFO imstat::image::statistics Determining stats for image testy.image 2023-11-08 07:12:23 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2023-11-08 07:12:23 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-08 07:12:23 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2023-11-08 07:12:23 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2023-11-08 07:12:23 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2023-11-08 07:12:23 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2023-11-08 07:12:23 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.569, -29.59.36.399, I, 2.25982e+11Hz 2023-11-08 07:12:23 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.591, -29.59.23.699, I, 2.25982e+11Hz 2023-11-08 07:12:23 INFO imstat:::: Values --- 2023-11-08 07:12:23 INFO imstat:::: -- flux density [flux]: 0.0113138 Jy 2023-11-08 07:12:23 INFO imstat:::: -- number of points [npts]: 16384 2023-11-08 07:12:23 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2023-11-08 07:12:23 INFO imstat:::: -- minimum value [min]: -0.185778 Jy/beam 2023-11-08 07:12:23 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2023-11-08 07:12:23 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2023-11-08 07:12:23 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.076, -29.59.29.999, I, 2.25982e+11Hz 2023-11-08 07:12:23 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.053, -29.59.34.499, I, 2.25982e+11Hz 2023-11-08 07:12:23 INFO imstat:::: -- Sum of pixel values [sum]: 5.87683 Jy/beam 2023-11-08 07:12:23 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.04 Jy/beam.Jy/beam 2023-11-08 07:12:23 INFO imstat:::: Statistics --- 2023-11-08 07:12:23 INFO imstat:::: -- Mean of the pixel values [mean]: 0.000358693 Jy/beam 2023-11-08 07:12:23 INFO imstat:::: -- Variance of the pixel values : 0.132945 Jy/beam 2023-11-08 07:12:23 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364616 Jy/beam 2023-11-08 07:12:23 INFO imstat:::: -- Root mean square [rms]: 0.364605 Jy/beam 2023-11-08 07:12:23 INFO imstat:::: -- Median of the pixel values [median]: -0.0988264 Jy/beam 2023-11-08 07:12:23 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.064913 Jy/beam 2023-11-08 07:12:23 INFO imstat:::: -- IQR [quartile]: 0.146844 Jy/beam 2023-11-08 07:12:23 INFO imstat:::: -- First quartile [q1]: -0.155414 Jy/beam 2023-11-08 07:12:23 INFO imstat:::: -- Third quartile [q3]: -0.00856965 Jy/beam 2023-11-08 07:12:23 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2023-11-08 07:12:23 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2023-11-08 07:12:23 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2023-11-08 07:12:23 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2023-11-08 07:12:23 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2023-11-08 07:12:23 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2023-11-08 07:12:23 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-08 07:12:23 INFO imstat::::casa Task imstat complete. Start time: 2023-11-08 07:12:23.335214 End time: 2023-11-08 07:12:23.351749 2023-11-08 07:12:23 INFO imstat::::casa ##### End Task: imstat ##### 2023-11-08 07:12:23 INFO imstat::::casa ########################################## . Using configuration file ~/.casa/config.py Using matplotlib backend: agg CASA 6.6.3.7a14239.dev1 -- Common Astronomy Software Applications [6.6.3.7a14239.dev1] Testing non-existing datacolumn - expected error *** Expected error *** ---------------------------------------------------------------------- Ran 15 tests in 26.616s OK