I think I have found a problem with LISTR and RLDIF application of calibration. I find this hard to believe since these are bread-and-butter programs, but here is the story. I am processing BW088A (Observed 27 January) which is 7mm VLBA observations of M87. I am working on my Mac PowerBook G4. I have a multisource file with extensive calibration tables and a split file made from it for calibrator OJ287. Here are the SPLIT inputs used to create the SPLIT file: AIPS 1: SPLIT Task to split multi-source uv data to single source AIPS 1: Adverbs Values Comments AIPS 1: ---------------------------------------------------------------- AIPS 1: also works on single files. AIPS 1: INNAME 'BW088A' Input UV file name (name) AIPS 1: INCLASS 'MULTI' Input UV file name (class) AIPS 1: INSEQ 1 Input UV file name (seq. #) AIPS 1: INDISK 4 Input UV file disk unit # AIPS 1: SOURCES 'OJ287' Source list AIPS 1: *rest ' ' AIPS 1: QUAL -1 Source qualifier -1=>all AIPS 1: CALCODE ' ' Calibrator code ' '=>all AIPS 1: TIMERANG *all 0 Time range to copy AIPS 1: STOKES ' ' Stokes type to pass. 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: BIF 0 Lowest IF number 0=>all AIPS 1: EIF 0 Highest IF number 0=>all AIPS 1: BCHAN 3 Lowest channel number 0=>all AIPS 1: ECHAN 28 Highest channel number AIPS 1: SUBARRAY 0 Subarray, 0=>all AIPS 1: DOCALIB 2 > 0 calibrate data & weights AIPS 1: > 99 do NOT calibrate weights AIPS 1: GAINUSE 15 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: DOBAND 1 If >0 apply bandpass cal. AIPS 1: Method used depends on value AIPS 1: of DOBAND (see HELP file). AIPS 1: BPVER 4 Bandpass table version AIPS 1: SMOOTH *all 0 Smoothing function. See AIPS 1: HELP SMOOTH for details. AIPS 1: OUTCLASS 'SPLITP' Output UV file name (class) AIPS 1: OUTSEQ 0 Output UV file name (seq. #) AIPS 1: OUTDISK 4 Output UV file disk unit #. AIPS 1: DOUVCOMP 1 1 (T) => compressed data AIPS 1: APARM 2 1.9661 Control information: AIPS 1: 0 1 1 = 1 => avg. freq. in IF AIPS 1: 0 1 multi-channel out AIPS 1: *rest 0 = 2 => avg. freq. in IF AIPS 1: single channel out AIPS 1: = 3 => avg IF's also AIPS 1: 2 = Input avg. time (sec) AIPS 1: 3 > 0 => Drop subarrays AIPS 1: 4 > 0 => calibrate weights AIPS 1: 5 = 0 pass only xc data AIPS 1: = 1 pass xc and ac data AIPS 1: = 2 pass only ac data AIPS 1: 6 > 0 add full source name AIPS 1: to header AIPS 1: NCHAV 1 Number of chan. to average. AIPS 1: (used if APARM(1) = 1) AIPS 1: <= 0 -> ALL AIPS 1: CHINC 1 Channel incr. between output AIPS 1: channels (used if APARM(1)=1) AIPS 1: ICHANSEL *all 0 Array of channel start, stop, AIPS 1: and increment numbers and IF AIPS 1: number to be used when AIPS 1: averaging in frequency. AIPS 1: (used if APARM(1) = 2, 3) AIPS 1: BADDISK *all 0 Disks to avoid for scratch When I use VPLOT to display the data from either the MULTI file with the same calibration as above or the SPLIT file using no calibration, I get the same plot. Here are the VPLOT inputs to look at the multi plot. For looking at the split file, I set DOCALIB 0; DOPOL=0; DOBAND=0; BCHAN=1; ECHAN=0 and otherwise make no changes. Staring at the screen as it flips from one to the other, they look the same. >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 'BW088A' Input UV file name (name) AIPS 1: INCLASS 'MULTI' Input UV file name (class) AIPS 1: INSEQ 1 Input UV file name (seq. #) AIPS 1: INDISK 4 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 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 'OJ287' Source list AIPS 1: *rest ' ' 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 3 Lowest channel number 0=>1 AIPS 1: ECHAN 28 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 2 Highest IF number 0=>high AIPS 1: AVGIF 0 > 0 => average IFs AIPS 1: DO3COL 1 Plot separate IFs in separate AIPS 1: colors. AIPS 1: XINC 1 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 'RL' Stokes' parameters AIPS 1: ANTENNAS *all 0 Antennas to be plotted AIPS 1: BASELINE *all 0 Baselines with ANTENNAS AIPS 1: See explanation. AIPS 1: DOCALIB 2 > 0 calibrate data & weights AIPS 1: > 99 do NOT calibrate weights AIPS 1: GAINUSE 15 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 4 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 .6 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 'AMP' '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: *rest 0 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 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: 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 23 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. When I run RLDIF on the SPLIT file with no calibration applied, I get numbers that look reasonable relative to the VPLOT numbers. Here is an example matrix: File = BW088AOJ287 .SPLITP. 1 Vol = 4 Userid = 501 Channel = 1 IF = 1 Freq= 43.126927500 GHz Ncor= 4 No. vis= 13663 Stokes = FULL Subarray = 1 Time = 0/08:37:41 to 0/08:42:58 Source = OJ287 : 0000 Flux = .0000 Jy, Calcode = , Freq = 43.126927500 GHz, IF = 1 Phases in degrees, averaging type = Vector RL in upper right, conjg(LR) in lower left Ant -- 1-- 2-- 3-- 4-- 5-- 6-- 7-- 8-- 9 1| -74 -82 -78 -33 -81 -94 -90 -93 2|-164 -150-113 14-112-139-114-118 3|-142 -93 -97 -24 -86-106-112 -75 4|-145 -96-107 5 -99-117-117 -74 5|-128 -89 -99 -92 -107 -92 -94 -96 6|-116 -76-106 -81 -31 -96 -84 -57 7|-130 -78-112 -88 -32 -76 -88 -92 8|-160-102-139-110 -4 -86-130 -99 9|-171-102-139-118 -41-127-146-126 Ampscalar average of matrix = -98.42( 1.864) sigma = 15.906 Vector average of upper data = -89.22 sigma =106.873 Vector average of lower data = -98.30 sigma =115.401 When I run on the multi file, applying calibration, I get numbers that just don't look right. Here is the same matrix File = BW088A .MULTI . 1 Vol = 4 Userid = 501 Channel = 3 IF = 1 Freq= 43.119490000 GHz Ncor= 4 No. vis= 17845 Stokes = FULL Subarray = 1 Applying calibration table 15 Applying polarization corrections Applying flag table 2 Applying bandpass table 4 Time = 0/08:37:39 to 0/08:42:58 Source = OJ287 : 0000 Flux = .0000 Jy, Calcode = V , Freq = 43.119490000 GHz, IF = 1 Phases in degrees, averaging type = Vector RL in upper right, conjg(LR) in lower left Ant -- 1-- 2-- 3-- 4-- 5-- 6-- 7-- 8-- 9 1| -60 -46 -99 26 35 -71-100-119 2|-130 -53-114 -40-112 -41 -91-156 3|-150-105 -114 -78-177-112 -91 -39 4|-114-116 135 27-131-156 174 -16 5|-101 -82 -71 -75 -91 -86-110 -80 6| 9 -97 128-118 -51 3 -94 7 7|-130-142 163 -53 -46 -91 -91-157 8| -94-102 -81 -83 6-102 -87 -162 9| 39-117-167 -49 27 26 79 -69 Ampscalar average of matrix = -88.72( 3.786) sigma = 33.801 Vector average of upper data = -89.82 sigma =535.460 Vector average of lower data = -94.68 sigma =561.085 Note that they don't agree at a rather serious level on at some baselines. Here are the RLDIF inputs: >inp rldif AIPS 1: RLDIF: Task to return Right - Left phase difference AIPS 1: Adverbs Values Comments AIPS 1: ---------------------------------------------------------------- AIPS 1: USERID 0 User number. AIPS 1: INNAME 'BW088A' UV data (name). AIPS 1: INCLASS 'MULTI' UV data (class). AIPS 1: INSEQ 1 UV data (seq. #). 0 => high AIPS 1: INDISK 4 Disk unit #. 0 => any AIPS 1: SOURCES 'OJ287' Source list AIPS 1: *rest ' ' AIPS 1: CALCODE ' ' Calibrator code ' '=>all AIPS 1: TIMERANG *all 0 Time range to list 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: None selected => 1. AIPS 1: BIF 1 Lowest IF number 0=1 AIPS 1: EIF 2 Highest IF number AIPS 1: BCHAN 3 Low channel number 0=>1 AIPS 1: ECHAN 28 High channel number AIPS 1: ANTENNAS *all 0 Antennas to list AIPS 1: UVRANGE 0 0 UV range in kilolambda AIPS 1: SUBARRAY 0 Subarray, 0=>1 AIPS 1: Cal. info for input: AIPS 1: DOCALIB 2 > 0 calibrate data & weights AIPS 1: > 99 do NOT calibrate weights AIPS 1: GAINUSE 15 CAL (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: DOBAND 1 If >0 apply bandpass cal. AIPS 1: Method used depends on value AIPS 1: of DOBAND (see HELP file). AIPS 1: BPVER 4 Bandpass table version AIPS 1: SMOOTH *all 0 Smoothing function. See AIPS 1: HELP SMOOTH for details. AIPS 1: DOCRT 132 > 0 -> use the terminal, AIPS 1: else use the line printer AIPS 1: > 72 => terminal width AIPS 1: = 0 -> no printing AIPS 1: OUTPRINT ' ' AIPS 1: Printer disk file to save AIPS 1: BADDISK *all 0 Disk to avoid for scratch. When I run LISTR MATX set up to look like RLDIF on the MULTI file, it gives the same answers. So whatever is causeing the problem seems to be shared between those two programs (I'll hazard a guess that RLDIF is a LISTR clone). I looked casually for a correlation between the difference in the two RLDIF results and the delays, thinking there might be some effect related to the reference frequency and application, or not, of delay, but I don't see it. While neither RLDIF run gives results of the quality one would like for an EVPA calibration, the one on the SPLIT file does look rather better. I sure hope I'm doing something stupid here, but I don't see it.