[daip] CALIB broken?
Andy Biggs
abiggs at eso.org
Fri Feb 22 16:12:28 EST 2013
Hi Eric. Here are those various things you asked for.
Andy
AIPS 2: Image=MULTI (UV) Filename=MULTI .VBGLU . 3
AIPS 2: Telescope=GMRT Receiver=GMRT
AIPS 2: Observer=ANDY BIG User #= 4000
AIPS 2: Observ. date=26-JAN-2008 Map date=15-NOV-2011
AIPS 2: # visibilities 4391274 Sort order TB
AIPS 2: Rand axes: UU-L-SIN VV-L-SIN WW-L-SIN BASELINE TIME1
AIPS 2: SOURCE
AIPS 2: ----------------------------------------------------------------
AIPS 2: Type Pixels Coord value at Pixel Coord incr Rotat
AIPS 2: COMPLEX 3 1.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 2: STOKES 2 -1.0000000E+00 1.00 -1.0000000E+00 0.00
AIPS 2: FREQ 128 5.9406250E+08 1.00 1.2500000E+05 0.00
AIPS 2: IF 2 1.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 2: RA 1 00 00 00.000 1.00 3600.000 0.00
AIPS 2: DEC 1 00 00 00.000 1.00 3600.000 0.00
AIPS 2: ----------------------------------------------------------------
AIPS 2: Coordinate equinox 2000.00
AIPS 2: Maximum version number of extension files of type FQ is 1
AIPS 2: Maximum version number of extension files of type NX is 1
AIPS 2: Maximum version number of extension files of type HI is 1
AIPS 2: Maximum version number of extension files of type CL is 3
AIPS 2: Maximum version number of extension files of type SU is 1
AIPS 2: Maximum version number of extension files of type AN is 1
AIPS 2: Maximum version number of extension files of type FG is 5
AIPS 2: Maximum version number of extension files of type SN is 3
AIPS 2: Maximum version number of extension files of type BP is 1
AIPS 2: Maximum version number of extension files of type PL is 7
AIPS 2: CALIB: Task to determine calibration for data.
AIPS 2: Adverbs Values Comments
AIPS 2: ----------------------------------------------------------------
AIPS 2: Input uv data.
AIPS 2: INNAME 'MULTI' UV file name (name)
AIPS 2: INCLASS 'VBGLU' UV file name (class)
AIPS 2: INSEQ 3 UV file name (seq. #)
AIPS 2: INDISK 2 UV file disk drive #
AIPS 2: Data selection (multisource):
AIPS 2: CALSOUR '3C147' Calibrator sources
AIPS 2: '1400+621' *rest ' '
AIPS 2: QUAL -1 Calibrator qualifier -1=>all
AIPS 2: CALCODE ' ' Calibrator code ' '=>all
AIPS 2: SELBAND -1 Bandwidth to select (kHz)
AIPS 2: SELFREQ -1 Frequency to select (MHz)
AIPS 2: FREQID -1 Freq. ID to select.
AIPS 2: TIMERANG *all 0 Time range to use.
AIPS 2: ICHANSEL 75 75 Array of start and stop chn
AIPS 2: 1 *rest 0 numbers, plus a channel
AIPS 2: increment and IF to be used
AIPS 2: for channel selection in the
AIPS 2: averaging. See HELP ICHANSEL.
AIPS 2: Default = center 75% of band.
AIPS 2: ANTENNAS *all 0 Antennas to select. 0=all
AIPS 2: DOFIT *all 0 Subset of ANTENNAS list for
AIPS 2: which solns are desired.
AIPS 2: ANTUSE *all 0 Mean gain is calculated
AIPS 2: (CPARM(2)>0) using only the
AIPS 2: listed antennas. See explain.
AIPS 2: SUBARRAY 0 Subarray, 0=>all
AIPS 2: UVRANGE 2 300 Range of uv distance for full
AIPS 2: weight
AIPS 2: WTUV 0 Weight outside UVRANGE 0=0.
AIPS 2: WEIGHTIT 3 Modify data weights function
AIPS 2: Cal. info for input:
AIPS 2: DOCALIB 1 > 0 calibrate data & weights
AIPS 2: > 99 do NOT calibrate weights
AIPS 2: NOT for single-source files
AIPS 2: GAINUSE 2 CL table to apply.
AIPS 2: DOPOL -1 If >0 correct polarization.
AIPS 2: PDVER 0 PD table to apply (DOPOL>0)
AIPS 2: BLVER -1 BL table to apply.
AIPS 2: FLAGVER 0 Flag table version
AIPS 2: DOBAND -1 If >0 apply bandpass cal.
AIPS 2: Method used depends on value
AIPS 2: of DOBAND (see HELP file).
AIPS 2: BPVER -1 Bandpass table version
AIPS 2: SMOOTH *all 0 Smoothing function. See
AIPS 2: HELP SMOOTH for details.
AIPS 2:
AIPS 2: CLEAN map. See HELP.
AIPS 2: IN2NAME ' ' Cleaned map name (name)
AIPS 2: IN2CLASS ' ' Cleaned map name (class)
AIPS 2: IN2SEQ 0 Cleaned map name (seq. #)
AIPS 2: IN2DISK 0 Cleaned map disk unit #
AIPS 2: INVERS 0 CC file version #.
AIPS 2: NCOMP *all 0 # comps to use for model.
AIPS 2: 1 value per field
AIPS 2: FLUX 0 Lowest CC component used.
AIPS 2: NMAPS 0 No. Clean map files
AIPS 2: CMETHOD ' ' Modeling method:
AIPS 2: 'DFT','GRID',' '
AIPS 2: CMODEL ' ' Model type: 'COMP','IMAG'
AIPS 2: 'SUBI' (see HELP re images)
AIPS 2: SMODEL *all 0 Source model, 1=flux,2=x,3=y
AIPS 2: See HELP SMODEL for models.
AIPS 2: Output uv data file.
AIPS 2: OUTNAME ' ' UV file name (name)
AIPS 2: OUTCLASS ' ' UV file name (class)
AIPS 2: OUTSEQ 0 UV file name (seq. #)
AIPS 2: OUTDISK 1 UV file disk drive #
AIPS 2: DOAPPLY 0 >= 0 write output if
AIPS 2: single source
AIPS 2: Solution control adverbs:
AIPS 2: REFANT 4 Reference antenna
AIPS 2: SOLINT 0.5 Solution interval (min)
AIPS 2: SOLSUB 0 Solution subinterval
AIPS 2: SOLMIN 0 Min solution interval
AIPS 2: APARM *all 0 General parameters
AIPS 2: 1=min. no. antennas
AIPS 2: 2 > 0 => data divided
AIPS 2: 3 > 0 => avg. RR,LL
AIPS 2: 5 > 0 => avg. IFs.
AIPS 2: 6=print level, 1=good,
AIPS 2: 2 closure, 3 SNR
AIPS 2: 7=SNR cutoff (0=>5)
AIPS 2: 8=max. ant. # (no AN)
AIPS 2: 9 > 0 => pass failed soln
AIPS 2: 10 < 99 cal output weights
AIPS 2: Phase-amplitude Parameters:
AIPS 2: DOFLAG 0 Flag on closure error?
AIPS 2: SOLTYPE 'L1R' Soln type,' ','L1','GCON',
AIPS 2: 'R', 'L1R', 'GCOR'
AIPS 2: SOLMODE 'A&P' Soln. mode: 'A&P','P','P!A',
AIPS 2: 'GCON',
AIPS 2: SOLCON 0 Gain constraint factor.
AIPS 2: MINAMPER 0 Amplitude closure error
AIPS 2: regarded as excessive in %
AIPS 2: MINPHSER 0 Phase closure error regarded
AIPS 2: as excessive in degrees
AIPS 2: NORMALIZ 0 >0 => normalize gain:
AIPS 2: 1 globally, 2 by subarray,
AIPS 2: 3 by subarray,IF
AIPS 2: 4 by subarray,IF,pol
AIPS 2: CPARM *all 0 Phase-amp. parameters
AIPS 2: 1 = Min el for gain
AIPS 2: normalization (deg)
AIPS 2: 2 > 0 normalize w median
AIPS 2: else use mean
AIPS 2: 3 avg. amp. closure err
AIPS 2: 4 avg. ph. closure err
AIPS 2: 5 = 1 vector average
AIPS 2: channels, scalar avg
AIPS 2: between times
AIPS 2: >= 2 scalar average
AIPS 2: 6 limit clipping in robust
AIPS 2: 7 limit display of closure
AIPS 2: errors
AIPS 2: SNVER 3 Output SN table, 0=>new table
AIPS 2: ANTWT *all 0 Ant. weights (0=>1.0)
AIPS 2: GAINERR *all 0 Std. Dev. of antenna gains.
AIPS 2: BADDISK *all 0 Disk no. not to use for
AIPS 2: scratch files.
CALIB2: Task CALIB (release of 31DEC13) begins
CALIB2: Using MULTI . VBGLU . 3 Disk= 2 Usid= 4000
CALIB2: L1 Solution type
CALIB2: UVGET: Using flag table version 5 to edit data
CALIB2: Selecting, editing and calibrating the data
AIPS 2: Resumes
CALIB2: This is not self-calibration
CALIB2: Dividing data by source flux densities
CALIB2: Determining solutions using full vector averaging
CALIB2: Writing SN table 3
CALIB2: RPOL, IF= 1 The average gain over these antennas is 1.102E+00
CALIB2: RPOL, IF= 2 The average gain over these antennas is 1.014E+00
CALIB2: LPOL, IF= 1 The average gain over these antennas is 9.611E-01
CALIB2: LPOL, IF= 2 The average gain over these antennas is 9.061E-01
CALIB2: Found 21044 good solutions
CALIB2: 1708 solutions had no data
CALIB2: Average closure rms = 0.05253 +- 0.02232
CALIB2: Fraction of times having data > 2.5 rms from solution
CALIB2: 0.40603 of the times had 0 - 2 percent outside 2.5 times rms
CALIB2: 0.19196 of the times had 2 - 4 percent outside 2.5 times rms
CALIB2: 0.01910 of the times had 4 - 6 percent outside 2.5 times rms
CALIB2: 0.03819 of the times had 6 - 8 percent outside 2.5 times rms
CALIB2: 0.04623 of the times had 8 - 10 percent outside 2.5 times rms
CALIB2: 0.05729 of the times had 10 - 12 percent outside 2.5 times rms
CALIB2: 0.04623 of the times had 12 - 14 percent outside 2.5 times rms
CALIB2: 0.00402 of the times had 14 - 16 percent outside 2.5 times rms
CALIB2: Appears to have ended successfully
CALIB2: pc016987 31DEC13 TST: Cpu= 17.3 Real= 153 IO= 4845
LISTR2: Task LISTR (release of 31DEC13) begins
pc016987 LISTR(31DEC13) 4000 22-FEB-2013 22:07:32 Page 1
File = MULTI .VBGLU . 3 Vol = 2 Userid = 4000
Freq = 0.594062500 GHz Ncor = 2 No. vis = 4391274
Scan summary listing
Scan Source Qual Calcode Sub Timerange FrqID
START VIS END VIS
1 3C147 : 0000 1 0/18:43:37 - 0/19:11:27
1 1 173565
2 1400+621 : 0000 1 0/19:16:41 - 0/19:21:35
1 173566 204450
3 GROTHN3 : 0000 1 0/19:22:38 - 0/19:53:53
1 204451 399330
4 1400+621 : 0000 1 0/19:54:47 - 0/19:59:41
1 399331 430215
5 GROTHN3 : 0000 1 0/20:00:31 - 0/20:31:46
1 430216 625095
6 1400+621 : 0000 1 0/20:32:53 - 0/20:37:47
1 625096 655980
7 GROTHN3 : 0000 1 0/20:38:41 - 0/21:10:09
1 655981 852165
8 1400+621 : 0000 1 0/21:11:07 - 0/21:16:01
1 852166 883050
9 GROTHN3 : 0000 1 0/21:17:00 - 0/21:48:19
1 883051 1078365
10 1400+621 : 0000 1 0/21:49:22 - 0/21:54:15
1 1078366 1109250
11 GROTHN3 : 0000 1 0/21:55:10 - 0/22:26:37
1 1109251 1305435
12 1400+621 : 0000 1 0/22:27:49 - 0/22:32:42
1 1305436 1336320
13 GROTHN3 : 0000 1 0/22:33:41 - 0/23:04:47
1 1336321 1525806
14 1400+621 : 0000 1 0/23:05:54 - 0/23:10:48
1 1525807 1556691
15 GROTHN3 : 0000 1 0/23:11:43 - 0/23:43:06
1 1556692 1752441
16 1400+621 : 0000 1 0/23:44:00 - 0/23:48:54
1 1752442 1783326
17 GROTHN3 : 0000 1 0/23:50:05 - 1/00:21:33
1 1783327 1982556
18 1400+621 : 0000 1 1/00:22:31 - 1/00:27:25
1 1982557 2013876
19 GROTHN3 : 0000 1 1/00:28:28 - 1/00:59:55
1 2013877 2216151
20 1400+621 : 0000 1 1/01:00:54 - 1/01:05:48
1 2216152 2247906
21 GROTHN3 : 0000 1 1/01:06:51 - 1/01:38:01
1 2247907 2444961
22 1400+621 : 0000 1 1/01:39:17 - 1/01:44:10
1 2444962 2475846
23 GROTHN3 : 0000 1 1/01:45:09 - 1/02:16:36
1 2475847 2675076
24 1400+621 : 0000 1 1/02:17:39 - 1/02:22:37
1 2675077 2706396
25 GROTHN3 : 0000 1 1/02:23:36 - 1/02:55:03
1 2706397 2904756
26 1400+621 : 0000 1 1/02:56:15 - 1/03:01:08
1 2904757 2936511
27 GROTHN3 : 0000 1 1/03:02:03 - 1/03:33:22
1 2936512 3133131
28 1400+621 : 0000 1 1/03:34:25 - 1/03:39:18
1 3133132 3164016
29 GROTHN3 : 0000 1 1/03:40:21 - 1/04:11:45
1 3164017 3359012
30 1400+621 : 0000 1 1/04:12:43 - 1/04:17:37
1 3359013 3389462
31 GROTHN3 : 0000 1 1/04:18:44 - 1/04:50:03
1 3389463 3574962
32 1400+621 : 0000 1 1/04:51:06 - 1/04:55:59
1 3574963 3602556
33 GROTHN3 : 0000 1 1/04:56:54 - 1/05:28:17
1 3602557 3788098
34 1400+621 : 0000 1 1/05:29:29 - 1/05:34:22
1 3788099 3817736
35 GROTHN3 : 0000 1 1/05:35:25 - 1/06:06:48
1 3817737 4021548
36 1400+621 : 0000 1 1/06:07:51 - 1/06:12:45
1 4021549 4056058
37 GROTHN3 : 0000 1 1/06:13:35 - 1/06:23:39
1 4056059 4118848
38 1400+621 : 0000 1 1/06:24:55 - 1/06:32:49
1 4118849 4174470
39 3C48 : 0000 1 1/06:54:16 - 1/07:30:50
1 4174471 4391274
Type Q to stop, just hit RETURN to continue
pc016987 LISTR(31DEC13) 4000 22-FEB-2013 22:07:34 Page 2
File = MULTI .VBGLU . 3 Vol = 2 Userid = 4000
Source summary
Velocity type = 'TOPOCENT' Definition = 'RADIO '
ID Source Qual Calcode RA(2000.0) Dec(2000.0) IFlux
QFlux UFlux VFlux No. vis
1 3C147 : 0000 05:42:36.1330 49:51:07.165 36.714
0.000 0.000 0.000 173565
IF( 2) 36.254
0.000 0.000 0.000
2 1400+621 : 0000 14:00:28.6536 62:10:38.465 6.302
0.000 0.000 0.000 612814
IF( 2) 6.396
0.000 0.000 0.000
3 GROTHN3 : 0000 14:22:24.7041 53:18:10.967 0.000
0.000 0.000 0.000 3388091
IF( 2) 0.000
0.000 0.000 0.000
4 3C48 : 0000 01:37:41.2777 33:09:35.683 29.414
0.000 0.000 0.000 216804
IF( 2) 28.932
0.000 0.000 0.000
ID Source Freq(GHz) Velocity(Km/s) Rest freq (GHz)
1 All Sources 0.5941 0.0000 0.0000
IF( 2) 0.6101 0.0000 0.0000
Frequency Table summary
FQID IF# Freq(GHz) BW(kHz) Ch.Sep(kHz) Sideband Bandcode
1 1 0.59406250 16000.0010 125.0000 1
2 0.61006250 16000.0010 125.0000 1
LISTR2: Appears to have ended successfully
LISTR2: pc016987 31DEC13 TST: Cpu= 0.0 Real= 2
AIPS 2: Resumes
On 22/02/2013 21:34, Eric Greisen wrote:
> Andrew Biggs wrote:
>> Hi Eric. The colours are different IFs/polarizations and the station IDs
>> are C for "central" (a square km full of antennas) and E means "East"
>> i.e.
>> for the East arm with numbers increasing according to how far away they
>> are from the centre. S therefore refers to the South arm.
>>
>> Anyway, I've tracked down the source of the problem. The clue is that, in
>> the plot I sent you, the first source looks fine. If I therefore run the
>> 31DEC13 CALIB with just the second source (I usually run it with both
>> simultaneously), the resulting solutions look absolutely fine. Does this
>> ring any bells?
>
> Given that I cannot find any difference between 31DEC13 and 31DEC12
> CALIB in either CALIB.FOR or any relevant subroutines, it is hard to
> come up with a clue. Is your 31DEC12 up to date or when was the last
> time you ran a MNJ on it? Can I see your inputs to CALIB when you try to
> do the 2 sources at once (and it fails)? Also an IMHEADER of the input
> file and perhaps a SCAN output of LISTR? It is as if some model is being
> applied to both sources but with a position error in the 2nd when you do
> both at once - perhaps the CALIB messages as well might indicate this.
>
> Thanks
>
> Eric
--
ESO, ALMA Regional Centre
Karl-Schwarzschild-Strasse 2
D-85748 Garching, Germany
tel. +49-89-32006471
fax. +49-89-32006898
More information about the Daip
mailing list