[daip] feature in VPLOT
Leonia Kogan
lkogan at aoc.nrao.edu
Mon Oct 31 11:20:33 EST 2005
Eric,
I have found VPLOT has a limit for number of visibilitties:
MSAMP = APARM(7)
IF (MSAMP.LE.100) MSAMP=CATBLK(KIGCN)
MSAMP = MIN (MSAMP, 1000000)
Leonia
Claire Chandler wrote:
> 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
>
> _______________________________________________
> Daip mailing list
> Daip at listmgr.cv.nrao.edu
> http://listmgr.cv.nrao.edu/mailman/listinfo/daip
>
More information about the Daip
mailing list