Source code for suncasa.eovsa.msUtils

from ..casa_compat import import_casatools

[docs]
tools = import_casatools(['tbtool'])



[docs]
tbtool = tools['tbtool']


[docs]
tb = tbtool()




[docs]
def getAntennaPosition(vis):
    tb.open(vis + '/ANTENNA')
    position = tb.getcol('POSITION')
    diameter = tb.getcol('DISH_DIAMETER')
    antenna = tb.getcol('NAME')
    tb.close()
    return position, diameter, antenna



[docs]
def getObservatoryName(ms):
    """
    Returns the observatory name in the specified ms, using the tb tool.
    -- Todd Hunter
    """
    antTable = ms + '/OBSERVATION'
    try:
        tb.open(antTable)
        myName = tb.getcell('TELESCOPE_NAME')
        tb.close()
    except:
        print("Could not open OBSERVATION table to get the telescope name: {}".format(antTable))
        myName = ''
    return (myName)




[docs]
def buildConfigurationFile(vis='', cfgfile=None):
    observatory = getObservatoryName(vis)
    position, diameter, antenna = getAntennaPosition(vis)
    if cfgfile == '':
        cfgfile = vis + '.cfg'
    cfg = open(cfgfile, 'w')
    coordsys = '# coordsys=XYZ\n'
    cfg.write('# observatory={}\n'.format(observatory))
    cfg.write(coordsys)
    for s in range(len(antenna)):
        x, y, z = position[:, s]
        line = '{:.6f} {:.6f} {:.6f} {:5.1f}  {}\n'.format(x, y, z, diameter[s], antenna[s])
        cfg.write(line)
    cfg.close()
    return cfgfile