[daip] UVCOP destroys an IM table keyword

Mark Claussen mclausse at aoc.nrao.edu
Wed Jul 12 12:40:50 EDT 2000


****************************************************************
* Ian Hoffman and I have run into a problem, which appears to
* be a problem with UVCOP.  It first manifested itself when
* FXPOL failed, having a problem with the IM table:
* 
****************************************************************


FXPOL2: Fixing FQ table version    1
FXPOL2: Copied CT file from vol/cno/vers  2   10   1 to  2   13   1
FXPOL2: Copied OB file from vol/cno/vers  2   10   1 to  2   13   1
FXPOL2: Copied WX file from vol/cno/vers  2   10   1 to  2   13   1
FXPOL2: Copied NX file from vol/cno/vers  2   10   1 to  2   13   1
FXPOL2: Fixing AN table version    1
FXPOL2: Fixing CL table version    1
FXPOL2: Fixing CQ table version    1
FXPOL2: Fixing FG table version    1
FXPOL2: Fixing GC table version    1
>FXPOL2: Fixing IM table version    1
FXPOL2: TABKEY: TABLE KEYWORD RDATE    NOT FOUND
FXPOL2: IMINIT: FAILED TO READ KEYWORD VALUES (ERROR   21)
FXPOL2: OIMINT: FAILED TO OPEN TABLE (ERROR   21)
FXPOL2: FXPIM3: FAILED TO OPEN INPUT TABLE (ERROR   21)
FXPOL2: Destroyed  1 extension files of type GC
FXPOL2: Destroyed  1 extension files of type FG
FXPOL2: Destroyed  1 extension files of type CQ
FXPOL2: Destroyed  1 extension files of type CL
FXPOL2: Destroyed  1 extension files of type AN
FXPOL2: Destroyed  1 extension files of type NX
FXPOL2: Destroyed  1 extension files of type WX
FXPOL2: Destroyed  1 extension files of type OB
FXPOL2: Destroyed  1 extension files of type CT
FXPOL2: Destroyed  1 extension files of type FQ
FXPOL2: Destroyed UV image file: catno=     13 disk= 2
FXPOL2: Purports to die of UNNATURAL causes
FXPOL2: pinenut      31DEC99 TST: Cpu=     422.9  Real=     590

********************************************************************
* 
* I looked a bit more carefully at what was going on.  It turns
* out that Ian had UVCOPied the data from the original.  An
* IMHEAD on the original data follows, and after that a PRTAB
* on the IM table in the original data:
* 
********************************************************************

#imhe
AIPS 2: Image=MULTI     (UV)         Filename=BC101       .MULTTB.   1
AIPS 2: Telescope=VLBA               Receiver=VLBA
AIPS 2: Observer=BC101               User #= 2314
AIPS 2: Observ. date=01-MAY-2000     Map date=06-JUL-2000
AIPS 2: # visibilities    149889     Sort order  TB
AIPS 2: Rand axes: UU-L-SIN  VV-L-SIN  WW-L-SIN  TIME1  BASELINE
AIPS 2:            SOURCE  FREQSEL  INTTIM  GATEID  CORR-ID  WEIGHT
AIPS 2:            SCALE
AIPS 2: ----------------------------------------------------------------
AIPS 2: Type    Pixels   Coord value     at Pixel     Coord incr   Rotat
AIPS 2: COMPLEX      1   1.0000000E+00       1.00  1.0000000E+00    0.00
AIPS 2: STOKES       1  -1.0000000E+00       1.00 -1.0000000E+00    0.00
AIPS 2: FREQ      1024   1.7198500E+09       1.00  9.7656250E+02    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       0.000000    0.00
AIPS 2: DEC          1    00 00 00.000       1.00       0.000000    0.00
AIPS 2: ----------------------------------------------------------------
AIPS 2: Coordinate equinox 2000.00
AIPS 2: Maximum version number of extension files of type HI is   1
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 AT is   1
AIPS 2: Maximum version number of extension files of type IM is   1
AIPS 2: Maximum version number of extension files of type CT is   1
AIPS 2: Maximum version number of extension files of type MC is   1
AIPS 2: Maximum version number of extension files of type OB is   1
AIPS 2: Maximum version number of extension files of type GC is   1
AIPS 2: Maximum version number of extension files of type TY is   1
AIPS 2: Maximum version number of extension files of type WX is   1
AIPS 2: Maximum version number of extension files of type PC is   1
AIPS 2: Maximum version number of extension files of type CL is   1
AIPS 2: Maximum version number of extension files of type CQ is   1
AIPS 2: Maximum version number of extension files of type FG is   1
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

