[daip] may be a problem with new ULVSF

[Miller Goss]

Eric- I have had success with your new UVLSF in tst
to fix up the EVLA - EVLA spectra. Thanks.

But I cannot flag within UVLSF using FLUX= some value .
I always find this handy.

I enclose a test with flux=0.1 - in this case all the
data should be flagged. And the tst version does not
do so.

Here is the NEW version and this flags as I would expect :

VLSF1: Create ORIONA      .BASFIT.   1 (UV)  on disk  3  cno   99
inp
AIPS 1: Found in Version=NEW
AIPS 1: UVLSF     Least squares fit baseline and subtracts from uv data.
AIPS 1: Adverbs     Values                 Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME     'ORIONA'                Input UV file name (name)
AIPS 1: INCLASS    'LIND4'                 Input UV file name (class)
AIPS 1: INSEQ         1                    Input UV file name (seq. #)
AIPS 1: INDISK        3                    Input UV file disk unit #
AIPS 1: OUTNAME    ' '                     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       3                    Output UV file disk unit #.
AIPS 1: BCHAN         1                    Lowest channel to write
AIPS 1: ECHAN         0                    Highest channel to write
AIPS 1: ICHANSEL     25          84        Select channels to fit: NOTE
AIPS 1:               1           1        this is start,end,increment
AIPS 1:              90         100        and IF for each region
AIPS 1:               1           1         148         192
AIPS 1:               1           1         199         235
AIPS 1:               1           1        *rest 0
AIPS 1: ORDER         1                    Order of fit line (0 -> DC)
AIPS 1: DOOUTPUT      1                    > 0 => write fit baseline as
AIPS 1:                                    a continuum uv data base
AIPS 1: CHANNEL       0                    Channel of fit to be written
AIPS 1:                                    as continuum: 0 -> reference
AIPS 1: FLUX          0.1                  Flag data if residual flux in
AIPS 1:                                    any fit channel is > FLUX
AIPS 1:                                    0 -> 10**20
AIPS 1: CUTOFF        0                    Flag data if RMS of residual
AIPS 1:                                    in fit channels > CUTOFF
AIPS 1:                                    0 -> 10**20
AIPS 1: SHIFT         0           0        Shift in asec to RA, DEC
> 
IPS 1: Resumes
>type version
AIPS 1: 'NEW'
>UVLSF1: Create ORIONA      .UVLSF .   1 (UV)  on disk  3  cno   98
UVLSF1: Create ORIONA      .BASFIT.   1 (UV)  on disk  3  cno   99
inp>UVLSF1: Flagged   721294 time/IF/Stokes samples for excess residuals
UVLSF1: Copied FQ file from vol/cno/vers  3   96   1 to  3   98   1
UVLSF1: Copied AN file from vol/cno/vers  3   96   1 to  3   98   1
UVLSF1: Copied WX file from vol/cno/vers  3   96   1 to  3   98   1
UVLSF1: Copied OF file from vol/cno/vers  3   96   1 to  3   98   1
UVLSF1: Copied FQ file from vol/cno/vers  3   96   1 to  3   99   1
UVLSF1: Copied AN file from vol/cno/vers  3   96   1 to  3   99   1
UVLSF1: Copied WX file from vol/cno/vers  3   96   1 to  3   99   1
UVLSF1: Copied OF file from vol/cno/vers  3   96   1 to  3   99   1
UVLSF1: Appears to have ended successfully
UVLSF1: parkes       31DEC06 NEW: Cpu=      32.0  Real=     362


Note that all the data is flagged as expected.



VLSF1: parkes       31DEC06 NEW: Cpu= 
32.0  Real=     362
qh
AIPS 1: Image=ORIONA    (UV)         Filename=ORIONA      .LIND4 .   1
AIPS 1: # visibilities    360647     Sort order  TB
AIPS 1: --------------------------------------------------------------
LOGFILE FOR USER  321 GETTING LARGE: USE PRTMSG AND CLRMSG NOW!
AIPS 1: Type    Pixels   Coord value  at Pixel    Coord incr   Rotat
AIPS 1: COMPLEX      3    1.000000E+00       1    1.0000E+00    0.00
AIPS 1: STOKES       2    RR                 1       -1.0000    0.00
AIPS 1: FREQ       255    1.419985E+09       1*   3.0518E+03    0.00
AIPS 1: IF           1    1.000000E+00       1    1.0000E+00    0.00
AIPS 1: RA           1    05 35 17.458       1      3600.000    0.00
AIPS 1: DEC          1   -05 23 36.785       1      3600.000    0.00
AIPS 1: --------------------------------------------------------------
AIPS 1: Max ext file #: FQ   1    HI   1    AN   1    WX   1    OF   1
AIPS 1: * => NOT the formal reference pixel on this axis
>


Next the exact same data and the TST version ( the one
we use for ORDER=3 , 4 etc ).

goUVLSF1: Create ORIONA      .UVLSF .   2 (UV)  on disk  3  cno  100

UVLSF1: Create ORIONA      .BASFIT.   2 (UV)  on disk  3  cno  101
POS


in TST is the version. and order is STILL one- just like I ran
in NEW>

Here only one point is flagged for flux =0.1 ( I would normally use s=16 
JY)

>UVLSF1: Flagged        1 time/IF/Stokes samples for excess residuals
UVLSF1: Copied AN file from vol/cno/vers  3   96   1 to  3  100   1
UVLSF1: Copied OF file from vol/cno/vers  3   96   1 to  3  100   1
UVLSF1: Updating tables for IF/FREQID/channel selection
UVLSF1: Copied NX file from vol/cno/vers  3  100   1 to  3  101   1
UVLSF1: Copied AN file from vol/cno/vers  3   96   1 to  3  101   1
UVLSF1: Copied OF file from vol/cno/vers  3   96   1 to  3  101   1
UVLSF1: Updating tables for IF/FREQID/channel selection
UVLSF1: Appears to have ended successfully
UVLSF1: parkes       31DEC07 TST: Cpu=      72.5  Real=     532




Thanks.
miller