[daip] Problem with CLCOR
Mark Reid
reid at cfa.harvard.edu
Wed Feb 27 11:38:35 EST 2002
Dear Eric or Leonia,
I'm doing an unusual thing (as usual) and hit a problem
with CLCOR.
BACKGROUND (in case you're interested):
I've written a program that takes phase-refernced VLBA data
and fits for position offsets and a vertical atmospheric error
in the processor model for each station. I recently added
multi-band delay and fringe rate data to the fitting program
in order to help constrain the atmospheric parameters.
At present I'm trying to verify all the sign conventions.
So I'm using CLCOR to shift positions and also to change
the atmospheric model. Then I see if my fitting program
reproduces the CLCOR changes.
For testing purposes, I want infinite SNR data. So I took
an typical data set and wanted to replace all visibility data
by amp=1 and phase=0. The plan was then to run CLCOR to
shift things and then process the data in my usual manner
and fit for position and/or atmospheric deviations.
PROBLEM:
First, UVMOD doesn't run on multi-source data, and CLCOR
doesn't run on single source data. This complicates
things. So I did the following:
SPLIT: make a single source data set
UVMOD: replace all visibilities with amp=1 and phase=0
MULTI: reproducing original scan structure
INDXR: make sure I have CL entries every scan
TACOP: copy CL1 -> CL2
CLCOR: shift atmosphere of one antenna (see parameters
used given below)
SNPLT: to look at updated CL2 table
The CL2 table is now totally messed up. SNPLT says it now has
entries from day 0 to day 13 and does not have the original structure.
LISTR (see below) gives other messed up table info.
Any ideas?
Mark
--------------------------
AIPS 1: CLCOR Task which applies various corrections to CL tables.
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME '1745-283 ' Input UV file name (name)
AIPS 1: INCLASS 'UVMOD ' Input UV file name (class)
AIPS 1: INSEQ 2 Input UV file name (seq. #)
AIPS 1: INDISK 1 Input UV file disk unit #
AIPS 1: SOURCES *all ' ' Source list ' '=>all.
AIPS 1: STOKES ' ' Stokes type to process
AIPS 1: SELBAND -1 Bandwidth to select (kHz)
AIPS 1: SELFREQ -1 Frequency to select (MHz)
AIPS 1: FREQID -1 Freq. ID to select, 0=>all
AIPS 1: BIF 1 Lowest IF number 0=>all
AIPS 1: EIF 0 Highest IF number 0=>all
AIPS 1: TIMERANG *all 0 Time range to use.
AIPS 1: ANTENNAS 8 *rest 0 Antennas to correct.
AIPS 1: SUBARRAY 0 Subarray; 0 => 1.
AIPS 1: GAINVER 2 CL table version to update
AIPS 1: OPCODE 'ADEL' Operation code.
AIPS 1: CLCORPRM 500 *rest 0 Parameters (see HELP CLCOR).
AIPS 1: BADDISK *all 0 Disks to aviod for scratch
FYI the header for the file looks like this:
AIPS 1: Image=1745-283 (UV) Filename=1745-283 .UVMOD . 2
AIPS 1: Telescope=VLBA Receiver=VLBA
AIPS 1: Observer=BR066 User #= 180
AIPS 1: Observ. date=18-OCT-1999 Map date=27-FEB-2002
AIPS 1: # visibilities 148688 Sort order TB
AIPS 1: Rand axes: UU-L-SIN VV-L-SIN WW-L-SIN TIME1 BASELINE
AIPS 1: INTTIM GATEID CORR-ID SOURCE FREQSEL WEIGHT
AIPS 1: SCALE
AIPS 1: ----------------------------------------------------------------
AIPS 1: Type Pixels Coord value at Pixel Coord incr Rotat
AIPS 1: COMPLEX 1 1.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: STOKES 1 1.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: FREQ 1 4.3205474E+10 1.00 7.0000000E+06 0.00
AIPS 1: IF 4 1.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: RA 1 17 45 52.497 1.00 3600.000 0.00
AIPS 1: DEC 1 -28 20 26.294 1.00 3600.000 0.00
AIPS 1: ----------------------------------------------------------------
AIPS 1: Coordinate equinox 2000.00
AIPS 1: Rest freq 0.000 Vel type: OPTICAL wrt YOU
AIPS 1: Alt ref. value 0.00000E+00 wrt pixel 0.00
AIPS 1: Maximum version number of extension files of type SU is 1
AIPS 1: Maximum version number of extension files of type CL is 2
AIPS 1: Maximum version number of extension files of type HI is 1
AIPS 1: Maximum version number of extension files of type FQ is 1
AIPS 1: Maximum version number of extension files of type OB is 1
AIPS 1: Maximum version number of extension files of type WX is 1
AIPS 1: Maximum version number of extension files of type AN is 1
AIPS 1: Maximum version number of extension files of type NX is 1
===========================================
A LISTR (opco 'gain') printout of CL1 is as follows:
File = 1745-283 .UVMOD . 2 Vol = 1 Userid = 180 IF = 1
Freq= 43.205474375 GHz Ncor= 1 No. vis= 148688
Polarization = R Subarray = 1
Listing CL table, version 1
Gain phases in degrees
Stokes = R IF = 1 Freq = 43.205474375 GHz
Time Source -- 1-- 2-- 3-- 4-- 5-- 6-- 7-- 8-- 9
Day # 0
20:25:38 1745-283 0 0 0 0 0
20:25:53 1745-283 0 0 0 0 0
20:26:38 1745-283 0 0 0 0
20:26:53 1745-283 0 0 0 0
20:27:38 1745-283 0 0 0 0 0
20:27:53 1745-283 0 0 0 0 0
20:28:38 1745-283 0 0 0 0 0
20:28:53 1745-283 0 0 0 0 0
20:29:38 1745-283 0 0 0 0 0
20:29:53 1745-283 0 0 0 0 0
20:30:38 1745-283 0 0 0 0 0
20:30:53 1745-283 0 0 0 0 0
20:31:38 1745-283 0 0 0 0 0
20:31:53 1745-283 0 0 0 0 0
20:32:38 1745-283 0 0 0 0 0
20:32:53 1745-283 0 0 0 0 0
20:33:38 1745-283 0 0 0 0 0
20:33:53 1745-283 0 0 0 0 0
20:34:38 1745-283 0 0 0 0 0
20:34:53 1745-283 0 0 0 0 0
20:35:38 1745-283 0 0 0 0 0
================================
After running CLCOR the CL2 table (which was a TACOP of CL1)
looks like this:
LISTR1: Task LISTR (release of 31DEC01) begins
LISTR1: Determining scaling parameters
LISTR1: Maximum absolute phase = 179.90112 degrees
mjr LISTR(31DEC01) 180 27-FEB-2002 11:35:27 Page 1
File = 1745-283 .UVMOD . 2 Vol = 1 Userid = 180 IF = 1
Freq= 43.205474375 GHz Ncor= 1 No. vis= 148688
Polarization = R Subarray = 1
Listing CL table, version 2
Gain phases in degrees
Stokes = R IF = 1 Freq = 43.205474375 GHz
Time Source -- 1-- 2-- 3-- 4-- 5-- 6-- 7-- 8-- 9
Day # 0
00:42:16 1745-283 99
20:25:38 1745-283 0 0 0 0
17:20:00 1745-283 0 0 0 0 108
20:26:38 1745-283 0
Day # 2
13:07:06 1745-283 0 0 0 156
06:12:51 1745-283 0 0 0 0 174
Day # 0
17:00:48 1745-283 0 0 0 0-100
15:39:43 1745-283 0 0 0 0 66
17:17:43 1745-283 0 0 0 0 25
21:28:13 1745-283 0 0 0 0-131
17:51:52 1745-283 0 0 0 0 -23
17:57:59 1745-283 0 0 0 0 76
17:23:03 1745-283 0 0 0 0 108
11:37:39 1745-283 0 0 0 0 127
20:09:44 1745-283 0 0 0 0 135
15:50:00 1745-283 0 0 0 0 -92
14:10:46 1745-283 0 0 0 0-140
Type Q to stop, just hit RETURN to continue
mjr LISTR(31DEC01) 180 27-FEB-2002 11:35:40 Page 2
21:20:58 1745-283 0 0 0-127
15:57:16 1745-283 0 0 0-144
LISTR1: Determining scaling parameters
LISTR1: Maximum absolute phase = 179.50562 degrees
Type Q to stop, just hit RETURN to continue
mjr LISTR(31DEC01) 180 27-FEB-2002 11:35:53 Page 3
File = 1745-283 .UVMOD . 2 Vol = 1 Userid = 180 IF = 2
Freq= 43.205474375 GHz Ncor= 1 No. vis= 148688
Polarization = R Subarray = 1
Listing CL table, version 2
Gain phases in degrees
Stokes = R IF = 2 Freq = 43.213474375 GHz
Time Source -- 1-- 2-- 3-- 4-- 5-- 6-- 7-- 8-- 9
Day # 0
00:42:16 1745-283 60
20:25:38 1745-283 0 0 0 0
17:20:00 1745-283 0 0 0 0 69
20:26:38 1745-283 0
Day # 2
13:07:06 1745-283 0 0 0 118
06:12:51 1745-283 0 0 0 0 136
Day # 0
17:00:48 1745-283 0 0 0 0-131
15:39:43 1745-283 0 0 0 0 38
17:17:43 1745-283 0 0 0 0 -2
21:28:13 1745-283 0 0 0 0-157
17:51:52 1745-283 0 0 0 0 -47
17:57:59 1745-283 0 0 0 0 52
17:23:03 1745-283 0 0 0 0 85
11:37:39 1745-283 0 0 0 0 103
20:09:44 1745-283 0 0 0 0 111
15:50:00 1745-283 0 0 0 0-116
14:10:46 1745-283 0 0 0 0-167
Type Q to stop, just hit RETURN to continue
-------------------------------------------------------------
Mark J. Reid Phone: 617-495-7470
Harvard-Smithsonian CfA Fax : 617-495-7345
60 Garden Street Email: reid at cfa.harvard.edu
Cambridge, MA 02138, USA Web : cfa-www.harvard.edu/~reid
-------------------------------------------------------------
More information about the Daip
mailing list