[daip] DOFARS problem

[Eric Greisen]

Sayali Avachat wrote:
> Hello Dr. Greisen, 
> 
> Thanks for your reply.
> I rectified my mistake of not using a cube in DOFARS, even after that I 
> am getting same error as before.
> 
> windu.> FARS 1: ZTXOP2: using translated file name =
> windu.> FARS 1: ZTXOP2: /usr/AIPS/FITS/DOFARS.WEIGHTS
> windu.> FARS 1: ZTXOP2: FILE TO BE READ DOES NOT EXIST
> windu.> FARS 1: ERROR
> 
> Here is my corrected RSPEC input and header if you want to have a look 
> at it again..
> 
> AIPS 1: RSPEC : Task to plot spectrum of the rms of an image
> AIPS 1: Adverbs     Values                 Comments
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: USERID        0                    User ID.  ignored
> AIPS 1: INNAME     'M87_03_AXUK'           Image name (name)
> AIPS 1: INCLASS    'FQBEQ'                 Image name (class)
> AIPS 1: INSEQ         1                    Image name (seq. #)
> AIPS 1: INDISK        1                    Disk drive #
> AIPS 1: BLC        *all 0                  Bottom left corner 0->use all
> AIPS 1:                                      (1,2) -> exclude area
> AIPS 1: TRC        *all 0                  Top right corner 0 -> use all
> AIPS 1:                                      (1,2) -> exclude area
> AIPS 1: OPTYPE     ' '                     'MEDI' median window
> AIPS 1:                                    'LOW' use tight limits in rms
> AIPS 1:                                       else use less tight limits
> AIPS 1: OPCODE     'WEIT'                  '    ' => RMS
> AIPS 1:                                    'WEIT' => WEIGHT=1/(RMS**2)
> AIPS 1: DOOUTPUT     -1                    > 0 => write S/N data cube
> AIPS 1: OUTNAME    'M87_03_AXUK'           Output name(name).
> AIPS 1: OUTCLASS   'FARSAM'                Output name(class).
> AIPS 1: OUTSEQ        0                    Output name(seq. #).
> AIPS 1:                                      0 => highest unique
> AIPS 1: OUTDISK       1                    Output image disk drive #
> AIPS 1:                                      0 => highest with room
> AIPS 1: PIXRANGE      0           0        Range of intensities to plot
> AIPS 1: ZINC          1                    Increment on freq axis.
> AIPS 1:                                    < 0 => invert Z axis in plot
> AIPS 1: SMOOTH     *all 0                  Frequency smoothing function
> AIPS 1: LTYPE         3                    Type of labeling: 1 border,
> AIPS 1:                                    2 no ticks, 3 standard, 4 rel
> AIPS 1:                                    to center, 5 rel to subim cen
> AIPS 1:                                    6 pixels, 7-10 as 3-6 with
> AIPS 1:                                    only tick labels
> AIPS 1:                                    <0 -> no date/time
> AIPS 1: DOCENTER      1                    > 0 => plot with box-like
> AIPS 1:                                    lines, else plot line to
> AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
> #
> AIPS 1:                                    points
> AIPS 1: DOTV         -1                    > 0 Do plot on the TV, else
> AIPS 1:                                    make a plot file
> AIPS 1: GRCHAN        0                    Graphics channel 0 => 1.
> AIPS 1: DOPRINT       0                    > 0 => use the terminal, 0 =>
> AIPS 1:                                    no output, < 0 => printer or
> AIPS 1:                                    file, > 72 => terminal width
> AIPS 1: OUTPRINT   'FITS:DOFARS.WEIGHTS'
> AIPS 1:                                    Printer disk file to save
> AIPS 1: DOSLICE      -1                    > 0 => save as SLice file
> AIPS 1: XYRATIO       0                    Stretch the X axis wrt Y axis
> AIPS 1:                                    0 -> 1.33
> AIPS 1: Image=M87       (MA)         Filename=M87_03_AXUK .FQBEQ .   1
> AIPS 1: Telescope=VLA                Receiver=VLA
> AIPS 1: Observer=AH822               User #= 1306
> AIPS 1: Observ. date=02-JUN-2003     Map date=18-OCT-2012
> AIPS 1: Minimum=-6.45825407E-03      Maximum= 2.40508169E-02 JY/BEAM
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: Type    Pixels   Coord value     at Pixel     Coord incr   Rotat
> AIPS 1: RA---SIN  2048    12 30 49.000     776.15      -0.025000    0.00
> AIPS 1: DEC--SIN  1024    12 23 30.000     434.98       0.025000    0.00
> AIPS 1: FQID         3   1.0000000E+00       1.00  1.0000000E+00    0.00
> AIPS 1: STOKES       1   2.0000000E+00       1.00  1.0000000E+00    0.00
> AIPS 1: FREQ         1   8.4601000E+09       1.00  6.4797998E+09    0.00
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: Coordinate equinox 2000.00
> AIPS 1: Map type=NORMAL              Number of iterations=    1824
> AIPS 1: Conv size= 0.22968 X 0.21635   Position angle= -26.62
> AIPS 1: Rest freq      0.000         Vel type: OPTICAL wrt YOU
> AIPS 1: Alt ref. value  0.00000E+00  wrt pixel    0.00
> AIPS 1: Maximum version number of extension files of type FQ is   1
> AIPS 1: Maximum version number of extension files of type CG 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 CC is   1
> AIPS 1: Maximum version number of extension files of type PL is   1
> 
> What is wrong this time?

RSPEC did not do what you want because DOPRINT is 0.  What version of 
aips are you using>  If 31DEc12, is it up to date (recent MNJ)?  Have 
you entered RUN DOFARS recently?  The current DOFARS proc does set 
DOPRINT appropriately.

Now I need to tell you that starting with only 3 frequency samples is 
not likely to work well with FARS.  There is a task called TARS which 
will let you study this - make a text file listing your 3 frequencies 
and calling for a single RM model and then run TARS with the adverbs you 
would give FARS.  It will write a text file but TARPL will make a plot 
if you prefer from that text file.

Eric Greisen