Commits

Takahiro Tsutsumi authored a1013edb630 
Removed the duplicated calls for cube imaging
48b797291ef
master

casa5/gcwrap/python/scripts/tests/test_stk_alma_pipeline_imaging.py

Modified
518 518 ''' Modified test_dict costructed by casatestutils add_to_dict to include only
519 519 the task commands executed and also add self.parallel value to the dictionary.
520 520 The cube imaging cases usually have if-else conditional based on parallel mode is on or not
521 521 to trigger different set of tclean commands.
522 522 Assumption: self.parallel is used to trigger different tclean commands at iter1 step.
523 523 For self.parallel=True, iter1 has two tclean commands (2nd and 3rd tclean commands within
524 524 each relevante test(cube) and so in test_dict['taskcall'], 1st(iter0) and 2nd and 3rd commands
525 525 are the ones acutually executed and should remove 4th (self.parallel=False) case.
526 526 '''
527 527 if testname in output:
528 - if 'taskcall' in output[testname] and len(output[testname]['taskcall'])==4:
528 + if 'taskcall' in output[testname] and len(output[testname]['taskcall'])==3:
529 529 if parallel:
530 530 # 0,1,2th in the list are used pop last one
531 531 output[testname]['taskcall'].pop()
532 532 else:
533 533 output[testname]['taskcall'].pop(1)
534 - output[testname]['taskcall'].pop(1)
534 + #output[testname]['taskcall'].pop(1)
535 535 output[testname]['self.parallel']=parallel
536 536 ​
537 537 def remove_prefix(self,string, prefix):
538 538 ''' Remove a specified prefix string from string '''
539 539 return string[string.startswith(prefix) and len(prefix):]
540 540 ​
541 541 ##############################################
542 542 # TESTS
543 543 ##############################################
544 544 test_dict = {}
577 577 fastnoise=False, savemodel='none', parallel=self.parallel,
578 578 verbose=True)
579 579 ​
580 580 # move files to iter1
581 581 print('Copying iter0 files to iter1')
582 582 self.copy_products(file_name+'0', file_name+'1')
583 583 ​
584 584 print("STARTING: iter1 routine")
585 585 ​
586 586 # iter1 (restart)
587 - if self.parallel:
588 - tclean(vis=self.msfile, imagename=file_name+'1', field='1', \
589 - spw=['0'], imsize=[80, 80], antenna=['0,1,2,3,4,5,6,7,8'], \
590 - scan=['8,12,16'], intent='OBSERVE_TARGET#ON_SOURCE', \
591 - datacolumn='data', cell=['1.1arcsec'], phasecenter='ICRS '
592 - '00:45:54.3836 -073.15.29.413', stokes='I', specmode='cube', \
593 - nchan=508, start='220.2526743594GHz', width='0.2441741MHz',\
594 - outframe='LSRK', perchanweightdensity=False, \
595 - usepointing=False, pblimit=0.2, nsigma=0.0, \
596 - gridder='standard', mosweight=False, \
597 - deconvolver='hogbom', restoringbeam='common', restoration=True, pbcor=True, \
598 - weighting='briggs', robust=0.5, npixels=0, niter=20000, \
599 - threshold='0.354Jy', interactive=0, usemask='auto'
600 - '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
601 - lownoisethreshold=2.0, negativethreshold=0.0, \
602 - minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
603 - minpercentchange=1.0, fastnoise=False, restart=True, \
604 - calcres=False, calcpsf=False, savemodel='none', \
605 - parallel=True, verbose=True)
606 - else:
607 - tclean(vis=self.msfile, imagename=file_name+'1', field='1', \
608 - spw=['0'], imsize=[80, 80], antenna=['0,1,2,3,4,5,6,7,8'], \
609 - scan=['8,12,16'], intent='OBSERVE_TARGET#ON_SOURCE', \
610 - datacolumn='data', cell=['1.1arcsec'], phasecenter='ICRS '
611 - '00:45:54.3836 -073.15.29.413', stokes='I', specmode='cube', \
612 - nchan=508, start='220.2526743594GHz', width='0.2441741MHz',\
613 - outframe='LSRK', perchanweightdensity=False, \
614 - usepointing=False, pblimit=0.2, nsigma=0.0, \
615 - gridder='standard', mosweight=False, \
616 - deconvolver='hogbom', restoration=True, pbcor=True, \
617 - weighting='briggs', robust=0.5, npixels=0, niter=20000, \
618 - threshold='0.354Jy', interactive=0, usemask='auto'
619 - '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
620 - lownoisethreshold=2.0, negativethreshold=0.0, \
621 - minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
622 - minpercentchange=1.0, fastnoise=False, restart=True, \
623 - calcres=False, calcpsf=False, savemodel='none', \
624 - restoringbeam='common', parallel=False, verbose=True)
625 -​
587 + tclean(vis=self.msfile, imagename=file_name+'1', field='1', \
588 + spw=['0'], imsize=[80, 80], antenna=['0,1,2,3,4,5,6,7,8'], \
589 + scan=['8,12,16'], intent='OBSERVE_TARGET#ON_SOURCE', \
590 + datacolumn='data', cell=['1.1arcsec'], phasecenter='ICRS '
591 + '00:45:54.3836 -073.15.29.413', stokes='I', specmode='cube', \
592 + nchan=508, start='220.2526743594GHz', width='0.2441741MHz',\
593 + outframe='LSRK', perchanweightdensity=False, \
594 + usepointing=False, pblimit=0.2, nsigma=0.0, \
595 + gridder='standard', mosweight=False, \
596 + deconvolver='hogbom', restoringbeam='common', restoration=True, pbcor=True, \
597 + weighting='briggs', robust=0.5, npixels=0, niter=20000, \
598 + threshold='0.354Jy', interactive=0, usemask='auto'
599 + '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
600 + lownoisethreshold=2.0, negativethreshold=0.0, \
601 + minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
602 + minpercentchange=1.0, fastnoise=False, restart=True, \
603 + calcres=False, calcpsf=False, savemodel='none', \
604 + parallel=self.parallel, verbose=True)
626 605 ​
627 606 # retrieve per-channel beam statistics
628 607 bmin_dict, bmaj_dict, pa_dict = \
629 608 self.cube_beam_stats(image=img+'.psf')
630 609 ​
631 610 report0 = th.checkall(imgexist = self.image_list(img, 'standard'))
632 611 ​
633 612 # .image report(test_standard_cube)
634 613 im_stats_dict = self.image_stats(image=img+'.image', fit_region = \
635 614 'ellipse[[11.47881897deg, -73.25881015deg], [9.0414arcsec, 8.4854arcsec], 90.00000000deg]')
921 900 fastnoise=False, savemodel='none', parallel=self.parallel,
922 901 verbose=True)
923 902 ​
924 903 # move files to iter1
925 904 print('Copying iter0 files to iter1')
926 905 self.copy_products(file_name+'0', file_name+'1')
927 906 ​
928 907 print("STARTING: iter1 routine")
929 908 ​
930 909 # iter1 (restart)
931 - if self.parallel:
932 - tclean(vis=self.msfile, imagename=file_name+'1', field='1', \
933 - spw=['0'], imsize=[80, 80], antenna=['0,1,2,3,4,5,6,7,8'], \
934 - scan=['8,12,16'], intent='OBSERVE_TARGET#ON_SOURCE', \
935 - datacolumn='data', cell=['1.1arcsec'], phasecenter='ICRS '
936 - '00:45:54.3836 -073.15.29.413', stokes='I', specmode='cube', \
937 - nchan=508, start='220.2526743594GHz', width='0.2441741MHz',\
938 - outframe='LSRK', perchanweightdensity=True, \
939 - usepointing=False, pblimit=0.2, nsigma=0.0, \
940 - gridder='standard', mosweight=False, \
941 - deconvolver='hogbom', restoration=True, restoringbeam='common', pbcor=True, \
942 - weighting='briggs', robust=0.5, npixels=0, niter=20000, \
943 - threshold='0.354Jy', interactive=0, usemask='auto'
944 - '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
945 - lownoisethreshold=2.0, negativethreshold=0.0, \
946 - minbeamfrac=0.08, growiterations=75, dogrowprune=True, \
947 - minpercentchange=1.0, fastnoise=False, restart=True, \
948 - calcres=False, calcpsf=False, savemodel='none', \
949 - parallel=True, verbose=True)
950 - else:
951 - tclean(vis=self.msfile, imagename=file_name+'1', field='1', \
952 - spw=['0'], imsize=[80, 80], antenna=['0,1,2,3,4,5,6,7,8'], \
953 - scan=['8,12,16'], intent='OBSERVE_TARGET#ON_SOURCE', \
954 - datacolumn='data', cell=['1.1arcsec'], phasecenter='ICRS '
955 - '00:45:54.3836 -073.15.29.413', stokes='I', specmode='cube', \
956 - nchan=508, start='220.2526743594GHz', width='0.2441741MHz',\
957 - outframe='LSRK', perchanweightdensity=True, \
958 - usepointing=False, pblimit=0.2, nsigma=0.0, \
959 - gridder='standard', mosweight=False, \
960 - deconvolver='hogbom', restoration=True, pbcor=True, \
961 - weighting='briggs', robust=0.5, npixels=0, niter=20000, \
962 - threshold='0.354Jy', interactive=0, usemask='auto'
963 - '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
964 - lownoisethreshold=2.0, negativethreshold=0.0, \
965 - minbeamfrac=0.08, growiterations=75, dogrowprune=True, \
966 - minpercentchange=1.0, fastnoise=False, restart=True, \
967 - calcres=False, calcpsf=False, savemodel='none', \
968 - restoringbeam='common', parallel=False, verbose=True)
910 + tclean(vis=self.msfile, imagename=file_name+'1', field='1', \
911 + spw=['0'], imsize=[80, 80], antenna=['0,1,2,3,4,5,6,7,8'], \
912 + scan=['8,12,16'], intent='OBSERVE_TARGET#ON_SOURCE', \
913 + datacolumn='data', cell=['1.1arcsec'], phasecenter='ICRS '
914 + '00:45:54.3836 -073.15.29.413', stokes='I', specmode='cube', \
915 + nchan=508, start='220.2526743594GHz', width='0.2441741MHz',\
916 + outframe='LSRK', perchanweightdensity=True, \
917 + usepointing=False, pblimit=0.2, nsigma=0.0, \
918 + gridder='standard', mosweight=False, \
919 + deconvolver='hogbom', restoration=True, restoringbeam='common', pbcor=True, \
920 + weighting='briggs', robust=0.5, npixels=0, niter=20000, \
921 + threshold='0.354Jy', interactive=0, usemask='auto'
922 + '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
923 + lownoisethreshold=2.0, negativethreshold=0.0, \
924 + minbeamfrac=0.08, growiterations=75, dogrowprune=True, \
925 + minpercentchange=1.0, fastnoise=False, restart=True, \
926 + calcres=False, calcpsf=False, savemodel='none', \
927 + parallel=self.parallel, verbose=True)
969 928 ​
970 929 # retrieve per-channel beam statistics
971 930 bmin_dict, bmaj_dict, pa_dict = \
972 931 self.cube_beam_stats(image=img+'.psf')
973 932 ​
974 933 ​
975 934 report0 = th.checkall(imgexist = self.image_list(img, 'standard'))
976 935 ​
977 936 # .image report(test_standard_cube_pcwdT)
978 937 im_stats_dict = self.image_stats(image=img+'.image', fit_region = \
1268 1227 fastnoise=False, savemodel='none', parallel=self.parallel,
1269 1228 verbose=True)
1270 1229 ​
1271 1230 # move files to iter1
1272 1231 print('Copying iter0 files to iter1')
1273 1232 self.copy_products(file_name+'0', file_name+'1')
1274 1233 ​
1275 1234 print("STARTING: iter1 routine")
1276 1235 ​
1277 1236 # iter1 (restart)
1278 - if self.parallel:
1279 - tclean(vis=self.msfile, imagename=file_name+'1', field='1', \
1280 - spw=['0'], imsize=[80, 80], antenna=['0,1,2,3,4,5,6,7,8'], \
1281 - scan=['8,12,16'], intent='OBSERVE_TARGET#ON_SOURCE', \
1282 - datacolumn='data', cell=['1.1arcsec'], phasecenter='ICRS '
1283 - '00:45:54.3836 -073.15.29.413', stokes='I', specmode='cube', \
1284 - nchan=508, start='220.2526743594GHz', width='0.2441741MHz',\
1285 - outframe='LSRK', perchanweightdensity=True, \
1286 - usepointing=False, pblimit=0.2, nsigma=0.0, \
1287 - gridder='standard', mosweight=False, \
1288 - deconvolver='hogbom', restoration=True, restoringbeam='common', pbcor=True, \
1289 - weighting='briggsbwtaper', robust=0.5, npixels=0, niter=20000, \
1290 - threshold='0.354Jy', interactive=0, usemask='auto'
1291 - '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
1292 - lownoisethreshold=2.0, negativethreshold=0.0, \
1293 - minbeamfrac=0.08, growiterations=75, dogrowprune=True, \
1294 - minpercentchange=1.0, fastnoise=False, restart=True, \
1295 - calcres=False, calcpsf=False, savemodel='none', \
1296 - parallel=True, verbose=True)
1297 - else:
1298 - tclean(vis=self.msfile, imagename=file_name+'1', field='1', \
1299 - spw=['0'], imsize=[80, 80], antenna=['0,1,2,3,4,5,6,7,8'], \
1300 - scan=['8,12,16'], intent='OBSERVE_TARGET#ON_SOURCE', \
1301 - datacolumn='data', cell=['1.1arcsec'], phasecenter='ICRS '
1302 - '00:45:54.3836 -073.15.29.413', stokes='I', specmode='cube', \
1303 - nchan=508, start='220.2526743594GHz', width='0.2441741MHz',\
1304 - outframe='LSRK', perchanweightdensity=True, \
1305 - usepointing=False, pblimit=0.2, nsigma=0.0, \
1306 - gridder='standard', mosweight=False, \
1307 - deconvolver='hogbom', restoration=True, pbcor=True, \
1308 - weighting='briggsbwtaper', robust=0.5, npixels=0, niter=20000, \
1309 - threshold='0.354Jy', interactive=0, usemask='auto'
1310 - '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
1311 - lownoisethreshold=2.0, negativethreshold=0.0, \
1312 - minbeamfrac=0.08, growiterations=75, dogrowprune=True, \
1313 - minpercentchange=1.0, fastnoise=False, restart=True, \
1314 - calcres=False, calcpsf=False, savemodel='none', \
1315 - restoringbeam='common', parallel=False, verbose=True)
1237 + tclean(vis=self.msfile, imagename=file_name+'1', field='1', \
1238 + spw=['0'], imsize=[80, 80], antenna=['0,1,2,3,4,5,6,7,8'], \
1239 + scan=['8,12,16'], intent='OBSERVE_TARGET#ON_SOURCE', \
1240 + datacolumn='data', cell=['1.1arcsec'], phasecenter='ICRS '
1241 + '00:45:54.3836 -073.15.29.413', stokes='I', specmode='cube', \
1242 + nchan=508, start='220.2526743594GHz', width='0.2441741MHz',\
1243 + outframe='LSRK', perchanweightdensity=True, \
1244 + usepointing=False, pblimit=0.2, nsigma=0.0, \
1245 + gridder='standard', mosweight=False, \
1246 + deconvolver='hogbom', restoration=True, restoringbeam='common', pbcor=True, \
1247 + weighting='briggsbwtaper', robust=0.5, npixels=0, niter=20000, \
1248 + threshold='0.354Jy', interactive=0, usemask='auto'
1249 + '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
1250 + lownoisethreshold=2.0, negativethreshold=0.0, \
1251 + minbeamfrac=0.08, growiterations=75, dogrowprune=True, \
1252 + minpercentchange=1.0, fastnoise=False, restart=True, \
1253 + calcres=False, calcpsf=False, savemodel='none', \
1254 + parallel=self.parallel, verbose=True)
1316 1255 ​
1317 1256 # retrieve per-channel beam statistics
1318 1257 bmin_dict, bmaj_dict, pa_dict = \
1319 1258 self.cube_beam_stats(image=img+'.psf')
1320 1259 ​
1321 1260 ​
1322 1261 report0 = th.checkall(imgexist = self.image_list(img, 'standard'))
1323 1262 ​
1324 1263 # .image report(test_standard_cube_briggsbwtaper)
1325 1264 im_stats_dict = self.image_stats(image=img+'.image', fit_region = \
4294 4233 fastnoise=False, savemodel='none', parallel=self.parallel,
4295 4234 verbose=True)
4296 4235 ​
4297 4236 # move files to iter1
4298 4237 print('Copying iter0 files to iter1')
4299 4238 self.copy_products(file_name+'0', file_name+'1')
4300 4239 ​
4301 4240 print("STARTING: iter1 routine")
4302 4241 ​
4303 4242 # iter1 (restart)
4304 - if self.parallel:
4305 - tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
4306 - antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
4307 - intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
4308 - imagename=file_name+'1', imsize=[108, 108], \
4309 - cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
4310 - ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
4311 - start='220.2526743594GHz', width='0.2441741MHz', \
4312 - outframe='LSRK', perchanweightdensity=False, \
4313 - gridder='mosaic', mosweight=True, \
4314 - usepointing=False, pblimit=0.2, deconvolver='hogbom', \
4315 - restoration=True, restoringbeam='common', pbcor=True, \
4316 - weighting='briggs', robust=0.5,\
4317 - npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
4318 - interactive=0, usemask='auto-multithresh', \
4319 - sidelobethreshold=1.25, noisethreshold=5.0, \
4320 - lownoisethreshold=2.0, negativethreshold=0.0, \
4321 - minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4322 - minpercentchange=1.0, fastnoise=False, restart=True, \
4323 - savemodel='none', calcres=False, calcpsf=False, \
4324 - parallel=True, verbose=True)
4325 - else:
4326 - tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
4327 - antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
4328 - intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
4329 - imagename=file_name+'1', imsize=[108, 108], \
4330 - cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
4331 - ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
4332 - start='220.2526743594GHz', width='0.2441741MHz', \
4333 - outframe='LSRK', perchanweightdensity=False, \
4334 - gridder='mosaic', mosweight=True, \
4335 - usepointing=False, pblimit=0.2, deconvolver='hogbom', \
4336 - restoration=True, pbcor=True, weighting='briggs', robust=0.5,\
4337 - npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
4338 - interactive=0, usemask='auto-multithresh', \
4339 - sidelobethreshold=1.25, noisethreshold=5.0, \
4340 - lownoisethreshold=2.0, negativethreshold=0.0, \
4341 - minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4342 - minpercentchange=1.0, fastnoise=False, restart=True, \
4343 - savemodel='none', calcres=False, calcpsf=False, \
4344 - restoringbeam='common', parallel=False, verbose=True)
4243 + tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
4244 + antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
4245 + intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
4246 + imagename=file_name+'1', imsize=[108, 108], \
4247 + cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
4248 + ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
4249 + start='220.2526743594GHz', width='0.2441741MHz', \
4250 + outframe='LSRK', perchanweightdensity=False, \
4251 + gridder='mosaic', mosweight=True, \
4252 + usepointing=False, pblimit=0.2, deconvolver='hogbom', \
4253 + restoration=True, restoringbeam='common', pbcor=True, \
4254 + weighting='briggs', robust=0.5,\
4255 + npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
4256 + interactive=0, usemask='auto-multithresh', \
4257 + sidelobethreshold=1.25, noisethreshold=5.0, \
4258 + lownoisethreshold=2.0, negativethreshold=0.0, \
4259 + minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4260 + minpercentchange=1.0, fastnoise=False, restart=True, \
4261 + savemodel='none', calcres=False, calcpsf=False, \
4262 + parallel=self.parallel, verbose=True)
4345 4263 ​
4346 4264 # retrieve per-channel beam statistics
4347 4265 bmin_dict, bmaj_dict, pa_dict = \
4348 4266 self.cube_beam_stats(img+'.psf')
4349 4267 ​
4350 4268 report0 = th.checkall(imgexist = self.image_list(img, 'mosaic'))
4351 4269 ​
4352 4270 # .image report (test_mosaic_cube)
4353 4271 im_stats_dict = self.image_stats(img+'.image', fit_region = \
4354 4272 'ellipse[[11.48661818deg, -73.26292371deg], [8.2211arcsec, 7.4698arcsec], 90.00000000deg]', field_regions = \
4644 4562 fastnoise=False, savemodel='none', parallel=self.parallel,
4645 4563 verbose=True)
4646 4564 ​
4647 4565 # move files to iter1
4648 4566 print('Copying iter0 files to iter1')
4649 4567 self.copy_products(file_name+'0', file_name+'1')
4650 4568 ​
4651 4569 print("STARTING: iter1 routine")
4652 4570 ​
4653 4571 # iter1 (restart)
4654 - if self.parallel:
4655 - tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
4656 - antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
4657 - intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
4658 - imagename=file_name+'1', imsize=[108, 108], \
4659 - cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
4660 - ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
4661 - start='220.2526743594GHz', width='0.2441741MHz', \
4662 - outframe='LSRK', perchanweightdensity=True, \
4663 - gridder='mosaic', mosweight=True, \
4664 - usepointing=False, pblimit=0.2, deconvolver='hogbom', \
4665 - restoration=True, restoringbeam='common', pbcor=True, \
4666 - weighting='briggs', robust=0.5,\
4667 - npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
4668 - interactive=0, usemask='auto-multithresh', \
4669 - sidelobethreshold=1.25, noisethreshold=5.0, \
4670 - lownoisethreshold=2.0, negativethreshold=0.0, \
4671 - minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4672 - minpercentchange=1.0, fastnoise=False, restart=True, \
4673 - savemodel='none', calcres=False, calcpsf=False, \
4674 - parallel=True, verbose=True)
4675 - else:
4676 - tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
4677 - antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
4678 - intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
4679 - imagename=file_name+'1', imsize=[108, 108], \
4680 - cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
4681 - ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
4682 - start='220.2526743594GHz', width='0.2441741MHz', \
4683 - outframe='LSRK', perchanweightdensity=True, \
4684 - gridder='mosaic', mosweight=True, \
4685 - usepointing=False, pblimit=0.2, deconvolver='hogbom', \
4686 - restoration=True, pbcor=True, weighting='briggs', robust=0.5,\
4687 - npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
4688 - interactive=0, usemask='auto-multithresh', \
4689 - sidelobethreshold=1.25, noisethreshold=5.0, \
4690 - lownoisethreshold=2.0, negativethreshold=0.0, \
4691 - minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4692 - minpercentchange=1.0, fastnoise=False, restart=True, \
4693 - savemodel='none', calcres=False, calcpsf=False, \
4694 - restoringbeam='common', parallel=False, verbose=True)
4572 + tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
4573 + antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
4574 + intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
4575 + imagename=file_name+'1', imsize=[108, 108], \
4576 + cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
4577 + ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
4578 + start='220.2526743594GHz', width='0.2441741MHz', \
4579 + outframe='LSRK', perchanweightdensity=True, \
4580 + gridder='mosaic', mosweight=True, \
4581 + usepointing=False, pblimit=0.2, deconvolver='hogbom', \
4582 + restoration=True, restoringbeam='common', pbcor=True, \
4583 + weighting='briggs', robust=0.5,\
4584 + npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
4585 + interactive=0, usemask='auto-multithresh', \
4586 + sidelobethreshold=1.25, noisethreshold=5.0, \
4587 + lownoisethreshold=2.0, negativethreshold=0.0, \
4588 + minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4589 + minpercentchange=1.0, fastnoise=False, restart=True, \
4590 + savemodel='none', calcres=False, calcpsf=False, \
4591 + parallel=self.parallel, verbose=True)
4695 4592 ​
4696 4593 # retrieve per-channel beam statistics
4697 4594 bmin_dict, bmaj_dict, pa_dict = \
4698 4595 self.cube_beam_stats(img+'.psf')
4699 4596 ​
4700 4597 report0 = th.checkall(imgexist = self.image_list(img, 'mosaic'))
4701 4598 ​
4702 4599 # .image report (test_mosaic_cube_pcwdT)
4703 4600 im_stats_dict = self.image_stats(img+'.image', fit_region = \
4704 4601 'ellipse[[11.48661818deg, -73.26292371deg], [8.2211arcsec, 7.4698arcsec], 90.00000000deg]', field_regions = \
5025 4922 fastnoise=False, savemodel='none', parallel=self.parallel,
5026 4923 verbose=True)
5027 4924 ​
5028 4925 # move files to iter1
5029 4926 print('Copying iter0 files to iter1')
5030 4927 self.copy_products(file_name+'0', file_name+'1')
5031 4928 ​
5032 4929 print("STARTING: iter1 routine")
5033 4930 ​
5034 4931 # iter1 (restart)
5035 - if self.parallel:
5036 - tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
5037 - antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
5038 - intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
5039 - imagename=file_name+'1', imsize=[108, 108], \
5040 - cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
5041 - ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
5042 - start='220.2526743594GHz', width='0.2441741MHz', \
5043 - outframe='LSRK', perchanweightdensity=True, \
5044 - gridder='mosaic', mosweight=True, \
5045 - usepointing=False, pblimit=0.2, deconvolver='hogbom', \
5046 - restoration=True, restoringbeam='common', \
5047 - pbcor=True, weighting='briggsbwtaper', robust=0.5,\
5048 - npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
5049 - interactive=0, usemask='auto-multithresh', \
5050 - sidelobethreshold=1.25, noisethreshold=5.0, \
5051 - lownoisethreshold=2.0, negativethreshold=0.0, \
5052 - minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
5053 - minpercentchange=1.0, fastnoise=False, restart=True, \
5054 - savemodel='none', calcres=False, calcpsf=False, \
5055 - parallel=True, verbose=True)
5056 - else:
5057 - tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
5058 - antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
5059 - intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
5060 - imagename=file_name+'1', imsize=[108, 108], \
5061 - cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
5062 - ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
5063 - start='220.2526743594GHz', width='0.2441741MHz', \
5064 - outframe='LSRK', perchanweightdensity=True, \
5065 - gridder='mosaic', mosweight=True, \
5066 - usepointing=False, pblimit=0.2, deconvolver='hogbom', \
5067 - restoration=True, pbcor=True, weighting='briggsbwtaper', robust=0.5,\
5068 - npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
5069 - interactive=0, usemask='auto-multithresh', \
5070 - sidelobethreshold=1.25, noisethreshold=5.0, \
5071 - lownoisethreshold=2.0, negativethreshold=0.0, \
5072 - minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
5073 - minpercentchange=1.0, fastnoise=False, restart=True, \
5074 - savemodel='none', calcres=False, calcpsf=False, \
5075 - restoringbeam='common', parallel=False, verbose=True)
4932 + tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
4933 + antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
4934 + intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
4935 + imagename=file_name+'1', imsize=[108, 108], \
4936 + cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
4937 + ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
4938 + start='220.2526743594GHz', width='0.2441741MHz', \
4939 + outframe='LSRK', perchanweightdensity=True, \
4940 + gridder='mosaic', mosweight=True, \
4941 + usepointing=False, pblimit=0.2, deconvolver='hogbom', \
4942 + restoration=True, restoringbeam='common', \
4943 + pbcor=True, weighting='briggsbwtaper', robust=0.5,\
4944 + npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
4945 + interactive=0, usemask='auto-multithresh', \
4946 + sidelobethreshold=1.25, noisethreshold=5.0, \
4947 + lownoisethreshold=2.0, negativethreshold=0.0, \
4948 + minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4949 + minpercentchange=1.0, fastnoise=False, restart=True, \
4950 + savemodel='none', calcres=False, calcpsf=False, \
4951 + parallel=self.parallel, verbose=True)
5076 4952 ​
5077 4953 # retrieve per-channel beam statistics
5078 4954 bmin_dict, bmaj_dict, pa_dict = \
5079 4955 self.cube_beam_stats(img+'.psf')
5080 4956 ​
5081 4957 report0 = th.checkall(imgexist = self.image_list(img, 'mosaic'))
5082 4958 ​
5083 4959 # .image report (test_mosaic_cube_briggsbwtaper)
5084 4960 im_stats_dict = self.image_stats(img+'.image', fit_region = \
5085 4961 'ellipse[[11.48661818deg, -73.26292371deg], [8.2211arcsec, 7.4698arcsec], 90.00000000deg]', field_regions = \

Everything looks good. We'll let you know here if there's anything you should know about.

Add shortcut