casa6

import os

import numpy

from taskinit import *

​

def wvrgcal(vis=None, caltable=None, toffset=None, segsource=None,

        sourceflag=None, tie=None, nsol=None, disperse=None, 

        wvrflag=None, statfield=None, statsource=None, smooth=None,

        scale=None, spw=None, wvrspw=None,

        reversespw=None,  cont=None, maxdistm=None,

        minnumants=None, mingoodfrac=None, usefieldtab=None, 

        refant=None, offsetstable=None):

    """

    Generate a gain table based on Water Vapour Radiometer data.

    Returns a dictionary containing the RMS of the path length variation

    with time towards that antenna (RMS) and the discrepency between the RMS

    path length calculated separately for different WVR channels (Disc.).

​

      vis -- Name of input visibility file                                          

                  default: none; example: vis='ngc5921.ms'                                  

      caltable -- Name of output gain calibration table                                 

                  default: none; example: caltable='ngc5921.wvr'                                

      toffset -- Time offset (sec) between interferometric and WVR data                         

                 default: 0 (ALMA default for cycle 1, for cycle 0 it was -1)                                           

      segsource -- Do a new coefficient calculation for each source                             

                 default: True                                              

      tie -- Prioritise tieing the phase of these sources as well as possible                       

             (requires segsource=True)                                          

                 default: [] example: ['3C273,NGC253', 'IC433,3C279']                           

      sourceflag -- Flag the WVR data for these source(s) as bad and do not produce corrections for it          

                   (requires segsource=True)                                        

                   default: [] (none) example: ['3C273']                                    

      nsol -- Number of solutions for phase correction coefficients during this observation.                

              By default only one set of coefficients is generated for the entire observation.              

              If more sets are requested, then they will be evenly distributed in time throughout               

              the observation. Values > 1 require segsource=False.                              

                 default: 1                                                 

      disperse -- Apply correction for dispersion                                       

                 default: False                                             

      wvrflag -- Regard the WVR data for these antenna(s) as bad and use interpolated values instead        

                   default: [] (none) example: ['DV03','DA05','PM02']                                   

      statfield -- Compute the statistics (Phase RMS, Disc) on this field only                      

                   default: '' (all)                                            

      statsource -- Compute the statistics (Phase RMS, Disc) on this source only                        

                   default: '' (all)                                                    

      smooth -- Smooth the calibration solution on the given timescale 

                 default: '' (no smoothing), example: '3s' smooth on a timescale of 3 seconds

      scale -- Scale the entire phase correction by this factor                             

                 default: 1. (no scaling)                                           

​

          spw -- List of the spectral window IDs for which solutions should be saved into the caltable

                 default: [] (all spectral windows), example [17,19,21,23]

​

          wvrspw -- List of the spectral window IDs from which the WVR data should be taken

                 default: [] (all WVR spectral windows), example [0]

      reversespw -- Reverse the sign of the correction for the listed SPWs                          

                    (only needed for early ALMA data before Cycle 0)                                

                 default: '' (none), example: reversespw='0~2,4'; spectral windows 0,1,2,4                  

      cont -- Estimate the continuum (e.g., due to clouds)                                  

                     default: False

​

          maxdistm -- maximum distance (m) an antenna may have to be considered for being part

                  of the <=3 antenna set for interpolation of a solution for a flagged antenna

              default: 500

​

          minnumants -- minimum number of near antennas required for interpolation

                    default: 2

​

          mingoodfrac -- If the fraction of unflagged data for an antenna is below this value (0. to 1.),

                     the antenna is flagged.

             default: 0.8

​

          usefieldtab -- derive the antenna AZ/EL values from the FIELD rather than the POINTING table

                     default: False

​

      refant -- use the WVR data from this antenna for calculating the dT/dL parameters (can give ranked list)

                  default: '' (use the first good or interpolatable antenna),