[daip] FLAGR not checking everything ?

Lorant Sjouwerman lsjouwer at nrao.edu
Mon Jun 20 17:52:50 EDT 2005


Eric,

I'm tying to get FLAGR working with mostly defaults.
It seems that FLAGR is not checking/testing (all) stuff in 
LL and/or IF2. Like it is not making a switch in LL to IF2.

Maybe it is because of all the parameter complications, but 
e.g. forcing the amplitude clip CPARM(4) does flag RR, but not 
LL/IF2, not even if I make an LL or IF2 only data set. It is not 
in the FG-table. WIPER clearly shows there are LL-leftovers over
CPARM(4).. I couldn't really find it in the code (where are e.g.
OPUT and OUVGET?)

Could you please have a look at this when you have time? 
Data is at /home/euro/ERIC.FITS


Regards,

Loránt Sjouwerman - Scientific Services - lsjouwerman at nrao.edu
--------------------------------------------------------------
c/o NRAO Array Operations Center       Phone:  +1-505-835-7332
    P.O. Box 0 (1003 Lopezville Rd)    Switch: +1-505-835-7000
    Socorro NM 87801                   Fax:    +1-505-835-7027



AIPS 1: FLAGR    Use rms and amplitudes in uv data set to edit the data
AIPS 1: Adverbs     Values                 Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME     '00427+41161 '          Input UV data
AIPS 1: INCLASS    'UVCOP '                Input UV data (class)
AIPS 1: INSEQ         1                    Input UV data (seq. #)
AIPS 1: INDISK        1                    Input UV data disk drive #
AIPS 1: SOURCES    *all ' '                Source (pointings) list
AIPS 1: QUAL         -1                    Source qualifier -1=>all
AIPS 1: CALCODE    '    '                  Calibrator code '    '=>all
AIPS 1: TIMERANG   *all 0                  Time range to process.
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           0                    Lowest IF number 0=>all
AIPS 1: EIF           0                    Highest IF number 0=>all
AIPS 1: BCHAN         1                    Lowest channel number 0=>all
AIPS 1: ECHAN         0                    Highest channel number 0=>all
AIPS 1: SUBARRAY      0                    Subarray, 0=>all
AIPS 1: DOCALIB      -1                    Calibrate data?
AIPS 1:                                    = 2 calibrate weights
AIPS 1: GAINUSE       0                    CL/SN table to apply
AIPS 1: DOPOL        -1                    If >0 correct polarization.
AIPS 1: BLVER        -1                    BL table to apply.
AIPS 1: FLAGVER       1                    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        -1                    Bandpass table version
AIPS 1: UVRANGE       0           0        Min & max baseline (klambda)
AIPS 1: SOLINT        9.8                  Time interval (sec)
AIPS 1: SCANLENG      0                    Scan length (sec)
AIPS 1: REFANT        0                    Reference antenna for VDIF
AIPS 1: OPTYPE     '    '                  Type of operation on the rms
AIPS 1:                                    'RFI ' 'HIST' 'ANTE' 'GAIN'
AIPS 1:                                    'VDIF' 'VRFI'
AIPS 1:                                    'TIME' is default
AIPS 1: DOCAT        -1                    > 0 => keep the XX file
AIPS 1: VECTOR       -1                    > 0 -> vector averaging
AIPS 1:                                    <= 0 -> use amplitudes only
AIPS 1: BPARM      *all 0                  (1) Fraction allowable ??
AIPS 1:                                        correlators
AIPS 1:                                    (2) Fraction allowable ??
AIPS 1:                                        baselines
AIPS 1:                                    (3) rms cutoff (Jy) for ???
AIPS 1:                                    (4) Amp. factor for cutoff
AIPS 1:                                    (5) rms cutoff (Jy) for bad
AIPS 1:                                    (6) Amp. factor for cutoff
AIPS 1:                                    (7) Min hist level Jy for ???
AIPS 1:                                    (8) Min hist level Jy for bad
AIPS 1:                                    (9) > 0 save 1-sample SOLINTs
AIPS 1:                                    (10) best left <= 0.
AIPS 1: CPARM         0           0.12     (1) Question if amp<CPARM(1)
AIPS 1:               0           0.15     (2) Question if amp>CPARM(2)
AIPS 1:               0           0        (3) Flag if amp < CPARM(3)
AIPS 1:               0           0        (4) Flag if amp > CPARM(4)
AIPS 1:               5          10        (5) min amp rms allowed
AIPS 1:                                    (6) min rms rms allowed
AIPS 1:                                    (7) Question if amp >
AIPS 1:                                        CPARM(7)*rms away from
AIPS 1:                                        mean amplitude
AIPS 1:                                    (8) Question if rms >
AIPS 1:                                        CPARM(8)*(rms of rms's)
AIPS 1:                                        away from mean rms
AIPS 1:                                    (9) Flag if amp >
AIPS 1:                                        CPARM(8)*rms away from
AIPS 1:                                        mean amplitude.
AIPS 1:                                    (10) Flag if rms >
AIPS 1:                                        CPARM(10)*(rms of rms's)
AIPS 1:                                        away from mean rms.
AIPS 1: DPARM      *all 0                  (1) Question if wt<DPARM(1)
AIPS 1:                                    (2) Question if wt>DPARM(2)
AIPS 1:                                    (3) Flag if wt < DPARM(3)
AIPS 1:                                    (4) Flag if wt > DPARM(4)
AIPS 1:                                    (5) Question if closure
AIPS 1:                                        fraction > DPARM(5)
AIPS 1:                                    (6) Flag if closure error
AIPS 1:                                        fraction > DPARM(6)
AIPS 1: DOSTOKES      1                    > 0 -> flag all Stokes if
AIPS 1:                                    one is bad
AIPS 1: DOIFS        -1                    > 0 -> flag all IFs if one
AIPS 1:                                    is bad
AIPS 1: DOALL         1                    > Flag all channels, else
AIPS 1:                                    flag BCHAN - ECHAN
AIPS 1: DOCRT         1                    > 0 -> use the terminal,
AIPS 1:                                    > 72 => terminal width
AIPS 1:                                    0 -> no printing
AIPS 1:                                    else use the line printer
AIPS 1: OUTPRINT   '                                                '
AIPS 1:                                    Printer disk file to save
AIPS 1: PRTLEV        0                    > 0 => debug messages
AIPS 1: BADDISK    *all 0                  Disks to avoid for scratch.


FLAGR1: Task FLAGR  (release of 31DEC05) begins
 localhos  FLAGR(31DEC05)   ****     20-JUN-2005  15:17:21    Page    1
File = 00427+41161 .UVCOP .   1   Vol = 1    Userid =15661
Source= 00427+41   RA       =  00 42 44.33   DEC      =  41 16  8.4
Freq=  4.885100000 GHz   Ncor=  8   No. vis=   2252525   Sort order= TB

NO Calibration applied
For 'questionable' data: amp clipping range.=  0.0000   0.1200
For 'questionable' data: closure error fraction 0.6000
Max. fract. questionable corrs.= 0.510  antennas= 0.250
For 'bad' data: amp clipping range.=  0.0000   0.1500
FLAGR1: Overall average amp= 3.474E-02 +- 4.878E-03 source= 00427+41
FLAGR1: Overall average rms= 1.608E-02 +- 1.807E-03 source= 00427+41
Data flagged independent of antenna
       Time range        Source           Qual  Corr flags
 0/11:42:50 - 0/11:42:57 00427+41        :00000   FFFT
 0/11:43:10 - 0/11:43:17 00427+41        :00000   FFFT
 0/11:43:20 - 0/11:43:27 00427+41        :00000   FFFT
 0/11:43:30 - 0/11:43:37 00427+41        :00000   FFFT
 0/11:43:40 - 0/11:43:47 00427+41        :00000   FFFT
 0/11:44:30 - 0/11:44:37 00427+41        :00000   FFFT
 0/11:44:40 - 0/11:44:47 00427+41        :00000   FFFT
 0/11:44:50 - 0/11:44:57 00427+41        :00000   FFFT
 0/11:45:30 - 0/11:45:37 00427+41        :00000   FFFT
 0/11:45:40 - 0/11:45:47 00427+41        :00000   FFFT
FLAGR1: Found    10 bad times,   151 total flags
Found    10 bad times,   151 total flags
FLAGR1: Wrote     0 antenna & IF independent flags
Wrote     0 antenna & IF independent flags
FLAGR1: Wrote    10 antenna independent, IF dependent flags
Wrote    10 antenna independent, IF dependent flags
FLAGR1: Wrote   141 antenna & IF dependent flags
Wrote   141 antenna & IF dependent flags
FLAGR1: Wrote     0 antenna & IF dependent, all-times flags
Wrote     0 antenna & IF dependent, all-times flags
FLAGR1: Appears to have ended successfully
FLAGR1: localhost    31DEC05 TST: Cpu=      62.8  Real=      66
AIPS 1: Resumes

>ine'fg';box 1 4 5 6 8 9
>go prtab
PRTAB1: Task PRTAB  (release of 31DEC05) begins
 localhos  PRTAB(31DEC05)   ****     20-JUN-2005  15:24:15    Page    1
00427+41161 .UVCOP .   1  Disk= 1    FG Table version   1
Title: AIPS UV DATA FILE FLAG TABLE
Created by      FLAGR on 20-JUN-2005 15:17:21
Last written by FLAGR on 20-JUN-2005 15:17:21
Ncol   9  Nrow     151    Sort cols:
   Table format incompatable with FITS ASCII tables

COL. NO.      1       4         5        6       8                 9
     ROW   SOURCE   ANTS   TIME RANGE   IFS   PFLAGS   REASON
  NUMBER                   D/HMS
       1      0       4    08:57:00.0    2     1111    FLAGR 20-JUN-05 15:17:21
       1              0    08:57:06.7    2
       2      0       2    08:57:40.0    1     1111    FLAGR 20-JUN-05 15:17:21
       2              0    08:57:46.7    1
       3      0      15    09:01:00.0    2     1111    FLAGR 20-JUN-05 15:17:21
       3              0    09:01:06.7    2
       4      0      15    09:04:50.0    2     1111    FLAGR 20-JUN-05 15:17:21
       4              0    09:04:56.7    2
       5      0      15    09:08:50.0    2     1111    FLAGR 20-JUN-05 15:17:21
       5              0    09:08:56.7    2
       6      0      26    09:11:53.3    2     1111    FLAGR 20-JUN-05 15:17:21
       6              0    09:12:00.0    2
       7      0       4    09:12:53.3    1     1111    FLAGR 20-JUN-05 15:17:21
       7              0    09:13:00.0    1
       8      0      15    09:15:23.3    1     1111    FLAGR 20-JUN-05 15:17:21
       8              0    09:15:30.0    1
       9      0      21    09:15:53.3    2     1111    FLAGR 20-JUN-05 15:17:21
       9              0    09:16:00.0    2
      10      0      15    09:20:03.3    2     1111    FLAGR 20-JUN-05 15:17:21
      10              0    09:20:10.0    2
      11      0       6    09:20:23.3    2     1111    FLAGR 20-JUN-05 15:17:21
      11              0    09:20:30.0    2
      12      0      21    09:29:46.6    2     1111    FLAGR 20-JUN-05 15:17:21
      12              0    09:29:53.4    2
      13      0      12    09:31:56.6    1     1111    FLAGR 20-JUN-05 15:17:21
      13              0    09:32:03.4    1
etc

Wiper on RR all behaves okay
WIPER2: PLOTUV: X axis in Kilo Wavlngth        0.000E+00  1.859E+02
WIPER2: PLOTUV: Y axis in MilliJanskys         0.000E+00  1.872E+02

Wiper on LL stokes NOT:
WIPER2: PLOTUV: X axis in Kilo Wavlngth        0.000E+00  1.859E+02
WIPER2: PLOTUV: Y axis in      Janskys         0.000E+00  5.712E+00
which is way more than CPARM(4)=0.150

bif=eif=1 is okay
WIPER2: PLOTUV: X axis in Kilo Wavlngth        0.000E+00  1.859E+02
WIPER2: PLOTUV: Y axis in MilliJanskys         0.000E+00  1.701E+02

bif2=eif=2 is NOT okay in LL
WIPER2: PLOTUV: X axis in Kilo Wavlngth        0.000E+00  1.840E+02
WIPER2: PLOTUV: Y axis in      Janskys         0.000E+00  5.712E+00

same if I copy only the second IF in full pol to a new file
WIPER2: PLOTUV: X axis in Kilo Wavlngth        0.000E+00  1.840E+02
WIPER2: PLOTUV: Y axis in      Janskys         0.000E+00  5.712E+00





More information about the Daip mailing list