[daip] Problem with FITS and TASAV file.

R. Craig Walker cwalker at nrao.edu
Sun Nov 17 02:33:20 EST 2013


First, my issues with UVCOP with an edited TY table and with USUBA are
fixed after the recent midnight job.  Thanks.

I still have the odd behavior when I try to write a FITS file of a TASAV
file and read it back in.  Below is the sequence with the inputs for each
task and the two image headers.

First the image header of the file being sent to FITS:

AIPS 2: Image=MULTI     (UV)         Filename=BW098A_TST  .TASAV .   1
AIPS 2: Telescope=VLBA               Receiver=VLBA
AIPS 2: Observer=BW098A              User #=  604
AIPS 2: Observ. date=12-JAN-2013     Map date=14-NOV-2013
AIPS 2: # visibilities         1     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       4  -1.0000000E+00       1.00 -1.0000000E+00    0.00
AIPS 2: FREQ        64   4.3024000E+10       0.75  5.0000000E+05    0.00
AIPS 2: IF           8   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 FQ 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 CT 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 SN 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   2
AIPS 2: Maximum version number of extension files of type GC 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 MC 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 SU is   1
AIPS 2: Maximum version number of extension files of type TY is   2
AIPS 2: Maximum version number of extension files of type WX is   1


The inputs to FITTP:

