1 GAUSS_1D Multi component one-dimension gaussian fitting
        GAUSS_1D

    This task makes a  multi-component  non-linear  gaussian  fit  into col-
    umns of a table.

    It fits a function Y = f(X) where f  is  the   sum   of   up   to   five
    gaussian,  and  X and Y data are taken in specified columns of the input
    table.  The method is identical to that used by  the  GAUSS  command  in
    CLASS (simplex followed by conjugate gradient).  The format for input of
    initial parameters  is  also  similar.   All  messages  are  sent  to  a
    specified  output  file.   The  fitted  profiles can be kept in separate
    columns of the input table.
2 IN$
        TASK\FILE "Input Table" IN$

    The name of the input table.
2 LIST$
        TASK\FILE "Output File" LIST$

    The name of the output formatted file.  All messages are written to this
    file, including the parameters of the fitted gaussian(s).
2 COLUMN_IN$
        TASK\INTEGER "Column Numbers" COLUMN_IN$[2]

    The number of the table column used as X and Y data.
2 NLINE$
        TASK\INTEGER "Number of gaussian components" NLINE$

    The number of gaussian components used in the fit.  NLINE$ should not be
    larger  than  5.   If  NLINE$  is  0,  the program will attempt to guess
    initial values for a single component from the moments of the  spectrum.
    If NLINE$ > 0, it will need some initial values.
2 LINE_1$
        TASK\REAL "Parameters for first component" LINE_1$[6]

    The initial parameters for the first gaussian component in  the  profile
    fit. These values are entered in free format as follows :
     Code, Value,     Code, Value,     Code, Value,     Code, Value
                       -Intensity-     -Position-        -Width-
    The code is an integer interpreted as follows:
         0       adjustable parameter
         1       fixed parameter
         2       adjustable parameter (head of group)
         3       parameter fixed with respect to parameter coded 2 or 4
         4       fixed parameter (head of group)
    Codes 2 3 and 4 are used to fit depending lines (e.g.   HCN,  for  which
    the  displacements  are  4.842 and -7.064 km/s, or -1.431 and 2.088 MHz,
    and line ratios 1:0.6:0.2).
2 LINE_2$
        TASK\REAL "Parameters for second component" LINE_2$[6]

    The initial parameters for the second gaussian component in the  profile
    fit.

2 LINE_3$
        TASK\REAL "Parameters for second component" LINE_3$[6]

    The initial parameters for the third gaussian component in  the  profile
    fit.
2 LINE_4$
        TASK\REAL "Parameters for second component" LINE_4$[6]

    The initial parameters for the fourth gaussian component in the  profile
    fit.
2 LINE_5$
        TASK\REAL "Parameters for second component" LINE_5$[6]

    The initial parameters for the fifth gaussian component in  the  profile
    fit.
2 COLUMN_OUT$
        TASK\INTEGER    "Column    numbers   for   fitted   profiles"   COL-
    UMN_OUT$['NLINE$']

    These are the numbers  of  columns  where values  of  the  fitted   com-
    ponents  will be kept.  Specify one empty (or useless) column of the in-
    put table for each component that you  want  to store, and  0  for  each
    component that you don't want to keep.
1 ENDOFHELP





















































