[daip] VPLOT model plot.
Craig Walker
cwalker at nrao.edu
Thu Jul 10 11:16:09 EDT 2008
Eric Greisen wrote:
> Craig Walker wrote:
>> When plotting models on the data with VPLOT with NMAPS < 0 for a
>> smooth plot, the model is there most, but not all, of the time. There
>> seem to be periods when it just stops being plotted. This is for
>> phase. The vanished model seems to be for periods when the absolute
>> value of the phase would be above about 100 deg (happens on both sides
>> of zero). When self scaling, it sometimes plots the data on a range up
>> to over 400 degrees and, in those cases, very little of the model
>> shows, so I think this is consistent with staying between -100 and +100.
>>
>> This is happening on my Linux box with data on Fred's box. Since this
>> is a fixed machine at NRAO, I presume the midnight job is up-to-date.
>>
>> The data are on SUM-FDAVIES_1 under user 542. This is plotting of a
>> CC table attached to a UV data set.
>>
>> You might hear more later about the amplitudes as the model from UVFIT
>> doesn't match the data at all. But I just noticed that and haven't
>> got it figured out yet.
>>
>> Cheers,
>>
>> Craig
>>
>> Here are the inputs:
>>
>> inp vplot
>> 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 'BW090LM84' Input UV file name (name)
>> AIPS 1: INCLASS 'SPLIT4' Input UV file name (class)
>> AIPS 1: INSEQ 1 Input UV file name (seq. #)
>> AIPS 1: INDISK 3 Input UV file disk unit #
>> AIPS 1: IN2NAME 'BW090LM84' Input CC file name (name)
>> AIPS 1: IN2CLASS 'SPLIT4' Input CC file name (class)
>> AIPS 1: IN2SEQ 1 Input CC file name (seq. #)
>> AIPS 1: IN2DISK 3 Input CC file disk unit #
>> AIPS 1: IN2VERS 1 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 -1 # 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 0 Freq. ID to select.
>> AIPS 1: BCHAN 1 Lowest channel number 0=>1
>> AIPS 1: ECHAN 1 Highest channel number
>> AIPS 1: AVGCHAN 1 > 0 => average BCHAN-ECHAN
>> AIPS 1: BIF 1 Lowest IF number 0=>1
>> AIPS 1: EIF 2 Highest IF number 0=>high
>> AIPS 1: AVGIF 1 > 0 => average IFs
>> AIPS 1: CROWDED 3 Plot chans/IFs together or
>> AIPS 1: separately? 0 separate,
>> AIPS 1: 3 all together, 1 & 2 partly
>> AIPS 1: DO3COL 1 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 'I' Stokes' parameters
>> AIPS 1: ANTENNAS *all 0 Antennas to be plotted
>> AIPS 1: BASELINE 1 2 Baselines with ANTENNAS
>> AIPS 1: 3 4 See explanation.
>> AIPS 1: 5 6 7 8
>> AIPS 1: 9 *rest 0
>> AIPS 1: DOCALIB 0 > 0 calibrate data & weights
>> AIPS 1: > 99 do NOT calibrate weights
>> AIPS 1: GAINUSE -1 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 -1 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 0 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 -1 If > 0, plot error bars,
>> AIPS 1: else don't
>> AIPS 1: SOLINT 0.4 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: else =>crosscorrelation
>> AIPS 1: OPCODE ' ' Type of plot:
>> AIPS 1: 'IFRA' => ratio BIF and EIF
>> AIPS 1: else => plot requested IF
>> 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 12 2 Control parameters
>> 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: REFANT 0 > 0 -> use REFANT for hour
>> AIPS 1: angle, elevation, et al.
>> AIPS 1: NPLOTS 9 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 0 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: ratios of polarizations:
>> AIPS 1: 'RL/RR', 'RL/LL', 'LR/RR',
>> AIPS 1: 'LR/LL', 'RR/LL' 'LL/RR';
>> AIPS 1: 'Q/I', 'U/I', 'V/I'.
>> AIPS 1: other = plot data directly
>> AIPS 1: SYMBOL 23 Symbol to use when plotting.
>> AIPS 1: 0 => Vertical line
>> AIPS 1: 1 - 23, standard types
>> AIPS 1: FACTOR 0 Scale symbol by FACTOR
>> AIPS 1: LTYPE 4 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.
>>
>>
>>
>>
>
> If you tell me that NMAPS > 0 is not working I will be interested. NMAPS
> < 0 is only expected to work once in a while.
>
> Eric Greisen
That's unfortunate. Our data is roughly equivalent to 4 or 5 scans
spread across the day - we were mainly imaging another source. When
NMAPS>0, the model is only calculated for the actual data points, but a
line connecting those points is plotted. True, that line seems to
always be there. But it does not reflect what the phases are actually
doing in the gaps. If there is a significant change between scans, as
there is in some of our data, the line will not reflect the right wrap
direction and the result is rather deceptive. What should be a simple
smooth model has some weird kinks in it. That is true of our M84 data.
It was also true of the data that caused me to find the scaling
problem with the case of a CC component on a UV data set. I finally
understood what the crazy looking model plot was trying to say when I
went to NMAPS<1 and saw the rapid phase winding. Clearly the program
knows how to calculate the model properly for nmaps<1. The difficulty
seems to be in putting the result on the plot.
Cheers,
Craig
---------------------------------------------------------------------
R. Craig Walker Array Operations Center
cwalker at nrao.edu National Radio Astronomy Observatory
Phone 505 835 7247 P. O. Box O
Fax 505 835 7027 Socorro NM 87801 USA
---------------------------------------------------------------------
More information about the Daip
mailing list