[daip] AIPS error
Eric Greisen
egreisen at nrao.edu
Fri Dec 14 14:13:15 EST 2012
Sayali Avachat wrote:
> Hello Dr. Greisen,
>
> I am sorry to bother you over and over again. But I just wanted to make
> sure whatever I am doing is correct.
> I could run task FARS successfully after I tweaked the parameters a
> little as your and Dr. Owen's advise. I set the APARM(4)=1 to get
> uncleaned output and APARM(5) =2 to get Amplitude and Phase cubes.
>
> Amplitude cube looks like -
>
> AIPS 1: Image=M87 (MA) Filename=M87_03_BUKQ .FARSam. 4
> AIPS 1: Telescope=VLA Receiver=VLA
> AIPS 1: Observer=AH822 User #= 1306
> AIPS 1: Observ. date=16-NOV-2003 Map date=13-DEC-2012
> AIPS 1: Minimum= 1.15295640E-08 Maximum= 3.59934296E+02 JY/BEAM
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: Type Pixels Coord value at Pixel Coord incr Rotat
> AIPS 1: FARADROT 201 -2.0000000E+04 1.00 2.0000000E+02 0.00
> ---------------------> I wonder why this is always at pixel 1.00?
> AIPS 1: RA---SIN 2048 12 30 49.000 776.54 -0.025000 0.00
> AIPS 1: DEC--SIN 1024 12 23 30.000 435.48 0.025000 0.00
> AIPS 1: FREQ 1 1.4939900E+10 1.00 7.5202002E+09 0.00
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: Coordinate equinox 2000.00
> AIPS 1: Map type=NORMAL Number of iterations= 788
> AIPS 1: Conv size= 0.46440 X 0.41353 Position angle= 31.76
> 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 HI 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 CG is 1
>
The FARADROT axis is an arbitrary and linear axis. So long as the
reference value is correct for the chosen reference pixel everything is
fine. One could have chosen ref value 0 and ref pixel = 101 (in this
case) but this choice is fine.
> and Phase cube looks like -
>
> AIPS 1: Image=M87 (MA) Filename=M87_03_BUKQ .FARSph. 4
> AIPS 1: Telescope=VLA Receiver=VLA
> AIPS 1: Observer=AH822 User #= 1306
> AIPS 1: Observ. date=16-NOV-2003 Map date=13-DEC-2012
> AIPS 1: Minimum=-1.80000000E+02 Maximum= 1.80000000E+02 DEGREES
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: Type Pixels Coord value at Pixel Coord incr Rotat
> AIPS 1: FARADROT 201 -2.0000000E+04 1.00 2.0000000E+02 0.00
> AIPS 1: RA---SIN 2048 12 30 49.000 776.54 -0.025000 0.00
> AIPS 1: DEC--SIN 1024 12 23 30.000 435.48 0.025000 0.00
> AIPS 1: FREQ 1 1.4939900E+10 1.00 7.5202002E+09 0.00
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: Coordinate equinox 2000.00
> AIPS 1: Map type=NORMAL Number of iterations= 788
> AIPS 1: Conv size= 0.46440 X 0.41353 Position angle= 31.76
> 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 HI 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 CG is 1
>
> I next want to inspect them in SLICE, here I am not sure of what BLC and
> TRC shall I use. I tried the following input to inspect the inner part
> of the jet.
>
SLICE makes 1 dimensional slices not two dimensional. You can look
at any plane in this cube (or in the transposed cubes) using many
display tasks/verbs like TVLOD, KNTR, etc.
> AIPS 1: SLICE: Task to make a slice file from an image
> AIPS 1: Adverbs Values Comments
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: INNAME 'M87_03_BUKQ' Image name (name)
> AIPS 1: INCLASS 'FARSam' Image name (class)
> AIPS 1: INSEQ 4 Image name (seq. #)
> AIPS 1: INDISK 1 Disk drive #
> AIPS 1: BLC 0 747 Beginning of slice
> AIPS 1: 413 *rest 0
> AIPS 1: TRC 0 827 End of slice
> AIPS 1: 469 *rest 0
> AIPS 1: OUTTEXT ' '
> AIPS 1: Disk file to write slice text
> AIPS 1: output. Default (blank) = do
> AIPS 1: not write disk file.
> AIPS 1: DOKEEP 0 > 0 keep simple slices
> AIPS 1: matching input (no interp)
>
> After I ran this, I got this message (so I am sure I have set something
> wrong for BLC-TRC)
>
> solo.a> SLICE1: Task SLICE (release of 31DEC12) begins
> solo.a> SLICE1: CANNOT DO TRUE 3+ DIMENSIONS
> solo.a> SLICE1: Purports to die of UNNATURAL causes
>
> The help file for SLICE is not very clear about what BLC and TRC values
> should be exactly. SLICE is not working unless I set the them equal, but
> I think setting then equal is meaning less, since it will be just one
> point instead of a 2 dimensional slice.
>
> Next, AFARS output RM file looks like this -
>
> AIPS 1: Image=M87 (MA) Filename=M87_03_BUKQ .AFARrm. 1
> AIPS 1: Telescope=VLA Receiver=VLA
> AIPS 1: Observer=AH822 User #= 1306
> AIPS 1: Observ. date=16-NOV-2003 Map date=13-DEC-2012
> AIPS 1: Minimum=-2.00000000E+04 Maximum= 2.00000000E+04 RAD/M/M
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: Type Pixels Coord value at Pixel Coord incr Rotat
> AIPS 1: RA---SIN 2048 12 30 49.000 776.54 -0.025000 0.00
> AIPS 1: DEC--SIN 1024 12 23 30.000 435.48 0.025000 0.00
> AIPS 1: FREQ 1 1.4939900E+10 1.00 7.5202002E+09 0.00
> AIPS 1: ----------------------------------------------------------------
> AIPS 1: Coordinate equinox 2000.00
> AIPS 1: Map type=NORMAL Number of iterations= 788
> AIPS 1: Conv size= 0.46440 X 0.41353 Position angle= 31.76
> 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 HI 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 CG is 1
> >
>
> Amplitude and Phase output files have the same 3 axes as RM.
,
Yes - but what is your question? These images may be displayed with the
same tasks/verbs I just mentioned.
Eric Greisen
More information about the Daip
mailing list