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

# hif_editimlist(parameter_file='parameter.13879.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_T24t16.J153544+523000_P13879v1_2018_01_18T22_25_04.462/working/VLASS1.1.sb34539903.eb34684984.58068.81223880787_split.ms',
            vis='/lustre/aoc/cluster/pipeline/vlass_auto/spool/quicklook/VLASS1.1_T24t16.J153544+523000_P13879v1_2018_01_18T22_25_04.462/rawdata/VLASS1.1.sb34539903.eb34684984.58068.81223880787.ms',
            datacolumn='corrected', reindex=False,
            spw='10,11,12,13,14,15,16,17,2,3,4,5,6,7,8,9',
            field='270,271,272,273,274,275,276,277,278,279,1018,1019,1014,1015,1016,1017,1010,1011,1012,1013,1234,1235,1236,1237,1230,1231,1232,1233,1238,1239,13860,108,109,102,103,100,101,106,107,104,105,641,640,1243,642,645,644,1247,1246,649,648,854,855,856,857,850,851,852,853,858,859,13843,1759,1758,1757,1756,1755,1754,1753,1752,1751,1750,13834,13835,13836,13837,13830,13831,13832,13833,13838,13839,873,130,498,499,494,495,496,497,490,491,492,493,1371,1370,1373,1372,1375,1374,1377,1376,1379,1378,13689,13688,13683,13682,13681,13680,13687,13686,13685,13684,1009,1008,1007,1006,1005,1004,1003,1002,1001,1000,1227,1226,1225,1224,1223,626,1221,1220,629,628,1229,1228,656,654,655,652,653,650,651,869,868,861,860,863,862,865,864,867,866,1601,1600,1603,1602,13706,13707,1607,1606,1616,1614,1615,1612,1613,1610,1611,1222,1762,1763,1760,1761,1764,1765,13841,13840,997,13842,13845,13844,13847,13846,13849,13848,999,998,120,121,122,123,124,125,126,127,128,129,1388,1389,1384,1385,1386,1387,1380,1381,1382,1383,13694,13695,13696,13697,13690,13691,13692,13693,13698,13699,258,259,252,253,256,257,254,255,1588,1589,1586,1587,504,505,502,503,500,501,630,631,632,633,634,635,636,637,638,639,1218,1219,1241,1240,876,877,874,875,872,643,870,871,1242,1245,1244,647,646,1248,1609,1608,13702,13703,13700,13701,1605,1604,13704,13705,281,280,282,263,262,261,260,267,266,265,264,269,268,13852,13853,13850,13851,13856,13857,13854,13855,13858,13859,115,114,117,116,111,110,113,112,119,118,1397,1396,1395,1394,1393,1392,1391,1390,847,849,848,489,1739,1738,13678,13679,1735,13677,1737,1736,1595,1594,1597,1596,1591,1590,1593,1592,1599,1598,1025,1024,1027,1026,1021,1020,1023,1022,1028,1748,1749,1744,1745,1746,1747,1740,1741,1742,1743,1217,627,488,487,486,485,484,483,482,481,480,1367,476,477,474,475,478,479,1368,1369',
            chanbin=1)

