[daip] Fittp/Fitld impact on indxr
Eric Greisen
egreisen at nrao.edu
Fri Jul 1 10:50:00 EDT 2011
Arnaud Collioud wrote:
> Dear AIPS developers,
>
> I have to AIPS installations: one DEC07 on the lab server and one DEC11 on my mac. I was processing a session on the DEC07 version and I wanted to transfer the uv-data to continue on the DEC11 (up-to-date version).
>
> On DEC07, I used:
>
>> imh
> AIPS 1: Image=MULTI (UV) Filename=IYA09_S .UV . 1
> AIPS 1: Telescope=VLBI Receiver=MKIV
> AIPS 1: Observer=Bonn Cor User #= 4
> AIPS 1: Observ. date=18-NOV-2009 Map date=30-JUN-2011
> AIPS 1: # visibilities 4842987 Sort order TB
> AIPS 1: Rand axes: UU-L-SIN VV-L-SIN WW-L-SIN TIME1 BASELINE
> AIPS 1: SOURCE FREQSEL INTTIM WEIGHT 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: FREQ 32 2.2329900E+09 1.00 2.5000000E+05 0.00
> AIPS 1: STOKES 1 -1.0000000E+00 1.00 -1.0000000E+00 0.00
> AIPS 1: IF 4 1.0000000E+00 1.00 1.0000000E+00 0.00
> AIPS 1: RA 1 00 00 00.000 1.00 3600.000 0.00
> AIPS 1: DEC 1 00 00 00.000 1.00 3600.000 0.00
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: Coordinate equinox 2000.00
> 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 BS 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 CL is 1
> AIPS 1: Maximum version number of extension files of type SU is 1
>
>> inp
> AIPS 1: FITTP: Task to write maps or uv data on a FITS-tape.
> AIPS 1: Adverbs Values Comments
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: DOALL -1 1 (T) do all, -1 (F) do 1st.
> AIPS 1: INNAME 'IYA09_S' Image name(name). Standard.
> AIPS 1: INCLASS 'UV' Image name(class). Standard.
> AIPS 1: INSEQ 1 Image name(seq. #).
> AIPS 1: INDISK 1 Disk drive #. 0 => any
> AIPS 1: INTYPE 'UV' File type. Blank => any
> AIPS 1: OUTTAPE 1 Tape drive #. 0 => 1
> AIPS 1: OUTFILE 'TEMP:IYA09_S.UV'
> AIPS 1: Disk file name
> AIPS 1: DOEOT 1 1 (T) go to end tape first.
> AIPS 1: DOSTOKES -1 > 0 => Convert to true
> AIPS 1: Stokes (Formal I)
> AIPS 1: ******** Changed 12/2/93 ***********************
> AIPS 1: DONEWTAB 1 >0 standard FITS binary
> AIPS 1: tables (15JAN94 or later)
> AIPS 1: <= 0 AIPS FITS binary tables
> AIPS 1: ************************************************
> AIPS 1: FORMAT 0 tape format: 0 => 3.
> AIPS 1: 1 for 16-bit integer
> AIPS 1: 2 for 32-bit integer
> AIPS 1: 3 for 32-bit IEEE
> AIPS 1: floating point
> AIPS 1: BLOCKING 10 tape blocking factor (<= 10)
> AIPS 1: CATNO 0 0 Limit range of catalog number
>
>> go fittp
>
> Then, on DEC11:
>
>> inp
> AIPS 1: FITLD: Task to store an image or UV data from a FITS tape
> AIPS 1: Adverbs Values Comments
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: INTAPE 1 Input tape drive # (0 => 1)
> AIPS 1: NFILES 0 # of files to advance on tape
> AIPS 1: DATAIN 'TEMP:IYA09_S.UV' Disk file name
> AIPS 1: OUTNAME 'IYA09_S' File name (name)
> AIPS 1: OUTCLASS 'UV' File name (class)
> AIPS 1: OUTSEQ 0 File name (seq. #)
> AIPS 1: 0 => highest unique number
> AIPS 1: => matching (on VLBA)
> AIPS 1: -1 => FITS tape value
> AIPS 1: OUTDISK 1 Disk drive # (0 => any)
> AIPS 1: OPTYPE ' ' Type of data to load,
> AIPS 1: ' ' => all types
> AIPS 1: 'UV' => UV data
> AIPS 1: 'IM' => images
> AIPS 1: NCOUNT 0 Number of files to load.
> AIPS 1: DOTABLE 1 True (1.0) means load tables
> AIPS 1: for images.
> AIPS 1: DOUVCOMP 0 >0 => compressed data (FITS)
> AIPS 1: DOCONCAT 0 >0 -> if VLBA correlator data
> AIPS 1: append data to existing
> AIPS 1: files, or if no appropriate
> AIPS 1: files exist create a new file
> AIPS 1: and append all data to that
> AIPS 1: file.
> AIPS 1: *****************************
> AIPS 1: Following adverbs are useful
> AIPS 1: only when reading VLBA
> AIPS 1: distribution tapes.
> AIPS 1: *****************************
> AIPS 1: CLINT 0.5 CL entry interval (min)
> AIPS 1: 0 => 1 minute.
> AIPS 1: SOURCES *all ' ' Source list to accept.
> AIPS 1: QUAL -1 Source qualifier -1=>all
> AIPS 1: TIMERANG *all 0 Timerange selected
> AIPS 1: BCHAN 1 Lowest spectral channel
> AIPS 1: number to select in each IF.
> AIPS 1: 0=>1
> AIPS 1: ECHAN 0 Highest spectral channel
> AIPS 1: number to select in each IF.
> AIPS 1: 0=>highest
> AIPS 1: BIF 1 Lowest IF number 0=>1
> AIPS 1: EIF 4 Highest IF number 0=>all
> AIPS 1: DIGICOR -1 Controls the application of
> AIPS 1: the VLBA correlator's digital
> AIPS 1: correction.
> AIPS 1: 0 or 1 => perform corrections
> AIPS 1: 2 => perform cross-power
> AIPS 1: corrections, do total-power
> AIPS 1: only if zero padding used.
> AIPS 1: -1 => do NOT perform
> AIPS 1: corrections.
> AIPS 1: SEE HELP for more details
> AIPS 1: SELBAND -1 Bandwidth to select (kHz)
> AIPS 1: SELFREQ -1 Frequency to select (MHz).
> AIPS 1: See HELP.
> AIPS 1: FQTOL -1 Frequency tolerance assigned
> AIPS 1: to SELFREQ, FITLD will select
> AIPS 1: data with freq. of SELFREQ
> AIPS 1: +/- FQTOL.
> AIPS 1: The unit of FQTOL is kHz.
> AIPS 1: <=0 => 10 kHz
> AIPS 1: WTTHRESH 0 Flagging threshold based on
> AIPS 1: weights. See HELP.
> AIPS 1: 0 => no flagging
> AIPS 1: 1 => will flag ALL data
> AIPS 1: OPCODE ' ' Specify if wish to keep VT
> AIPS 1: (VLBA Tape Statistics) table
> AIPS 1: and other tables. See HELP.
> AIPS 1: ' ' => do not save tables
> AIPS 1: NPIECE 0 Maximum uv table piece to
> AIPS 1: load (ignored for tape unless
> AIPS 1: NCOUNT = 1)
> AIPS 1: ERROR -1 >= 2 -> do not use AIPS
> AIPS 1: history in the FITS file
> AIPS 1: ANTNAME *all ' ' List antenna station names in
> AIPS 1: desired order (IDI data only)
>
> Corresponding AIPS messages:
>
> FITLD1: localhost 31DEC11 TST: Cpu= 0.0 Real= 0
> FITLD1: Task FITLD (release of 31DEC11) begins
> FITLD1: Found MULTI observed on 18-NOV-2009
> FITLD1: Create IYA09_S .UV . 1 (UV) on disk 1 cno 180
> FITLD1: UVFDAT: counts for cross-correlations
> FITLD1: UVFDAT: flagged 4598 valid spectra samples which were pure 0
> FITLD1: UVFDAT: flagged 0.024 percent valid samples which were pure 0
> FITLD1: UVFDAT: found 2011630 spectra samples already flagged
> FITLD1: UVFDAT: found 10.384 percent samples already flagged
> FITLD1: Image=MULTI (UV) Filename=IYA09_S .UV . 1
> FITLD1: Telescope=VLBI Receiver=MKIV
> FITLD1: Observer=Bonn Cor User #= 4
> FITLD1: Observ. date=18-NOV-2009 Map date=01-JUL-2011
> FITLD1: # visibilities 4842987 Sort order TB
> FITLD1: Rand axes: UU-L-SIN VV-L-SIN WW-L-SIN TIME1 BASELINE
> FITLD1: SOURCE FREQSEL INTTIM
> FITLD1: ----------------------------------------------------------------
> FITLD1: Type Pixels Coord value at Pixel Coord incr Rotat
> FITLD1: COMPLEX 3 1.0000000E+00 1.00 1.0000000E+00 0.00
> FITLD1: FREQ 32 2.2329900E+09 1.00 2.5000000E+05 0.00
> FITLD1: STOKES 1 -1.0000000E+00 1.00 -1.0000000E+00 0.00
> FITLD1: IF 4 1.0000000E+00 1.00 1.0000000E+00 0.00
> FITLD1: RA 1 00 00 00.000 1.00 3600.000 0.00
> FITLD1: DEC 1 00 00 00.000 1.00 3600.000 0.00
> FITLD1: ----------------------------------------------------------------
> FITLD1: Coordinate equinox 2000.00
> FITLD1: Maximum version number of extension files of type HI is 1
> FITLD1: Maximum version number of extension files of type FQ is 1
> FITLD1: Maximum version number of extension files of type BS is 1
> FITLD1: Maximum version number of extension files of type AN is 1
> FITLD1: Maximum version number of extension files of type CL is 1
> FITLD1: Maximum version number of extension files of type SU is 1
> FITLD1: Appears to have ended successfully
> FITLD1: localhos 31DEC11 TST: Cpu= 73.0 Real= 476 IO= 7555
>
> First question: Why this uv-data file has a COMPLEX value of 3 ? The original value is 1.
>
>
> My uv-data file is TB-sorted, but with no index table. If I use INDXR on DEC07, everything runs smoothly:
>
> VANOIS> INDXR1: Task INDXR (release of 31DEC07) begins
> VANOIS> INDXR1: A CL table already exists: a new one will not be created
> VANOIS> INDXR1: Appears to have ended successfully
> VANOIS> INDXR1: vanoise 31DEC07 TST: Cpu= 72.8 Real= 8
>
> But if I do it on DEC11, something seems to go wrong:
>
>> inp
> AIPS 1: INDXR Task to index a uv data base.
> AIPS 1: Adverbs Values Comments
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: INNAME 'IYA09_S' Input UV file name (name)
> AIPS 1: INCLASS 'UV' Input UV file name (class)
> AIPS 1: INSEQ 1 Input UV file name (seq. #)
> AIPS 1: INDISK 1 Input UV file disk unit #
> AIPS 1: INFILE *all ' ' Input file for forced scan
> AIPS 1: breaks
> AIPS 1: PRTLEV 0 Print level
> AIPS 1: CPARM 0 0 1=> max. time gap (min).
> AIPS 1: 0.5 1 0 => use self-adaptive
> AIPS 1: 1 *rest 0 method, like other tasks
> AIPS 1: 2=> max scan length (min).
> AIPS 1: 0 => 60 min
> AIPS 1: 3=> CL/CS entry interval in
> AIPS 1: minutes. 0 => 5 min,
> AIPS 1: < 0 => don't create a
> AIPS 1: new table.
> AIPS 1: 4=> VLBA only: recalculate
> AIPS 1: CL entry group delays
> AIPS 1: using IM table data.
> AIPS 1: 0 => No recalculation,
> AIPS 1: 1 => Recalculate delays.
> AIPS 1: 5=> VLBA only: recalculate
> AIPS 1: CL entry atmospheric
> AIPS 1: group delays and clock
> AIPS 1: offsets using MC table
> AIPS 1: data.
> AIPS 1: 0 => No recalculation
> AIPS 1: 1 => Recalculate delays
> AIPS 1: 6=> single-dish only: maximum
> AIPS 1: "antenna"/"beam" number
> AIPS 1: in data set (if no AN)
> AIPS 1: BPARM *all 0 VLA and EVLA ONLY: Opacity
> AIPS 1: and Gain-curve control
> AIPS 1: (see help)
> AIPS 1: CALIN ' '
> AIPS 1: Antenna gains file
>
> macpco> INDXR1: Task INDXR (release of 31DEC11) begins
> macpco> INDXR1: ANREFM: FOUND 1 LEFT AND 2 RIGHT-HAND By VALUES
> macpco> INDXR1: TABINI: SOME ERROR IN LIST OF TABLE DATA TYPES
> macpco> INDXR1: ANINI: ERROR 1 FROM TABINI OPCODE = WRIT
> macpco> INDXR1: ANINI: ERROR 1 INITIALIZING ANTENNA TABLE FOR WRIT
> macpco> INDXR1: ANTINI: ERROR 1 RETURNED BY ANREFM
> macpco> INDXR1: PROBLEM WITH TABLE OBJECT
> macpco> INDXR1: Purports to die of UNNATURAL causes
> macpco> INDXR1: localhost 31DEC11 TST: Cpu= 0.0 Real= 1
>
> Second question: Do you know where this error could come from ? Maybe something related to the antenna table (just by looking at the function names) ?
> The 2 AN tables seems to be identical...
>
> Best regards,
> Arnaud
>
>
> _______________________________________________
> Daip mailing list
> Daip at listmgr.cv.nrao.edu
> http://listmgr.cv.nrao.edu/mailman/listinfo/daip
Your first question:
1. FITLD has an adverb DOUVCOMP which you set to zero so it did not do
compressed data. That should not be a problem except the file is larger.
2. The antenna file issue is more serious and INDXR may be more
vulnerable to issues with it than most tasks. Can you do a PRTAB on the
AN file(s) to a text file and e-mail them to me? It is clear that the
(new) code thinks it recognizes 3 of the antennas in the file and thinks
that they are using different conventions for the handedness of the
coordinates. I will need to see the full details which PRTAB will show
me. Perhaps the code does not really know about the array or antennas
you are using....
ERic Greisen
More information about the Daip
mailing list