; TARS ;-------------------------------------------------------------------- ;! Simulation of Faraday rotation synthesis (mainly task FARS) ;# TASK ANALYSIS POLARIZATION ;-------------------------------------------------------------------- ;; Copyright (C) 2009-2012 ;; Associated Universities, Inc. Washington DC, USA. ;; ;; This program is free software; you can redistribute it and/or ;; modify it under the terms of the GNU General Public License as ;; published by the Free Software Foundation; either version 2 of ;; the License, or (at your option) any later version. ;; ;; This program is distributed in the hope that it will be useful, ;; but WITHOUT ANY WARRANTY; without even the implied warranty of ;; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the ;; GNU General Public License for more details. ;; ;; You should have received a copy of the GNU General Public ;; License along with this program; if not, write to the Free ;; Software Foundation, Inc., 675 Massachusetts Ave, Cambridge, ;; MA 02139, USA. ;; ;; Correspondence concerning AIPS should be addressed as follows: ;; Internet email: aipsmail@nrao.edu. ;; Postal address: AIPS Project Office ;; National Radio Astronomy Observatory ;; 520 Edgemont Road ;; Charlottesville, VA 22903-2475 USA ;----------------------------------------------------------------------- TARS LLLLLLLLLLLLUUUUUUUUUUUU CCCCCCCCCCCCCCCCCCCCCCCCCCCCC TARS: Faraday rotation synthesis on simulated input data INFILE Input file of the U,Q dependence on frequency. OUTFILE Output file of the depedence on RM. BLANK=>OUTFILE is not created APARM Parameters for algorithm: 1 number of pixels at half of the Fourier transform output The whole number is 2*APARM(1)+1 2 cell size in 1/m^2 0 => PI/(4*(Lmax^2-Lmin^2)) Lmax,Lmin-max,min lambda at the data 3 0 => regular output 1 => output is RMTF 4 0=> CLEANed Fourier transform 1=> unCLEANed Fourier transform 5 0=>original(shifted back) RE/IM are sent out 1=>the shifted RE/IM are sent out 2=>amplitude and phase of the data are sent out 6 Number of rows to use in INFILE 0=> all rows in INFILE 7 0=> convolve the clean components 1=> no convolve 8 0=> use the Gaussian as the convolve function 1=> use the Re of RMTF as the convolve function 9 full width of Gaussian convolve function, at 0.5 level, in 1/m^2 0 => fit to real RMTF 10 what send to output? 0 => sum of CLEAN and residual 1 => CLEAN result 2 => residual GAIN Gain in the CLEAN NITER Maximum number of clean components FLUX Minimum flux of clean component (Jy) OPCODE 'CMPL' - new method Clean else peak amplitude Clean ---------------------------------------------------------------- TARS Type: Task The task is for testing the result of Faraday rotation Synthesis (FARS) Adverbs: INFILE......User-supplied text file giving Q (second column) and U (third column) as a function of frequency (first column). A weight may be given in the fourth column. These cards determine the frequecies that will be used and their weights. Default values for Q and U are 0 while the default weight is 1. Model components (up to 20) may be added to the data contained in the cards described above. If the first symbol is "M", then the following 3 numbers give the parameters of a model component: RM, in 1/m^2; AMP, and Phase in degrees. The models are computed and added to the data given in the main data cards --- if you want no data added to the model, specify Q and U as 0.0 in the data cards. If the first symbol in a row is a semicolon, the row is skipped and acts as a comment. The format of columns in INFILE is free with any number of blank characters between columns. The following is example of an INFILE with 2 components in a model and 3 of the N rows of data ; Input file for TARS With the top rows for model. ; The model rows are started with symbol M ; and includes RM in 1/m^2, AMP and PHASE in Degrees ; RM, 1/m^2 AMP Phase, deg M -600 5 60 M -800 3 0 ; FREQ, Hz Q U weight 1266000000 -3.6542873 -2.3332453 0.9 1276000000 -5.0204344 -1.8241444 1.0 1286000000 -4.3795457 0.393933 0.7 OUTFILE.....The result of calculation of Faraday rotation synthesis is written to this file consist of 4 columns 1. the row number; 2. RM in 1/m^2 3. Amplitude (or Real, depending on APARM(5)) 4. Phase (or Imaginary, depending on APARM(5)) ' ' => results are directed to the message file and message display APARM.......Parameters: the same as in FARS APARM(1) number of pixels at half of the Fourier transform output. The whole number is 2*APARM(1)+1 with zero at the center. The value of APARM(1) should be chosen in accordance of expected range of the Faraday rotation measures and value of CELL (APARM(2)) APARM(2) cell size of the outputs, in 1/m^2. The cell is recomended to be less or around the default value (pi)/4/(lambda^2max - lambda^2min) APARM(3) 0 => regular outputs 1 => outputs are RMTF APARM(4) 0=> CLEANed Fourier transform, using inputs: NITER, FLUX, GAIN. The CLEAN uses the shifed (at lambda^2) data but the cleaned components correspond to the original lambda^2 RE is recorded to the first output IM is recorded to the second output 1=> the uncleaned Fourier transform is recorded RE is recorded to the first output IM is recorded to the second output APARM(5) 0=>original data (RE and IM) are used at the Fourier transform 1=> the shifted (to the center) data (RE, IM) are used at the Fourier transform. This option allows better discrimination of different features at the Fourier transform. 2=> amplitude and phase of the data are sent out. APARM(6) Number of rows to use in INFILE 0=> all rows in INFILE If the first symbol at the infile row .EQ. semicolumn, this row is skipped. THIS SEICOLUMN CONCEPT CAN BE USED FOR SELECTION THE ROWS! APARM(7) 0=> convolve the clean components 1 => no convolve; So just the set of the clean components is sent to the output files No convolution is forced if: uncleaned Fourier is sent out (APARM(4)=1) APARM(8) 0=> use the Gaussian as the convolve function 1=> use the Re of RMTF as the convolve function APARM(9) full width of the Gaussian convolve function, at the 0.5 level, in 1/m^2. 0 -> fit the real part of the RMTF. APARM(10) What send to the output? 0 => sum of CLEAN and residual 1 => CLEAN result 2 => residual GAIN........Gain in the CLEAN NITER.......Maximum number of clean components 0 => 1 FLUX........Minimum Clean component (Jy) The task can subtract the given number of complex components. On each iteration, the maximum (and its position) of the spectrum amplitude is determined. The complex function RMTF is multiplied by the complex value of the spectrum at the position of the found amplitude maximum and by GAIN. The found function(production) is put by its median on the position of found amplitude maximum, and is subtracted from the having evaluated spectrum. The process of the subtraction is terminated having achieved number of iterations NITER or the flux FLUX. OPCODE......'CMPL' invokes a (slower) true complex Clean rather than the one that looks for the peak in the amplitude. other - standard peak in the amplitude Clean. ----------------------------------------------------------------