Opened 15 years ago
Closed 15 years ago
#88 closed defect (fixed)
CABB: fits output of UV data does not read into AIPS properly
| Reported by: | Neeraj Gupta | Owned by: | MarkCalabretta |
|---|---|---|---|
| Priority: | major | Milestone: | |
| Component: | MIRIAD - main line | Version: | |
| Keywords: | fits, cabb | Cc: | |
| Estimated Number of Hours: | 0 | Add Hours to Ticket: | 0 |
| Billable?: | yes | Total Hours: | 0 |
Description
When I try to load CABB data into AIPS I get strange behaviour.
AIPS has problems with the RPFITS files, so I tried to use miriad 'atlod' and write the uvdata with 'fits'. I am selecting a subset of channels because of another problem within AIPS (UV data buffer sizes).
When I try to load the fits file into AIPS I get a UV data set with no SU and no FQ table.
The source data was loaded from RPFITS as follows
Task: atlod in = 2009-07-20_1440.C2098 out = 1.uv ifsel = restfreq = options = nfiles = nscans =
I wrote a fits file as follows:
fits% pwd /data/SACHIN_1/gup017/GOJKO fits% inp Task: fits in = 1.uv op = uvout out = ATCA100CH.FITS line = channel,100,1,1,1 region = select = stokes = options = nocal velocity =
I load the file into AIPS like this:
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 'G:ATCA100CH.FITS' Disk file name AIPS 1: OUTNAME 'ATCA100CH' File name (name) AIPS 1: OUTCLASS ' ' 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 1 >0 => compressed data (FITS) AIPS 1: DOCONCAT -1 >0 -> if VLBA or uv table, 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 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 0 Lowest IF number 0=>1 AIPS 1: EIF 0 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 0 >= 2 -> do not use AIPS AIPS 1: history in the FITS file
and then IMHEAD gives me:
AIPS 1: IMHEADER: Verb to list the image header AIPS 1: Adverbs Values Comments AIPS 1: ---------------------------------------------------------------- AIPS 1: USERID 0 User ID. 0 => current user AIPS 1: 32000 => any user AIPS 1: INNAME 'ATCA100CH' Image name(name). AIPS 1: INCLASS 'UVDATA' Image name(class). AIPS 1: INSEQ 1 Image name(seq. #). AIPS 1: INDISK 1 Disk drive #. 0 => any go imhead AIPS 1: Image=1934-638 (UV) Filename=ATCA100CH .UVDATA. 1 AIPS 1: Telescope=ATCA Receiver=ATCA AIPS 1: Observer=observer User #= 76 AIPS 1: Observ. date=20-JUL-2009 Map date=23-JUL-2009 AIPS 1: # visibilities 1944 Sort order AIPS 1: Rand axes: UU-L-NCP VV-L-NCP WW-L-NCP BASELINE TIME1 AIPS 1: 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: STOKES 4 -5.0000000E+00 1.00 -1.0000000E+00 0.00 AIPS 1: FREQ 100 2.7739999E+09 1.00 -1.0000000E+06 0.00 AIPS 1: RA 1 19 39 25.026 1.00 3600.000 0.00 AIPS 1: DEC 1 -63 42 45.630 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 AN is 1
What I should be seeing is the following tables:
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 SU 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 AIPS 1: Maximum version number of extension files of type CL is 1
Also I am not seeing the 'IF' line in the header. This should be shown even if there is only one IF in the data set.
Note: See
TracTickets for help on using
tickets.
This ticket relates to processing of CABB data using AIPS. AIPS is not supported for CABB data processing, please use Miriad.