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.sb34451818.eb34599612.58040.92692017361.ms'], session=['session_1'])
#
# If required, ASDMs are converted to measurement sets.
#

# hif_editimlist(parameter_file='parameter.8343.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_T20t24.J184131+373000_P8343v1_2017_11_02T23_34_08.229/working/VLASS1.1.sb34451818.eb34599612.58040.92692017361_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T20t24.J184131+373000_P8343v1_2017_11_02T23_34_08.229/rawdata/VLASS1.1.sb34451818.eb34599612.58040.92692017361.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='18125,16073,16072,16071,16070,16075,16074,16578,16579,16576,16577,16574,16575,16572,16573,16571,18141,18140,17229,17228,17225,17224,17227,17226,17221,17220,17223,17222,17982,17983,17984,17985,17986,17987,17988,17989,16957,16958,16959,16819,16818,16817,16816,16216,16211,16210,17346,17347,17345,16198,16199,17618,17619,16215,16214,16213,16212,17348,17349,17610,17611,17612,17613,17614,17615,17616,17617,16587,16586,16585,16584,16583,16582,16581,16580,16194,16195,16196,16197,16190,16191,16192,16193,16589,16588,18149,18148,16448,16449,16442,16443,16440,16441,16446,16447,16444,16445,17760,17761,17762,17763,17764,17765,16846,16968,16844,16845,16842,16843,16840,16841,16961,16960,16963,16962,16965,16964,16967,16966,16064,16065,16066,16067,16060,16061,16062,16063,16068,16069,17621,17620,17622,17234,17232,17233,17230,17231,18134,18135,18136,16433,16432,16431,16430,16437,16436,16435,16434,16439,16438,18143,18142,17351,18145,18137,18144,18147,18146,18130,17759,17758,17751,17750,17753,17752,17755,17754,17757,17756,18131,16059,16058,18132,16055,16054,16057,16056,16051,16050,16053,16052,18133,17355,17354,17605,16200,16201,16202,17356,16204,16205,16206,16207,17373,17372,17371,17370,17376,17375,17602,16590,16591,16592,16593,16594,17353,16596,16597,16598,16599,17600,18012,18011,18010,18152,18153,18150,18151,18154,18155,16978,16979,16972,16973,16970,16971,16976,16977,16974,16975,16969,17209,17208,17203,17207,17206,17205,17204,16837,16836,16835,16834,16833,16832,16831,16830,18129,18128,16839,16838,16459,16458,16988,18138,18139,16828,17748,17749,16450,17742,17743,17740,17741,17746,17747,17744,17745,17593,17592,16452,17597,17596,17595,17594,16455,17599,17598,16454,16457,16456,16048,16049,16046,16047,16045,16600,16601,16460,18008,18009,18004,18005,18006,18007,18000,18001,18002,18003,17997,17996,17995,17994,17993,17992,17991,17990,17999,17998,17364,17365,17366,17367,17360,17361,17362,17363,17368,17369,16451,16203,17374,16208,16209,16453,17737,17736,17735,17739,17738,16595,17359,17358,17609,17608,17607,17606,17357,17604,17603,17350,17601,17352,16187,16186,16189,16188,16820,16821,16822,16823,16824,16825,16826,16827,16983,16829,16981,16980,16987,16986,16985,16984,16982,17214,17215,17216,17217,17210,17211,17212,17213,17218,17219,18127,18126',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 18:41:31.794 +37.30.0.0000',
       scan=['290,293,296,299,306,303,274,277,280,309,284,287'],
       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_T20t24.J184131+373000_P8343v1_2017_11_02T23_34_08.229/working/VLASS1.1.sb34451818.eb34599612.58040.92692017361_split.ms'],
       imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='16045,16046,16047,16048,16049,16050,16051,16052,16053,16054,16055,16056,16057,16058,16059,16060,16061,16062,16063,16064,16065,16066,16067,16068,16069,16070,16071,16072,16073,16074,16075,16186,16187,16188,16189,16190,16191,16192,16193,16194,16195,16196,16197,16198,16199,16200,16201,16202,16203,16204,16205,16206,16207,16208,16209,16210,16211,16212,16213,16214,16215,16216,16430,16431,16432,16433,16434,16435,16436,16437,16438,16439,16440,16441,16442,16443,16444,16445,16446,16447,16448,16449,16450,16451,16452,16453,16454,16455,16456,16457,16458,16459,16460,16571,16572,16573,16574,16575,16576,16577,16578,16579,16580,16581,16582,16583,16584,16585,16586,16587,16588,16589,16590,16591,16592,16593,16594,16595,16596,16597,16598,16599,16600,16601,16816,16817,16818,16819,16820,16821,16822,16823,16824,16825,16826,16827,16828,16829,16830,16831,16832,16833,16834,16835,16836,16837,16838,16839,16840,16841,16842,16843,16844,16845,16846,16957,16958,16959,16960,16961,16962,16963,16964,16965,16966,16967,16968,16969,16970,16971,16972,16973,16974,16975,16976,16977,16978,16979,16980,16981,16982,16983,16984,16985,16986,16987,16988,17203,17204,17205,17206,17207,17208,17209,17210,17211,17212,17213,17214,17215,17216,17217,17218,17219,17220,17221,17222,17223,17224,17225,17226,17227,17228,17229,17230,17231,17232,17233,17234,17345,17346,17347,17348,17349,17350,17351,17352,17353,17354,17355,17356,17357,17358,17359,17360,17361,17362,17363,17364,17365,17366,17367,17368,17369,17370,17371,17372,17373,17374,17375,17376,17592,17593,17594,17595,17596,17597,17598,17599,17600,17601,17602,17603,17604,17605,17606,17607,17608,17609,17610,17611,17612,17613,17614,17615,17616,17617,17618,17619,17620,17621,17622,17735,17736,17737,17738,17739,17740,17741,17742,17743,17744,17745,17746,17747,17748,17749,17750,17751,17752,17753,17754,17755,17756,17757,17758,17759,17760,17761,17762,17763,17764,17765,17982,17983,17984,17985,17986,17987,17988,17989,17990,17991,17992,17993,17994,17995,17996,17997,17998,17999,18000,18001,18002,18003,18004,18005,18006,18007,18008,18009,18010,18011,18012,18125,18126,18127,18128,18129,18130,18131,18132,18133,18134,18135,18136,18137,18138,18139,18140,18141,18142,18143,18144,18145,18146,18147,18148,18149,18150,18151,18152,18153,18154,18155',
       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 18:41:31.794 +37.30.0.0000',
       scan=['290,293,296,299,306,303,274,277,280,309,284,287'],
       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_T20t24.J184131+373000_P8343v1_2017_11_02T23_34_08.229/working/VLASS1.1.sb34451818.eb34599612.58040.92692017361_split.ms'],
       imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00052253988458505769, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='16045,16046,16047,16048,16049,16050,16051,16052,16053,16054,16055,16056,16057,16058,16059,16060,16061,16062,16063,16064,16065,16066,16067,16068,16069,16070,16071,16072,16073,16074,16075,16186,16187,16188,16189,16190,16191,16192,16193,16194,16195,16196,16197,16198,16199,16200,16201,16202,16203,16204,16205,16206,16207,16208,16209,16210,16211,16212,16213,16214,16215,16216,16430,16431,16432,16433,16434,16435,16436,16437,16438,16439,16440,16441,16442,16443,16444,16445,16446,16447,16448,16449,16450,16451,16452,16453,16454,16455,16456,16457,16458,16459,16460,16571,16572,16573,16574,16575,16576,16577,16578,16579,16580,16581,16582,16583,16584,16585,16586,16587,16588,16589,16590,16591,16592,16593,16594,16595,16596,16597,16598,16599,16600,16601,16816,16817,16818,16819,16820,16821,16822,16823,16824,16825,16826,16827,16828,16829,16830,16831,16832,16833,16834,16835,16836,16837,16838,16839,16840,16841,16842,16843,16844,16845,16846,16957,16958,16959,16960,16961,16962,16963,16964,16965,16966,16967,16968,16969,16970,16971,16972,16973,16974,16975,16976,16977,16978,16979,16980,16981,16982,16983,16984,16985,16986,16987,16988,17203,17204,17205,17206,17207,17208,17209,17210,17211,17212,17213,17214,17215,17216,17217,17218,17219,17220,17221,17222,17223,17224,17225,17226,17227,17228,17229,17230,17231,17232,17233,17234,17345,17346,17347,17348,17349,17350,17351,17352,17353,17354,17355,17356,17357,17358,17359,17360,17361,17362,17363,17364,17365,17366,17367,17368,17369,17370,17371,17372,17373,17374,17375,17376,17592,17593,17594,17595,17596,17597,17598,17599,17600,17601,17602,17603,17604,17605,17606,17607,17608,17609,17610,17611,17612,17613,17614,17615,17616,17617,17618,17619,17620,17621,17622,17735,17736,17737,17738,17739,17740,17741,17742,17743,17744,17745,17746,17747,17748,17749,17750,17751,17752,17753,17754,17755,17756,17757,17758,17759,17760,17761,17762,17763,17764,17765,17982,17983,17984,17985,17986,17987,17988,17989,17990,17991,17992,17993,17994,17995,17996,17997,17998,17999,18000,18001,18002,18003,18004,18005,18006,18007,18008,18009,18010,18011,18012,18125,18126,18127,18128,18129,18130,18131,18132,18133,18134,18135,18136,18137,18138,18139,18140,18141,18142,18143,18144,18145,18146,18147,18148,18149,18150,18151,18152,18153,18154,18155',
       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.T20t24.J184131+373000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.pb.tt0')

# hif_makermsimages(pipelinemode="automatic")
#
imdev(zscore=-1, ylength='60arcsec', stretch=False, interp='cubic',
      outfile='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.image.pbcor.tt0.rms',
      imagename='VLASS1.1.ql.T20t24.J184131+373000.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.T20t24.J184131+373000.10.2048.v1.I.iter1.image.pbcor.tt0')
imhead(imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.image.pbcor.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.image.pbcor.tt0.subim',
           imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.image.pbcor.tt0')
imstat(imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imhead(imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.image.pbcor.tt0.rms.subim',
           imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.image.pbcor.tt0.rms')
imstat(imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.psf.tt0')
imhead(imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.psf.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.psf.tt0.subim',
           imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.psf.tt0')
imstat(imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.residual.tt0')
imhead(imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.residual.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.residual.tt0.subim',
           imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.residual.tt0')
imstat(imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.pb.tt0')
imhead(imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.pb.tt0')
imsubimage(box='1785,1785,5506,5506',
           outfile='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.pb.tt0.subim',
           imagename='VLASS1.1.ql.T20t24.J184131+373000.10.2048.v1.I.iter1.pb.tt0')

# hifv_exportvlassdata(pipelinemode="automatic")
#