; MTARS
;---------------------------------------------------------------
;! Compute model inputs to TARS
;# TASK ANALYSIS POLARIZATION
;-----------------------------------------------------------------------
;;  Copyright (C) 2020
;;  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
;-----------------------------------------------------------------------
MTARS     LLLLLLLLLLLLUUUUUUUUUUUU CCCCCCCCCCCCCCCCCCCCCCCCCCCCC
MTARS:  Generate model inputs for TARS
INTEXT
                                   Text file with frequencies
FLUX                               Noise level to add
INLIST
                                   List of sources up to 9999
OUTFILE
                                   Output text file
----------------------------------------------------------------
MTARS
Type: Task
Use:  MTARS will take models of source flux, rotation measure, and
      "RM thickness" and compute input files for TARS.
Adverbs:
  INTEXT.....Text file containing a sample input file for TARS or
             simply a list of frequencies (in Hz) one per line.
  FLUX.......Noise level to be added.
  OPTYPE.....
  INLIST.....Text file containing one line per source, giving all 5 of

             Q  U  Spix  RM  RMthick  Type

             blank separated free format.  The 6th value, type has a
             default of 1.  Q, and U are those that would be
             measured at 1 GHz wavelength.  The resulting number of
             components will be the number of lines in INLIST
             neglecting comments (which start with # or ; in column 1)
             and lines with format errors.  Units are Jy/beam,
             Jy/beam, unitless, radians/m^2, and radians/m^2.  The
             type are
             The thickness is applied with a function specified by:
             1. 'SLAB' - sin(x)/x         with x = beta * lambda^2
             2. 'GAUS' - exp(-f * x * x)  with f = ln(2) / (1.8954)^2
             3. 'EXP ' - exp(-f * abs(x)) with f = ln(2) / 1.8954
             (Note 1.8954 radians is the point where sin(x)/x = 0.5
             and beta is called RMthick below in radians/m/m)
  OUTFILE....Output text file suitable for input to tars
----------------------------------------------------------------