pinenut   PRTAB(31DEC99)   2314     12-JUL-2000  10:22:16    Page    1
BC101       .MULTTB.   1  Disk= 1    IM Table version   1
Title: INTERFEROMETER_MODEL
Created by      FITLD on 05-JUL-2000 17:01:17
Last written by FITLD on 06-JUL-2000 15:52:26
Ncol  14  Nrow     561    Sort cols:   1 (ASCEND)   4 (ASCEND)
    Table has    17 keyword-value pairs:
   OBSCODE  =  BC101
   RDATE    =  20000501
   NO_STKD  =            1
   STK_1    =           -1
   NO_BAND  =            2
   NO_CHAN  =         1024
   REF_FREQ =  1.7198500000000E+09
   CHAN_BW  =  9.7656250000000E+02
   REF_PIXL =  1.0000000000000E+00
   TABREV   =            2
   NO_POL   =            1
   GSTIA0   =  0.0000000000000E+00
   DEGPDY   =  0.0000000000000E+00
   RDATE    =  20000501
   CDATE    =  20000522
   NPOLY    =            6
   REVISION =  1.0000000000000E+00
   Table format incompatable with FITS ASCII tables


***********************************************************************
*
* One can see that the original IM table has the RDATE keyword, and
* has the proper string in it.  I then ran UVCOP; the inputs, messages,
* from UVCOP are below, as well as running PRTAB on the IM table that
* UVCOP made:
*
***********************************************************************


