[daip] feature in VPLOT
Claire Chandler
cchandle at aoc.nrao.edu
Sun Oct 30 19:53:10 EST 2005
I've found a bizarre feature in VPLOT. I have a continuum dataset with a
total number of visibilities of 1194364 and 23 antennas, so ~4720
visibilities per baseline. When I use VPLOT with the following inputs to
plot all sources but limiting the antennas to plot, it says it's plotting
sensible numbers of points:
habane> VPLOT1: PLTVB: 4878 points plotted: baseline 1 - 2
habane> VPLOT1: PLTVB: 4560 points plotted: baseline 1 - 3
habane> VPLOT1: PLTVB: 4746 points plotted: baseline 1 - 4
habane> VPLOT1: PLTVB: 4881 points plotted: baseline 1 - 5
habane> VPLOT1: PLTVB: 4883 points plotted: baseline 1 - 6
habane> VPLOT1: PLTVB: 4766 points plotted: baseline 1 - 7
habane> VPLOT1: PLTVB: 4773 points plotted: baseline 1 - 8
habane> VPLOT1: PLTVB: 4866 points plotted: baseline 1 - 9
plotted with inputs:
AIPS 1: VPLOT Plots data from a u,v data base, model from a CC file
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: USERID 0 Data base owner number
AIPS 1: INNAME '20051017' Input UV file name (name)
AIPS 1: INCLASS 'CH 0' Input UV file name (class)
AIPS 1: INSEQ 1 Input UV file name (seq. #)
AIPS 1: INDISK 1 Input UV file disk unit #
AIPS 1: IN2NAME ' ' Input CC file name (name)
AIPS 1: IN2CLASS ' ' Input CC file name (class)
AIPS 1: IN2SEQ 0 Input CC file name (seq. #)
AIPS 1: IN2DISK 0 Input CC file disk unit #
AIPS 1: IN2VERS 0 Input CC file version #
AIPS 1: NCOMP *all 0 # comps to use for model.
AIPS 1: 1 value per field
AIPS 1: total 25000 CCs maximum
AIPS 1: FLUX 0 Lowest CC component used.
AIPS 1: NMAPS 0 # clean maps (fields).
AIPS 1: 0 => no model
AIPS 1: >0 => model at samples
AIPS 1: <0 => model smoothly
AIPS 1: Data selection:
AIPS 1: SOURCES *all ' ' Source list
AIPS 1: QUAL -1 Source qualifier -1=>all
AIPS 1: CALCODE ' ' Calibrator code ' '=>all
AIPS 1: SELBAND -1 Bandwidth to select (kHz)
AIPS 1: SELFREQ -1 Frequency to select (MHz)
AIPS 1: FREQID -1 Freq. ID to select.
AIPS 1: BCHAN 1 Lowest channel number 0=>1
AIPS 1: ECHAN 0 Highest channel number
AIPS 1: 0=>BCHAN
AIPS 1: If BCHAN and ECHAN are
AIPS 1: different, then data will be
AIPS 1: averaged across specified
AIPS 1: channel range before plotting
AIPS 1: BIF 1 Lowest IF number 0=>1
AIPS 1: EIF 1 Highest IF number 0=>high
AIPS 1: AVGIF 1 > 0 => average IFs
AIPS 1: DO3COL 0 Plot separate IFs in separate
AIPS 1: colors.
AIPS 1: XINC 0 Plot every XINC'th visibility
AIPS 1: 0 => 1
AIPS 1: UVRANGE 0 0 Range of projected spacings
AIPS 1: (thousands of wavelengths)
AIPS 1: TIMERANG *all 0 Time: start day,hr,min,sec
AIPS 1: stop day,hr,min,sec.
AIPS 1: STOKES ' ' Stokes' parameters
AIPS 1: ANTENNAS 1 2 Antennas to be plotted
AIPS 1: *rest 0
AIPS 1: BASELINE *all 0 Baselines with ANTENNAS
AIPS 1: See explanation.
AIPS 1: DOCALIB 2 If >0 calibrate data
AIPS 1: = 2 calibrate weights
AIPS 1: GAINUSE 3 CL (or SN) table to apply
AIPS 1: DOPOL -1 If >0 correct polarization.
AIPS 1: BLVER -1 BL table to apply.
AIPS 1: FLAGVER 2 Flag table version
AIPS 1: See explanation.
AIPS 1: DOBAND -1 If >0 apply bandpass cal.
AIPS 1: Method used depends on value
AIPS 1: of DOBAND (see HELP file).
AIPS 1: BPVER 1 Bandpass table version
AIPS 1: SMOOTH *all 0 Smoothing function. See
AIPS 1: SUBARRAY 0 Subarray number 0=>1.
AIPS 1: -1 => all, but beware!
AIPS 1: DOEBAR 0 If > 0, plot error bars,
AIPS 1: else don't
AIPS 1: SOLINT 0 Time over which to average
AIPS 1: data before plotting (mins)
AIPS 1: =0 => no averaging
AIPS 1: >0 SOLINT=MIN(SOLINT, SCAN)
AIPS 1: <0 SOLINT=ABS(SOLINT)
AIPS 1: average through scans
AIPS 1: OPTYPE ' ' 'AUTO' -> autocorrelation
AIPS 1: only, else crosscorrelations
AIPS 1: APARM *all 0 Control parameters
AIPS 1: 1 : > 0 => edit bad
AIPS 1: amplitude points.
AIPS 1: BPARM(2)=0, or 1
AIPS 1: = 0 => do not edit
AIPS 1: bad amplitude points.
AIPS 1: Number of averaged points
AIPS 1: for each baseline has to
AIPS 1: be more than 10.
AIPS 1: 2 : number of sigmas outside
AIPS 1: of which bad points are
AIPS 1: clipped (See HELP); 0=>3.
AIPS 1: 3 : > 0 plot min/max range
AIPS 1: with average
AIPS 1: = 0 do not plot min/max
AIPS 1: lines
AIPS 1: 4 : = 0 min/max over selected
AIPS 1: IFs after averaging
AIPS 1: selected CHs
AIPS 1: > 0 min/max over all
AIPS 1: selected CH/IFs
AIPS 1: 5 : > 0 scalar average
AIPS 1: = 0 vector average
AIPS 1: 6 : used with APARM(1)>0
AIPS 1: >0 flag the bad points
AIPS 1: for all channels
AIPS 1: =0 flag the bad points
AIPS 1: only for the selected
AIPS 1: channels
AIPS 1: 7 : est total plot points
AIPS 1: BPARM *all 0 Control parameters
AIPS 1: 1 : X-axis type code
AIPS 1: 2 : Y-axis type code;
AIPS 1: if <0 => AMP & PHASE
AIPS 1: plotted together
AIPS 1: 3 : >0.0 => fixed scale
AIPS 1: <0.0 => fixed range
AIPS 1: 0.0 => sep. scales
AIPS 1: 4 : Xmin (fixed scale)
AIPS 1: 5 : Xmax (fixed scale)
AIPS 1: 6 : Ymin (fixed scale)
AIPS 1: 7 : Ymax (fixed scale)
AIPS 1: 8 : PHmin (fixed scale)
AIPS 1: 9 : PHmax (fixed scale)
AIPS 1: 10: > 0 -> plot models for
AIPS 1: baselines with no data
AIPS 1: ALIAS *all 0 Antenna alias selection
AIPS 1: If ALIAS(1) is non-zero, any
AIPS 1: antennas matching ALIAS(2-30)
AIPS 1: are treated as if they were
AIPS 1: actually ALIAS(1)
AIPS 1: NPLOTS 12 Max number of plots per page.
AIPS 1: 0 => 3; If bparm(2)<0 =>
AIPS 1: Max number of pairs of plots.
AIPS 1: DOARRAY -1 Array type : 0 = VLBI
AIPS 1: 1 = VLA and ATCA;
AIPS 1: 2 = MERLIN; 3 = WSRT
AIPS 1: POLPLOT ' ' Option to display various
AIPS 1: combinations of polzns:
AIPS 1: 'RL/RR', 'RL/LL',
AIPS 1: 'LR/RR', 'LR/LL', 'RR/LL'
AIPS 1: 'LL/RR'; other = don't use
AIPS 1: this option.
AIPS 1: SYMBOL 1 Symbol to use when plotting.
AIPS 1: 0 => Vertical line
AIPS 1: 1 - 23, standard types
AIPS 1: LTYPE 3 Type of labeling: 1 border,
AIPS 1: 2 no ticks, 3 - 6 standard,
AIPS 1: 7 - 10 only tick labels
AIPS 1: <0 -> no date/time
AIPS 1: DOTV 1 > 0 Do plot on the TV, else
AIPS 1: make a plot file
AIPS 1: GRCHAN 0 Graphics channel 0 => 1.
However, when I set ANTENNAS=0 with all other inputs the same, it (1) runs
very slowly (not unexpected), (2) claims to be plotting ~150,000
visibilities per baseline, and (3) does not include the full timerange on
X axis of the plot:
habane> VPLOT1: PLTVB: 149926 points plotted: baseline 1 - 3
habane> VPLOT1: PLTVB: 150019 points plotted: baseline 1 - 4
habane> VPLOT1: PLTVB: 150258 points plotted: baseline 1 - 5
habane> VPLOT1: PLTVB: 150254 points plotted: baseline 1 - 6
habane> VPLOT1: PLTVB: 150147 points plotted: baseline 1 - 7
habane> VPLOT1: PLTVB: 150130 points plotted: baseline 1 - 8
habane> VPLOT1: PLTVB: 150212 points plotted: baseline 1 - 9
I assume this is something to do with coping with being asked to plot too
many visibilities, but these numbers make no sense.
Claire
More information about the Daip
mailing list