[daip] 31DEC04 CALIB giving peculiar closure errors
Michael Rupen
mrupen at xeros.gsfc.nasa.gov
Wed Nov 10 10:19:20 EST 2004
Dear DAip,
the new version of CALIB is giving me some odd results, even for the old
SOLTYPEs ('' or 'L1'). Specifically, I get slightly fewer solutions and
dramatically higher reported closure errors in TST than in NEW (on alcor in
the AOC). The actual solutions are identical apart from the ones missing
from TST, as found by diff'ing PRTABs of the NEW and TST SN tables; but
the weights given in those tables are often (but not always!) different.
This can also be seen with SNPLT, where the SNRs for some antennas, at some
times, differ significantly from one another.
The data set is available on alcor, user 213, and is the only one present
(indisk 1 ; inname '20041109' ; incla 'C band '; inseq 1). It's a simple
continuum VLA data set taken in the current A configuration. See below for
details of CALIB setup and output. The details here are for soltype 'L1',
but '' and 'L1R' give similar results (I've not tried any others). I get
the same results from the MNJ'd version of AIPS at GSFC, for whatever that's
worth.
Any thoughts?
Michael
================================================================================
alcor AIPS NEW:
vers 'new' ; default CALIB
indisk 1 ; inname '20041109' ; incla 'C band '; inseq 1
calsour '17444-31166',''
docal 1 ; gainuse 1
uvra 10, 0
wtuv 0.0
ante 0
flagver 1
solint 0 ; solmode ''
soltype 'l1'
cparm 0 ; cparm(5) 0
aparm 0 ; aparm(6) 2
minamper 10 ; minphser 10 ; cparm(3) 1 ; cparm(4) 1
refant 10
timer 1 0 0 0 10 0 0 0
inp
AIPS 1: Found in Version=NEW
AIPS 1: CALIB: Task to determine calibration for data.
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: Input uv data.
AIPS 1: INNAME '20041109 ' UV file name (name)
AIPS 1: INCLASS 'C BAND' UV file name (class)
AIPS 1: INSEQ 1 UV file name (seq. #)
AIPS 1: INDISK 1 UV file disk drive #
AIPS 1: Data selection (multisource):
AIPS 1: CALSOUR '17444-31166 ' Calibrator sources
AIPS 1: *rest ' '
AIPS 1: QUAL -1 Calibrator 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: TIMERANG *all 0 Time range to use.
AIPS 1: BCHAN 1 Lowest channel number 0=>all
AIPS 1: ECHAN 0 Highest channel number
AIPS 1: ANTENNAS *all 0 Antennas to select. 0=all
AIPS 1: DOFIT *all 0 Subset of ANTENNAS list for
AIPS 1: which solns are desired. 0
AIPS 1: => all in ANTENNAS, < 0 all
AIPS 1: but those in DOFIT
AIPS 1: ANTUSE *all 0 Mean gain is calculated
AIPS 1: (CPARM(2)>0) using only the
AIPS 1: listed antennas. See explain.
AIPS 1: SUBARRAY 0 Subarray, 0=>all
AIPS 1: UVRANGE 10 0 Range of uv distance for full
AIPS 1: weight
AIPS 1: WTUV 0 Weight outside UVRANGE 0=0.
AIPS 1:
AIPS 1: Cal. info for input:
AIPS 1: DOCALIB 1 If >0 calibrate data
AIPS 1: = 2 calibrate weights
AIPS 1: GAINUSE 1 CL 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: 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: HELP SMOOTH for details.
AIPS 1:
AIPS 1: CLEAN map. See HELP.
AIPS 1: IN2NAME ' ' Cleaned map name (name)
AIPS 1: IN2CLASS ' ' Cleaned map name (class)
AIPS 1: IN2SEQ 0 Cleaned map name (seq. #)
AIPS 1: IN2DISK 0 Cleaned map disk unit #
AIPS 1: INVERS 0 CC file version #.
AIPS 1: NCOMP *all 0 # comps to use for model.
AIPS 1: 1 value per field
AIPS 1: FLUX 0 Lowest CC component used.
AIPS 1: NMAPS 0 No. Clean map files
AIPS 1: CMETHOD ' ' Modeling method:
AIPS 1: 'DFT','GRID',' '
AIPS 1: CMODEL ' ' Model type: 'COMP','IMAG'
AIPS 1: SMODEL *all 0 Source model, 1=flux,2=x,3=y
AIPS 1: See HELP SMODEL for models.
AIPS 1:
AIPS 1: Output uv data file.
AIPS 1: OUTNAME ' ' UV file name (name)
AIPS 1: OUTCLASS ' ' UV file name (class)
AIPS 1: OUTSEQ 0 UV file name (seq. #)
AIPS 1: OUTDISK 1 UV file disk drive #
AIPS 1:
AIPS 1: Solution control adverbs:
AIPS 1: REFANT 10 Reference antenna
AIPS 1: SOLINT 0 Solution interval (min)
AIPS 1: APARM 0 0 General parameters
AIPS 1: 0 0 1=min. no. antennas
AIPS 1: 0 2 2 > 0 => data divided
AIPS 1: *rest 0 3 > 0 => avg. RR,LL
AIPS 1: 5 > 0 => avg. IFs.
AIPS 1: 6=print level, 1=good,
AIPS 1: 2 closure, 3 SNR
AIPS 1: 7=SNR cutoff (0=>5)
AIPS 1: 8=max. ant. # (no AN)
AIPS 1: 9 > 0 => pass failed soln
AIPS 1: Phase-amplitude Parameters:
AIPS 1: SOLTYPE 'L1 ' Soln type,' ','L1','GCON'
AIPS 1: SOLMODE ' ' Soln. mode: 'A&P','P','P!A',
AIPS 1: 'GCON'
AIPS 1: SOLCON 0 Gain constraint factor.
AIPS 1: MINAMPER 10 Amplitude closure error
AIPS 1: regarded as excessive in %
AIPS 1: MINPHSER 10 Phase closure error regarded
AIPS 1: as excessive in degrees
AIPS 1: CPARM 0 0 Phase-amp. parameters
AIPS 1: 1 1 1 = Min el for gain
AIPS 1: *rest 0 normalization (deg)
AIPS 1: 2 >0 => normalize gain
AIPS 1: 3 avg. amp. closure err
AIPS 1: 4 avg. ph. closure err
AIPS 1: 5 >0 => scalar average
AIPS 1:
AIPS 1: SNVER 0 Output SN table, 0=>new table
AIPS 1: ANTWT *all 0 Ant. weights (0=>1.0)
AIPS 1: GAINERR *all 0 Std. Dev. of antenna gains.
AIPS 1: BADDISK *all 0 Disk no. not to use for
AIPS 1: scratch files.
>
> go calib
AIPS 1: Found in Version=NEW
CALIB1: Task CALIB (release of 31DEC03) begins
CALIB1: CALIB USING 20041109 . C BAND . 1 DISK= 1 USID= 213
CALIB1: L1 Solution type
CALIB1: Selecting, editing and calibrating the data
CALIB1: Doing cal transfer mode with point model for each source
CALIB1: This is not self-calibration
CALIB1: Dividing data by source flux densities
CALIB1: Determining solutions
CALIB1: Writing SN table 1
CALIB1: Closure errors at 1/00:00:26 17444-31166 IF no. 1 Rpol
CALIB1: Average closure error 1.315% 0.73d
CALIB1: Closure errors at 1/00:00:26 17444-31166 IF no. 2 Rpol
CALIB1: Average closure error 1.505% 0.81d
CALIB1: Closure errors at 1/00:00:26 17444-31166 IF no. 1 Lpol
CALIB1: Average closure error 1.204% 0.71d
CALIB1: Closure errors at 1/00:00:26 17444-31166 IF no. 2 Lpol
CALIB1: Average closure error 1.290% 0.79d
CALIB1: Closure errors at 1/00:03:26 17444-31166 IF no. 1 Rpol
CALIB1: Average closure error 1.414% 0.74d
CALIB1: Closure errors at 1/00:03:26 17444-31166 IF no. 2 Rpol
CALIB1: Average closure error 1.419% 0.73d
CALIB1: Closure errors at 1/00:03:26 17444-31166 IF no. 1 Lpol
CALIB1: Average closure error 1.328% 0.66d
CALIB1: Closure errors at 1/00:03:26 17444-31166 IF no. 2 Lpol
CALIB1: Average closure error 1.251% 0.63d
CALIB1: Closure errors at 1/00:06:26 17444-31166 IF no. 1 Rpol
CALIB1: Average closure error 1.286% 0.71d
CALIB1: Closure errors at 1/00:06:26 17444-31166 IF no. 2 Rpol
CALIB1: Average closure error 1.474% 0.77d
CALIB1: Closure errors at 1/00:06:26 17444-31166 IF no. 1 Lpol
CALIB1: Average closure error 1.265% 0.67d
CALIB1: Closure errors at 1/00:06:26 17444-31166 IF no. 2 Lpol
CALIB1: Average closure error 1.300% 0.69d
CALIB1: Closure errors at 1/00:09:26 17444-31166 IF no. 1 Rpol
CALIB1: Average closure error 1.460% 0.75d
CALIB1: Closure errors at 1/00:09:26 17444-31166 IF no. 2 Rpol
CALIB1: Average closure error 1.410% 0.77d
CALIB1: Closure errors at 1/00:09:26 17444-31166 IF no. 1 Lpol
CALIB1: Average closure error 1.339% 0.70d
CALIB1: Closure errors at 1/00:09:26 17444-31166 IF no. 2 Lpol
CALIB1: Average closure error 1.328% 0.62d
CALIB1: Found 398 good solutions
CALIB1: Failed on 1 solutions
CALIB1: 1 solutions had insufficient data
CALIB1: Appears to have ended successfully
CALIB1: alcor 31DEC03 NEW: Cpu= 0.2 Real= 0
>
================================================================================
alcor AIPS TST:
> vers ''
> inp
AIPS 1: CALIB: Task to determine calibration for data.
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: Input uv data.
AIPS 1: INNAME '20041109 ' UV file name (name)
AIPS 1: INCLASS 'C BAND' UV file name (class)
AIPS 1: INSEQ 1 UV file name (seq. #)
AIPS 1: INDISK 1 UV file disk drive #
AIPS 1: Data selection (multisource):
AIPS 1: CALSOUR '17444-31166 ' Calibrator sources
AIPS 1: *rest ' '
AIPS 1: QUAL -1 Calibrator 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: TIMERANG 1 0 Time range to use.
AIPS 1: 0 0 10 *rest 0
AIPS 1: BCHAN 1 Lowest channel number 0=>all
AIPS 1: ECHAN 0 Highest channel number
AIPS 1: ANTENNAS *all 0 Antennas to select. 0=all
AIPS 1: DOFIT *all 0 Subset of ANTENNAS list for
AIPS 1: which solns are desired.
AIPS 1: ANTUSE *all 0 Mean gain is calculated
AIPS 1: (CPARM(2)>0) using only the
AIPS 1: listed antennas. See explain.
AIPS 1: SUBARRAY 0 Subarray, 0=>all
AIPS 1: UVRANGE 10 0 Range of uv distance for full
AIPS 1: weight
AIPS 1: WTUV 0 Weight outside UVRANGE 0=0.
AIPS 1: WEIGHTIT 0 Modify data weights function
AIPS 1: Cal. info for input:
AIPS 1: DOCALIB 1 If >0 calibrate data
AIPS 1: = 2 calibrate weights
AIPS 1: GAINUSE 1 CL 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: 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: HELP SMOOTH for details.
AIPS 1:
AIPS 1: CLEAN map. See HELP.
AIPS 1: IN2NAME ' ' Cleaned map name (name)
AIPS 1: IN2CLASS ' ' Cleaned map name (class)
AIPS 1: IN2SEQ 0 Cleaned map name (seq. #)
AIPS 1: IN2DISK 0 Cleaned map disk unit #
AIPS 1: INVERS 0 CC file version #.
AIPS 1: NCOMP *all 0 # comps to use for model.
AIPS 1: 1 value per field
AIPS 1: FLUX 0 Lowest CC component used.
AIPS 1: NMAPS 0 No. Clean map files
AIPS 1: CMETHOD ' ' Modeling method:
AIPS 1: 'DFT','GRID',' '
AIPS 1: CMODEL ' ' Model type: 'COMP','IMAG'
AIPS 1: SMODEL *all 0 Source model, 1=flux,2=x,3=y
AIPS 1: See HELP SMODEL for models.
AIPS 1:
AIPS 1: Output uv data file.
AIPS 1: OUTNAME ' ' UV file name (name)
AIPS 1: OUTCLASS ' ' UV file name (class)
AIPS 1: OUTSEQ 0 UV file name (seq. #)
AIPS 1: OUTDISK 1 UV file disk drive #
AIPS 1:
AIPS 1: Solution control adverbs:
AIPS 1: REFANT 10 Reference antenna
AIPS 1: SOLINT 0 Solution interval (min)
AIPS 1: SOLSUB 0 Solution subinterval
AIPS 1: SOLMIN 0 Min solution interval
AIPS 1: APARM 0 0 General parameters
AIPS 1: 0 0 1=min. no. antennas
AIPS 1: 0 2 2 > 0 => data divided
AIPS 1: *rest 0 3 > 0 => avg. RR,LL
AIPS 1: 5 > 0 => avg. IFs.
AIPS 1: 6=print level, 1=good,
AIPS 1: 2 closure, 3 SNR
AIPS 1: 7=SNR cutoff (0=>5)
AIPS 1: 8=max. ant. # (no AN)
AIPS 1: 9 > 0 => pass failed soln
AIPS 1: Phase-amplitude Parameters:
AIPS 1: DOFLAG 0 Flag on closure error?
AIPS 1: SOLTYPE 'L1 ' Soln type,' ','L1','GCON',
AIPS 1: 'R', 'L1R', 'GCOR'
AIPS 1: SOLMODE ' ' Soln. mode: 'A&P','P','P!A',
AIPS 1: 'GCON',
AIPS 1: SOLCON 0 Gain constraint factor.
AIPS 1: MINAMPER 10 Amplitude closure error
AIPS 1: regarded as excessive in %
AIPS 1: MINPHSER 10 Phase closure error regarded
AIPS 1: as excessive in degrees
AIPS 1: CPARM 0 0 Phase-amp. parameters
AIPS 1: 1 1 1 = Min el for gain
AIPS 1: *rest 0 normalization (deg)
AIPS 1: 2 >0 => normalize gain
AIPS 1: 3 avg. amp. closure err
AIPS 1: 4 avg. ph. closure err
AIPS 1: 5 >0 => scalar average
AIPS 1:
AIPS 1: SNVER 0 Output SN table, 0=>new table
AIPS 1: ANTWT *all 0 Ant. weights (0=>1.0)
AIPS 1: GAINERR *all 0 Std. Dev. of antenna gains.
AIPS 1: BADDISK *all 0 Disk no. not to use for
AIPS 1: scratch files.
>
> go calib
CALIB1: Task CALIB (release of 31DEC04) begins
CALIB1: CALIB USING 20041109 . C BAND . 1 DISK= 1 USID= 213
CALIB1: L1 Solution type
CALIB1: Selecting, editing and calibrating the data
CALIB1: Doing cal transfer mode with point model for each source
CALIB1: This is not self-calibration
CALIB1: Dividing data by source flux densities
CALIB1: Determining solutions
CALIB1: Writing SN table 2
CALIB1: Closure errors at 1/00:00:26 17444-31166 IF no. 1 Rpol
CALIB1: Average closure error 1.313% 0.72d
CALIB1: Closure errors at 1/00:00:26 17444-31166 IF no. 2 Rpol
CALIB1: Average closure error 1.510% 0.80d
CALIB1: Closure errors at 1/00:00:26 17444-31166 IF no. 1 Lpol
CALIB1: Average closure error 1.202% 0.70d
CALIB1: Closure errors at 1/00:00:26 17444-31166 IF no. 2 Lpol
CALIB1: Average closure error 1.290% 0.79d
CALIB1: Closure errors at 1/00:03:26 17444-31166 IF no. 1 Rpol
CALIB1: Average closure error 1.382% 0.72d
CALIB1: Closure errors at 1/00:03:26 17444-31166 IF no. 2 Rpol
CALIB1: Average closure error 1.419% 0.73d
CALIB1: Closure errors at 1/00:03:26 17444-31166 IF no. 1 Lpol
CALIB1: Average closure error 1.328% 0.66d
CALIB1: Closure errors at 1/00:03:26 17444-31166 IF no. 2 Lpol
CALIB1: Average closure error 1.251% 0.62d
CALIB1: Closure errors at 1/00:06:26 17444-31166 IF no. 1 Rpol
CALIB1: Average closure error 26.584% 8.32d
CALIB1: 21-26 1576.0% 160d 21-27 1453.5% 123d 21-28 1486.4% 41d
CALIB1: 22-23 1475.8% 81d 22-24 1489.2% 125d 22-25 1480.9% 108d
CALIB1: 22-26 1508.9% 109d 22-27 1527.0% 147d 22-28 1455.3% 47d
CALIB1: 23-24 1480.7% 64d 23-25 1463.5% 64d 23-26 1485.0% 61d
CALIB1: 23-27 1444.2% 25d 23-28 1462.8% 142d 24-25 1473.2% 142d
CALIB1: 24-26 1485.2% 145d 24-27 1502.0% 179d 24-28 1489.0% 13d
CALIB1: 25-26 1555.1% 87d 25-27 1447.6% 51d 25-28 1547.1% 116d
CALIB1: 26-27 1496.8% 54d 26-28 1488.7% 112d 27-28 1480.7% 76d
CALIB1: Antenna 21 had 3 excess closure errors
CALIB1: Antenna 22 had 6 excess closure errors
CALIB1: Antenna 23 had 6 excess closure errors
CALIB1: Antenna 24 had 6 excess closure errors
CALIB1: Antenna 25 had 6 excess closure errors
CALIB1: Antenna 26 had 7 excess closure errors
CALIB1: Antenna 27 had 7 excess closure errors
CALIB1: Antenna 28 had 7 excess closure errors
CALIB1: Closure errors at 1/00:06:26 17444-31166 IF no. 2 Rpol
CALIB1: Average closure error 1.474% 0.77d
CALIB1: Closure errors at 1/00:06:26 17444-31166 IF no. 1 Lpol
CALIB1: Average closure error 26.408% 8.27d
CALIB1: 21-26 1518.4% 161d 21-27 1507.8% 123d 21-28 1453.1% 43d
CALIB1: 22-23 1455.0% 81d 22-24 1429.8% 123d 22-25 1450.5% 107d
CALIB1: 22-26 1492.7% 111d 22-27 1514.2% 149d 22-28 1490.7% 46d
CALIB1: 23-24 1469.5% 67d 23-25 1437.9% 67d 23-26 1431.5% 62d
CALIB1: 23-27 1457.5% 25d 23-28 1411.7% 143d 24-25 1442.0% 139d
CALIB1: 24-26 1493.9% 144d 24-27 1517.8% 180d 24-28 1455.7% 15d
CALIB1: 25-26 1499.0% 84d 25-27 1493.2% 50d 25-28 1496.7% 117d
CALIB1: 26-27 1528.1% 53d 26-28 1496.0% 112d 27-28 1510.4% 75d
CALIB1: Antenna 21 had 3 excess closure errors
CALIB1: Antenna 22 had 6 excess closure errors
CALIB1: Antenna 23 had 6 excess closure errors
CALIB1: Antenna 24 had 6 excess closure errors
CALIB1: Antenna 25 had 6 excess closure errors
CALIB1: Antenna 26 had 7 excess closure errors
CALIB1: Antenna 27 had 7 excess closure errors
CALIB1: Antenna 28 had 7 excess closure errors
CALIB1: Closure errors at 1/00:06:26 17444-31166 IF no. 2 Lpol
CALIB1: Average closure error 1.289% 0.67d
CALIB1: Closure errors at 1/00:09:26 17444-31166 IF no. 1 Rpol
CALIB1: Average closure error 2.404% 0.99d
CALIB1: 27-28 1499.4% 74d
CALIB1: Antenna 27 had 1 excess closure errors
CALIB1: Antenna 28 had 1 excess closure errors
CALIB1: Closure errors at 1/00:09:26 17444-31166 IF no. 2 Rpol
CALIB1: Average closure error 3.293% 0.99d
CALIB1: 27-28 9999.0% 71d
CALIB1: Antenna 27 had 1 excess closure errors
CALIB1: Antenna 28 had 1 excess closure errors
CALIB1: Closure errors at 1/00:09:26 17444-31166 IF no. 1 Lpol
CALIB1: Average closure error 4.194% 0.94d
CALIB1: 27-28 9999.0% 72d
CALIB1: Antenna 27 had 1 excess closure errors
CALIB1: Antenna 28 had 1 excess closure errors
CALIB1: Closure errors at 1/00:09:26 17444-31166 IF no. 2 Lpol
CALIB1: Average closure error 5.156% 0.86d
CALIB1: 27-28 9999.0% 75d
CALIB1: Antenna 27 had 1 excess closure errors
CALIB1: Antenna 28 had 1 excess closure errors
CALIB1: Found 397 good solutions
CALIB1: Failed on 3 solutions
CALIB1: Fraction of times having data > 2.5 rms from solution
CALIB1: 0.45455 of the times had 0 - 2 percent outside 2.5 times rms
CALIB1: Appears to have ended successfully
CALIB1: alcor 31DEC04 TST: Cpu= 0.3 Real= 0
================================================================================
More information about the Daip
mailing list