Ejemplo n.º 1
0
# from miriad import uvspec, imspec
import miriad

if __name__ == '__main__':
    for sb in ['usb', 'lsb']:
        miriad.uvspec({
            'vis':
            'UVDATA/irc+10216_{0}.vis.uvcal,UVoffsetCORRECT/irc+10216.{0}.corrected.slfc'
            .format(sb),
            'device':
            '{0}.ps/cps'.format(sb),
            'interval':
            9999,
            'options':
            'avall,nobase',
            'nxy':
            '1,2',
            'stokes':
            'v',
            'axis':
            'freq,amp',
            'line':
            'chan,96,1,64.0',
        })

    lines = ['cnt', 'co3-2', 'cs7-6', 'h13cn4-3', 'sis19-18']

    for line in lines:
        miriad.uvspec({
            'vis':
Ejemplo n.º 2
0

def plotPhase(vis, device='2/xs'):
    miriad.smauvplt({
        'vis': vis,
        'device': device,
        'interval': '1e3',
        'stokes': 'i,v',
        'axis': 'time,pha',
        'nxy': '2,2'
    })


# vis = 'UVDATA/NGC7538S-s4.usb'
vis = 'UVOffsetCorrect/NGC7538S-s4.usb.corrected.slfc'
numChannels = miriad.getNumChannels(vis)
velrange = miriad.getVelocityRange(vis)

miriad.uvspec({
    'vis':
    'UVDATA/NGC7538S-s4.usb,UVOffsetCorrect/NGC7538S-s4.usb.corrected.slfc',
    'device': '1/xs',
    'interval': 9999,
    'options': 'avall,nobase',
    'nxy': '1,2',
    'stokes': 'v',
    'axis': 'freq,phase',
    # 'yrange': '0,0.5',
    'line': miriad.averageVelocityLine(vis, factor=20),
})
Ejemplo n.º 3
0
#!/usr/bin/python3

import miriad


def plotPhase(vis, device='2/xs'):
    miriad.smauvplt({
        'vis': vis,
        'device': device,
        'interval': '1e3',
        'stokes': 'i,v',
        'axis': 'time,pha',
        'nxy': '2,2'
    })


vis = 'UVDATA/HH211.usb'

miriad.uvspec({
    'vis': vis,
    'device': '1/xs',
    'interval': 9999,
    'options': 'avall,nobase',
    'nxy': '1,1',
    'stokes': 'i',
    'axis': 'freq,amp',
    # 'line': miriad.averageVelocityLine(vis, factor=1),
})
Ejemplo n.º 4
0
#!/usr/bin/python3

import miriad

if __name__ == '__main__':
    miriad.showChannels('tmp.3')
    # exit()
    miriad.uvspec({
        'vis':
        'UVDATA/orkl_080106.usb,UVOffsetCorrect/orkl_080106.usb.corrected.slfc',
        'device': '2/xs',
        'interval': 9999,
        'options': 'avall,nobase',
        'nxy': '1,2',
        'stokes': 'v',
        'axis': 'chan,amp',
        'line': 'chan,48,1,64.0',
    })