casa6

##########################################################################

# test_task_wvrgcal.py

#

# Copyright (C) 2018, 2024

# Associated Universities, Inc. Washington DC, USA.

# European Sourthern Observatory, Garching, Germany

#

# This script is free software; you can redistribute it and/or modify it

# under the terms of the GNU Library General Public License as published by

# the Free Software Foundation; either version 2 of the License, or (at your

# option) any later version.

#

# This library is distributed in the hope that it will be useful, but WITHOUT

# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or

# FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Library General Public

# License for more details.

#

# [Add the link to the JIRA ticket here once it exists]

#

#

##########################################################################

import os

import sys

import shutil

import numpy as np

import matplotlib.pyplot as plt

import math

​

from casatools import ctsys, table

from casatasks import flagdata, smoothcal, split

from casatasks import wvrgcal

​

import unittest

​

from casatestutils import testhelper as th

​

tb = table()

​

def comp_gainphase(oldctab, newctab, myant=0, spw=None, tol_deg=0.2, figfile='', symsize=12, x_is_time=True):

    """

    Test the phaseangle dfference of the two gaintables oldctab and newctab for antenna myant

    together with their difference and write plot to "figfile".png .

​

    If myant==-1, check all antennas.

​

    tol_deg is the max permitted difference in phase angle (degrees) for return value True.

​

    If figfile is not an empty string, dignostic plots of phase angle difference between

    the result and the reference are produced with the filenames starting with the value of figfile.

​

    symsize specifies the plot symbol size in pixels. The minimum size is limited to 5 pixels.

​

    x_is_time specifies whether the plot x-axis is in units of time or table rows.

​

    """

​

    rval = True

    mytb = table()

    mytb.open(oldctab)

    old = mytb.getcol('CPARAM')

    ants = mytb.getcol('ANTENNA1')