Program:      CONREM

Purpose:      Removes continuum from channel maps by fitting a
              polynomial to the continuum channels.

Category:     ANALYSIS, COMBINATION

File:         conrem.c

Author:       K.G. Begeman

Keywords:

   INSET=     Set and subsets from which to subtract the continuum.
              Maximum number of subsets is 2048.

   FITSET=    Set and subsets for which to fit the continuum.
              Maximum number of subsets is 2048.

   NPOLY=     Degree of polynomial to be fitted [1].
              Maximum degree is 5.

   OUTSET=    Set and subsets for the continuum subtracted data
              [no continuum will be subtracted].

   CONPOS=    Coordinate for which to calculate continuum [No continuum
              is calculated]. If the interpolated continuum is wanted,
              enter here the coordinate for which it should be interpolated.

   CONSET=    Set and subset for the continuum. The interpolated continuum
              is calculated for the coordinate entered at CONPOS=.

Description:  A polynomial of degree NPOLY is fitted to the subsets
              selected with FITSET for each grid position in the subset
              separately. This fitted polynomial can then be subtracted
              from the subsets selected with INSET and/or can be used
              to calculate an interpolated continuum in CONSET at user
              selected coordinate CONPOS. The program is especially
              useful for removing the continuum from a series of line
              maps.

Example:      <USER> CONREM
              CONREM  Version 1.0  (Jul 30 1991)
              <USER> CONREM INSET=NGC4214 3:58
              Set NGC4214 has 3 axes
              RA-NCP             from  -127 to   128
              DEC-NCP            from  -127 to   128
              FREQ-OHEL          from     1 to    63
              <USER> CONREM FITSET=NGC4214 3:16 42:58
              Set NGC4214 has 3 axes
              RA-NCP             from  -127 to   128
              DEC-NCP            from  -127 to   128
              FREQ-OHEL          from     1 to    63
              <USER> CONREM NPOLY=1
              <USER> CONREM OUTSET=NGC4214_SUB
              Set NGC4214_SUB has 3 axes
              RA-NCP             from  -127 to   128
              DEC-NCP            from  -127 to   128
              FREQ-OHEL          from     3 to    58
              <USER> CONREM CONPOS=300 km/s
              Continuum calculated at 300.000000 KM/S
              <USER> CONREM CONSET=NGC4214_CON
              Set NGC4214_CON has 3 axes
              RA-NCP             from  -127 to   128
              DEC-NCP            from  -127 to   128
              FREQ-OHEL          from    32 to    32
              <STATUS>  CONREM   +++ FINISHED +++

Updates:      May 13, 1990: KGB Document created.
              Dec 11, 1991: KGB Cancel replace by reject.