[daip] LINIMAGE output

[Eric Greisen]

Emmanuel Momjian wrote:
> Below are:
> 1) The RSPEC output when it cannot produce a spectrum
> 2) The RSPEC inputs
> 3) The FQ table where the FRQSEL entries of the first channels are 
> marked as 0. Could this be the issue? (note that the first 5 channels 
> are not flagged; they have images in the cube).
> 
> 
> 
> In RSPEC, I get  (similar messages show up in ISPEC)
> 
> RSPEC2: Task RSPEC  (release of 31DEC12) begins
> RSPEC2: DISK/SLOT  1    27 FQ REC  1 FQID POINT   0 OUTSIDE AXIS
> RSPEC2: FQID AXIS NOT CHANGED TO FREQUENCIES
> RSPEC2: Purports to die of UNNATURAL causes
> RSPEC2: evla         31DEC12 TST: Cpu=       0.0  Real=       0
> AIPS 2: Resumes
> 
> The inputs I have in RSPEC are:
> AIPS 2: RSPEC : Task to plot spectrum of the rms of an image
> AIPS 2: Adverbs     Values                 Comments
> AIPS 2: ----------------------------------------------------------------
> AIPS 2: USERID        0                    User ID.  ignored
> AIPS 2: INNAME     'TEST'                  Image name (name)
> AIPS 2: INCLASS    'FQUBE'                 Image name (class)
> AIPS 2: INSEQ         1                    Image name (seq. #)
> AIPS 2: INDISK        1                    Disk drive #
> AIPS 2: BLC         231         324        Bottom left corner 0->use all
> AIPS 2:               1           1          (1,2) -> exclude area
> AIPS 2:               1           1           1
> AIPS 2: TRC         699         745        Top right corner 0 -> use all
> AIPS 2:             512           1          (1,2) -> exclude area
> AIPS 2:               1           1           1
> AIPS 2: OPTYPE     ' '                     'MEDI' median window
> AIPS 2:                                    'LOW' use tight limits in rms
> AIPS 2:                                       else use less tight limits
> AIPS 2: OPCODE     ' '                     '    ' => RMS
> AIPS 2:                                    'WEIT' => WEIGHT=1/(RMS**2)
> AIPS 2: DOOUTPUT     -1                    > 0 => write S/N data cube
> AIPS 2: OUTNAME    ' '                     Output name(name).
> AIPS 2: OUTCLASS   ' '                     Output name(class).
> AIPS 2: OUTSEQ        0                    Output name(seq. #).
> AIPS 2:                                      0 => highest unique
> AIPS 2: OUTDISK       1                    Output image disk drive #
> AIPS 2:                                      0 => highest with room
> AIPS 2: PIXRANGE      0           0        Range of intensities to plot
> AIPS 2: ZINC          1                    Increment on freq axis.
> AIPS 2:                                    < 0 => invert Z axis in plot
> AIPS 2: SMOOTH     *all 0                  Frequency smoothing function
> AIPS 2: LTYPE         3                    Type of labeling: 1 border,
> AIPS 2:                                    2 no ticks, 3 standard, 4 rel
> AIPS 2:                                    to center, 5 rel to subim cen
> AIPS 2:                                    6 pixels, 7-10 as 3-6 with
> AIPS 2:                                    only tick labels
> AIPS 2:                                    <0 -> no date/time
> AIPS 2: DOCENTER      1                    > 0 => plot with box-like
> AIPS 2:                                    lines, else plot line to
> AIPS 2:                                    points
> AIPS 2: DOTV          1                    > 0 Do plot on the TV, else
> AIPS 2:                                    make a plot file
> AIPS 2: GRCHAN        0                    Graphics channel 0 => 1.
> AIPS 2: DOPRINT       0                    > 0 => use the terminal, 0 =>
> AIPS 2:                                    no output, < 0 => printer or
> AIPS 2:                                    file, > 72 => terminal width
> AIPS 2: OUTPRINT   ' '
> AIPS 2:                                    Printer disk file to save
> AIPS 2: DOSLICE       1                    > 0 => save as SLice file
> AIPS 2: XYRATIO       0                    Stretch the X axis wrt Y axis
> AIPS 2:                                    0 -> 1.33
> 
> 
> 
> The FQ table has 0 values in the first few channels (these channels were 
> flagged). The first 5 channels in the cube are not flagged.
> 
> COL. NO.      1            2           3          4          5          6
>      ROW  FRQSEL     IF FREQ       CH WIDTH   TOTAL BAN  SIDEBAND   BANDCODE
>   NUMBER             HZ            HZ         HZ
>        1          0  6.013470-154  1.562E+04  1.562E+04          1
>        2          0  6.013469-154  1.356E-19  1.562E+04          1
>        3          0  3.125000D+04  1.356E-19  1.356E-19          1
>        4          0  4.687500D+04  1.562E+04  1.356E-19          1
>        5          0  6.250000D+04  1.562E+04  1.562E+04          1
>        6  538976288  7.812500D+04  1.562E+04  1.562E+04          1
>        7  538976288  9.375000D+04  1.562E+04  1.562E+04  538976288
>        8  538976288  1.093750D+05  1.562E+04  1.562E+04  538976288
>        9  538976288  1.250000D+05  1.562E+04  1.562E+04          1
>       10  538976288  1.406250D+05  1.562E+04  1.562E+04          1
>       11  538976288  1.562500D+05  1.562E+04  1.562E+04          1
>       12  538976288  1.718750D+05  1.562E+04  1.562E+04          1

This FQ table is garbage and unusable so I am glad that the ISPEC and 
RSPEC reject it.  The issue is how did this arise - I suspect FQUBE of 
being blasted in some way.

The problem is that the individual FLATNed cubes were zapped.

Do the following

cp $TST/RUN/LINIMAGE.001 $HOME
edit $HOME/LINIMAGE.001 to put a $ sign in front of the line very near 
the end that says for inseq=cubeseq:in2seq; zap; end

Then in aips say
version = 'home'
run linimage
you should see this line commented out

Then do LINIMAGE on your data

Then send me the FQ table (more than 12 lines) and the headers of each 
of the FLATNed cubes and the final LINIMAGE output cube.

Eric