Usage for plotWVRSolutions
Return to directory of Todd's CASA extensions
Although plotcal can be used to plot WVR solutions on an antenna-basis, it is not very insightful because it includes the common mode changes in delay which makes can be large enough that it becomes hard to see temporal differences between antennas. The function plotWVRSolutions takes pairwise difference between antennas and produces plots of the effective phase correction on a baseline basis. The plots can be sorted by antenna number pair or by baseline length. Fields are plotted in different colors. All spws are displayed sequentially, so be sure to limit them to (one of) the science spws, because they should only differ by the frequency ratio. The antenna parameter accepts the standard '!' modifier, but not '&'.
CASA <2>: au.plotWVRSolutions(help=True)
Usage: plotWVRSolutions(caltable='', spw=, field=, subplot=22, sort='number',
antenna='', xrange=[0,0], yrange=[0,0], figfile='', xaxis='ut',
interactive=True, ms='', showFlaggedSolutions=False, usems=True,
- antenna: can be ID or name; single value, or list: 0,'0','0,1','DV04,DV05'. All baselines including this antenna will be plotted.
- field: must be field ID, or list of IDs, or comma-delimited string of IDs
- figfile: will append .0.png, .1.png, etc. for successive pages
- interactive: True or False (True will require user input after showing each page)
- ms: specify the ms name rather than look for it in the caltable
- showFlaggedSolutions: Boolean, if False, then don't show points for which the solution is flagged on either antenna
- sort: by baseline 'number', or by baseline 'length'
- subplot: valid values: 11,12,21,22,23,32,42 (any 3rd digit is ignored)
- unwrap: attempt to further unwrap the 360-deg phase wraps after Ed's algorithm
- usems: set to False if you don't want to risk disturbing someone else's use of the measurement set by causing a write-lock on the ms
- xaxis: 'seconds' or 'ut' (default)
- xrange: min/max time to show: specify in seconds if xaxis='seconds' or floating point UT hours if 'ut'. The latter will probably not work well if the data cross the 0/24 hour mark.
- yrange: min/max phase to show
CASA <6>: au.plotWVRSolutions('X279.wvr',spw='19',yrange=[-300,300],figfile=T)