[daip] Fittp/Fitld impact on indxr

Arnaud Collioud Arnaud.Collioud at obs.u-bordeaux1.fr
Fri Jul 1 09:42:30 EDT 2011


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





More information about the Daip mailing list