[daip] FARS and DOFARS script
Lawrence Rudnick
larry at umn.edu
Wed Jan 7 21:39:45 EST 2015
Trying to run FARS after install 31DEC15 and just running MNJ. Two problems
1.
>tget fars
AIPS 1: CHARACTER 1 BAD IN
AIPS 1: BAD CHARACTERS IN ADVERB INFILE
AIPS 1: CHARACTER?
>inp
AIPS 1: FARS: Task to create Faraday rotation measure image
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME '2SIG+NOISB' First image name
AIPS 1: INCLASS 'TRANS1' First image class
AIPS 1: INSEQ 0 First image seq. #
AIPS 1: INDISK 1 First image disk drive #
AIPS 1: IN2NAME '2SIG+NOISB' Second image name
AIPS 1: IN2CLASS 'TRANS2' Second image class
AIPS 1: IN2SEQ 0 Second image seq. #
AIPS 1: IN2DISK 1 Second image disk drive #
AIPS 1: OUTNAME '2SIG+NOISB' First output image name
AIPS 1: OUTSEQ 9.030E+11 First output image seq. #
AIPS 1: ****** VALUE =************ EXCEEDS LIMIT 9999.0000
AIPS 1: OUTDISK 1.358E-19 First output image disk #
AIPS 1: CHARACTER 1 BAD IN
AIPS 1: BAD CHARACTERS IN ADVERB INFILE
AIPS 1: CHARACTER?
2. After doing restore 0 and setting all the inputs in FARS, then run
DOFARS. Following is
a) input variables, then
b) what happens when invoking DOFARS
c) the DOFARS script itself
a) ##############
AIPS 1: FARS: Task to create Faraday rotation measure image
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME 'A2256_X4' First image name
AIPS 1: INCLASS 'QFLATC' First image class
AIPS 1: INSEQ 1 First image seq. #
AIPS 1: INDISK 1 First image disk drive #
AIPS 1: IN2NAME 'A2256_X4' Second image name
AIPS 1: IN2CLASS 'UFLATC' Second image class
AIPS 1: IN2SEQ 1 Second image seq. #
AIPS 1: IN2DISK 1 Second image disk drive #
AIPS 1: OUTNAME 'A2256_X4' First output image name
AIPS 1: OUTSEQ 0 First output image seq. #
AIPS 1: OUTDISK 1 First output image disk #
AIPS 1: INFILE *all ' ' Input file of weights
AIPS 1: DOALIGN 1 Should images be coincident?
AIPS 1: (See HELP.)
AIPS 1: BLC *all 0 Bottom left corner
AIPS 1: TRC *all 0 Top right corner
AIPS 1: See help for the 1st axis
AIPS 1: (must be frequencies!!!)
AIPS 1: APARM 400 10 Parameters for algorithm:
AIPS 1: 0 1 1 number of pixels at
AIPS 1: 2 0 half of the Fourier
AIPS 1: 1 *rest 0 transform output
AIPS 1: The whole number is
AIPS 1: 2*APARM(1)+1
AIPS 1: 2 cell size in 1/m^2
AIPS 1: 0 =>
AIPS 1: PI/(4*(Lmax^2-Lmin^2))
AIPS 1: Lmax,Lmin-max,min lambda
AIPS 1: at the data
AIPS 1: 3 0 => regular output
AIPS 1: 1 => output is RMTF
AIPS 1: 4 0=> CLEANed Fourier
AIPS 1: transform
AIPS 1: 1=> unCLEANed Fourier
AIPS 1: transform
AIPS 1: 5 0=>original(shifted back)
AIPS 1: RE/IM are sent out
AIPS 1: 1=>the shifted RE/IM are
AIPS 1: sent out
AIPS 1: 2=>amplitude and phase of
AIPS 1: the data are sent out
AIPS 1: 6 is not used
AIPS 1: 7 0=> convolve the clean
AIPS 1: components
AIPS 1: 1=> no convolve
AIPS 1: 8 0=> use the Gaussian as
AIPS 1: the convolve function
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
#
AIPS 1: 1=> use the Re of RMTF as
AIPS 1: the convolve function
AIPS 1: 9 full width of Gaussian
AIPS 1: convolve function, at 0.5
AIPS 1: level, in 1/m^2
AIPS 1: 0 => fit to real RMTF
AIPS 1: 10 what send to output?
AIPS 1: 0 => sum of CLEAN and
AIPS 1: residual
AIPS 1: 1 => CLEAN result
AIPS 1: 2 => residual
AIPS 1: GAIN 0.1 Gain in the CLEAN
AIPS 1: NITER 0 Maximum number of clean
AIPS 1: components
AIPS 1: FLUX 0 Minimum flux of clean
AIPS 1: component (Jy)
AIPS 1: DOHIST -1 -2 => copy 1st HI only
AIPS 1: -3 => copy no HI files
AIPS 1: PRTLEV 0 Number of Cleans to print
b) ##############
>
>dofars
AIPS 1: STRING SIZE
>
[No messages in Message window]
c) ###############
>list dofars
AIPS 1: 1 PROC DOFARS
AIPS 1: 2 $ defaults
AIPS 1: 3 IF (OUTNAME=' ') THEN; OUTNAME = INNAME; END
AIPS 1: 4 IF (APARM(5)>=1.5) THEN
AIPS 1: 5 IF (OUTCLASS=' ') THEN; OUTCLASS = 'FARSam; end;
AIPS 1: 6 ELSE
AIPS 1: 7 IF (OUTCLASS=' ') THEN; OUTCLASS = 'FARSre; end;
AIPS 1: 8 END;
AIPS 1: 9 IF (OUTDISK=0) THEN; OUTDISK = INDISK; END
AIPS 1: 10 FBLC = BLC(3); FTRC=TRC(3); BLC(3)=0; TRC(3)=0
AIPS 1: 11 DOWINAIT=DOWAIT; DOWAIT=1
AIPS 1: 12 IF (INFILE<>' ') THEN; FINFIL=' ';
AIPS 1: 13 ELSE; FINFIL='FITS:DOFARS.WEIGHTS'; END
AIPS 1: 14 VNUM = 33 - GETPOPSN; VPUT dofars
AIPS 1: 15 $ rspec weights
AIPS 1: 16 IF (FINFIL<>' ') THEN;
AIPS 1: 17 FILEZAP FINFIL
AIPS 1: 18 DEFAULT rspec; VGET dofars; TASK 'RSPEC'
AIPS 1: 19 BLC 0; TRC 0; OPCODE 'WEIT'; DOSLICE=-1; DOPRINT=-3;
AIPS 1: 20 OUTPRINT=FINFIL; GO;
AIPS 1: 21 END
AIPS 1: 22 $ trans 1
AIPS 1: 23 DEFAULT trans; VGET dofars; TASK 'TRANS'
AIPS 1: 24 OUTNAME=INNAME; OUTCLASS='TRANS1'; OUTSEQ=0;
AIPS 1: 25 OUTDISK=INDISK; TRANSCOD='312'; GO
AIPS 1: 26 $ trans 2
AIPS 1: 27 DEFAULT trans; VGET dofars; TASK 'TRANS'
AIPS 1: 28 INNAME=IN2NAME; INCLASS=IN2CLASS; INSEQ=IN2SEQ;
AIPS 1: 29 INDISK=IN2DISK; OUTNAME=INNAME; OUTCLASS='TRANS2';
AIPS 1: 30 OUTSEQ=0; OUTDISK=INDISK; TRANSCOD='312'; GO
AIPS 1: 31 $ fars
AIPS 1: 32 DEFAULT fars; VGET dofars; TASK 'FARS'
AIPS 1: 33 INCLASS='TRANS1'; IN2CLASS='TRANS2'; INSEQ=0; IN2SEQ=0;
AIPS 1: 34 IF (INFILE=' ') THEN; INFILE=FINFIL; END
AIPS 1: 35 IF (INFILE='1') THEN; INFILE=' '; END
AIPS 1: 36 BLC=FBLC,0; TRC=FTRC,0; GO
AIPS 1: 37 $ clean up
AIPS 1: 38 ZAP; INNAME=IN2NAME; INCLASS=IN2CLASS; INSEQ=IN2SEQ;
AIPS 1: 39 INDISK=IN2DISK; ZAP
AIPS 1: 40 $ more trans
AIPS 1: 41 TASK = 'TRANS'; INNAME=OUTNAME; INCLASS='FARSre
AIPS 1: 42 INDISK=OUTDISK; INSEQ=OUTSEQ; BLC=0; TRC=0; TRANSC='231';
AIPS 1: 43 IF (APARM(5)>=1.5) THEN; OUTCLASS='TRANam
AIPS 1: 44 INCLASS='FARSam
AIPS 1: 45 ELSE; OUTCLASS='TRANre
AIPS 1: 46 END; GO
AIPS 1: 47 INNAME=OUTNAME; INDISK=OUTDISK; INCLASS='FARSim
AIPS 1: 48 INSEQ=OUTSEQ; OUTNAME=OUTNAME; OUTSEQ=OUTSEQ;
AIPS 1: 49 OUTDISK=OUTDISK;
AIPS 1: 50 IF (APARM(5)>=1.5) THEN; OUTCLASS='TRANph
AIPS 1: 51 INCLASS='FARSph
AIPS 1: 52 ELSE; OUTCLASS='TRANim
AIPS 1: 53 END; GO
AIPS 1: 54 $ clean return
AIPS 1: 55 DOWAIT=DOWINAIT; VGET dofars; BLC(3)=FBLC; TRC(3)=FTRC
AIPS 1: 56 RETURN
--
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/20150107/45797bc9/attachment.html>
More information about the Daip
mailing list