; RIRMS ;--------------------------------------------------------------- ;! computes rms of real/imag of a selected subset of a uv data set ;# Task UV PLOT ;----------------------------------------------------------------------- ;; Copyright (C) 2020 ;; Associated Universities, Inc. Washington DC, USA. ;; ;; This program is free software; you can redistribute it and/or ;; modify it under the terms of the GNU General Public License as ;; published by the Free Software Foundation; either version 2 of ;; the License, or (at your option) any later version. ;; ;; This program 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 General Public License for more details. ;; ;; You should have received a copy of the GNU General Public ;; License along with this program; if not, write to the Free ;; Software Foundation, Inc., 675 Massachusetts Ave, Cambridge, ;; MA 02139, USA. ;; ;; Correspondence concerning AIPS should be addressed as follows: ;; Internet email: aipsmail@nrao.edu. ;; Postal address: AIPS Project Office ;; National Radio Astronomy Observatory ;; 520 Edgemont Road ;; Charlottesville, VA 22903-2475 USA ;----------------------------------------------------------------------- UVRMS LLLLLLLLLLLLUUUUUUUUUUUU CCCCCCCCCCCCCCCCCCCCCCCCCCCCC UVRMS Computes statistics of a subset of a uv data set INNAME Input UV file name (name) INCLASS Input UV file name (class) INSEQ 0.0 9999.0 Input UV file name (seq. #) INDISK 0.0 9.0 Input UV file disk unit # SOURCES Source name QUAL -10.0 Calibrator qualifier -1=>all CALCODE Calibrator code ' '=>all STOKES Stokes of output TIMERANG Time range to use SELBAND Bandwidth to select (kHz) SELFREQ Frequency to select (MHz) FREQID Freq. ID to select. ANTENNAS Antennas to plot 0=>all BASELINE Baselines with ANTENNAS UVRANGE UV range in kilolambda. SUBARRAY 0.0 1000.0 Sub-array, 0=>all BIF Low IF number included EIF Highest IF number included BCHAN 0.0 First channel included ECHAN 0.0 last channel included DOCALIB -1.0 101.0 > 0 calibrate data & weights > 99 do NOT calibrate weights GAINUSE CL (or SN) table to apply DOPOL -1.0 10.0 If >0.5 correct polarization. PDVER PD table to apply (DOPOL>0) BLVER BL table to apply. FLAGVER Flag table version DOBAND -1.0 10.0 If >0.5 apply bandpass cal. Method used depends on value of DOBAND (see HELP file). BPVER Bandpass table version SMOOTH Smoothing function. See HELP SMOOTH for details. SOLINT Averaging time (min) NDIG 5.0 Number digits in matrix print. 0 -> fit >= 2 DOCRT -3.0 132.0 = 0 -> no matrix print > 0 -> use the terminal, else use the line printer > 72 => terminal width OUTPRINT Printer disk file to save DOPLOT 7.0 1,3,5,7 => plot the histogram of RMSes 2,3,6,7 => plot vs time 4,5,6,7 => plot the histogram of MEANs FACTOR Scale plotted points NBOXES Number boxes in histogram CUTOFF Histogram +- CUTOFF*rms FUNCTYPE 'LG' use LOG10 in hist. plot APARM (1,2) Plot range of time plot DOTV > 0 -> plot on TV, else make a plot file GRCHAN Graphics channel to use LTYPE -410.0 410.0 Type of labeling: 1 border, 2 no ticks, 3 - 6 standard, 7 - 10 only tick labels <0 -> no date/time BADDISK Disks to avoid for scratch ---------------------------------------------------------------- UVRMS Task: This task applies calibration to a selected subset of a uv data file and then averages all included data over the included sectral channels and IFs and the stated averaging time. The resulting list of values for mean and rms as a function of time and baseline are analyzed for there mean and rms. These are then printed in a matrix form. Following that plots of these values as a functions of time and as histograms may be generated. Note that this task depends on having a number of spectral channels and IFs and times in each group in order to form a reasonable rms and mean. Adverbs: INNAME.....Input UV file name (name). Standard defaults. INCLASS....Input UV file name (class). Standard defaults. INSEQ......Input UV file name (seq. #). 0 => highest. INDISK.....Disk drive # of input UV file. 0 => any. SOURCES....Source to be included. ' '=> all; if any starts with a '-' then all except ANY source named. QUAL.......Qualifier of source to be included. -1 => all. CALCODE....Calibrator code of sources to include. ' '=> all. STOKES.....Specifies which 2 STOKES parameters are included. Values allowed are 'HALF', 'RLLR', and 'XYYX' TIMERANG...Time range of the data to be included. In order: Start day, hour, min. sec, end day, hour, min. sec. Days relative to ref. date. SELBAND....Bandwidth of data to be selected. If more than one IF is present SELBAND is the width of the first IF required. Units = kHz. For data which contain multiple bandwidths/frequencies the task will insist that some form of selection be made by frequency or bandwidth. SELFREQ....Frequency of data to be selected. If more than one IF is present SELFREQ is the frequency of the first IF required. Units = MHz. FREQID.....Frequency identifier to select (you may determine which is applicable from the OPTYPE='SCAN' listing produced by LISTR). If either SELBAND or SELFREQ are set, their values override that of FREQID. However, setting SELBAND and SELFREQ may result in an ambiguity. In that case, the task will request that you use FREQID. ANTENNAS...A list of the antennas to be plotted. All 0 => all. Note that the values are computed for all baselines. ANTENNAS only controls which baselines are plotted. If any number is negative then all antennas listed are NOT iplotted and all others are. BASELINE...Baselines between antennas named in ANTENNAS and those named in BASELINE are selected.. There are four possible combinations of ANTENNAS and BASELINE: 1. ANTENNAS = 0; BASELINE = 0. All possible baselines are selected. 2. ANTENNAS <>0; BASELINE = 0. a)All ANTENNAS > 0 Baselines including an antenna in the ANTENNAS list are selected; b)Some ANTENNAS < 0 All baselines NOT including an antenna in the ANTENNAS list are selected; 3. ANTENNAS = 0; BASELINE <> 0. a)All BASELINE > 0 Baselines including an antenna in the BASELINE list are selected; b)Some BASELINE < 0 All baselines NOT including an antenna in the BASELINE list are selected; 4. ANTENNAS <> 0; BASELINE <> 0. a)All ANTENNAS>0 and all BASELINE>0 Baselines between antennas named in ANTENNAS and those named in BASELINE are selected. b)Some ANTENNAS<0 .OR. some BASELINE<0 Baselines between antennas named in ANTENNAS and those named in BASELINE are DE-selected, all others are selected. UVRANGE....Range of projected spacings to be plotted in 1000's of wavelengths. 0 => 1, 1.E10 SUBARRAY...Sub-array number to include. 0=>all. BIF........First IF to include. 0 -> 1. EIF........Last IF to include. 0 -> max. BCHAN......First channel to include. 0=>all. ECHAN......Highest channel to include. 0=>all higher than BCHAN DOCALIB....If true (>0), calibrate the data using information in the specified Cal (CL) table for multi-source or SN table for single-source data. Also calibrate the weights unless DOCALIB > 99 (use this for old non-physical weights). GAINUSE....version number of the CL table to apply to multi-source files or the SN table for single source files. 0 => highest. DOPOL......If > 0 then correct data for instrumental polarization as represented in the AN or PD table. This correction is only useful if PCAL has been run or feed polarization parameters have been otherwise obtained. See HELP DOPOL for available correction modes: 1 is normal, 2 and 3 are for VLBI. 1-3 use a PD table if available; 6, 7, 8 are the same but use the AN (continuum solution) even if a PD table is present. PDVER......PD table to apply if PCAL was run with SPECTRAL true and 0 < DOPOL < 6. <= 0 => highest. BLVER......Version number of the baseline based calibration (BL) table to apply. <0 => apply no BL table, 0 => highest. FLAGVER....specifies the version of the flagging table to be applied. 0 => highest numbered table. <0 => no flagging to be applied. DOBAND.....If true (>0) then correct the data for the shape of the antenna bandpasses using the BP table specified by BPVER. The correction has five modes: (a) if DOBAND=1 all entries for an antenna in the table are averaged together before correcting the data. (b) if DOBAND=2 the entry nearest in time (including solution weights) is used to correct the data. (c) if DOBAND=3 the table entries are interpolated in time (using solution weights) and the data are then corrected. (d) if DOBAND=4 the entry nearest in time (ignoring solution weights) is used to correct the data. (e) if DOBAND=5 the table entries are interpolated in time (ignoring solution weights) and the data are then corrected. IMAGR uses DOBAND as the nearest integer; 0.1 is therefore "false". BPVER......Specifies the version of the BP table to be applied 0 => highest numbered table. <0 => no bandpass correction to be applied. SMOOTH.....Specifies the type of spectral smoothing to be applied to a uv database . The default is not to apply any smoothing. The elements of SMOOTH are as follows: SMOOTH(1) = type of smoothing to apply: 0 => no smoothing To smooth before applying bandpass calibration 1 => Hanning, 2 => Gaussian, 3 => Boxcar, 4 => Sinc To smooth after applying bandpass calibration 5 => Hanning, 6 => Gaussian, 7 => Boxcar, 8 => Sinc SMOOTH(2) = the "diameter" of the function, i.e. width between first nulls of Hanning triangle and sinc function, FWHM of Gaussian, width of Boxcar. Defaults (if < 0.1) are 4, 2, 2 and 3 channels for SMOOTH(1) = 1 - 4 and 5 - 8, resp. SMOOTH(3) = the diameter over which the convolving function has value - in channels. Defaults: 1,3,1,4 times SMOOTH(2) used when input SMOOTH(3) < net SMOOTH(2). SOLINT.....Averaging time in minutes: 0 -> 0.007 NDIG.......Number of digits to use in matrix print. Will actually use 1 more for a blank or minus sign. DOCRT......= 0 -> do not print the matrix average values False (< 0) use the line printer if OUTPRINT = ' ' else write named OUTPRINT file only. When OUTPRINT is not blank, DOCRT=-2 suppresses the page-feed character on page headers and DOCRT=-3 suppresses page headers and most other header information. When OUTPRINT is blank, the line printer will be used. PRTSY now counts the lines to be printed before printing any and asks for permission to continue if the count > 500. True (> 0) use the terminal interactively. The task will use the actual terminal width as a display limit unless 72 < DOCRT < width. In that case, the display limit will be DOCRT characters. OUTPRINT...Disk file name in which to save the line printer output. ' ' => use scratch and print immediately for interactive jobs - batch jobs use OUTPRINT = 'PRTFIL:BATCHjjj.nnn' (jjj= job #, nnn = user #). When OUTPRINT is not blank, multiple outputs are concatenated, and the file is not actually printed. DOPLOT.....1,3,5,7=> plot the histogram of the RMS samples using the following adverbs 2,3,6,7 => plot the averaged values versus time 4,5,6,7=> plot the histogram of the MEAN samples using the following adverbs 0 => no plots, only the printed matrices FACTOR.....Scale plotted points by FACTOR (in the time plots). 0 -> 1 NBOXES.....Number boxes in histogram: 6 - 1024. < 6 -> 128. A smaller number may be used depending on the number of time samples. Histograms with too few samples per bin are not very useful. CUTOFF.....Each histogram is computed in the range +- CUTOFF * rms about the mean of that parameter. 0 -> 3. FUNCTYPE...'LG' use log of histogram value in plot APARM......Plot range of values in time plot. APARM(2) <= APARM(1) -> self scale. DOTV.......> 0 => plot directly on the TV device, otherwise make a plot file for later display on one or more devices (including the TV if desired). GRCHAN.....Graphics channel to use. 0 -> multiple LTYPE......Labelling type, see HELP LTYPE for details: 1 = border, 2 = no ticks, 3 or 7 = standard, 4 or 8 = relative to ref. pixel, 5 or 9 = relative to subimage (BLC, TRC) center, 6 or 10 = pixels. 7-10 all labels other than tick numbers and axis type are omitted. Less than 0 is the same except that the plot file version number and create time are omitted. Add n * 100 to alter the metric scaling. BADDISK....The disk numbers to avoid for scratch files (sorting tables mostly). ----------------------------------------------------------------