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