>tget fittp
>inp
AIPS 2: FITTP:  Task to write maps or uv data on a FITS-tape.
AIPS 2: Adverbs     Values                 Comments
AIPS 2: ----------------------------------------------------------------
AIPS 2: DOALL        -1                    1 (T) do all, -1 (F) do 1st.
AIPS 2: INNAME     'BW098A_TST'            Image name(name).  Standard.
AIPS 2: INCLASS    'TASAV'                 Image name(class). Standard.
AIPS 2: INSEQ         1                    Image name(seq. #).
AIPS 2: INDISK        6                    Disk drive #.      0 => any
AIPS 2: INTYPE     'UV'                    File type.     Blank => any
AIPS 2: OUTTAPE       1                    Tape drive #.  0 => 1
AIPS 2: DATAOUT    'FITS:FITSTST2.FITS'
AIPS 2:                                    Disk file name
AIPS 2: DOEOT         1                    1 (T) go to end tape first.
AIPS 2: FORMAT        0                    tape format:   0 => 3.
AIPS 2:                                        1 for 16-bit integer
AIPS 2:                                        2 for 32-bit integer
AIPS 2:                                        3 for 32-bit IEEE
AIPS 2:                                          floating point
AIPS 2:                                        4 for 8-bit integer maps
AIPS 2: BLOCKING     10                    tape blocking factor (<= 10)
AIPS 2: CATNO         0           0        Limit range of catalog number

>tget fitld
>inp
AIPS 2: FITLD:  Task to store an image or UV data from a FITS tape
AIPS 2: Adverbs     Values                 Comments
AIPS 2: ----------------------------------------------------------------
AIPS 2: INTAPE        1                    Input tape drive # (0 => 1)
AIPS 2: NFILES        1                    # of files to advance on tape
AIPS 2: DATAIN     'FITS:FITSTST2.FITS'    Disk file name
AIPS 2: OUTNAME    'FITSTEST'              File name (name)
AIPS 2: OUTCLASS   'TASAV'                 File name (class)
AIPS 2: OUTSEQ        4                    File name (seq. #)
AIPS 2:                                       0 => highest unique number
AIPS 2:                                         => matching (on VLBA)
AIPS 2:                                      -1 => FITS tape value
AIPS 2: OUTDISK       6                    Disk drive # (0 => any)
AIPS 2: OPTYPE     ' '                     Type of data to load,
AIPS 2:                                    '  ' => all types
AIPS 2:                                    'UV' => UV data
AIPS 2:                                    'IM' => images
AIPS 2: NCOUNT        0                    Number of files to load.
AIPS 2: DOTABLE       1                    True (1.0) means load tables
AIPS 2:                                    for images.
AIPS 2: DOUVCOMP     -1                    >0 => compressed data (FITS)
AIPS 2: DOCONCAT     -1                    >0 -> if VLBA correlator data
AIPS 2:                                    append data to existing
AIPS 2:                                    files, or if no appropriate
AIPS 2:                                    files exist create a new file
AIPS 2:                                    and append all data to that
AIPS 2:                                    file.
AIPS 2: DOKEEP        0                    > 0 -> keep fully flagged
AIPS 2:                                    data records
AIPS 2: NPIECE        0                    Maximum uv table piece to
AIPS 2:                                    load (ignored for tape unless
AIPS 2:                                    NCOUNT = 1)
AIPS 2: ERROR        -1                    >= 2 -> do not use AIPS
AIPS 2:                                    history in the FITS file
AIPS 2:                                    *****************************

IDIFITS parameters not shown.


Header of loaded file:
>imh
AIPS 2: Image=MULTI     (UV)         Filename=FITSTEST    .TASAV .   4
AIPS 2: Telescope=VLBA               Receiver=VLBA
AIPS 2: Observer=BW098A              User #=  604
AIPS 2: Observ. date=12-JAN-2013     Map date=16-NOV-2013
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       4  -1.0000000E+00       1.00 -1.0000000E+00    0.00
AIPS 2: FREQ        64   4.3024000E+10       0.75  5.0000000E+05    0.00
AIPS 2: IF           8   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 HI 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 AT 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 AN is   1
AIPS 2: Maximum version number of extension files of type SN 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   2
AIPS 2: Maximum version number of extension files of type GC 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 MC 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 SU is   1
AIPS 2: Maximum version number of extension files of type TY is   2
AIPS 2: Maximum version number of extension files of type WX is   1


Note missing number of records and missing list of random parameters.


Now try UVCOP:

>inp
AIPS 2: UVCOP     Task to copy selected times and channels of uv data.
AIPS 2: Adverbs     Values                 Comments
AIPS 2: ----------------------------------------------------------------
AIPS 2: INNAME     'FITSTEST'              Input UV file name (name)
AIPS 2: INCLASS    'TASAV'                 Input UV file name (class)
AIPS 2: INSEQ         4                    Input UV file name (seq. #)
AIPS 2: INDISK        6                    Input UV file disk unit #
AIPS 2: OUTNAME    'FITSTEST'              Output UV file name (name)
AIPS 2: OUTCLASS   'FITCOP'                Output UV file name (class)
AIPS 2: OUTSEQ        4                    Output UV file name (seq. #)
AIPS 2: OUTDISK       6                    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         1                    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       0                    FG table to apply to data
AIPS 2:                                    <= 0 => do not apply
AIPS 2: UVCOPPRM   *all 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 => flag data with weights
AIPS 2:                                          below this limit
AIPS 2:                                     6 = 1,3 -> do NOT flag TY/SY
AIPS 2:                                       = 2,3 -> do NOT flag SN
AIPS 2:                                     7 = extend copy of tables
AIPS 2:                                         by UVCO(7) min, 0 -> 15
AIPS 2: NCOUNT        0                    If UVCOPPRM(3) = 2 then
AIPS 2:                                    NCOUNT is the approximate
AIPS 2:                                    ratio of AC to XC in file.

Reproduce the UVCOP error:

>go uvcop
UVCOP2: Task UVCOP  (release of 31DEC13) begins
UVCOP2: Create FITSTEST    .FITCOP.   4 (UV)  on disk  6  cno    5
UVCOP2: No flag table is applied this time
UVCOP2: UVINIT: LREC, NPIO =   6154      0 ILLEGAL
UVCOP2: UVINIT: OR NVIS, BUFSZ, OR BO =         02390016      1 ILLEGAL
UVCOP2: COPYUV: ERROR   2 INITING OUTPUT VIS FILE
UVCOP2: Destroyed  1 extension files of type FQ
UVCOP2: Destroyed UV image file: catno=      5 disk= 6
UVCOP2: Purports to die of UNNATURAL causes
UVCOP2: localhost    31DEC13 TST: Cpu=       0.0  Real=       0
AIPS 2: Resumes
AIPS 2: RETURN CODE      2 RECEIVED: STOPPING


I have put the FITS file on NOATAK at
/DATA/NOATAK_1/FITS/FITSTST2.FITS

Cheers,

Craig


---------------------------------------------------------------------
    R. Craig Walker            Array Operations Center
    cwalker at nrao.edu           National Radio Astronomy Observatory
    Phone  575 835 7247        P. O. Box O
    Fax    575 835 7027        Socorro NM 87801   USA
--------------------------------------------------------------web





More information about the Daip mailing list