Source of task_appendantab.py - casa6

 
  'NAME': 'StandardStMan',
  'SEQNR': 0,
  'SPEC': {'BUCKETSIZE': 1920,
   'IndexLength': 142,
   'MaxCacheSize': 2,
   'PERSCACHESIZE': 2},
  'TYPE': 'StandardStMan'}}
​
#################################################################################################
# Main task
def appendantab(vis=None, outvis=None, antab=None, overwrite=False):
def appendantab(vis=None, outvis=None, antab=None, 
                overwrite=False, append_tsys=True, append_gc=True):
​
    # get info from vis
    msmd.open(vis)
    ant_names = msmd.antennastations()
    n_band = len(msmd.bandwidths())
    spws = msmd.spwfordatadesc()
​
    # is gc_begin even needed?
    gc_begin = msmd.timerangeforobs(0)['begin']['m0']['value'] * 86400
    gc_end = msmd.timerangeforobs(0)['end']['m0']['value'] * 86400
Show more
        tm_hour = int(ts.split(':')[0])
        tm_min = int(ts.split(':')[1])
        tm_sec = int(ts.split(':')[2])
    elif re.match(r"[0-9]{2}:[0-9]{2}:[0-9]{2}\.[0-9]+", ts):
        # hh:mm:ss.ss
        tm_hour = int(ts.split(':')[0])
        tm_min = int(ts.split(':')[1])
        tm_sec = float(ts.split(':')[2])
    return tm_hour, tm_min, tm_sec
​
def process(infp, keys, pols, data, n_band, spws, first_time, last_time, ant_names):
def process_tsys(infp, keys, pols, data, n_band, spws, first_time, last_time, ant_names):
    # Get the antenna name
    antenna_name = find_antenna(keys[0], ['SRC/SYS'])
    # Skip if no antenna name found
    if not antenna_name:
        print('ANTENNA missing from TSYS group')
        skip_values(infp)
        return
    try:
        antenna = get_antenna_index(antenna_name, ant_names)
    except:
Show more
        if line.startswith('!'):
            continue
        keys.write(line)
        if line.strip().endswith('/'):
            keys.seek(0)
            try:
                tsys = read_keyfile(keys)
            except RuntimeError:
                pass
            if tsys and tsys[0] and tsys[0][0][0] == 'TSYS':
                process(fp, tsys, pols, data, n_band, spws, first_time, last_time, ant_names)
                process_tsys(fp, tsys, pols, data, n_band, spws, first_time, last_time, ant_names)
                pass
            elif tsys and tsys[0] and tsys[0][0][0] == 'GAIN':
                process_gc(fp, tsys, pols, data_gc, ant_names)
            keys = StringIO()
            continue
        continue
    
    # Create the subtable for syscal
    create_syscal_subtable(data, outvis, spws)
    create_gaincurve_subtable(data_gc, outvis, spws, gc_interval, gc_time)
    # Create the subtable for syscal if selected in params
    if append_tsys: create_syscal_subtable(data, outvis, spws)
    if append_gc: create_gaincurve_subtable(data_gc, outvis, spws, gc_interval, gc_time)
​
def create_gaincurve_subtable(data, outvis, spws, gc_interval, gc_time):
    # Create the subtable for gain curve
    tb.create(f'{outvis}/GAIN_CURVE', desc_gc, dminfo=dminfo_gc)
    tb.putkeyword('GAIN_CURVE', f'Table: {outvis}/GAIN_CURVE')
​
    # assemble columns
    ANTENNA_ID = []
    FEED_ID = []
    SPECTRAL_WINDOW_ID = []
casa6

Source

Add shortcut
