[daip] Problem with WSRT data
Eric Greisen
egreisen at nrao.edu
Fri Apr 26 11:02:53 EDT 2013
Mpati Ramatsoku wrote:
> Dear daip ,
>
>
> I have a problem with the WSRT data, please see the details below
>
> Through one observation we alternated between two sources.
>
> Which i have renamed to my prefered name to make it easy to identify
> them please see a copy of the SU table below:
>
> Scan Source Qual Calcode Sub Timerange FrqID
> 1 3C48 : 0000 1 0/22:12:00 - 0/22:41:40 1
> 2 PP08A : 0000 1 0/22:46:30 - 0/22:46:30 1
> 3 PP09B : 0000 1 0/22:46:55 - 0/23:01:30 1
> 4 PP08A : 0000 1 0/23:01:55 - 0/23:16:30 1
> 5 PP09B : 0000 1 0/23:16:55 - 0/23:31:30 1
> 6 PP08A : 0000 1 0/23:31:55 - 0/23:46:30 1
> 7 PP09B : 0000 1 0/23:46:55 - 1/00:01:30 1
> 8 PP08A : 0000 1 1/00:01:55 - 1/00:16:30 1
> 9 PP09B : 0000 1 1/00:16:55 - 1/00:31:30 1
> 10 PP08A : 0000 1 1/00:31:55 - 1/00:46:30 1
> 11 PP09B : 0000 1 1/00:46:55 - 1/01:01:30 1
> 12 PP08A : 0000 1 1/01:01:55 - 1/01:16:30 1
> 13 PP09B : 0000 1 1/01:16:55 - 1/01:31:30 1
> 14 PP08A : 0000 1 1/01:31:55 - 1/01:46:30 1
> 15 PP09B : 0000 1 1/01:46:55 - 1/02:01:30 1
> 16 PP08A : 0000 1 1/02:01:55 - 1/02:16:30 1
> 17 PP09B : 0000 1 1/02:16:55 - 1/02:31:30 1
> 18 PP08A : 0000 1 1/02:31:55 - 1/02:46:30 1
> 19 PP09B : 0000 1 1/02:46:55 - 1/03:01:30 1
> 20 PP08A : 0000 1 1/03:01:55 - 1/03:16:30 1
> 21 PP09B : 0000 1 1/03:16:55 - 1/03:31:30 1
> 22 PP08A : 0000 1 1/03:31:55 - 1/03:46:30 1
> 23 PP09B : 0000 1 1/03:46:55 - 1/04:01:30 1
> 24 PP08A : 0000 1 1/04:01:55 - 1/04:16:30 1
> 25 PP09B : 0000 1 1/04:16:55 - 1/04:31:30 1
> .....
>
>
> I then ran "SPLAT" to separate the PP08As from the PP09B
>
> Below are my splat inputs
>
> AIPS 1: SPLAT Split/assemble the sources in single/multi source file
> AIPS 1: Adverbs Values Comments
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: INNAME 'PPOBS08_SPLT' Input UV file name (name)
> AIPS 1: INCLASS 'UVDATA' 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: SOURCES 'PP08A' Source list
> AIPS 1: *rest ' '
> AIPS 1: QUAL -1 Source qualifier -1=>all
> AIPS 1: CALCODE ' ' Calibrator code ' '=>all
> AIPS 1: TIMERANG *all 0 Time range to copy
> AIPS 1: STOKES ' ' Stokes type to pass.
> AIPS 1: Look HELP.
> AIPS 1: SELBAND -1 Bandwidth to select (kHz)
> AIPS 1: SELFREQ -1 Frequency to select (MHz)
> AIPS 1: FREQID -1 Freq. ID to select.
> AIPS 1: BIF 1 Lowest IF number 0=>all
> AIPS 1: EIF 8 Highest IF number 0=>all
> AIPS 1: BCHAN 0 Lowest channel number 0=>all
> AIPS 1: ECHAN 0 Highest channel number
> AIPS 1: SUBARRAY 0 Subarray, 0=>all
> AIPS 1: DOCALIB 1 > 0 calibrate data & weights
> AIPS 1: > 99 do NOT calibrate weights
> AIPS 1: GAINUSE 2 CL (or SN) table to apply
> AIPS 1: DOPOL -1 If >0 correct polarization.
> AIPS 1: PDVER 0 PD table to apply (DOPOL>0)
> AIPS 1: BLVER -1 BL table to apply.
> AIPS 1: FLAGVER 0 Flag table version.
> AIPS 1: 0 => highest numbered table
> AIPS 1: <0 => no flagging
> AIPS 1: DOBAND 3 If >0 apply bandpass cal.
> AIPS 1: Method used depends on value
> AIPS 1: of DOBAND (see HELP file).
> AIPS 1: BPVER 4 Bandpass table version
> AIPS 1: SMOOTH *all 0 Smoothing function. See
> AIPS 1: HELP SMOOTH for details.
> AIPS 1: OUTNAME 'PP08A' Output UV file name (name)
> AIPS 1: OUTCLASS ' ' Output UV file name (class)
> AIPS 1: OUTSEQ 0 Output UV file name (seq. #)
> AIPS 1: OUTDISK 1 Output UV file disk unit #.
> AIPS 1: DOUVCOMP -1 1 (T) => compressed data
> AIPS 1: APARM *all 0 Control information:
> AIPS 1: 1 = 1 => avg. freq. in IF
> AIPS 1: Use ICHANSEL.
> AIPS 1: CHANNEL is not used
> AIPS 1: = 2 => avg IF's also
> AIPS 1: Use ICHANSEL.
> AIPS 1: CHANNEL is not used
> AIPS 1: = 3 => average each
> AIPS 1: CHANNEL channels.
> AIPS 1: ICHANSEL is not used.
> AIPS 1: 2 = Input avg. time (sec)
> AIPS 1: 3 > 0 => Drop subarrays
> AIPS 1: 5 = 0 pass only xc data
> AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
> #
> AIPS 1: = 1 pass xc and ac data
> AIPS 1: = 2 pass only ac data
> AIPS 1: 6 > 0 add full source name
> AIPS 1: to header
> AIPS 1: 7 = 0 assemble all selected
> AIPS 1: sources in one
> AIPS 1: multiple source file.
> AIPS 1: > 0 split for single
> AIPS 1: source files.
> AIPS 1: ICHANSEL *all 0 Array of channel start, stop,
> AIPS 1: and increment numbers and IF
> AIPS 1: number to be used when
> AIPS 1: averaging in frequency.
> AIPS 1: It's used if APARM(1)=1 or 2
> AIPS 1: CHANNEL 0 Number of chans to average
> AIPS 1: together; if APARM(1) = 3
> AIPS 1: CHINC 1 Increment in input chans for
> AIPS 1: output 0 -> CHANNEL;
> AIPS 1: if APARM(1) = 3
> AIPS 1: SOLINT 0 Time of averaging in min.
> AIPS 1: 0 => no averaging
> AIPS 1: BADDISK *all 0
>
> This gives me an SU table that looks like this:
>
> kreiken LISTR(31DEC11) 129 26-APR-2013 16:06:25 Page 1
> File = PP08A .SPLAT . 1 Vol = 1 Userid = 129
> Freq = 1.408916275 GHz Ncor = 2 No. vis = 82392
> Scan summary listing
>
> Scan Source Qual Calcode Sub Timerange FrqID
> 1 PP08A : 0000 1 0/22:46:30 - 0/22:46:30 1
> 2 PP08A : 0000 1 0/23:01:55 - 0/23:16:30 1
> 3 PP08A : 0000 1 0/23:31:55 - 0/23:46:30 1
> 4 PP08A : 0000 1 1/00:01:55 - 1/00:16:30 1
> 5 PP08A : 0000 1 1/00:31:55 - 1/00:46:30 1
> 6 PP08A : 0000 1 1/01:01:55 - 1/01:16:30 1
> 7 PP08A : 0000 1 1/01:31:55 - 1/01:46:30 1
> 8 PP08A : 0000 1 1/02:01:55 - 1/02:16:30 1
> 9 PP08A : 0000 1 1/02:31:55 - 1/02:46:30 1
> 10 PP08A : 0000 1 1/03:01:55 - 1/03:16:30 1
> 11 PP08A : 0000 1 1/03:31:55 - 1/03:46:30 1
> 12 PP08A : 0000 1 1/04:01:55 - 1/04:16:30 1
> 13 PP08A : 0000 1 1/04:31:55 - 1/04:46:30 1
> 14 PP08A : 0000 1 1/05:01:55 - 1/05:16:30 1
> 15 PP08A : 0000 1 1/05:31:55 - 1/05:46:30 1
> 16 PP08A : 0000 1 1/06:01:55 - 1/06:16:30 1
> 17 PP08A : 0000 1 1/06:31:55 - 1/06:46:30 1
> 18 PP08A : 0000 1 1/07:01:55 - 1/07:16:30 1
> 19 PP08A : 0000 1 1/07:31:55 - 1/07:46:30 1
> 20 PP08A : 0000 1 1/08:01:55 - 1/08:16:30 1
> 21 PP08A : 0000 1 1/08:31:55 - 1/08:46:30 1
> 22 PP08A : 0000 1 1/09:01:55 - 1/09:16:30 1
> 23 PP08A : 0000 1 1/09:31:55 - 1/09:46:30 1
> 24 PP08A : 0000 1 1/10:01:55 - 1/10:16:30 1
> 25 PP08A : 0000 1 1/10:31:55 - 1/10:46:10 1
> Type Q to stop, just hit RETURN to continue
>
> kreiken LISTR(31DEC11) 129 26-APR-2013 16:06:26 Page 2
> File = PP08A .SPLAT . 1 Vol = 1 Userid = 129
>
> Source summary
> Velocity type = ' ' Definition = ' '
>
> ID Source Qual Calcode RA(2000.0) Dec(2000.0) No. vis
> 2 PP08A : 0000 04:53:58.3099 45:40:43.676 21
> 4 PP08A : 0000 04:53:58.3099 45:40:43.676 3510
> 6 PP08A : 0000 04:53:58.3099 45:40:43.676 3510
> 8 PP08A : 0000 04:53:58.3099 45:40:43.676 3510
> 10 PP08A : 0000 04:53:58.3099 45:40:43.676 3498
> 12 PP08A : 0000 04:53:58.3099 45:40:43.676 3498
> 14 PP08A : 0000 04:53:58.3099 45:40:43.676 3510
> 16 PP08A : 0000 04:53:58.3099 45:40:43.676 3510
> 18 PP08A : 0000 04:53:58.3099 45:40:43.676 3510
> 20 PP08A : 0000 04:53:58.3099 45:40:43.676 3378
> 22 PP08A : 0000 04:53:58.3099 45:40:43.676 3510
> 24 PP08A : 0000 04:53:58.3099 45:40:43.676 3486
> 26 PP08A : 0000 04:53:58.3099 45:40:43.676 3474
> 28 PP08A : 0000 04:53:58.3099 45:40:43.676 3510
> 30 PP08A : 0000 04:53:58.3099 45:40:43.676 3498
> 32 PP08A : 0000 04:53:58.3099 45:40:43.676 3498
> 34 PP08A : 0000 04:53:58.3099 45:40:43.676 3510
> 36 PP08A : 0000 04:53:58.3099 45:40:43.676 3330
> 38 PP08A : 0000 04:53:58.3099 45:40:43.676 3310
> 40 PP08A : 0000 04:53:58.3099 45:40:43.676 3330
> 42 PP08A : 0000 04:53:58.3099 45:40:43.676 3320
> 44 PP08A : 0000 04:53:58.3099 45:40:43.676 3277
> 46 PP08A : 0000 04:53:58.3099 45:40:43.676 3330
> 48 PP08A : 0000 04:53:58.3099 45:40:43.676 3320
> 50 PP08A : 0000 04:53:58.3099 45:40:43.676 3234
>
>
>
>
> What I noticed is that scan numbers no longer correspond to the ID
> numbers. This is a problem because I cannot run "IMAGR" as it says that
> it requires "ONE and ONLY source". Is there a way to fix this? or a way
> to create a single source file?
>
> I should mention that I also used tried "SPLIT" with the following inputs:
>
> AIPS 1: SPLIT Task to split multi-source uv data to single source
> AIPS 1: Adverbs Values Comments
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: also works on single files.
> AIPS 1: INNAME 'PPOBS08_SPLT' Input UV file name (name)
> AIPS 1: INCLASS 'UVDATA' 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: SOURCES 'PP08A' Source list
> AIPS 1: *rest ' '
> AIPS 1: QUAL -1 Source qualifier -1=>all
> AIPS 1: CALCODE ' ' Calibrator code ' '=>all
> AIPS 1: TIMERANG *all 0 Time range to copy
> AIPS 1: STOKES ' ' Stokes type to pass.
> AIPS 1: SELBAND -1 Bandwidth to select (kHz)
> AIPS 1: SELFREQ -1 Frequency to select (MHz)
> AIPS 1: FREQID -1 Freq. ID to select.
> AIPS 1: BIF 1 Lowest IF number 0=>all
> AIPS 1: EIF 8 Highest IF number 0=>all
> AIPS 1: BCHAN 0 Lowest channel number 0=>all
> AIPS 1: ECHAN 0 Highest channel number
> AIPS 1: SUBARRAY 0 Subarray, 0=>all
> AIPS 1: DOCALIB -1 > 0 calibrate data & weights
> AIPS 1: > 99 do NOT calibrate weights
> AIPS 1: GAINUSE 0 CL (or SN) table to apply
> AIPS 1: DOPOL -1 If >0 correct polarization.
> AIPS 1: PDVER 0 PD table to apply (DOPOL>0)
> AIPS 1: BLVER -1 BL table to apply.
> AIPS 1: FLAGVER 0 Flag table version
> AIPS 1: DOBAND -1 If >0 apply bandpass cal.
> AIPS 1: Method used depends on value
> AIPS 1: of DOBAND (see HELP file).
> AIPS 1: BPVER 0 Bandpass table version
> AIPS 1: SMOOTH *all 0 Smoothing function. See
> AIPS 1: HELP SMOOTH for details.
> AIPS 1: OUTCLASS ' ' Output UV file name (class)
> AIPS 1: OUTSEQ 1 Output UV file name (seq. #)
> AIPS 1: OUTDISK 1 Output UV file disk unit #.
> AIPS 1: DOUVCOMP -1 1 (T) => compressed data
> AIPS 1: APARM *all 0 Control information:
> AIPS 1: 1 = 1 => avg. freq. in IF
> AIPS 1: multi-channel out
> AIPS 1: = 2 => avg. freq. in IF
> AIPS 1: single channel out
> AIPS 1: = 3 => avg IF's also
> AIPS 1: 2 = Input avg. time (sec)
> AIPS 1: 3 > 0 => Drop subarrays
> AIPS 1: 4 > 0 => calibrate weights
> AIPS 1: 5 = 0 pass only xc data
> AIPS 1: = 1 pass xc and ac data
> AIPS 1: = 2 pass only ac data
> AIPS 1: 6 > 0 add full source name
> AIPS 1: to header
> AIPS 1: NCHAV 1 Number of chan. to average.
> AIPS 1: (used if APARM(1) = 1)
> AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
> #
> AIPS 1: <= 0 -> ALL
> AIPS 1: CHINC 1 Channel incr. between output
> AIPS 1: channels (used if APARM(1)=1)
> AIPS 1: ICHANSEL *all 0 Array of channel start, stop,
> AIPS 1: and increment numbers and IF
> AIPS 1: number to be used when
> AIPS 1: averaging in frequency.
> AIPS 1: (used if APARM(1) = 2, 3)
> AIPS 1: BADDISK *all 0 Disks to avoid for scratch
>
>
> It does not seem to create single source file, instead gives me mutliple
> versions of the scans with the following tables:
>
> AIPS 1: Image=MULTI (UV) Filename=PP08A -CAL.SPLAT . 1
> AIPS 1: Telescope=WSRT Receiver=WSRT
> AIPS 1: Observer=ProcessD User #= 129
> AIPS 1: Observ. date=15-SEP-2012 Map date=26-APR-2013
> AIPS 1: # visibilities 82392 Sort order TB
> AIPS 1: Rand axes: UU-L-SIN VV-L-SIN WW-L-SIN TIME1 BASELINE
> AIPS 1: SOURCE 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: STOKES 2 -1.0000000E+00 1.00 -1.0000000E+00 0.00
> AIPS 1: FREQ 256 1.4089163E+09 1.00 -3.9062500E+04 0.00
> AIPS 1: IF 8 1.0000000E+00 1.00 1.0000000E+00 0.00
> AIPS 1: RA 1 01 37 41.300 1.00 3600.000 0.00
> AIPS 1: DEC 1 33 09 35.082 1.00 3600.000 0.00
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: Coordinate equinox 2000.00
> AIPS 1: Rest freq 1420.406 Vel type: OPTICAL wrt YOU
> AIPS 1: Alt ref. value 3.53591E+06 wrt pixel 1.00
> AIPS 1: Maximum version number of extension files of type NX 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 HI is 1
> AIPS 1: Maximum version number of extension files of type AN is 1
>
> Is there a way to possibly fix all of this?
It would be better to show me the PRTAB output on the SU table. LISTR
interprets things and does not print them exactly like the SU or NX
tables. It would appear that every scan was assigned a different source
number in the input multi-source data set. This means the you do not
have 1 source named PP08A at 04:53:58.3099 45:40:43.676 but 25
"sources" all named the same and at the same location. That is not how
the data set should look - a source of the same pointing position should
have the same source number throughout the data set. Fortunately,
we have seen this issue before and have written a task DSORC which you
can use. It only does 10 source numbers at a time so you will have to
run it multiple times (start aparm = 4 6 8 10 12 14 16 18 20 22;
bparm=2). IMAGR will still be confused since you will still have an SU
table with 25 sources named PP08A. TABPUT (a verb) can be used to
rename either source 2 - the one to which you renumber everything - or
sources 4, 6, 8, 10, ... 50. Then a source name select in IMAGR will
find only 1 source and it will have all of the data.
Eric Greisen
More information about the Daip
mailing list