casa6

import collections

import glob

import itertools

import numpy

import os

​

try:

    from casatools import msmetadata

    from casatools import table

    from casatools import measures

    from casatools import quanta

    from casatools import ms as mstool

    from casatools import image

    from casatasks import casalog

except Exception:

    from __casac__.msmetadata import msmetadata

    from __casac__.table import table

    from __casac__.measures import measures

    from __casac__.quanta import quanta

    from __casac__.ms import ms as mstool

    from __casac__.image import image

    from taskinit import casalog

​

MetaDataSet = collections.namedtuple(

    'MetaDataSet',

    ['msmeta', 'ephemmeta']

)

​

MsMeta = collections.namedtuple(

    'MsMeta',

    ['positions', 'times', 'freqmin', 'freqmax']

)

​

EphemMeta = collections.namedtuple(

    'EphemMeta',

    ['table', 'times', 'unit_time', 'velocities', 'unit_vel', 'frame_vel']

)

​

FrequencyRange = collections.namedtuple(

    'FrequencyRange',

    ['min', 'max', 'ref']

)

​

COMMON_VELOCITY_UNIT = 'km/s'

​

DEBUG = False

​

​

def debug_print(msg):

    if DEBUG:

        for m in msg.split('\n'):

            casalog.post('DEBUG: {}'.format(m))

​

​

def inspect_ms(vis, fieldid, spwid, chanstart=0, nchan=-1):

    """Inspect MS

​

    Arguments:

        vis {str} -- name of MS

        fieldid {int} -- FIELD_ID for target

        spwid {int} -- SPW_ID for target spw

        chanstart {int} -- start channel (default: 0)

        nchan {int} -- number of channels