raise Error('The casa commands log is not executable!')

# # This file contains CASA commands run by the pipeline. Although all
# commands required to calibrate the data are included here, this file cannot
# be executed, nor does it contain heuristic and flagging calculations
# performed by pipeline code. This file is useful to understand which CASA
# commands are being run by each pipeline task. If one wishes to re-run the
# pipeline, one should use the pipeline script linked on the front page or By
# Task page of the weblog. Some stages may not have any commands listed here,
# e.g. hifa_importdata if conversion from ASDM to MS is not required.

# hifv_importdata(nocopy=True, vis=['VLASS1.1.sb34516741.eb34573529.58031.90910875.ms'], session=['session_1'])
#
# If required, ASDMs are converted to measurement sets.
#

# hif_editimlist(parameter_file='parameter.12482.list')
#

# hif_transformimagedata(datacolumn='corrected', modify_weights=False, clear_pointing=True)
#
mstransform(timebin='0s',
            outputvis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T19t21.J163121+343000_P12482v1_2017_11_10T21_39_09.970/working/VLASS1.1.sb34516741.eb34573529.58031.90910875_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T19t21.J163121+343000_P12482v1_2017_11_10T21_39_09.970/rawdata/VLASS1.1.sb34516741.eb34573529.58031.90910875.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='10299,10298,10295,10294,10297,10296,10291,10290,10293,10292,10451,10450,11911,11910,11913,11912,11915,11914,11917,11916,11918,10709,10708,10703,10702,10701,10700,10707,10706,10705,10704,10422,11229,10034,10026,10027,10024,10025,10022,10023,10020,10021,10421,10426,10427,10424,10425,10840,10841,10842,10843,10844,10845,10846,10847,10848,10849,11228,12051,12050,12053,12052,12055,12054,12057,12056,10835,10834,10837,10836,10831,10830,10833,10832,10839,10838,11652,11653,11650,11651,11656,11657,11654,11655,10699,10698,10691,10690,10693,10692,10695,10694,10697,10696,10288,10289,10286,10287,10284,10285,10282,10283,10280,10281,10448,10449,10444,10445,10446,10447,10440,10441,10442,10443,9906,9904,9905,9902,9903,9900,9901,11902,11903,11900,11901,11906,11907,11904,11905,11908,11909,11889,11888,11089,11088,11087,11086,10851,10850,10309,10308,10307,10306,10305,10304,10303,10302,10301,10300,9878,9879,9876,9877,11627,11629,11628,11238,11239,10019,10018,11230,11231,11232,11233,10017,11235,11236,11237,10714,10710,10711,10712,10713,11489,11488,11249,11248,11245,11244,11247,11246,11241,11240,11243,11242,9894,9895,9896,9897,9890,9891,9892,9893,9898,9899,12042,12043,12040,12041,12046,12047,12044,12045,12048,12049,10826,10827,10824,10825,10828,10829,11111,11110,11113,11112,11115,11114,11117,11116,10688,10689,10686,10687,10684,10685,10439,10438,10431,10430,10433,10432,10435,10434,10437,10436,11898,11899,11894,11895,11896,11897,11890,11891,11892,11893,11518,11515,11514,11517,11516,11511,11510,11513,11512,11094,11095,11096,11097,11090,11091,11092,11093,11098,11099,12039,12038,12033,12032,12031,12030,12037,12036,12035,12034,10310,11638,11639,11634,11635,11636,11637,11630,11631,11632,11633,11498,11499,11490,11491,11492,11493,11494,11495,11496,11497,11256,11257,11254,11255,11252,11253,11250,11251,9887,9886,9885,9884,9883,9882,9881,9880,9889,9888,11108,11109,11102,11103,11100,11101,11106,11107,11104,11105,10028,10040,10041,10042,10029,10044,10045,10046,10047,11227,11226,10428,10429,10039,10038,10035,10423,10037,10036,10031,10030,10033,10032,10855,10854,10853,10852,11508,11509,11506,11507,11504,11505,11502,11503,11500,11501,10043,12027,12028,12029,11234,11649,11648,11641,11640,11643,11642,11645,11644,11647,11646',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 16:31:21.700 +34.30.0.0000',
       scan=['193,162,196,165,168,171,174,177,181,184,187,190'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T19t21.J163121+343000_P12482v1_2017_11_10T21_39_09.970/working/VLASS1.1.sb34516741.eb34573529.58031.90910875_split.ms'],
       imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='9876,9877,9878,9879,9880,9881,9882,9883,9884,9885,9886,9887,9888,9889,9890,9891,9892,9893,9894,9895,9896,9897,9898,9899,9900,9901,9902,9903,9904,9905,9906,10017,10018,10019,10020,10021,10022,10023,10024,10025,10026,10027,10028,10029,10030,10031,10032,10033,10034,10035,10036,10037,10038,10039,10040,10041,10042,10043,10044,10045,10046,10047,10280,10281,10282,10283,10284,10285,10286,10287,10288,10289,10290,10291,10292,10293,10294,10295,10296,10297,10298,10299,10300,10301,10302,10303,10304,10305,10306,10307,10308,10309,10310,10421,10422,10423,10424,10425,10426,10427,10428,10429,10430,10431,10432,10433,10434,10435,10436,10437,10438,10439,10440,10441,10442,10443,10444,10445,10446,10447,10448,10449,10450,10451,10684,10685,10686,10687,10688,10689,10690,10691,10692,10693,10694,10695,10696,10697,10698,10699,10700,10701,10702,10703,10704,10705,10706,10707,10708,10709,10710,10711,10712,10713,10714,10824,10825,10826,10827,10828,10829,10830,10831,10832,10833,10834,10835,10836,10837,10838,10839,10840,10841,10842,10843,10844,10845,10846,10847,10848,10849,10850,10851,10852,10853,10854,10855,11086,11087,11088,11089,11090,11091,11092,11093,11094,11095,11096,11097,11098,11099,11100,11101,11102,11103,11104,11105,11106,11107,11108,11109,11110,11111,11112,11113,11114,11115,11116,11117,11226,11227,11228,11229,11230,11231,11232,11233,11234,11235,11236,11237,11238,11239,11240,11241,11242,11243,11244,11245,11246,11247,11248,11249,11250,11251,11252,11253,11254,11255,11256,11257,11488,11489,11490,11491,11492,11493,11494,11495,11496,11497,11498,11499,11500,11501,11502,11503,11504,11505,11506,11507,11508,11509,11510,11511,11512,11513,11514,11515,11516,11517,11518,11627,11628,11629,11630,11631,11632,11633,11634,11635,11636,11637,11638,11639,11640,11641,11642,11643,11644,11645,11646,11647,11648,11649,11650,11651,11652,11653,11654,11655,11656,11657,11888,11889,11890,11891,11892,11893,11894,11895,11896,11897,11898,11899,11900,11901,11902,11903,11904,11905,11906,11907,11908,11909,11910,11911,11912,11913,11914,11915,11916,11917,11918,12027,12028,12029,12030,12031,12032,12033,12034,12035,12036,12037,12038,12039,12040,12041,12042,12043,12044,12045,12046,12047,12048,12049,12050,12051,12052,12053,12054,12055,12056,12057',
       outframe='LSRK', gridder='mosaic', cycleniter=500, niter=0,
       datacolumn='data', savemodel='none', reffreq='3.0GHz', restoration=False,
       intent='OBSERVE_TARGET#UNSPECIFIED', cyclefactor=2.0, robust=1.0,
       parallel=False, stokes='I', nchan=-1, conjbeams=False,
       deconvolver='mtmfs', weighting='briggs', pbcor=False, pblimit=0.2,
       restoringbeam='common', specmode='mfs', chanchunks=-1, interactive=0)
tclean(phasecenter='J2000 16:31:21.700 +34.30.0.0000',
       scan=['193,162,196,165,168,171,174,177,181,184,187,190'],
       spw=['2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17'],
       vis=['/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T19t21.J163121+343000_P12482v1_2017_11_10T21_39_09.970/working/VLASS1.1.sb34516741.eb34573529.58031.90910875_split.ms'],
       imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00050459947354138416, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='9876,9877,9878,9879,9880,9881,9882,9883,9884,9885,9886,9887,9888,9889,9890,9891,9892,9893,9894,9895,9896,9897,9898,9899,9900,9901,9902,9903,9904,9905,9906,10017,10018,10019,10020,10021,10022,10023,10024,10025,10026,10027,10028,10029,10030,10031,10032,10033,10034,10035,10036,10037,10038,10039,10040,10041,10042,10043,10044,10045,10046,10047,10280,10281,10282,10283,10284,10285,10286,10287,10288,10289,10290,10291,10292,10293,10294,10295,10296,10297,10298,10299,10300,10301,10302,10303,10304,10305,10306,10307,10308,10309,10310,10421,10422,10423,10424,10425,10426,10427,10428,10429,10430,10431,10432,10433,10434,10435,10436,10437,10438,10439,10440,10441,10442,10443,10444,10445,10446,10447,10448,10449,10450,10451,10684,10685,10686,10687,10688,10689,10690,10691,10692,10693,10694,10695,10696,10697,10698,10699,10700,10701,10702,10703,10704,10705,10706,10707,10708,10709,10710,10711,10712,10713,10714,10824,10825,10826,10827,10828,10829,10830,10831,10832,10833,10834,10835,10836,10837,10838,10839,10840,10841,10842,10843,10844,10845,10846,10847,10848,10849,10850,10851,10852,10853,10854,10855,11086,11087,11088,11089,11090,11091,11092,11093,11094,11095,11096,11097,11098,11099,11100,11101,11102,11103,11104,11105,11106,11107,11108,11109,11110,11111,11112,11113,11114,11115,11116,11117,11226,11227,11228,11229,11230,11231,11232,11233,11234,11235,11236,11237,11238,11239,11240,11241,11242,11243,11244,11245,11246,11247,11248,11249,11250,11251,11252,11253,11254,11255,11256,11257,11488,11489,11490,11491,11492,11493,11494,11495,11496,11497,11498,11499,11500,11501,11502,11503,11504,11505,11506,11507,11508,11509,11510,11511,11512,11513,11514,11515,11516,11517,11518,11627,11628,11629,11630,11631,11632,11633,11634,11635,11636,11637,11638,11639,11640,11641,11642,11643,11644,11645,11646,11647,11648,11649,11650,11651,11652,11653,11654,11655,11656,11657,11888,11889,11890,11891,11892,11893,11894,11895,11896,11897,11898,11899,11900,11901,11902,11903,11904,11905,11906,11907,11908,11909,11910,11911,11912,11913,11914,11915,11916,11917,11918,12027,12028,12029,12030,12031,12032,12033,12034,12035,12036,12037,12038,12039,12040,12041,12042,12043,12044,12045,12046,12047,12048,12049,12050,12051,12052,12053,12054,12055,12056,12057',
       outframe='LSRK', gridder='mosaic', cycleniter=500, niter=20000,
       datacolumn='data', savemodel='none', reffreq='3.0GHz', restoration=True,
       intent='OBSERVE_TARGET#UNSPECIFIED', cyclefactor=2.0, robust=1.0,
       calcpsf=False, restart=True, stokes='I', nchan=-1, conjbeams=False,
       deconvolver='mtmfs', weighting='briggs', pbcor=False, calcres=False,
       pblimit=0.2, restoringbeam='common', specmode='mfs', chanchunks=-1,
       interactive=0)

# hifv_pbcor(pipelinemode="automatic")
#
impbcor(cutoff=-1.0, stretch=False,
        outfile='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.pb.tt0')

# hif_makermsimages(pipelinemode="automatic")
#
imdev(zscore=-1, ylength='60arcsec', stretch=False, interp='cubic',
      outfile='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.image.pbcor.tt0',
      grid=[10, 10], overwrite=True, maxiter=-1, xlength='60arcsec',
      anchor='ref', stattype='xmadm', statalg='chauvenet')

# hif_makecutoutimages(pipelinemode="automatic")
#
imstat(imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imhead(imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imstat(imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.image.pbcor.tt0')
imhead(imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.image.pbcor.tt0')
imstat(imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.psf.tt0')
imhead(imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.psf.tt0.subim',
           imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.psf.tt0')
imstat(imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.residual.tt0')
imhead(imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.residual.tt0.subim',
           imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.residual.tt0')
imstat(imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.pb.tt0')
imhead(imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.pb.tt0.subim',
           imagename='VLASS1.1.ql.T19t21.J163121+343000.10.2048.v1.I.iter1.pb.tt0')

# hifv_exportvlassdata(pipelinemode="automatic")
#