AIPS 2: UVCOP     Task to copy selected times and channels of uv data.
AIPS 2: Adverbs         Values            Comments
AIPS 2: ----------------------------------------------------------------
AIPS 2: INNAME     'BC101       '          Input UV file name (name)
AIPS 2: INCLASS    'MULTTB'                Input UV file name (class)
AIPS 2: INSEQ         1                    Input UV file name (seq. #)
AIPS 2: INDISK        1                    Input UV file disk unit #
AIPS 2: OUTNAME    'BC101       '          Output UV file name (name)
AIPS 2: OUTCLASS   'TSTIM '                Output UV file name (class)
AIPS 2: OUTSEQ        0                    Output UV file name (seq. #)
AIPS 2: OUTDISK       1                    Output UV file disk unit #.
AIPS 2: SELBAND      -1                    Bandwidth to select (kHz)
AIPS 2: SELFREQ      -1                    Frequency to select (MHz)
AIPS 2: FREQID       -1                    Freq. ID to select. 0=> all
AIPS 2: BCHAN         0                    First channel selected.
AIPS 2: ECHAN         0                    Last channel selected.
AIPS 2: BIF           0                    Lowest IF number 0=>all
AIPS 2: EIF           0                    Highest IF number 0=>all
AIPS 2: SOURCES    *all ' '                Source list
AIPS 2: QUAL         -1                    Source qualifier -1=>all
AIPS 2: UVRANGE       0           0        Annulus in uv plane selected
AIPS 2:                                     (thousands of wavelengths)
AIPS 2: TIMERANG   *all 0                  Timerange:
AIPS 2:                                     1-4 = start day,hr,min,sec
AIPS 2:                                     5-8 = end   day,hr,min,sec
AIPS 2: ANTENNAS   *all 0                  Antennas to copy 0=>all
AIPS 2: BASELINE   *all 0                  Baselines with ANTENNAS
AIPS 2: SUBARRAY      0                    Limit copy to one subarray
AIPS 2:                                    0 => all.
AIPS 2: FLAGVER      -1                    FG table to apply to data
AIPS 2:                                    <= 0 => do not apply
AIPS 2: UVCOPPRM      1        *rest 0     Flags:
AIPS 2:                                     1 > 0 => keep flagged data
AIPS 2:                                     2 > 0 => drop subarray info
AIPS 2:                                     3 = 0 => copy xc and ac data
AIPS 2: SELBAND      -1                    Bandwidth to select (kHz)
AIPS 2: SELFREQ      -1                    Frequency to select (MHz)
AIPS 2: FREQID       -1                    Freq. ID to select. 0=> all
AIPS 2: BCHAN         0                    First channel selected.
AIPS 2: ECHAN         0                    Last channel selected.
AIPS 2: BIF           0                    Lowest IF number 0=>all
AIPS 2: EIF           0                    Highest IF number 0=>all
AIPS 2: SOURCES    *all ' '                Source list
AIPS 2: QUAL         -1                    Source qualifier -1=>all
AIPS 2: UVRANGE       0           0        Annulus in uv plane selected
AIPS 2:                                     (thousands of wavelengths)
AIPS 2: TIMERANG   *all 0                  Timerange:
AIPS 2:                                     1-4 = start day,hr,min,sec
AIPS 2:                                     5-8 = end   day,hr,min,sec
AIPS 2: ANTENNAS   *all 0                  Antennas to copy 0=>all
AIPS 2: BASELINE   *all 0                  Baselines with ANTENNAS
AIPS 2: SUBARRAY      0                    Limit copy to one subarray
AIPS 2:                                    0 => all.
AIPS 2: FLAGVER      -1                    FG table to apply to data
AIPS 2:                                    <= 0 => do not apply
AIPS 2: UVCOPPRM      1        *rest 0     Flags:
AIPS 2:                                     1 > 0 => keep flagged data
AIPS 2:                                     2 > 0 => drop subarray info
AIPS 2:                                     3 = 0 => copy xc and ac data
AIPS 2:                                       = 1 => copy only xc data
AIPS 2:                                       = 2 => copy only ac data
AIPS 2:                                          of antennas, named in
AIPS 2:                                          antenna list. Baseline
AIPS 2:                                          list must be zero.
AIPS 2:                                     4 > 0 => report progress
AIPS 2:                                     5 => clip data with weights
AIPS 2:                                          below this limit
AIPS 2: NCOUNT        0                    If UVCOPPRM(3) = 2 then
AIPS 2:                                    NCOUNT is the approximate
AIPS 2:                                    ratio of AC to XC in file.
>

>go
UVCOP2: Task UVCOP  (release of 31DEC99) begins

AIPS 2: Resumes
>
>UVCOP2: Create BC101       .TSTIM .   1 (UV)  on disk  1  cno   16

>UVCOP2: Copied from the beginning
UVCOP2:         to  the end
UVCOP2: Copied   149889. vis records
UVCOP2: Can include fully flagged samples
UVCOP2: Copied AT file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied CT file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied OB file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied WX file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied NX file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Updating tables for IF/FREQID/channel selection
UVCOP2: Copied AN file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied CL file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied CQ file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied FG file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied GC file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied IM file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied MC file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied PC file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied SU file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Copied TY file from vol/cno/vers  1   14   1 to  1   16   1
UVCOP2: Appears to have ended successfully
UVCOP2: pinenut      31DEC99 TST: Cpu=     400.4  Real=     576


tget prtab

>indi 1;getn 16
AIPS 2: Got(1)   disk= 1  user=2314   type=UV   BC101.TSTIM.1
>docrt 1  
>go

PRTAB2: Task PRTAB  (release of 31DEC99) begins
 pinenut   PRTAB(31DEC99)   2314     12-JUL-2000  10:17:15    Page    1
BC101       .TSTIM .   1  Disk= 1    IM Table version   1
Title: INTERFEROMETER MODEL TABLE
Created by      UVCOP on 12-JUL-2000 10:12:04
Last written by UVCOP on 12-JUL-2000 10:12:04
Ncol  14  Nrow     561    Sort cols:
    Table has    12 keyword-value pairs:
   OBSCODE  =  BC101
   NO_STKD  =            1
   STK_1    =           -1
   NO_BAND  =            2
   NO_CHAN  =         1024
   REF_FREQ =  1.7198500000000E+09
   CHAN_BW  =  9.7656250000000E+02
   REF_PIXL =  1.0000000000000E+00
   NO_POL   =            1
   NPOLY    =            6
   REVISION =  1.0000000000000E+00
   TABREV   =            2
   Table format incompatable with FITS ASCII tables
 
PRTAB2: Appears to have ended successfully
PRTAB2: pinenut      31DEC99 TST: Cpu=       0.0  Real=       7
AIPS 2: Resumes



*******************************************************************
*                                                                 
* So, somehow UVCOP destroys the RDATE keyword, either by itself
* or user error (though I don't see what mistake I made), and leads
* to the original FXPOL mistake.  We got around this by copying the
* original IM table to the UVCOPied data, and then ran FXPOL.
*
*******************************************************************


Mark Claussen



More information about the Daip mailing list