# hif_makeimages(hm_masking='none', hm_cleaning='manual')
#
# A list of target sources is cleaned.
#
tclean(phasecenter='J2000 15:35:44.035 +52.30.0.0000',
       scan=['34,37,9,12,242,15,18,21,24,239,28,31'],
       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_T24t16.J153544+523000_P13879v1_2018_01_18T22_25_04.462/working/VLASS1.1.sb34539903.eb34684984.58068.81223880787_split.ms'],
       imagename='VLASS1.1.ql.T24t16.J153544+523000.10.2048.v1.I.iter0',
       nterms=2, threshold='0.0mJy', imsize=[7290, 7290], npixels=0, scales=[0],
       cell='1.0arcsec',
       field='100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126,127,128,129,130,252,253,254,255,256,257,258,259,260,261,262,263,264,265,266,267,268,269,270,271,272,273,274,275,276,277,278,279,280,281,282,474,475,476,477,478,479,480,481,482,483,484,485,486,487,488,489,490,491,492,493,494,495,496,497,498,499,500,501,502,503,504,505,626,627,628,629,630,631,632,633,634,635,636,637,638,639,640,641,642,643,644,645,646,647,648,649,650,651,652,653,654,655,656,847,848,849,850,851,852,853,854,855,856,857,858,859,860,861,862,863,864,865,866,867,868,869,870,871,872,873,874,875,876,877,997,998,999,1000,1001,1002,1003,1004,1005,1006,1007,1008,1009,1010,1011,1012,1013,1014,1015,1016,1017,1018,1019,1020,1021,1022,1023,1024,1025,1026,1027,1028,1217,1218,1219,1220,1221,1222,1223,1224,1225,1226,1227,1228,1229,1230,1231,1232,1233,1234,1235,1236,1237,1238,1239,1240,1241,1242,1243,1244,1245,1246,1247,1248,1367,1368,1369,1370,1371,1372,1373,1374,1375,1376,1377,1378,1379,1380,1381,1382,1383,1384,1385,1386,1387,1388,1389,1390,1391,1392,1393,1394,1395,1396,1397,1586,1587,1588,1589,1590,1591,1592,1593,1594,1595,1596,1597,1598,1599,1600,1601,1602,1603,1604,1605,1606,1607,1608,1609,1610,1611,1612,1613,1614,1615,1616,1735,1736,1737,1738,1739,1740,1741,1742,1743,1744,1745,1746,1747,1748,1749,1750,1751,1752,1753,1754,1755,1756,1757,1758,1759,1760,1761,1762,1763,1764,1765,13677,13678,13679,13680,13681,13682,13683,13684,13685,13686,13687,13688,13689,13690,13691,13692,13693,13694,13695,13696,13697,13698,13699,13700,13701,13702,13703,13704,13705,13706,13707,13830,13831,13832,13833,13834,13835,13836,13837,13838,13839,13840,13841,13842,13843,13844,13845,13846,13847,13848,13849,13850,13851,13852,13853,13854,13855,13856,13857,13858,13859,13860',
       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 15:35:44.035 +52.30.0.0000',
       scan=['34,37,9,12,242,15,18,21,24,239,28,31'],
       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_T24t16.J153544+523000_P13879v1_2018_01_18T22_25_04.462/working/VLASS1.1.sb34539903.eb34684984.58068.81223880787_split.ms'],
       imagename='VLASS1.1.ql.T24t16.J153544+523000.10.2048.v1.I.iter1',
       nterms=2, threshold=0.00048091168110840952, imsize=[7290, 7290],
       npixels=0, scales=[0], parallel=False, cell='1.0arcsec',
       field='100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126,127,128,129,130,252,253,254,255,256,257,258,259,260,261,262,263,264,265,266,267,268,269,270,271,272,273,274,275,276,277,278,279,280,281,282,474,475,476,477,478,479,480,481,482,483,484,485,486,487,488,489,490,491,492,493,494,495,496,497,498,499,500,501,502,503,504,505,626,627,628,629,630,631,632,633,634,635,636,637,638,639,640,641,642,643,644,645,646,647,648,649,650,651,652,653,654,655,656,847,848,849,850,851,852,853,854,855,856,857,858,859,860,861,862,863,864,865,866,867,868,869,870,871,872,873,874,875,876,877,997,998,999,1000,1001,1002,1003,1004,1005,1006,1007,1008,1009,1010,1011,1012,1013,1014,1015,1016,1017,1018,1019,1020,1021,1022,1023,1024,1025,1026,1027,1028,1217,1218,1219,1220,1221,1222,1223,1224,1225,1226,1227,1228,1229,1230,1231,1232,1233,1234,1235,1236,1237,1238,1239,1240,1241,1242,1243,1244,1245,1246,1247,1248,1367,1368,1369,1370,1371,1372,1373,1374,1375,1376,1377,1378,1379,1380,1381,1382,1383,1384,1385,1386,1387,1388,1389,1390,1391,1392,1393,1394,1395,1396,1397,1586,1587,1588,1589,1590,1591,1592,1593,1594,1595,1596,1597,1598,1599,1600,1601,1602,1603,1604,1605,1606,1607,1608,1609,1610,1611,1612,1613,1614,1615,1616,1735,1736,1737,1738,1739,1740,1741,1742,1743,1744,1745,1746,1747,1748,1749,1750,1751,1752,1753,1754,1755,1756,1757,1758,1759,1760,1761,1762,1763,1764,1765,13677,13678,13679,13680,13681,13682,13683,13684,13685,13686,13687,13688,13689,13690,13691,13692,13693,13694,13695,13696,13697,13698,13699,13700,13701,13702,13703,13704,13705,13706,13707,13830,13831,13832,13833,13834,13835,13836,13837,13838,13839,13840,13841,13842,13843,13844,13845,13846,13847,13848,13849,13850,13851,13852,13853,13854,13855,13856,13857,13858,13859,13860',
       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.T24t16.J153544+523000.10.2048.v1.I.iter1.image.pbcor.tt0',
        imagename='VLASS1.1.ql.T24t16.J153544+523000.10.2048.v1.I.iter1.image.tt0',
        mode='divide',
        pbimage='VLASS1.1.ql.T24t16.J153544+523000.10.2048.v1.I.iter1.pb.tt0')

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

# hifv_exportvlassdata(pipelinemode="automatic")
#