# 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':
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), })
#!/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), })
#!/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', })