2024-05-02 08:58:17 INFO ::casa Using configuration file ~/.casa/config.py 2024-05-02 08:58:17 INFO ::casa 2024-05-02 08:58:17 INFO ::casa Checking Measures tables in data repository sub-directory /wheeldirectory/casa-CAS-14347-1-py3.8.el8/lib/py/lib/python3.8/site-packages/casadata/__data__/geodetic 2024-05-02 08:58:17 INFO ::casa IERSeop2000 (version date, last date in table (UTC)): 2024/04/13/15:00, 2024/03/14/00:00:00 2024-05-02 08:58:17 INFO ::casa IERSeop97 (version date, last date in table (UTC)): 2024/04/13/15:00, 2024/03/14/00:00:00 2024-05-02 08:58:17 INFO ::casa IERSpredict (version date, last date in table (UTC)): 2024/04/28/15:00, 2024/07/27/00:00:00 2024-05-02 08:58:17 INFO ::casa TAI_UTC (version date, last date in table (UTC)): 2024/04/22/15:00, 2017/01/01/00:00:00 2024-05-02 08:58:17 INFO fixvis::::casa ########################################## 2024-05-02 08:58:17 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:17 INFO fixvis::::casa fixvis( vis='0420+417.ms', outputvis='output.ms', field=[], refcode='J2000', reuse=True, phasecenter='', distances=[], datacolumn='all' ) 2024-05-02 08:58:17 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:17 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:17 INFO fixvis::::casa Will leave phase centers unchanged. 2024-05-02 08:58:17 INFO fixvis::::casa Recalculating the UVW coordinates ... 2024-05-02 08:58:18 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/output.ms 2024-05-02 08:58:18 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2024-05-02 08:58:18 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:18 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:17.442924 End time: 2024-05-02 08:58:17.625007 2024-05-02 08:58:18 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:18 INFO fixvis::::casa ########################################## 2024-05-02 08:58:18 INFO split::::casa ########################################## 2024-05-02 08:58:18 INFO split::::casa ##### Begin Task: split ##### 2024-05-02 08:58:18 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='' ) 2024-05-02 08:58:18 INFO MSTransformManager::parseMsSpecParams Input file name is testx.ms 2024-05-02 08:58:18 INFO MSTransformManager::parseMsSpecParams Data column is CORRECTED 2024-05-02 08:58:18 INFO MSTransformManager::parseMsSpecParams Output file name is test9tmp.ms 2024-05-02 08:58:18 INFO MSTransformManager::colCheckInfo Adding DATA column to output MS from input CORRECTED_DATA column 2024-05-02 08:58:18 INFO MSTransformManager::open Select data 2024-05-02 08:58:18 INFO MSTransformManager::createOutputMSStructure Create output MS structure 2024-05-02 08:58:18 INFO split::::casa Task split complete. Start time: 2024-05-02 08:58:17.748745 End time: 2024-05-02 08:58:18.377928 2024-05-02 08:58:18 INFO split::::casa ##### End Task: split ##### 2024-05-02 08:58:18 INFO split::::casa ########################################## 2024-05-02 08:58:18 INFO fixvis::::casa ########################################## 2024-05-02 08:58:18 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:18 INFO fixvis::::casa fixvis( vis='test9tmp.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='-0.2333266794430661s 0deg', distances=[], datacolumn='CORRECTED' ) 2024-05-02 08:58:18 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:18 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:18 SEVERE fixvis::::casa Exception Reported: Error in fixvis: Datacolumn CORRECTED_DATA not present 2024-05-02 08:58:18 INFO fixvis::::casa Traceback (most recent call last): 2024-05-02 08:58:18 INFO fixvis::::casa+ File "/wheeldirectory/casa-CAS-14347-1-py3.8.el8/lib/py/lib/python3.8/site-packages/casashell/private/fixvis.py", line 472, in __call__ 2024-05-02 08:58:18 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'] ) 2024-05-02 08:58:18 INFO fixvis::::casa+ File "/wheeldirectory/casa-CAS-14347-1-py3.8.el8/lib/py/lib/python3.8/site-packages/casatasks/private/task_fixvis.py", line 219, in fixvis 2024-05-02 08:58:18 INFO fixvis::::casa+ raise RuntimeError(msg) 2024-05-02 08:58:18 INFO fixvis::::casa+ RuntimeError: Datacolumn CORRECTED_DATA not present 2024-05-02 08:58:18 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:18.384280 End time: 2024-05-02 08:58:18.418030 2024-05-02 08:58:18 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:18 INFO fixvis::::casa ########################################## 2024-05-02 08:58:19 INFO split::::casa ########################################## 2024-05-02 08:58:19 INFO split::::casa ##### Begin Task: split ##### 2024-05-02 08:58:19 INFO split::::casa split( vis='testx.ms', outputvis='test9tmp.ms', keepmms=True, field='', spw='', scan='', antenna='', correlation='', timerange='', intent='', array='', uvrange='', observation='', feed='', datacolumn='corrected', keepflags=True, width=1, timebin='0s', combine='' ) 2024-05-02 08:58:19 INFO MSTransformManager::parseMsSpecParams Input file name is testx.ms 2024-05-02 08:58:19 INFO MSTransformManager::parseMsSpecParams Data column is CORRECTED 2024-05-02 08:58:19 INFO MSTransformManager::parseMsSpecParams Output file name is test9tmp.ms 2024-05-02 08:58:19 INFO MSTransformManager::colCheckInfo Adding DATA column to output MS from input CORRECTED_DATA column 2024-05-02 08:58:19 INFO MSTransformManager::open Select data 2024-05-02 08:58:19 INFO MSTransformManager::createOutputMSStructure Create output MS structure 2024-05-02 08:58:19 INFO split::::casa Task split complete. Start time: 2024-05-02 08:58:18.532453 End time: 2024-05-02 08:58:19.113186 2024-05-02 08:58:19 INFO split::::casa ##### End Task: split ##### 2024-05-02 08:58:19 INFO split::::casa ########################################## 2024-05-02 08:58:19 INFO fixvis::::casa ########################################## 2024-05-02 08:58:19 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:19 INFO fixvis::::casa fixvis( vis='test9tmp.ms', outputvis='testx.ms', field='1', refcode='J2000', reuse=True, phasecenter='-0.2333266794430661s 0deg', distances=[], datacolumn='DATA' ) 2024-05-02 08:58:19 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:19 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:19 INFO fixvis::::casa Will only modify the visibilities in the columns DATA 2024-05-02 08:58:19 INFO fixvis::::casa field: non-blanktest2_0 2024-05-02 08:58:19 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3331 -030.59.29.9987 2024-05-02 08:58:19 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.57062666] [-0.54090662] 2024-05-02 08:58:19 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:19 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2024-05-02 08:58:19 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2024-05-02 08:58:19 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0997 -030.59.29.9987 2024-05-02 08:58:19 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706436299021347 -0.5409066175416467 2024-05-02 08:58:19 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 1. 2024-05-02 08:58:19 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:19 INFO SubMS::parseColumnNames() Using DATA column. 2024-05-02 08:58:19 INFO FixVis::setFields Selecting fields 1 2024-05-02 08:58:19 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:19 INFO FixVis::makeSelection() Processing field 1 2024-05-02 08:58:19 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:19.122961 End time: 2024-05-02 08:58:19.368112 2024-05-02 08:58:19 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:19 INFO fixvis::::casa ########################################## 2024-05-02 08:58:19 INFO tclean::::casa ########################################## 2024-05-02 08:58:19 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:19 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 ) 2024-05-02 08:58:19 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:20 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:20 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:20 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:20 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:20 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:20 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:20 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:20 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:20 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:20 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:20 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:20 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.012305 2024-05-02 08:58:20 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:20 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:20 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:20 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58: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. 2024-05-02 08:58:20 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:20 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:20 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:19.491758 End time: 2024-05-02 08:58:20.370741 2024-05-02 08:58:20 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:20 INFO tclean::::casa ########################################## 2024-05-02 08:58:20 INFO imstat::::casa ########################################## 2024-05-02 08:58:20 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58: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 ) 2024-05-02 08:58:20 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:20 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:20 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:20 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:20 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58: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) 2024-05-02 08:58:20 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:20 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:20 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:20 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:20 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:20 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:20 INFO imstat:::: Values --- 2024-05-02 08:58:20 INFO imstat:::: -- flux density [flux]: 0.0112617 Jy 2024-05-02 08:58:20 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:20 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2024-05-02 08:58:20 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2024-05-02 08:58:20 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2024-05-02 08:58:20 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2024-05-02 08:58:20 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:20 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2024-05-02 08:58:20 INFO imstat:::: -- Sum of pixel values [sum]: 5.84977 Jy/beam 2024-05-02 08:58:20 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2024-05-02 08:58:20 INFO imstat:::: Statistics --- 2024-05-02 08:58:20 INFO imstat:::: -- Mean of the pixel values [mean]: 0.000357042 Jy/beam 2024-05-02 08:58:20 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2024-05-02 08:58:20 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2024-05-02 08:58:20 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2024-05-02 08:58:20 INFO imstat:::: -- Median of the pixel values [median]: -0.0988178 Jy/beam 2024-05-02 08:58:20 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649163 Jy/beam 2024-05-02 08:58:20 INFO imstat:::: -- IQR [quartile]: 0.146872 Jy/beam 2024-05-02 08:58:20 INFO imstat:::: -- First quartile [q1]: -0.155428 Jy/beam 2024-05-02 08:58:20 INFO imstat:::: -- Third quartile [q3]: -0.0085559 Jy/beam 2024-05-02 08:58:20 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:20 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:20 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:20 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:20 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:20 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58: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 2024-05-02 08:58:20 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:20.395216 End time: 2024-05-02 08:58:20.423562 2024-05-02 08:58:20 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:20 INFO imstat::::casa ########################################## 2024-05-02 08:58:20 INFO tclean::::casa ########################################## 2024-05-02 08:58:20 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58: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 ) 2024-05-02 08:58:20 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:20 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 1 | [Opened in readonly mode] 2024-05-02 08:58:20 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:21 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testz] : 2024-05-02 08:58:21 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:21 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56667e+07] 2024-05-02 08:58:21 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testz] with ftmachine : gridft 2024-05-02 08:58:21 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:21 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:21 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testz] : clark 2024-05-02 08:58:21 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:21 INFO SIImageStore::calcSensitivity [testz] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:21 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010901 2024-05-02 08:58:21 INFO SIImageStore::printBeamSet Beam : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2024-05-02 08:58:21 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:21 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:21 INFO task_tclean::SDAlgorithmBase::restore [testz] : Restoring model image. 2024-05-02 08:58:21 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. 2024-05-02 08:58:21 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2024-05-02 08:58:21 INFO tclean::::casa Searching for images with prefix 'testz'... Found these, writing history into them: ['testz.residual', 'testz.sumwt', 'testz.psf', 'testz.image', 'testz.model', 'testz.pb'] 2024-05-02 08:58:21 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:20.468280 End time: 2024-05-02 08:58:21.478284 2024-05-02 08:58:21 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:21 INFO tclean::::casa ########################################## 2024-05-02 08:58:21 INFO imstat::::casa ########################################## 2024-05-02 08:58:21 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58: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 ) 2024-05-02 08:58:21 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:21 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:21 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:21 INFO imstat::image::statistics Determining stats for image testz.image 2024-05-02 08:58:21 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58: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) 2024-05-02 08:58:22 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:22 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:22 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:22 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:22 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.597, -30.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:22 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.610, -30.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:22 INFO imstat:::: Values --- 2024-05-02 08:58:22 INFO imstat:::: -- flux density [flux]: 0.187451 Jy 2024-05-02 08:58:22 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:22 INFO imstat:::: -- maximum value [max]: 2.76828 Jy/beam 2024-05-02 08:58:22 INFO imstat:::: -- minimum value [min]: -0.185749 Jy/beam 2024-05-02 08:58:22 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2024-05-02 08:58:22 INFO imstat:::: -- position of min value (pixel) [minpos]: [31, 109, 0, 0] 2024-05-02 08:58:22 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.333, -30.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:22 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.356, -30.59.25.499, I, 2.25982e+11Hz 2024-05-02 08:58:22 INFO imstat:::: -- Sum of pixel values [sum]: 97.5737 Jy/beam 2024-05-02 08:58:22 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2189.16 Jy/beam.Jy/beam 2024-05-02 08:58:22 INFO imstat:::: Statistics --- 2024-05-02 08:58:22 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00595543 Jy/beam 2024-05-02 08:58:22 INFO imstat:::: -- Variance of the pixel values : 0.133588 Jy/beam 2024-05-02 08:58:22 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365497 Jy/beam 2024-05-02 08:58:22 INFO imstat:::: -- Root mean square [rms]: 0.365535 Jy/beam 2024-05-02 08:58:22 INFO imstat:::: -- Median of the pixel values [median]: -0.0939644 Jy/beam 2024-05-02 08:58:22 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0676308 Jy/beam 2024-05-02 08:58:22 INFO imstat:::: -- IQR [quartile]: 0.155104 Jy/beam 2024-05-02 08:58:22 INFO imstat:::: -- First quartile [q1]: -0.153538 Jy/beam 2024-05-02 08:58:22 INFO imstat:::: -- Third quartile [q3]: 0.00156572 Jy/beam 2024-05-02 08:58:22 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:22 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:22 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:22 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:22 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:22 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:22 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 2024-05-02 08:58:22 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:21.481953 End time: 2024-05-02 08:58:21.509656 2024-05-02 08:58:22 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:22 INFO imstat::::casa ########################################## 2024-05-02 08:58:22 INFO split::::casa ########################################## 2024-05-02 08:58:22 INFO split::::casa ##### Begin Task: split ##### 2024-05-02 08:58:22 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='' ) 2024-05-02 08:58:22 INFO MSTransformManager::parseMsSpecParams Input file name is testx.ms 2024-05-02 08:58:22 INFO MSTransformManager::parseMsSpecParams Data column is CORRECTED 2024-05-02 08:58:22 INFO MSTransformManager::parseMsSpecParams Output file name is test9tmp.ms 2024-05-02 08:58:22 INFO MSTransformManager::colCheckInfo Adding DATA column to output MS from input CORRECTED_DATA column 2024-05-02 08:58:22 INFO MSTransformManager::open Select data 2024-05-02 08:58:22 INFO MSTransformManager::createOutputMSStructure Create output MS structure 2024-05-02 08:58:22 INFO split::::casa Task split complete. Start time: 2024-05-02 08:58:21.673803 End time: 2024-05-02 08:58:22.358112 2024-05-02 08:58:22 INFO split::::casa ##### End Task: split ##### 2024-05-02 08:58:22 INFO split::::casa ########################################## 2024-05-02 08:58:22 INFO fixvis::::casa ########################################## 2024-05-02 08:58:22 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:22 INFO fixvis::::casa fixvis( vis='test9tmp.ms', outputvis='testx.ms', field='1', refcode='J2000', reuse=True, phasecenter='-0.2333266794430661s 0deg', distances=[], datacolumn='all' ) 2024-05-02 08:58:22 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:22 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:22 INFO fixvis::::casa Will modify the visibilities in all columns 2024-05-02 08:58:22 INFO fixvis::::casa field: non-blanktest2_0 2024-05-02 08:58:22 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3331 -030.59.29.9987 2024-05-02 08:58:22 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.57062666] [-0.54090662] 2024-05-02 08:58:22 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:22 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2024-05-02 08:58:22 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2024-05-02 08:58:22 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0997 -030.59.29.9987 2024-05-02 08:58:22 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706436299021347 -0.5409066175416467 2024-05-02 08:58:22 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 1. 2024-05-02 08:58:22 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:22 INFO SubMS::parseColumnNames() Using DATA column. 2024-05-02 08:58:22 INFO FixVis::setFields Selecting fields 1 2024-05-02 08:58:22 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:22 INFO FixVis::makeSelection() Processing field 1 2024-05-02 08:58:23 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:22.363565 End time: 2024-05-02 08:58:22.526426 2024-05-02 08:58:23 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:23 INFO fixvis::::casa ########################################## 2024-05-02 08:58:23 INFO tclean::::casa ########################################## 2024-05-02 08:58:23 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:23 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 ) 2024-05-02 08:58:23 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:23 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:23 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:23 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:23 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:23 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:23 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:23 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:23 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:23 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:23 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:23 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:23 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.02576 2024-05-02 08:58:23 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:23 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:23 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:23 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58: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. 2024-05-02 08:58:23 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:23 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:24 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:22.560839 End time: 2024-05-02 08:58:23.513505 2024-05-02 08:58:24 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:24 INFO tclean::::casa ########################################## 2024-05-02 08:58:24 INFO imstat::::casa ########################################## 2024-05-02 08:58:24 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:24 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 ) 2024-05-02 08:58:24 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:24 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:24 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:24 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:24 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:24 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) 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:24 INFO imstat:::: Values --- 2024-05-02 08:58:24 INFO imstat:::: -- flux density [flux]: 0.0112617 Jy 2024-05-02 08:58:24 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:24 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2024-05-02 08:58:24 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2024-05-02 08:58:24 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:24 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2024-05-02 08:58:24 INFO imstat:::: -- Sum of pixel values [sum]: 5.84977 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2024-05-02 08:58:24 INFO imstat:::: Statistics --- 2024-05-02 08:58:24 INFO imstat:::: -- Mean of the pixel values [mean]: 0.000357042 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Median of the pixel values [median]: -0.0988178 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649163 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- IQR [quartile]: 0.146872 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- First quartile [q1]: -0.155428 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Third quartile [q3]: -0.0085559 Jy/beam 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:24 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 2024-05-02 08:58:24 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:23.518126 End time: 2024-05-02 08:58:23.550073 2024-05-02 08:58:24 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:24 INFO imstat::::casa ########################################## 2024-05-02 08:58:24 INFO tclean::::casa ########################################## 2024-05-02 08:58:24 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:24 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 ) 2024-05-02 08:58:24 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:24 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 1 | [Opened in readonly mode] 2024-05-02 08:58:24 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:24 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testz] : 2024-05-02 08:58:24 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:24 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56667e+07] 2024-05-02 08:58:24 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testz] with ftmachine : gridft 2024-05-02 08:58:24 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:24 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:24 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testz] : clark 2024-05-02 08:58:24 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:24 INFO SIImageStore::calcSensitivity [testz] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:24 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010709 2024-05-02 08:58:24 INFO SIImageStore::printBeamSet Beam : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2024-05-02 08:58:24 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:24 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:24 INFO task_tclean::SDAlgorithmBase::restore [testz] : Restoring model image. 2024-05-02 08:58:24 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. 2024-05-02 08:58:24 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2024-05-02 08:58:24 INFO tclean::::casa Searching for images with prefix 'testz'... Found these, writing history into them: ['testz.residual', 'testz.sumwt', 'testz.psf', 'testz.image', 'testz.model', 'testz.pb'] 2024-05-02 08:58:24 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:23.598365 End time: 2024-05-02 08:58:24.427800 2024-05-02 08:58:24 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:24 INFO tclean::::casa ########################################## 2024-05-02 08:58:24 INFO imstat::::casa ########################################## 2024-05-02 08:58:24 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:24 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 ) 2024-05-02 08:58:24 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:24 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:24 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:24 INFO imstat::image::statistics Determining stats for image testz.image 2024-05-02 08:58:24 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:24 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) 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.597, -30.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.610, -30.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:24 INFO imstat:::: Values --- 2024-05-02 08:58:24 INFO imstat:::: -- flux density [flux]: 0.187451 Jy 2024-05-02 08:58:24 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:24 INFO imstat:::: -- maximum value [max]: 2.76828 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- minimum value [min]: -0.185749 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2024-05-02 08:58:24 INFO imstat:::: -- position of min value (pixel) [minpos]: [31, 109, 0, 0] 2024-05-02 08:58:24 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.333, -30.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:24 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.356, -30.59.25.499, I, 2.25982e+11Hz 2024-05-02 08:58:24 INFO imstat:::: -- Sum of pixel values [sum]: 97.5737 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2189.16 Jy/beam.Jy/beam 2024-05-02 08:58:24 INFO imstat:::: Statistics --- 2024-05-02 08:58:24 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00595543 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Variance of the pixel values : 0.133588 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365497 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Root mean square [rms]: 0.365535 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Median of the pixel values [median]: -0.0939644 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0676308 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- IQR [quartile]: 0.155104 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- First quartile [q1]: -0.153538 Jy/beam 2024-05-02 08:58:24 INFO imstat:::: -- Third quartile [q3]: 0.00156572 Jy/beam 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:24 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:24 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 2024-05-02 08:58:24 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:24.432035 End time: 2024-05-02 08:58:24.463443 2024-05-02 08:58:24 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:24 INFO imstat::::casa ########################################## 2024-05-02 08:58:25 INFO fixvis::::casa ########################################## 2024-05-02 08:58:25 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:25 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' ) 2024-05-02 08:58:25 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:25 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:25 INFO fixvis::::casa Will modify the visibilities in all columns 2024-05-02 08:58:25 INFO fixvis::::casa field: non-blanktest_0 2024-05-02 08:58:25 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:25 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2024-05-02 08:58:25 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:25 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.49.29.9987 2024-05-02 08:58:25 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5205444431047309 2024-05-02 08:58:25 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2024-05-02 08:58:25 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/outms3.ms 2024-05-02 08:58:25 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2024-05-02 08:58:25 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:25 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:25 INFO FixVis::makeSelection() Processing field 0 2024-05-02 08:58:25 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:24.580378 End time: 2024-05-02 08:58:24.812990 2024-05-02 08:58:25 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:25 INFO fixvis::::casa ########################################## 2024-05-02 08:58:25 INFO fixvis::::casa ########################################## 2024-05-02 08:58:25 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:25 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' ) 2024-05-02 08:58:25 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:25 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:25 INFO fixvis::::casa Will modify the visibilities in all columns 2024-05-02 08:58:25 INFO fixvis::::casa field: non-blanktest_0 2024-05-02 08:58:25 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.49.29.9987 2024-05-02 08:58:25 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52054444] 2024-05-02 08:58:25 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:25 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:25 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2024-05-02 08:58:25 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2024-05-02 08:58:25 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:25 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2024-05-02 08:58:25 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:25 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:25 INFO FixVis::makeSelection() Processing field 0 2024-05-02 08:58:25 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:24.814180 End time: 2024-05-02 08:58:25.034414 2024-05-02 08:58:25 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:25 INFO fixvis::::casa ########################################## 2024-05-02 08:58:25 INFO tclean::::casa ########################################## 2024-05-02 08:58:25 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:25 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 ) 2024-05-02 08:58:25 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:25 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:25 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:25 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:25 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:25 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:25 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:25 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:25 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:25 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:25 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:26 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:26 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.011391 2024-05-02 08:58:26 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:26 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:26 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:26 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58:26 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. 2024-05-02 08:58:26 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:26 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:26 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:25.081191 End time: 2024-05-02 08:58:26.142861 2024-05-02 08:58:26 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:26 INFO tclean::::casa ########################################## 2024-05-02 08:58:26 INFO imstat::::casa ########################################## 2024-05-02 08:58:26 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:26 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 ) 2024-05-02 08:58:26 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:26 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:26 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:26 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:26 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:26 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) 2024-05-02 08:58:26 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:26 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:26 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:26 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:26 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:26 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:26 INFO imstat:::: Values --- 2024-05-02 08:58:26 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2024-05-02 08:58:26 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:26 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2024-05-02 08:58:26 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2024-05-02 08:58:26 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2024-05-02 08:58:26 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2024-05-02 08:58:26 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:26 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2024-05-02 08:58:26 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2024-05-02 08:58:26 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2024-05-02 08:58:26 INFO imstat:::: Statistics --- 2024-05-02 08:58:26 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2024-05-02 08:58:26 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2024-05-02 08:58:26 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2024-05-02 08:58:26 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2024-05-02 08:58:26 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2024-05-02 08:58:26 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2024-05-02 08:58:26 INFO imstat:::: -- IQR [quartile]: 0.14686 Jy/beam 2024-05-02 08:58:26 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2024-05-02 08:58:26 INFO imstat:::: -- Third quartile [q3]: -0.00856261 Jy/beam 2024-05-02 08:58:26 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:26 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:26 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:26 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:26 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:26 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:26 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 2024-05-02 08:58:26 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:26.146603 End time: 2024-05-02 08:58:26.173515 2024-05-02 08:58:26 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:26 INFO imstat::::casa ########################################## 2024-05-02 08:58:26 INFO fixvis::::casa ########################################## 2024-05-02 08:58:26 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:26 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' ) 2024-05-02 08:58:26 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:26 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:26 INFO fixvis::::casa Will modify the visibilities in all columns 2024-05-02 08:58:26 INFO fixvis::::casa field: non-blanktest_0 2024-05-02 08:58:26 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:26 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2024-05-02 08:58:26 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:26 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:26 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2024-05-02 08:58:26 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2024-05-02 08:58:26 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:26 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2024-05-02 08:58:26 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:26 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:26 INFO FixVis::makeSelection() Processing field 0 2024-05-02 08:58:26 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:26.188838 End time: 2024-05-02 08:58:26.417214 2024-05-02 08:58:26 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:26 INFO fixvis::::casa ########################################## 2024-05-02 08:58:26 INFO tclean::::casa ########################################## 2024-05-02 08:58:26 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:26 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, 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 ) 2024-05-02 08:58:26 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:26 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:26 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:27 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:27 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:27 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:27 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:27 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:27 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:27 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:27 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:27 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:27 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010136 2024-05-02 08:58:27 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:27 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:27 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:28 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58:28 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. 2024-05-02 08:58:28 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:28 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:28 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:26.454340 End time: 2024-05-02 08:58:27.629324 2024-05-02 08:58:28 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:28 INFO tclean::::casa ########################################## 2024-05-02 08:58:28 INFO imstat::::casa ########################################## 2024-05-02 08:58:28 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:28 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2024-05-02 08:58:28 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:28 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:28 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:28 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:28 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:28 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) 2024-05-02 08:58:28 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:28 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:28 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:28 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:28 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:28 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:28 INFO imstat:::: Values --- 2024-05-02 08:58:28 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2024-05-02 08:58:28 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:28 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2024-05-02 08:58:28 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2024-05-02 08:58:28 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2024-05-02 08:58:28 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2024-05-02 08:58:28 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:28 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2024-05-02 08:58:28 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2024-05-02 08:58:28 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2024-05-02 08:58:28 INFO imstat:::: Statistics --- 2024-05-02 08:58:28 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2024-05-02 08:58:28 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2024-05-02 08:58:28 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2024-05-02 08:58:28 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2024-05-02 08:58:28 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2024-05-02 08:58:28 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2024-05-02 08:58:28 INFO imstat:::: -- IQR [quartile]: 0.146861 Jy/beam 2024-05-02 08:58:28 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2024-05-02 08:58:28 INFO imstat:::: -- Third quartile [q3]: -0.0085626 Jy/beam 2024-05-02 08:58:28 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:28 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:28 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:28 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:28 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:28 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:28 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 2024-05-02 08:58:28 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:27.631571 End time: 2024-05-02 08:58:27.648518 2024-05-02 08:58:28 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:28 INFO imstat::::casa ########################################## 2024-05-02 08:58:28 INFO fixvis::::casa ########################################## 2024-05-02 08:58:28 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:28 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' ) 2024-05-02 08:58:28 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:28 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:28 INFO fixvis::::casa Will refocus to the given distances: ['100AU'] 2024-05-02 08:58:28 INFO fixvis::::casa Will modify the visibilities in all columns 2024-05-02 08:58:28 INFO fixvis::::casa field: non-blanktest_0 2024-05-02 08:58:28 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:28 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2024-05-02 08:58:28 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:28 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:28 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2024-05-02 08:58:28 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2024-05-02 08:58:28 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:28 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2024-05-02 08:58:28 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:28 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:28 INFO FixVis::makeSelection() Processing field 0 2024-05-02 08:58:28 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:27.725099 End time: 2024-05-02 08:58:27.942715 2024-05-02 08:58:28 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:28 INFO fixvis::::casa ########################################## 2024-05-02 08:58:28 INFO tclean::::casa ########################################## 2024-05-02 08:58:28 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:28 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, 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 ) 2024-05-02 08:58:28 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:28 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:28 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:28 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:28 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:28 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:28 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:28 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:28 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:28 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:28 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:28 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:28 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010872 2024-05-02 08:58:28 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:28 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:29 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:29 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58:29 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. 2024-05-02 08:58:29 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:29 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:29 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:27.979694 End time: 2024-05-02 08:58:28.882069 2024-05-02 08:58:29 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:29 INFO tclean::::casa ########################################## 2024-05-02 08:58:29 INFO imstat::::casa ########################################## 2024-05-02 08:58:29 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:29 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 ) 2024-05-02 08:58:29 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:29 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:29 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:29 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:29 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:29 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz) 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:29 INFO imstat:::: Values --- 2024-05-02 08:58:29 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2024-05-02 08:58:29 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:29 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2024-05-02 08:58:29 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2024-05-02 08:58:29 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2024-05-02 08:58:29 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2024-05-02 08:58:29 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:29 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2024-05-02 08:58:29 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2024-05-02 08:58:29 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2024-05-02 08:58:29 INFO imstat:::: Statistics --- 2024-05-02 08:58:29 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2024-05-02 08:58:29 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2024-05-02 08:58:29 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2024-05-02 08:58:29 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2024-05-02 08:58:29 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2024-05-02 08:58:29 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2024-05-02 08:58:29 INFO imstat:::: -- IQR [quartile]: 0.146861 Jy/beam 2024-05-02 08:58:29 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2024-05-02 08:58:29 INFO imstat:::: -- Third quartile [q3]: -0.0085626 Jy/beam 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:29 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 5.849419e+00 3.570202e-04 3.646062e-01 3.646172e-01 -1.857702e-01 2.768178e+00 2024-05-02 08:58:29 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:28.884563 End time: 2024-05-02 08:58:28.902768 2024-05-02 08:58:29 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:29 INFO imstat::::casa ########################################## 2024-05-02 08:58:29 INFO fixvis::::casa ########################################## 2024-05-02 08:58:29 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:29 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='J2000 18h00m02.3092s -29d59m29.9987s', distances='100AU', datacolumn='all' ) 2024-05-02 08:58:29 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:29 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:29 INFO fixvis::::casa Will refocus to the given distances: 100AU 2024-05-02 08:58:29 INFO fixvis::::casa Will modify the visibilities in all columns 2024-05-02 08:58:29 INFO fixvis::::casa field: non-blanktest_0 2024-05-02 08:58:29 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:29 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2024-05-02 08:58:29 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:29 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:29 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2024-05-02 08:58:29 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2024-05-02 08:58:29 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:29 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2024-05-02 08:58:29 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:29 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:29 INFO FixVis::makeSelection() Processing field 0 2024-05-02 08:58:29 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:28.976278 End time: 2024-05-02 08:58:29.193134 2024-05-02 08:58:29 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:29 INFO fixvis::::casa ########################################## 2024-05-02 08:58:29 INFO tclean::::casa ########################################## 2024-05-02 08:58:29 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:29 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, 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 ) 2024-05-02 08:58:29 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:29 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:29 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:29 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:29 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:29 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:29 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:29 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:29 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:29 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:29 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:30 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:30 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010186 2024-05-02 08:58:30 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:30 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:30 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:30 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58:30 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. 2024-05-02 08:58:30 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:30 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:30 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:29.233208 End time: 2024-05-02 08:58:30.132154 2024-05-02 08:58:30 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:30 INFO tclean::::casa ########################################## 2024-05-02 08:58:30 INFO imstat::::casa ########################################## 2024-05-02 08:58:30 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:30 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2024-05-02 08:58:30 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:30 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:30 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:30 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:30 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:30 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz) 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:30 INFO imstat:::: Values --- 2024-05-02 08:58:30 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2024-05-02 08:58:30 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:30 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2024-05-02 08:58:30 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2024-05-02 08:58:30 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2024-05-02 08:58:30 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2024-05-02 08:58:30 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:30 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2024-05-02 08:58:30 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2024-05-02 08:58:30 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2024-05-02 08:58:30 INFO imstat:::: Statistics --- 2024-05-02 08:58:30 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2024-05-02 08:58:30 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2024-05-02 08:58:30 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2024-05-02 08:58:30 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2024-05-02 08:58:30 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2024-05-02 08:58:30 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2024-05-02 08:58:30 INFO imstat:::: -- IQR [quartile]: 0.146861 Jy/beam 2024-05-02 08:58:30 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2024-05-02 08:58:30 INFO imstat:::: -- Third quartile [q3]: -0.0085626 Jy/beam 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:30 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 5.849419e+00 3.570202e-04 3.646062e-01 3.646172e-01 -1.857702e-01 2.768178e+00 2024-05-02 08:58:30 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:30.134474 End time: 2024-05-02 08:58:30.151454 2024-05-02 08:58:30 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:30 INFO imstat::::casa ########################################## 2024-05-02 08:58:30 INFO fixvis::::casa ########################################## 2024-05-02 08:58:30 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:30 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' ) 2024-05-02 08:58:30 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:30 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:30 INFO fixvis::::casa Will modify the visibilities in all columns 2024-05-02 08:58:30 INFO fixvis::::casa field: non-blanktest_0 2024-05-02 08:58:30 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:30 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2024-05-02 08:58:30 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:30 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:30 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234533251913881 2024-05-02 08:58:30 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2024-05-02 08:58:30 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:30 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2024-05-02 08:58:30 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:30 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:30 INFO FixVis::makeSelection() Processing field 0 2024-05-02 08:58:30 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:30.229759 End time: 2024-05-02 08:58:30.441862 2024-05-02 08:58:30 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:30 INFO fixvis::::casa ########################################## 2024-05-02 08:58:30 INFO tclean::::casa ########################################## 2024-05-02 08:58:30 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:30 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, 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 ) 2024-05-02 08:58:30 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:30 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:30 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:31 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:31 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:31 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:31 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:31 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:31 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:31 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:31 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:31 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:31 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010184 2024-05-02 08:58:31 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:31 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:31 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:31 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58:31 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. 2024-05-02 08:58:31 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:31 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:31 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:30.480897 End time: 2024-05-02 08:58:31.418353 2024-05-02 08:58:31 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:31 INFO tclean::::casa ########################################## 2024-05-02 08:58:31 INFO imstat::::casa ########################################## 2024-05-02 08:58:31 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:31 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2024-05-02 08:58:31 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:31 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:31 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:31 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:31 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:31 INFO imstat::ImageRegion::fromRecord+ [0, 0, 0, 0] to [127, 127, 0, 0] (18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz to 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz) 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:31 INFO imstat:::: Values --- 2024-05-02 08:58:31 INFO imstat:::: -- flux density [flux]: 0.0112611 Jy 2024-05-02 08:58:31 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:31 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2024-05-02 08:58:31 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2024-05-02 08:58:31 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2024-05-02 08:58:31 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2024-05-02 08:58:31 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:31 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2024-05-02 08:58:31 INFO imstat:::: -- Sum of pixel values [sum]: 5.84942 Jy/beam 2024-05-02 08:58:31 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2024-05-02 08:58:31 INFO imstat:::: Statistics --- 2024-05-02 08:58:31 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00035702 Jy/beam 2024-05-02 08:58:31 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2024-05-02 08:58:31 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2024-05-02 08:58:31 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2024-05-02 08:58:31 INFO imstat:::: -- Median of the pixel values [median]: -0.0988091 Jy/beam 2024-05-02 08:58:31 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649207 Jy/beam 2024-05-02 08:58:31 INFO imstat:::: -- IQR [quartile]: 0.146861 Jy/beam 2024-05-02 08:58:31 INFO imstat:::: -- First quartile [q1]: -0.155423 Jy/beam 2024-05-02 08:58:31 INFO imstat:::: -- Third quartile [q3]: -0.0085626 Jy/beam 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:31 INFO imstat::ImageStatsCalculator::statistics 1.638400e+04 5.849419e+00 3.570202e-04 3.646062e-01 3.646172e-01 -1.857702e-01 2.768178e+00 2024-05-02 08:58:31 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:31.421774 End time: 2024-05-02 08:58:31.447286 2024-05-02 08:58:31 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:31 INFO imstat::::casa ########################################## 2024-05-02 08:58:32 INFO fixvis::::casa ########################################## 2024-05-02 08:58:32 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:32 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' ) 2024-05-02 08:58:32 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:32 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:32 INFO fixvis::::casa Will modify the visibilities in all columns 2024-05-02 08:58:32 INFO fixvis::::casa field: non-blanktest_0 2024-05-02 08:58:32 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:32 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2024-05-02 08:58:32 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:32 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.26.9987 2024-05-02 08:58:32 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283970320343 -0.5234387807809548 2024-05-02 08:58:32 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2024-05-02 08:58:32 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:32 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2024-05-02 08:58:32 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:32 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:32 INFO FixVis::makeSelection() Processing field 0 2024-05-02 08:58:32 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:31.548179 End time: 2024-05-02 08:58:31.777714 2024-05-02 08:58:32 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:32 INFO fixvis::::casa ########################################## 2024-05-02 08:58:32 INFO tclean::::casa ########################################## 2024-05-02 08:58:32 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:32 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, 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 ) 2024-05-02 08:58:32 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:32 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:32 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:32 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:32 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:32 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:32 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:32 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:32 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:32 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:32 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:32 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:32 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010315 2024-05-02 08:58:32 INFO SIImageStore::printBeamSet Beam : 2.2794 arcsec, 2.01116 arcsec, 30.5337 deg 2024-05-02 08:58:32 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:32 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:33 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58:33 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. 2024-05-02 08:58:33 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.2794 arcsec, 2.01116 arcsec, 30.5337 deg 2024-05-02 08:58:33 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:33 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:31.812653 End time: 2024-05-02 08:58:32.779949 2024-05-02 08:58:33 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:33 INFO tclean::::casa ########################################## 2024-05-02 08:58:33 INFO imstat::::casa ########################################## 2024-05-02 08:58:33 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:33 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2024-05-02 08:58:33 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:33 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:33 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:33 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:33 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:33 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) 2024-05-02 08:58:33 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:33 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:33 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:33 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:33 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.33.399, I, 2.25982e+11Hz 2024-05-02 08:58:33 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.20.699, I, 2.25982e+11Hz 2024-05-02 08:58:33 INFO imstat:::: Values --- 2024-05-02 08:58:33 INFO imstat:::: -- flux density [flux]: 0.249791 Jy 2024-05-02 08:58:33 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:33 INFO imstat:::: -- maximum value [max]: 2.76821 Jy/beam 2024-05-02 08:58:33 INFO imstat:::: -- minimum value [min]: -0.185766 Jy/beam 2024-05-02 08:58:33 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 34, 0, 0] 2024-05-02 08:58:33 INFO imstat:::: -- position of min value (pixel) [minpos]: [61, 79, 0, 0] 2024-05-02 08:58:33 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:33 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.332, -29.59.25.499, I, 2.25982e+11Hz 2024-05-02 08:58:33 INFO imstat:::: -- Sum of pixel values [sum]: 129.75 Jy/beam 2024-05-02 08:58:33 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2161.57 Jy/beam.Jy/beam 2024-05-02 08:58:33 INFO imstat:::: Statistics --- 2024-05-02 08:58:33 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00791933 Jy/beam 2024-05-02 08:58:33 INFO imstat:::: -- Variance of the pixel values : 0.131877 Jy/beam 2024-05-02 08:58:33 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.363149 Jy/beam 2024-05-02 08:58:33 INFO imstat:::: -- Root mean square [rms]: 0.363224 Jy/beam 2024-05-02 08:58:33 INFO imstat:::: -- Median of the pixel values [median]: -0.0875305 Jy/beam 2024-05-02 08:58:33 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0698687 Jy/beam 2024-05-02 08:58:33 INFO imstat:::: -- IQR [quartile]: 0.153301 Jy/beam 2024-05-02 08:58:33 INFO imstat:::: -- First quartile [q1]: -0.149306 Jy/beam 2024-05-02 08:58:33 INFO imstat:::: -- Third quartile [q3]: 0.00399534 Jy/beam 2024-05-02 08:58:33 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:33 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:33 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:33 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:33 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:33 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:33 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 2024-05-02 08:58:33 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:32.782331 End time: 2024-05-02 08:58:32.799711 2024-05-02 08:58:33 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:33 INFO imstat::::casa ########################################## 2024-05-02 08:58:33 INFO fixvis::::casa ########################################## 2024-05-02 08:58:33 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:33 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='0h -0d0m3s', distances='', datacolumn='all' ) 2024-05-02 08:58:33 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:33 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:33 INFO fixvis::::casa Will modify the visibilities in all columns 2024-05-02 08:58:33 INFO fixvis::::casa field: non-blanktest_0 2024-05-02 08:58:33 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:33 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2024-05-02 08:58:33 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:33 INFO fixvis::::casa No equinox given in parameter 'phasecenter': 0h -0d0m3s 2024-05-02 08:58:33 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2024-05-02 08:58:33 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.3092 -029.59.32.9987 2024-05-02 08:58:33 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706283963048133 -0.5234678696793916 2024-05-02 08:58:33 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2024-05-02 08:58:33 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:33 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2024-05-02 08:58:33 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:33 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:33 INFO FixVis::makeSelection() Processing field 0 2024-05-02 08:58:33 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:32.875117 End time: 2024-05-02 08:58:33.092937 2024-05-02 08:58:33 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:33 INFO fixvis::::casa ########################################## 2024-05-02 08:58:33 INFO tclean::::casa ########################################## 2024-05-02 08:58:33 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:33 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, 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 ) 2024-05-02 08:58:33 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:33 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:33 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:33 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:33 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:33 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:33 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:33 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:33 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:33 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:33 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:34 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:34 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010779 2024-05-02 08:58:34 INFO SIImageStore::printBeamSet Beam : 2.27942 arcsec, 2.01116 arcsec, 30.5324 deg 2024-05-02 08:58:34 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:34 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:34 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58:34 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. 2024-05-02 08:58:34 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27942 arcsec, 2.01116 arcsec, 30.5324 deg 2024-05-02 08:58:34 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:34 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:33.124752 End time: 2024-05-02 08:58:34.023772 2024-05-02 08:58:34 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:34 INFO tclean::::casa ########################################## 2024-05-02 08:58:34 INFO imstat::::casa ########################################## 2024-05-02 08:58:34 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:34 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2024-05-02 08:58:34 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:34 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:34 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:34 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:34 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:34 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) 2024-05-02 08:58:34 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:34 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:34 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:34 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:34 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.39.399, I, 2.25982e+11Hz 2024-05-02 08:58:34 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.26.699, I, 2.25982e+11Hz 2024-05-02 08:58:34 INFO imstat:::: Values --- 2024-05-02 08:58:34 INFO imstat:::: -- flux density [flux]: 0.274462 Jy 2024-05-02 08:58:34 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:34 INFO imstat:::: -- maximum value [max]: 2.76821 Jy/beam 2024-05-02 08:58:34 INFO imstat:::: -- minimum value [min]: -0.185768 Jy/beam 2024-05-02 08:58:34 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 94, 0, 0] 2024-05-02 08:58:34 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 49, 0, 0] 2024-05-02 08:58:34 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:34 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2024-05-02 08:58:34 INFO imstat:::: -- Sum of pixel values [sum]: 142.567 Jy/beam 2024-05-02 08:58:34 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2160.11 Jy/beam.Jy/beam 2024-05-02 08:58:34 INFO imstat:::: Statistics --- 2024-05-02 08:58:34 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00870157 Jy/beam 2024-05-02 08:58:34 INFO imstat:::: -- Variance of the pixel values : 0.131775 Jy/beam 2024-05-02 08:58:34 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.363008 Jy/beam 2024-05-02 08:58:34 INFO imstat:::: -- Root mean square [rms]: 0.363102 Jy/beam 2024-05-02 08:58:34 INFO imstat:::: -- Median of the pixel values [median]: -0.0857827 Jy/beam 2024-05-02 08:58:34 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0706394 Jy/beam 2024-05-02 08:58:34 INFO imstat:::: -- IQR [quartile]: 0.154155 Jy/beam 2024-05-02 08:58:34 INFO imstat:::: -- First quartile [q1]: -0.148851 Jy/beam 2024-05-02 08:58:34 INFO imstat:::: -- Third quartile [q3]: 0.00530347 Jy/beam 2024-05-02 08:58:34 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:34 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:34 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:34 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:34 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:34 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:34 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 2024-05-02 08:58:34 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:34.026099 End time: 2024-05-02 08:58:34.043582 2024-05-02 08:58:34 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:34 INFO imstat::::casa ########################################## 2024-05-02 08:58:34 INFO fixvis::::casa ########################################## 2024-05-02 08:58:34 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:34 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' ) 2024-05-02 08:58:34 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:34 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:34 INFO fixvis::::casa Will modify the visibilities in all columns 2024-05-02 08:58:34 INFO fixvis::::casa field: non-blanktest_0 2024-05-02 08:58:34 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:34 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2024-05-02 08:58:34 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:34 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.5401 -029.59.29.9987 2024-05-02 08:58:34 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706116055101893 -0.5234533251913881 2024-05-02 08:58:34 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2024-05-02 08:58:34 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:34 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2024-05-02 08:58:34 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:34 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:34 INFO FixVis::makeSelection() Processing field 0 2024-05-02 08:58:34 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:34.120970 End time: 2024-05-02 08:58:34.338887 2024-05-02 08:58:34 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:34 INFO fixvis::::casa ########################################## 2024-05-02 08:58:34 INFO tclean::::casa ########################################## 2024-05-02 08:58:34 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:34 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 ) 2024-05-02 08:58:34 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:34 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:34 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:34 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:34 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:34 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:34 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:34 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:34 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:34 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:34 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:35 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:35 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010585 2024-05-02 08:58:35 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01117 arcsec, 30.5334 deg 2024-05-02 08:58:35 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:35 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:35 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58:35 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. 2024-05-02 08:58:35 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01117 arcsec, 30.5334 deg 2024-05-02 08:58:35 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:35 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:34.371930 End time: 2024-05-02 08:58:35.310273 2024-05-02 08:58:35 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:35 INFO tclean::::casa ########################################## 2024-05-02 08:58:35 INFO imstat::::casa ########################################## 2024-05-02 08:58:35 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:35 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 ) 2024-05-02 08:58:35 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:35 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:35 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:35 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:35 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:35 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) 2024-05-02 08:58:35 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:35 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:35 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:35 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:35 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:03.033, -29.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:35 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:02.055, -29.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:35 INFO imstat:::: Values --- 2024-05-02 08:58:35 INFO imstat:::: -- flux density [flux]: 0.216271 Jy 2024-05-02 08:58:35 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:35 INFO imstat:::: -- maximum value [max]: 2.76821 Jy/beam 2024-05-02 08:58:35 INFO imstat:::: -- minimum value [min]: -0.185787 Jy/beam 2024-05-02 08:58:35 INFO imstat:::: -- position of max value (pixel) [maxpos]: [94, 64, 0, 0] 2024-05-02 08:58:35 INFO imstat:::: -- position of min value (pixel) [minpos]: [97, 19, 0, 0] 2024-05-02 08:58:35 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:35 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2024-05-02 08:58:35 INFO imstat:::: -- Sum of pixel values [sum]: 112.34 Jy/beam 2024-05-02 08:58:35 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2183.54 Jy/beam.Jy/beam 2024-05-02 08:58:35 INFO imstat:::: Statistics --- 2024-05-02 08:58:35 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00685667 Jy/beam 2024-05-02 08:58:35 INFO imstat:::: -- Variance of the pixel values : 0.133234 Jy/beam 2024-05-02 08:58:35 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365012 Jy/beam 2024-05-02 08:58:35 INFO imstat:::: -- Root mean square [rms]: 0.365066 Jy/beam 2024-05-02 08:58:35 INFO imstat:::: -- Median of the pixel values [median]: -0.092606 Jy/beam 2024-05-02 08:58:35 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0686032 Jy/beam 2024-05-02 08:58:35 INFO imstat:::: -- IQR [quartile]: 0.156016 Jy/beam 2024-05-02 08:58:35 INFO imstat:::: -- First quartile [q1]: -0.152775 Jy/beam 2024-05-02 08:58:35 INFO imstat:::: -- Third quartile [q3]: 0.00324101 Jy/beam 2024-05-02 08:58:35 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:35 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:35 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:35 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:35 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:35 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:35 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 2024-05-02 08:58:35 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:35.312632 End time: 2024-05-02 08:58:35.330485 2024-05-02 08:58:35 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:35 INFO imstat::::casa ########################################## 2024-05-02 08:58:35 INFO fixvis::::casa ########################################## 2024-05-02 08:58:35 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:35 INFO fixvis::::casa fixvis( vis='twocenteredpointsources.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='-0d0m3s 0deg', distances='', datacolumn='all' ) 2024-05-02 08:58:35 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:35 INFO fixvis::::casa Copying original MS to outputvis ... 2024-05-02 08:58:35 INFO fixvis::::casa Will modify the visibilities in all columns 2024-05-02 08:58:35 INFO fixvis::::casa field: non-blanktest_0 2024-05-02 08:58:35 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:35 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2024-05-02 08:58:35 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:35 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0d0m3s 0deg 2024-05-02 08:58:35 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2024-05-02 08:58:35 INFO fixvis::::casa RA offset is an angle (not a time). Will divide by cos(DEC) to compute time offset. 2024-05-02 08:58:35 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0783 -029.59.29.9987 2024-05-02 08:58:35 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706451893538216 -0.5234533252689584 2024-05-02 08:58:35 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2024-05-02 08:58:35 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:35 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2024-05-02 08:58:35 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:35 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:35 INFO FixVis::makeSelection() Processing field 0 2024-05-02 08:58:36 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:35.414664 End time: 2024-05-02 08:58:35.632750 2024-05-02 08:58:36 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:36 INFO fixvis::::casa ########################################## 2024-05-02 08:58:36 INFO tclean::::casa ########################################## 2024-05-02 08:58:36 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:36 INFO tclean::::casa tclean( vis='testx.ms', selectdata=True, field='0', spw='', timerange='', uvrange='', antenna='', scan='', observation='', intent='', datacolumn='corrected', imagename='testy', imsize=[128, 128], cell=['0.10000008arcsec', '0.10000008arcsec'], phasecenter='', stokes='I', projection='SIN', startmodel='', specmode='mfs', reffreq='', nchan=-1, start='', width='', outframe='LSRK', veltype='radio', restfreq=[], interpolation='linear', perchanweightdensity=True, gridder='standard', facets=1, psfphasecenter='', wprojplanes=1, vptable='', mosweight=True, aterm=True, psterm=False, wbawp=True, conjbeams=False, cfcache='', usepointing=False, computepastep=360.0, rotatepastep=360.0, pointingoffsetsigdev=[], pblimit=0.2, normtype='flatnoise', deconvolver='clark', scales=[], nterms=2, smallscalebias=0.0, fusedthreshold=0.0, largestscale=-1, restoration=True, restoringbeam=[], pbcor=False, outlierfile='', weighting='natural', robust=0.5, noise='1.0Jy', npixels=0, uvtaper=[], niter=0, gain=0.1, threshold='0.1mJy', nsigma=0.0, cycleniter=-1, cyclefactor=1.0, minpsffraction=0.05, maxpsffraction=0.8, interactive=False, 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 ) 2024-05-02 08:58:36 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:36 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:36 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:36 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:36 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:36 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:36 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:36 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:36 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:36 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:36 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:36 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:36 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.01023 2024-05-02 08:58:36 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2024-05-02 08:58:36 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:36 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:37 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58:37 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. 2024-05-02 08:58:37 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2024-05-02 08:58:37 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:37 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:35.666527 End time: 2024-05-02 08:58:36.605355 2024-05-02 08:58:37 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:37 INFO tclean::::casa ########################################## 2024-05-02 08:58:37 INFO imstat::::casa ########################################## 2024-05-02 08:58:37 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:37 INFO imstat::::casa imstat( imagename='testy.image', axes=[], region='', box='', chans='', stokes='', listit=True, verbose=True, mask='', stretch=False, logfile='', append=True, algorithm='classic', fence=-1.0, center='mean', lside=True, zscore=-1.0, maxiter=-1, clmethod='auto', niter=3 ) 2024-05-02 08:58:37 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:37 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:37 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:37 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:37 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:37 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) 2024-05-02 08:58:37 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:37 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:37 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:37 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:37 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.571, -29.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:37 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.593, -29.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:37 INFO imstat:::: Values --- 2024-05-02 08:58:37 INFO imstat:::: -- flux density [flux]: 0.179787 Jy 2024-05-02 08:58:37 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:37 INFO imstat:::: -- maximum value [max]: 2.76821 Jy/beam 2024-05-02 08:58:37 INFO imstat:::: -- minimum value [min]: -0.185809 Jy/beam 2024-05-02 08:58:37 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2024-05-02 08:58:37 INFO imstat:::: -- position of min value (pixel) [minpos]: [31, 109, 0, 0] 2024-05-02 08:58:37 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:37 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.332, -29.59.25.499, I, 2.25982e+11Hz 2024-05-02 08:58:37 INFO imstat:::: -- Sum of pixel values [sum]: 93.3882 Jy/beam 2024-05-02 08:58:37 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2184.71 Jy/beam.Jy/beam 2024-05-02 08:58:37 INFO imstat:::: Statistics --- 2024-05-02 08:58:37 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00569996 Jy/beam 2024-05-02 08:58:37 INFO imstat:::: -- Variance of the pixel values : 0.13332 Jy/beam 2024-05-02 08:58:37 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.36513 Jy/beam 2024-05-02 08:58:37 INFO imstat:::: -- Root mean square [rms]: 0.365163 Jy/beam 2024-05-02 08:58:37 INFO imstat:::: -- Median of the pixel values [median]: -0.0940252 Jy/beam 2024-05-02 08:58:37 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.067585 Jy/beam 2024-05-02 08:58:37 INFO imstat:::: -- IQR [quartile]: 0.154968 Jy/beam 2024-05-02 08:58:37 INFO imstat:::: -- First quartile [q1]: -0.153589 Jy/beam 2024-05-02 08:58:37 INFO imstat:::: -- Third quartile [q3]: 0.00137894 Jy/beam 2024-05-02 08:58:37 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:37 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:37 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:37 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:37 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:37 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:37 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 2024-05-02 08:58:37 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:36.607513 End time: 2024-05-02 08:58:36.623847 2024-05-02 08:58:37 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:37 INFO imstat::::casa ########################################## 2024-05-02 08:58:37 INFO fixvis::::casa ########################################## 2024-05-02 08:58:37 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:37 INFO fixvis::::casa fixvis( vis='testx.ms', outputvis='testx.ms', field='1', refcode='J2000', reuse=True, phasecenter='-0.2333266794430661s 0deg', distances=[], datacolumn='all' ) 2024-05-02 08:58:37 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:37 INFO fixvis::::casa Will overwrite original MS ... 2024-05-02 08:58:37 INFO fixvis::::casa Will modify the visibilities in all columns 2024-05-02 08:58:37 INFO fixvis::::casa field: non-blanktest2_0 2024-05-02 08:58:37 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3331 -030.59.29.9987 2024-05-02 08:58:37 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.57062666] [-0.54090662] 2024-05-02 08:58:37 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:37 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2024-05-02 08:58:37 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2024-05-02 08:58:37 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0997 -030.59.29.9987 2024-05-02 08:58:37 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706436299021347 -0.5409066175416467 2024-05-02 08:58:37 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 1. 2024-05-02 08:58:37 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:37 INFO SubMS::parseColumnNames() Using DATA MODEL_DATA CORRECTED_DATA columns. 2024-05-02 08:58:37 INFO FixVis::setFields Selecting fields 1 2024-05-02 08:58:37 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:37 INFO FixVis::makeSelection() Processing field 1 2024-05-02 08:58:37 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:36.741628 End time: 2024-05-02 08:58:37.018840 2024-05-02 08:58:37 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:37 INFO fixvis::::casa ########################################## 2024-05-02 08:58:37 INFO tclean::::casa ########################################## 2024-05-02 08:58:37 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:37 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 ) 2024-05-02 08:58:37 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:37 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:37 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:37 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:37 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:37 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:37 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:37 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:37 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:37 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:37 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:38 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:38 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010898 2024-05-02 08:58:38 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:38 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:38 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:38 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58:38 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. 2024-05-02 08:58:38 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01116 arcsec, 30.5332 deg 2024-05-02 08:58:38 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:38 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:37.057068 End time: 2024-05-02 08:58:38.059851 2024-05-02 08:58:38 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:38 INFO tclean::::casa ########################################## 2024-05-02 08:58:38 INFO imstat::::casa ########################################## 2024-05-02 08:58:38 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:38 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 ) 2024-05-02 08:58:38 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:38 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:38 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:38 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:38 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:38 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) 2024-05-02 08:58:38 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:38 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:38 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:38 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:38 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.802, -29.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:38 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.824, -29.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:38 INFO imstat:::: Values --- 2024-05-02 08:58:38 INFO imstat:::: -- flux density [flux]: 0.0112617 Jy 2024-05-02 08:58:38 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:38 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2024-05-02 08:58:38 INFO imstat:::: -- minimum value [min]: -0.18577 Jy/beam 2024-05-02 08:58:38 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2024-05-02 08:58:38 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2024-05-02 08:58:38 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.309, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:38 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.286, -29.59.34.499, I, 2.25982e+11Hz 2024-05-02 08:58:38 INFO imstat:::: -- Sum of pixel values [sum]: 5.84977 Jy/beam 2024-05-02 08:58:38 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.05 Jy/beam.Jy/beam 2024-05-02 08:58:38 INFO imstat:::: Statistics --- 2024-05-02 08:58:38 INFO imstat:::: -- Mean of the pixel values [mean]: 0.000357042 Jy/beam 2024-05-02 08:58:38 INFO imstat:::: -- Variance of the pixel values : 0.132946 Jy/beam 2024-05-02 08:58:38 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364617 Jy/beam 2024-05-02 08:58:38 INFO imstat:::: -- Root mean square [rms]: 0.364606 Jy/beam 2024-05-02 08:58:38 INFO imstat:::: -- Median of the pixel values [median]: -0.0988178 Jy/beam 2024-05-02 08:58:38 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0649163 Jy/beam 2024-05-02 08:58:38 INFO imstat:::: -- IQR [quartile]: 0.146872 Jy/beam 2024-05-02 08:58:38 INFO imstat:::: -- First quartile [q1]: -0.155428 Jy/beam 2024-05-02 08:58:38 INFO imstat:::: -- Third quartile [q3]: -0.0085559 Jy/beam 2024-05-02 08:58:38 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:38 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:38 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:38 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:38 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:38 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:38 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 2024-05-02 08:58:38 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:38.063474 End time: 2024-05-02 08:58:38.090347 2024-05-02 08:58:38 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:38 INFO imstat::::casa ########################################## 2024-05-02 08:58:38 INFO tclean::::casa ########################################## 2024-05-02 08:58:38 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:38 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 ) 2024-05-02 08:58:38 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:38 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 1 | [Opened in readonly mode] 2024-05-02 08:58:38 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:38 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testz] : 2024-05-02 08:58:38 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:38 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56667e+07] 2024-05-02 08:58:38 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testz] with ftmachine : gridft 2024-05-02 08:58:38 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:38 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:38 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testz] : clark 2024-05-02 08:58:38 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:39 INFO SIImageStore::calcSensitivity [testz] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:39 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.010463 2024-05-02 08:58:39 INFO SIImageStore::printBeamSet Beam : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2024-05-02 08:58:39 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:39 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:39 INFO task_tclean::SDAlgorithmBase::restore [testz] : Restoring model image. 2024-05-02 08:58:39 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. 2024-05-02 08:58:39 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.28347 arcsec, 2.0118 arcsec, 30.205 deg 2024-05-02 08:58:39 INFO tclean::::casa Searching for images with prefix 'testz'... Found these, writing history into them: ['testz.residual', 'testz.sumwt', 'testz.psf', 'testz.image', 'testz.model', 'testz.pb'] 2024-05-02 08:58:39 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:38.132975 End time: 2024-05-02 08:58:39.174258 2024-05-02 08:58:39 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:39 INFO tclean::::casa ########################################## 2024-05-02 08:58:39 INFO imstat::::casa ########################################## 2024-05-02 08:58:39 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:39 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 ) 2024-05-02 08:58:39 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:39 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:39 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:39 INFO imstat::image::statistics Determining stats for image testz.image 2024-05-02 08:58:39 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:39 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) 2024-05-02 08:58:39 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:39 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:39 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:39 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:39 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.597, -30.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:39 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.610, -30.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:39 INFO imstat:::: Values --- 2024-05-02 08:58:39 INFO imstat:::: -- flux density [flux]: 0.187451 Jy 2024-05-02 08:58:39 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:39 INFO imstat:::: -- maximum value [max]: 2.76828 Jy/beam 2024-05-02 08:58:39 INFO imstat:::: -- minimum value [min]: -0.185749 Jy/beam 2024-05-02 08:58:39 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2024-05-02 08:58:39 INFO imstat:::: -- position of min value (pixel) [minpos]: [31, 109, 0, 0] 2024-05-02 08:58:39 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.333, -30.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:39 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.356, -30.59.25.499, I, 2.25982e+11Hz 2024-05-02 08:58:39 INFO imstat:::: -- Sum of pixel values [sum]: 97.5737 Jy/beam 2024-05-02 08:58:39 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2189.16 Jy/beam.Jy/beam 2024-05-02 08:58:39 INFO imstat:::: Statistics --- 2024-05-02 08:58:39 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00595543 Jy/beam 2024-05-02 08:58:39 INFO imstat:::: -- Variance of the pixel values : 0.133588 Jy/beam 2024-05-02 08:58:39 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365497 Jy/beam 2024-05-02 08:58:39 INFO imstat:::: -- Root mean square [rms]: 0.365535 Jy/beam 2024-05-02 08:58:39 INFO imstat:::: -- Median of the pixel values [median]: -0.0939644 Jy/beam 2024-05-02 08:58:39 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0676308 Jy/beam 2024-05-02 08:58:39 INFO imstat:::: -- IQR [quartile]: 0.155104 Jy/beam 2024-05-02 08:58:39 INFO imstat:::: -- First quartile [q1]: -0.153538 Jy/beam 2024-05-02 08:58:39 INFO imstat:::: -- Third quartile [q3]: 0.00156572 Jy/beam 2024-05-02 08:58:39 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:39 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:39 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:39 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:39 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:39 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:39 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 2024-05-02 08:58:39 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:39.179141 End time: 2024-05-02 08:58:39.212754 2024-05-02 08:58:39 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:39 INFO imstat::::casa ########################################## 2024-05-02 08:58:39 INFO fixvis::::casa ########################################## 2024-05-02 08:58:39 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:39 INFO fixvis::::casa fixvis( vis='testx.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='-0.2333266794430661s 0deg', distances=[], datacolumn='DATA,CORRECTED' ) 2024-05-02 08:58:39 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:39 INFO fixvis::::casa Will overwrite original MS ... 2024-05-02 08:58:39 INFO fixvis::::casa Will only modify the visibilities in the columns DATA,CORRECTED 2024-05-02 08:58:39 INFO fixvis::::casa field: non-blanktest_0 2024-05-02 08:58:39 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:39 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2024-05-02 08:58:39 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:39 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2024-05-02 08:58:39 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2024-05-02 08:58:39 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0759 -029.59.29.9987 2024-05-02 08:58:39 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706453643230793 -0.5234533252689584 2024-05-02 08:58:39 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2024-05-02 08:58:39 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:39 INFO SubMS::parseColumnNames() Using DATA CORRECTED_DATA columns. 2024-05-02 08:58:39 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:39 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:39 INFO FixVis::makeSelection() Processing field 0 2024-05-02 08:58:40 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:39.360357 End time: 2024-05-02 08:58:39.572291 2024-05-02 08:58:40 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:40 INFO fixvis::::casa ########################################## 2024-05-02 08:58:40 INFO tclean::::casa ########################################## 2024-05-02 08:58:40 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:40 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 ) 2024-05-02 08:58:40 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:40 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:40 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:40 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:40 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:40 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:40 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:40 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:40 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:40 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:40 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:40 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:40 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.01066 2024-05-02 08:58:40 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2024-05-02 08:58:40 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:40 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:40 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58:40 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. 2024-05-02 08:58:41 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2024-05-02 08:58:41 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:41 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:39.612395 End time: 2024-05-02 08:58:40.584946 2024-05-02 08:58:41 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:41 INFO tclean::::casa ########################################## 2024-05-02 08:58:41 INFO imstat::::casa ########################################## 2024-05-02 08:58:41 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:41 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 ) 2024-05-02 08:58:41 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:41 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:41 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:41 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:41 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:41 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) 2024-05-02 08:58:41 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:41 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:41 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:41 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:41 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.569, -29.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:41 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.591, -29.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:41 INFO imstat:::: Values --- 2024-05-02 08:58:41 INFO imstat:::: -- flux density [flux]: 0.186776 Jy 2024-05-02 08:58:41 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:41 INFO imstat:::: -- maximum value [max]: 2.76655 Jy/beam 2024-05-02 08:58:41 INFO imstat:::: -- minimum value [min]: -0.185788 Jy/beam 2024-05-02 08:58:41 INFO imstat:::: -- position of max value (pixel) [maxpos]: [34, 64, 0, 0] 2024-05-02 08:58:41 INFO imstat:::: -- position of min value (pixel) [minpos]: [30, 110, 0, 0] 2024-05-02 08:58:41 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.307, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:41 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.338, -29.59.25.399, I, 2.25982e+11Hz 2024-05-02 08:58:41 INFO imstat:::: -- Sum of pixel values [sum]: 97.0186 Jy/beam 2024-05-02 08:58:41 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2184.63 Jy/beam.Jy/beam 2024-05-02 08:58:41 INFO imstat:::: Statistics --- 2024-05-02 08:58:41 INFO imstat:::: -- Mean of the pixel values [mean]: 0.00592155 Jy/beam 2024-05-02 08:58:41 INFO imstat:::: -- Variance of the pixel values : 0.133312 Jy/beam 2024-05-02 08:58:41 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.365119 Jy/beam 2024-05-02 08:58:41 INFO imstat:::: -- Root mean square [rms]: 0.365156 Jy/beam 2024-05-02 08:58:41 INFO imstat:::: -- Median of the pixel values [median]: -0.0935363 Jy/beam 2024-05-02 08:58:41 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.0680586 Jy/beam 2024-05-02 08:58:41 INFO imstat:::: -- IQR [quartile]: 0.155433 Jy/beam 2024-05-02 08:58:41 INFO imstat:::: -- First quartile [q1]: -0.153311 Jy/beam 2024-05-02 08:58:41 INFO imstat:::: -- Third quartile [q3]: 0.00212144 Jy/beam 2024-05-02 08:58:41 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:41 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:41 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:41 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:41 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:41 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:41 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 2024-05-02 08:58:41 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:40.587314 End time: 2024-05-02 08:58:40.605106 2024-05-02 08:58:41 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:41 INFO imstat::::casa ########################################## 2024-05-02 08:58:41 INFO fixvis::::casa ########################################## 2024-05-02 08:58:41 INFO fixvis::::casa ##### Begin Task: fixvis ##### 2024-05-02 08:58:41 INFO fixvis::::casa fixvis( vis='testx.ms', outputvis='testx.ms', field='0', refcode='J2000', reuse=True, phasecenter='-0.2333266794430661s 0deg', distances=[], datacolumn='DATA' ) 2024-05-02 08:58:41 WARN fixvis::::casa **ALERT:** Deprecated in CASA 5.9/6.3. Please use task **phaseshift** instead. 2024-05-02 08:58:41 INFO fixvis::::casa Will overwrite original MS ... 2024-05-02 08:58:41 INFO fixvis::::casa Will only modify the visibilities in the columns DATA 2024-05-02 08:58:41 INFO fixvis::::casa field: non-blanktest_0 2024-05-02 08:58:41 INFO fixvis::::casa old phasecenter RA, DEC J2000 18:00:02.3092 -029.59.29.9987 2024-05-02 08:58:41 INFO fixvis::::casa RA, DEC (rad) J2000 [-1.5706284] [-0.52345333] 2024-05-02 08:58:41 INFO fixvis::::casa When changing phase center, UVW coordinates will be recalculated. 2024-05-02 08:58:41 INFO fixvis::::casa No equinox given in parameter 'phasecenter': -0.2333266794430661s 0deg 2024-05-02 08:58:41 INFO fixvis::::casa Interpreting it as pair of offsets in (RA,DEC) ... 2024-05-02 08:58:41 INFO fixvis::::casa new phasecenter RA, DEC J2000 18:00:02.0759 -029.59.29.9987 2024-05-02 08:58:41 INFO fixvis::::casa RA, DEC (rad) J2000 -1.5706453643230793 -0.5234533252689584 2024-05-02 08:58:41 INFO fixvis::::casa FIELD table PHASE_DIR column changed for field 0. 2024-05-02 08:58:41 INFO Imager::open() Opening MeasurementSet /home/casatest/casa6/casatasks/tests/casashell_tests/work/test_task_fixvis_casashell/testx.ms 2024-05-02 08:58:41 INFO SubMS::parseColumnNames() Using DATA column. 2024-05-02 08:58:41 INFO FixVis::setFields Selecting fields 0 2024-05-02 08:58:41 INFO FixVis::makeSelection() 21600 rows selected out of 43200. 2024-05-02 08:58:41 INFO FixVis::makeSelection() Processing field 0 2024-05-02 08:58:41 INFO fixvis::::casa Task fixvis complete. Start time: 2024-05-02 08:58:40.727837 End time: 2024-05-02 08:58:40.950447 2024-05-02 08:58:41 INFO fixvis::::casa ##### End Task: fixvis ##### 2024-05-02 08:58:41 INFO fixvis::::casa ########################################## 2024-05-02 08:58:41 INFO tclean::::casa ########################################## 2024-05-02 08:58:41 INFO tclean::::casa ##### Begin Task: tclean ##### 2024-05-02 08:58:41 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 ) 2024-05-02 08:58:41 INFO tclean::::casa Verifying Input Parameters 2024-05-02 08:58:41 INFO SynthesisImagerVi2::selectData MS : testx.ms | Selecting on fields : 0 | [Opened in readonly mode] 2024-05-02 08:58:41 INFO SynthesisImagerVi2::selectData NRows selected : 21600 2024-05-02 08:58:41 INFO SynthesisImagerVi2::defineImage Define image coordinates for [testy] : 2024-05-02 08:58:41 INFO SynthesisImagerVi2::defineImage Impars : start 2024-05-02 08:58:41 INFO SynthesisImagerVi2::defineImage Shape : [128, 128, 1, 1]Spectral : [2.25982e+11] at [0] with increment [1.56671e+07] 2024-05-02 08:58:41 INFO SynthesisImagerVi2::defineImage Set Gridding options for [testy] with ftmachine : gridft 2024-05-02 08:58:41 INFO SynthesisImagerVi2::nSubCubeFitInMemory Required memory: 0.0008583 GB. Available mem.: 10.83 GB (rc, mem. fraction: 80%, memory: -) => Subcubes: 1. Processes on node: 1. 2024-05-02 08:58:41 INFO SynthesisImagerVi2::weight() Set imaging weights : Natural weighting 2024-05-02 08:58:41 INFO task_tclean::SynthesisDeconvolver::setupDeconvolution Set Deconvolution Options for [testy] : clark 2024-05-02 08:58:41 INFO SynthesisImager::makePSF ----------------------------------------------------------- Make PSF --------------------------------------------- 2024-05-02 08:58:42 INFO SIImageStore::calcSensitivity [testy] Theoretical sensitivity (Jy/bm):0.00481125 2024-05-02 08:58:42 INFO SIImageStore::getPSFGaussian Time to fit Gaussian to PSF 0.021847 2024-05-02 08:58:42 INFO SIImageStore::printBeamSet Beam : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2024-05-02 08:58:42 INFO task_tclean::SynthesisImagerVi2::makePrimaryBeam vi2 : Evaluating Primary Beam model onto image grid(s) 2024-05-02 08:58:42 INFO task_tclean::SynthesisImager::executeMajorCycle ----------------------------------------------------------- Run Major Cycle 1 ------------------------------------- 2024-05-02 08:58:42 INFO task_tclean::SDAlgorithmBase::restore [testy] : Restoring model image. 2024-05-02 08:58:42 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. 2024-05-02 08:58:42 INFO task_tclean::SIImageStore::restore Beam for chan : 0 : 2.27941 arcsec, 2.01117 arcsec, 30.5326 deg 2024-05-02 08:58:42 INFO tclean::::casa Searching for images with prefix 'testy'... Found these, writing history into them: ['testy.model', 'testy.residual', 'testy.psf', 'testy.image', 'testy.sumwt', 'testy.pb'] 2024-05-02 08:58:42 INFO tclean::::casa Task tclean complete. Start time: 2024-05-02 08:58:40.985278 End time: 2024-05-02 08:58:42.135666 2024-05-02 08:58:42 INFO tclean::::casa ##### End Task: tclean ##### 2024-05-02 08:58:42 INFO tclean::::casa ########################################## 2024-05-02 08:58:42 INFO imstat::::casa ########################################## 2024-05-02 08:58:42 INFO imstat::::casa ##### Begin Task: imstat ##### 2024-05-02 08:58:42 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 ) 2024-05-02 08:58:42 INFO imstat::CasacRegionManager::fromBCS No directional region specified. Using full positional plane. 2024-05-02 08:58:42 INFO imstat::CasacRegionManager::fromBCS Using all spectral channels. 2024-05-02 08:58:42 INFO imstat::CasacRegionManager::fromBCS Using polarizations ALL 2024-05-02 08:58:42 INFO imstat::image::statistics Determining stats for image testy.image 2024-05-02 08:58:42 INFO imstat::ImageRegion::fromRecord Selected bounding box : 2024-05-02 08:58:42 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) 2024-05-02 08:58:42 INFO imstat::ImageStatsCalculator::_logStartup Statistics calculated using Classic algorithm 2024-05-02 08:58:42 INFO imstat::ImageStatsCalculator::_logStartup Regions --- 2024-05-02 08:58:42 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (pixel) [blc]: [0, 0, 0, 0] 2024-05-02 08:58:42 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (pixel) [trc]: [127, 127, 0, 0] 2024-05-02 08:58:42 INFO imstat::ImageStatsCalculator::_logStartup -- bottom-left corner (world) [blcf]: 18:00:02.569, -29.59.36.399, I, 2.25982e+11Hz 2024-05-02 08:58:42 INFO imstat::ImageStatsCalculator::_logStartup -- top-right corner (world) [trcf]: 18:00:01.591, -29.59.23.699, I, 2.25982e+11Hz 2024-05-02 08:58:42 INFO imstat:::: Values --- 2024-05-02 08:58:42 INFO imstat:::: -- flux density [flux]: 0.0113138 Jy 2024-05-02 08:58:42 INFO imstat:::: -- number of points [npts]: 16384 2024-05-02 08:58:42 INFO imstat:::: -- maximum value [max]: 2.76818 Jy/beam 2024-05-02 08:58:42 INFO imstat:::: -- minimum value [min]: -0.185778 Jy/beam 2024-05-02 08:58:42 INFO imstat:::: -- position of max value (pixel) [maxpos]: [64, 64, 0, 0] 2024-05-02 08:58:42 INFO imstat:::: -- position of min value (pixel) [minpos]: [67, 19, 0, 0] 2024-05-02 08:58:42 INFO imstat:::: -- position of max value (world) [maxposf]: 18:00:02.076, -29.59.29.999, I, 2.25982e+11Hz 2024-05-02 08:58:42 INFO imstat:::: -- position of min value (world) [minposf]: 18:00:02.053, -29.59.34.499, I, 2.25982e+11Hz 2024-05-02 08:58:42 INFO imstat:::: -- Sum of pixel values [sum]: 5.87683 Jy/beam 2024-05-02 08:58:42 INFO imstat:::: -- Sum of squared pixel values [sumsq]: 2178.04 Jy/beam.Jy/beam 2024-05-02 08:58:42 INFO imstat:::: Statistics --- 2024-05-02 08:58:42 INFO imstat:::: -- Mean of the pixel values [mean]: 0.000358693 Jy/beam 2024-05-02 08:58:42 INFO imstat:::: -- Variance of the pixel values : 0.132945 Jy/beam 2024-05-02 08:58:42 INFO imstat:::: -- Standard deviation of the Mean [sigma]: 0.364616 Jy/beam 2024-05-02 08:58:42 INFO imstat:::: -- Root mean square [rms]: 0.364605 Jy/beam 2024-05-02 08:58:42 INFO imstat:::: -- Median of the pixel values [median]: -0.0988264 Jy/beam 2024-05-02 08:58:42 INFO imstat:::: -- Median of the deviations [medabsdevmed]: 0.064913 Jy/beam 2024-05-02 08:58:42 INFO imstat:::: -- IQR [quartile]: 0.146844 Jy/beam 2024-05-02 08:58:42 INFO imstat:::: -- First quartile [q1]: -0.155414 Jy/beam 2024-05-02 08:58:42 INFO imstat:::: -- Third quartile [q3]: -0.00856965 Jy/beam 2024-05-02 08:58:42 INFO imstat::ImageStatsCalculator::statistics Sum column unit = Jy/beam 2024-05-02 08:58:42 INFO imstat::ImageStatsCalculator::statistics+ Mean column unit = Jy/beam 2024-05-02 08:58:42 INFO imstat::ImageStatsCalculator::statistics+ Std_dev column unit = Jy/beam 2024-05-02 08:58:42 INFO imstat::ImageStatsCalculator::statistics+ Minimum column unit = Jy/beam 2024-05-02 08:58:42 INFO imstat::ImageStatsCalculator::statistics+ Maximum column unit = Jy/beam 2024-05-02 08:58:42 INFO imstat::ImageStatsCalculator::statistics Npts Sum Mean Rms Std_dev Minimum Maximum 2024-05-02 08:58:42 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 2024-05-02 08:58:42 INFO imstat::::casa Task imstat complete. Start time: 2024-05-02 08:58:42.138604 End time: 2024-05-02 08:58:42.161611 2024-05-02 08:58:42 INFO imstat::::casa ##### End Task: imstat ##### 2024-05-02 08:58:42 INFO imstat::::casa ########################################## Using configuration file ~/.casa/config.py Using matplotlib backend: agg CASA 6.6.4.20a14347.dev1 -- Common Astronomy Software Applications [6.6.4.20a14347.dev1] Testing non-existing datacolumn - expected error *** Expected error *** Running tests... ---------------------------------------------------------------------- ............... ---------------------------------------------------------------------- Ran 15 tests in 25.418s OK Generating XML reports...