'out': 'tmp.2',
'hdvar': 'restfreq',
'varval': freq
})
miriad.uvredo({'vis': 'tmp.2', 'out': 'tmp.3', 'options': 'velocity'})
# this is because we average 5 velocity channels into 1
miriad.uvflag({'vis': 'tmp.3', 'flagval': 'f', 'edge': '5,5'})
miriad.uvaver({'vis': 'tmp.3', 'out': 'tmp.4', 'line': ch})
miriad.uvlist({'vis': 'tmp.3', 'options': 'spec'})
miriad.uvlist({'vis': 'tmp.4', 'options': 'spec'})
miriad.smauvspec({
'vis': 'tmp.4',
'device': '1/xw',
'interval': '1e3',
'stokes': 'i',
'axis': 'freq,amp',
'nxy': '2,3'
})
input("Return to continue")
miriad.showChannels('tmp.3', subtitle='SplitCont')
miriad.uvlin({
'vis': 'tmp.3',
'out': '{}.{}'.format(vis, lab),
'chans': free,
'mode': 'chan0', # this says 'store a single average value only'
'order': 0
})
miriad.smauvplt({
miriad.uvputhd({
'vis': 'tmp.1',
'out': 'tmp.2',
'hdvar': 'restfreq',
'varval': freq
})
miriad.uvredo({'vis': 'tmp.2', 'out': 'tmp.3', 'options': 'velocity'})
miriad.uvlist({'vis': 'tmp.3', 'options': 'spec'})
miriad.smauvspec({
'vis': 'tmp.3',
'device': '1/xs',
'interval': '1e3',
'stokes': 'i',
'axis': 'chan,amp',
'nxy': '2,3'
})
input("Return to continue")
miriad.uvlin({
'vis': 'tmp.3',
'out': '{}.{}'.format(vis, lab),
'chans': free,
'mode': 'line',
'order': 0
})
for stk in ['i', 'v']: