<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<font size="-1">Eric - I just tried to use FARS in a different way
than before, specifically doing 1 clean component, 99% gain, and
asking for output of residuals. The amplitude and phase of
residuals was pm 1e30. I attach the headers for the input files
for DOFARS, the headers of output maps, the inputs to FARS used
by DOFARS and the prtmsg. I'm going to try changing one param at
a time and do a little diagnosis. note the additional little
problem of not finding the file at end to TRANS the output.<br>
<br>
Happy New year! lr<br>
<br>
INPUT HEADERS<br>
<br>
AIPS 2: Image=A2256 (MA) Filename=A2256ALL .QFTC6
. 1<br>
AIPS 2: Magic value blanking<br>
AIPS 2: Minimum= -1.3245 Maximum= 1.3876
MilliJY/BEAM<br>
AIPS 2:
--------------------------------------------------------------<br>
AIPS 2: Type Pixels Coord value at Pixel Coord incr
Rotat<br>
AIPS 2: RA---SIN 8000 17 03 43.500 4000 -0.400000
0.00<br>
AIPS 2: DEC--SIN 8000 78 39 59.600 4000* 0.400000
0.00<br>
AIPS 2: FQID 94 4.700000E+01 47* 1.0000E+00
0.00<br>
AIPS 2: STOKES 1 QPOL 1 1.0000
0.00<br>
AIPS 2: FREQ 1 1.016000E+09 1 1.0000E+07
0.00<br>
AIPS 2:
--------------------------------------------------------------<br>
AIPS 2: Map type=NORMAL Number of iterations= 24333<br>
AIPS 2: Conv size= 6.00 X 6.00 Position angle= 0.00<br>
AIPS 2: Max ext file #: HI 1 FQ 1<br>
AIPS 2: * => NOT the formal reference pixel on this axis<br>
<br>
AIPS 2: Image=A2256 (MA) Filename=A2256ALL .QFTC6
. 1<br>
AIPS 2: Magic value blanking<br>
AIPS 2: Minimum= -1.3245 Maximum= 1.3876
MilliJY/BEAM<br>
AIPS 2:
--------------------------------------------------------------<br>
AIPS 2: Type Pixels Coord value at Pixel Coord incr
Rotat<br>
AIPS 2: RA---SIN 8000 17 03 43.500 4000 -0.400000
0.00<br>
AIPS 2: DEC--SIN 8000 78 39 59.600 4000* 0.400000
0.00<br>
AIPS 2: FQID 94 4.700000E+01 47* 1.0000E+00
0.00<br>
AIPS 2: STOKES 1 QPOL 1 1.0000
0.00<br>
AIPS 2: FREQ 1 1.016000E+09 1 1.0000E+07
0.00<br>
AIPS 2:
--------------------------------------------------------------<br>
AIPS 2: Map type=NORMAL Number of iterations= 24333<br>
AIPS 2: Conv size= 6.00 X 6.00 Position angle= 0.00<br>
AIPS 2: Max ext file #: HI 1 FQ 1<br>
AIPS 2: * => NOT the formal reference pixel on this axis<br>
<br>
OUTPUT FILES<br>
AIPS 2: 35 2256 A2256RELCL .FARSam. 1 MA 01-JAN-2015
12:19:32<br>
AIPS 2: 36 2256 A2256RELCL .FARSph. 1 MA 01-JAN-2015
12:19:32<br>
><br>
>getn 35<br>
AIPS 2: Got(1) disk= 3 user=2256 type=MA
A2256RELCL.FARSam.1<br>
>imh<br>
AIPS 2: Image=A2256 (MA) Filename=A2256RELCL
.FARSam. 1<br>
AIPS 2: Telescope=EVLA Receiver=EVLA<br>
AIPS 2: Observer=Frazer O User #= 2256<br>
AIPS 2: Observ. date=29-JUL-2010 Map date=01-JAN-2015<br>
AIPS 2: Pixel type: FLOATING Magic value blanking<br>
AIPS 2: Minimum= 1.00000002E+30 Maximum=-1.00000002E+30
JY/BEAM<br>
AIPS 2:
----------------------------------------------------------------<br>
AIPS 2: Type Pixels Coord value at Pixel Coord incr
Rotat<br>
AIPS 2: FARADROT 301 -4.5000000E+02 1.00
3.0000000E+00 0.00<br>
AIPS 2: RA---SIN 2081 17 03 43.500 660.00
-0.400000 0.00<br>
AIPS 2: DEC--SIN 1941 78 39 60.000 341.00
0.400000 0.00<br>
AIPS 2: FREQ 1 1.0160000E+09 1.00
1.0000000E+07 0.00<br>
AIPS 2:
----------------------------------------------------------------<br>
AIPS 2: Coordinate equinox 2000.00<br>
AIPS 2: Map type=NORMAL Number of iterations= 24333<br>
AIPS 2: Conv size= 6.00 X 6.00 Position angle= 0.00<br>
AIPS 2: Rest freq 0.000 Vel type: OPTICAL wrt YOU<br>
AIPS 2: Alt ref. value 0.00000E+00 wrt pixel 0.00<br>
AIPS 2: Maximum version number of extension files of type HI is
1<br>
AIPS 2: Maximum version number of extension files of type FQ is
1<br>
>getn 36<br>
AIPS 2: Got(1) disk= 3 user=2256 type=MA
A2256RELCL.FARSph.1<br>
>imh<br>
AIPS 2: Image=A2256 (MA) Filename=A2256RELCL
.FARSph. 1<br>
AIPS 2: Telescope=EVLA Receiver=EVLA<br>
AIPS 2: Observer=Frazer O User #= 2256<br>
AIPS 2: Observ. date=29-JUL-2010 Map date=01-JAN-2015<br>
AIPS 2: Pixel type: FLOATING Magic value blanking<br>
AIPS 2: Minimum= 1.00000002E+30 Maximum=-1.00000002E+30
DEGREES<br>
AIPS 2:
----------------------------------------------------------------<br>
AIPS 2: Type Pixels Coord value at Pixel Coord incr
Rotat<br>
AIPS 2: FARADROT 301 -4.5000000E+02 1.00
3.0000000E+00 0.00<br>
AIPS 2: RA---SIN 2081 17 03 43.500 660.00
-0.400000 0.00<br>
AIPS 2: DEC--SIN 1941 78 39 60.000 341.00
0.400000 0.00<br>
AIPS 2: FREQ 1 1.0160000E+09 1.00
1.0000000E+07 0.00<br>
AIPS 2:
----------------------------------------------------------------<br>
AIPS 2: Coordinate equinox 2000.00<br>
AIPS 2: Map type=NORMAL Number of iterations= 24333<br>
AIPS 2: Conv size= 6.00 X 6.00 Position angle= 0.00<br>
AIPS 2: Rest freq 0.000 Vel type: OPTICAL wrt YOU<br>
AIPS 2: Alt ref. value 0.00000E+00 wrt pixel 0.00<br>
AIPS 2: Maximum version number of extension files of type HI is
1<br>
AIPS 2: Maximum version number of extension files of type FQ is
1<br>
<br>
INPUTS TO FARS, RUN WITHIN DOFARS<br>
>tget fars<br>
>inp<br>
AIPS 2: FARS: Task to create Faraday rotation measure image<br>
AIPS 2: Adverbs Values Comments<br>
AIPS 2:
----------------------------------------------------------------<br>
AIPS 2: INNAME 'A2256ALL' First image name<br>
AIPS 2: INCLASS 'TRANS1' First image class<br>
AIPS 2: INSEQ 0 First image seq. #<br>
AIPS 2: INDISK 2 First image disk drive
#<br>
AIPS 2: IN2NAME 'A2256ALL' Second image name<br>
AIPS 2: IN2CLASS 'TRANS2' Second image class<br>
AIPS 2: IN2SEQ 0 Second image seq. #<br>
AIPS 2: IN2DISK 2 Second image disk drive
#<br>
AIPS 2: OUTNAME 'A2256RELCL' First output image name<br>
AIPS 2: OUTSEQ 0 First output image seq.
#<br>
AIPS 2: OUTDISK 3 First output image
disk #<br>
AIPS 2: INFILE 'FITS:LR_A2256/DOFARS6. Input file of weights<br>
AIPS 2: WEIGHTS'<br>
AIPS 2: DOALIGN 1 Should images be
coincident?<br>
AIPS 2: (See HELP.)<br>
AIPS 2: BLC 1 *rest 0 Bottom left corner<br>
AIPS 2: TRC 1 *rest 0 Top right corner<br>
AIPS 2: See help for the 1st
axis<br>
AIPS 2: (must be
frequencies!!!)<br>
AIPS 2: APARM 150 3 Parameters for
algorithm:<br>
AIPS 2: 0 0 1 number of pixels at<br>
AIPS 2: 2 0 half of the Fourier<br>
AIPS 2: 1 0 transform output<br>
AIPS 2: 0 2 The whole number is<br>
AIPS 2: 2*APARM(1)+1<br>
AIPS 2: 2 cell size in 1/m^2<br>
AIPS 2: 0 =><br>
AIPS 2:
PI/(4*(Lmax^2-Lmin^2))<br>
AIPS 2: Lmax,Lmin-max,min
lambda<br>
AIPS 2: at the data<br>
AIPS 2: 3 0 => regular
output<br>
AIPS 2: 1 => output is
RMTF<br>
AIPS 2: 4 0=> CLEANed
Fourier<br>
AIPS 2: transform<br>
AIPS 2: 1=> unCLEANed
Fourier<br>
AIPS 2: transform<br>
AIPS 2: 5
0=>original(shifted back)<br>
AIPS 2: RE/IM are sent
out<br>
AIPS 2: ** press RETURN for more, enter Q or next line to quit
print **<br>
#<br>
AIPS 2: 1=>the shifted
RE/IM are<br>
AIPS 2: sent out<br>
AIPS 2: 2=>amplitude and
phase of<br>
AIPS 2: the data are
sent out<br>
AIPS 2: 6 is not used<br>
AIPS 2: 7 0=> convolve the
clean<br>
AIPS 2: components<br>
AIPS 2: 1=> no convolve<br>
AIPS 2: 8 0=> use the
Gaussian as<br>
AIPS 2: the convolve
function<br>
AIPS 2: 1=> use the Re
of RMTF as<br>
AIPS 2: the convolve
function<br>
AIPS 2: 9 full width of
Gaussian<br>
AIPS 2: convolve function,
at 0.5<br>
AIPS 2: level, in 1/m^2<br>
AIPS 2: 0 => fit to real
RMTF<br>
AIPS 2: 10 what send to
output?<br>
AIPS 2: 0 => sum of
CLEAN and<br>
AIPS 2: residual<br>
AIPS 2: 1 => CLEAN
result<br>
AIPS 2: 2 => residual<br>
AIPS 2: GAIN 0.999 Gain in the CLEAN<br>
AIPS 2: NITER 1 Maximum number of clean<br>
AIPS 2: components<br>
AIPS 2: FLUX 0 Minimum flux of clean<br>
AIPS 2: component (Jy)<br>
AIPS 2: DOHIST -2 -2 => copy 1st HI
only<br>
AIPS 2: -3 => copy no HI
files<br>
<br>
PRTMSG<br>
AIPS 2 11:48:48 tget fars<br>
AIPS 2 11:48:49 inp<br>
AIPS 2 11:49:07 dofars<br>
TRANS2 11:49:14 Task TRANS (release of 31DEC13) begins<br>
TRANS2 11:49:26 Create A2256ALL .TRANS1. 1 (MA) on disk 2
cno 20<br>
TRANS2 11:49:27 Copied FQ file from vol/cno/vers 2 1 1 to
2 20 1<br>
TRANS2 11:50:05 Begin row swap with subimaging<br>
TRANS2 11:55:35 Begin in-core x-y transpose<br>
TRANS2 11:58:28 Appears to have ended successfully<br>
TRANS2 11:58:28 mazel 31DEC13 TST: Cpu= 12.6 Real= 554
IO= 10865<br>
AIPS 2 11:58:29 Resumes<br>
TRANS2 11:58:29 Task TRANS (release of 31DEC13) begins<br>
TRANS2 11:58:41 Create A2256ALL .TRANS2. 1 (MA) on disk 2
cno 21<br>
TRANS2 11:58:44 Copied FQ file from vol/cno/vers 2 2 1 to
2 21 1<br>
TRANS2 11:59:20 Begin row swap with subimaging<br>
TRANS2 12:04:58 Begin in-core x-y transpose<br>
TRANS2 12:07:48 Appears to have ended successfully<br>
TRANS2 12:07:48 mazel 31DEC13 TST: Cpu= 12.1 Real= 559
IO= 10865<br>
AIPS 2 12:07:48 Resumes<br>
FARS 2 12:07:48 Task FARS (release of 31DEC13) begins<br>
FARS 2 12:07:48 ZTXOP2: using translated file name =<br>
FARS 2 12:07:48 ZTXOP2:
/local/astro/aips/FITS/LR_A2256/DOFARS6.WEIGHTS<br>
FARS 2 12:07:48 Line number= 1 Weight number= 1
Weight= 11.<br>
FARS 2 12:07:48 L2MEAN = 0.08707 m^2<br>
FARS 2 12:07:49 FARSIN: using restoring Gaussian of FWHM 2.000
1/m^2<br>
FARS 2 12:09:23 Create A2256RELCL .FARSam. 1 (MA) on disk 3
cno 35<br>
FARS 2 12:11:19 Create A2256RELCL .FARSph. 1 (MA) on disk 3
cno 36<br>
FARS 2 12:11:19 PROCESSING ROW 1<br>
FARS 2 12:11:25 PROCESSING ROW 2<br>
FARS 2 12:11:34 PROCESSING ROW 3<br>
<br>
... ### INTERVENING MESSAGES DELETED BY LR ###...<br>
<br>
FARS 2 12:19:32 PROCESSING ROW 1940<br>
FARS 2 12:19:32 PROCESSING ROW 1941<br>
FARS 2 12:19:32 OUTPUT Phase FILE COMPLETELY BLANKED -- saving
it anyway<br>
FARS 2 12:19:32 History file created and written for Phase
file<br>
FARS 2 12:19:32 Copied FQ file from vol/cno/vers 2 20 1 to
3 36 1<br>
FARS 2 12:19:32 OUTPUT Amplitude FILE COMPLETELY BLANKED -- saving
it anyway<br>
FARS 2 12:19:32 History file created and written for Amplitude
file<br>
FARS 2 12:19:32 Copied FQ file from vol/cno/vers 2 20 1 to
3 35 1<br>
FARS 2 12:19:32 Appears to have ended successfully<br>
FARS 2 12:19:32 mazel 31DEC13 TST: Cpu= 29.4 Real= 704
IO= 19734<br>
AIPS 2 12:19:32 Resumes<br>
AIPS 2 12:19:32 Destroyed 1 extension files of type FQ<br>
AIPS 2 12:19:32 Destroyed 1 extension files of type HI<br>
AIPS 2 12:19:32 Destroyed MA image file: catno= 20 disk= 2<br>
AIPS 2 12:19:32 Destroyed 1 extension files of type FQ<br>
AIPS 2 12:19:32 Destroyed 1 extension files of type HI<br>
AIPS 2 12:19:32 Destroyed MA image file: catno= 21 disk= 2<br>
TRANS2 12:19:33 Task TRANS (release of 31DEC13) begins<br>
TRANS2 12:19:33 A2256RELCL .TRANam. 0 MA 2256 NOT FOUND<br>
TRANS2 12:19:33 FAILS TO FIND AND OPEN INPUT IMAGE<br>
TRANS2 12:19:33 Purports to die of UNNATURAL causes<br>
TRANS2 12:19:33 mazel 31DEC13 TST: Cpu= 0.0
Real= 0<br>
AIPS 2 12:19:34 Resumes<br>
AIPS 2 12:19:34 RETURN CODE 8 RECEIVED: STOPPING<br>
<br>
<br>
<br>
</font>
<pre class="moz-signature" cols="72">--
Lawrence Rudnick
Distinguished Teaching Professor
Minnesota Institute for Astrophysics
School of Physics and Astronomy
University of Minnesota
<a class="moz-txt-link-freetext" href="http://umn.edu/%7Elarry">http://umn.edu/~larry</a> <a class="moz-txt-link-abbreviated" href="mailto:larry@umn.edu">larry@umn.edu</a></pre>
</body>
</html>