[daip] Problem with WSRT data

Mpati Ramatsoku mpatir8 at gmail.com
Fri Apr 26 10:29:26 EDT 2013


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?


Thank you


Kind regards
Mpati



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://listmgr.nrao.edu/pipermail/daip/attachments/20130426/64a39598/attachment.html>


More information about the Daip mailing list