1 FLOW
        FLOW

    This is a dedicated routine to  produce  bipolar  outflow   maps,   with
    optimum signal to noise ratio and little bias.  It takes as input a cube
    (N by NX by NY) with the velocity along the first axis, and  produces  a
    pseudo-cube  cube  (3 by NX by NY) containing the maps of the line width
    and of the red and blue lobes integrated intensities.

    For each spectrum, the algorithm  determines  the  peak  channel,   then
    determines  on  each  side  of  this  channel  the  velocity  at a given
    threshold, given as  ratio  to  peak  value.   From  this  velocity,  it
    integrates   out   to  a  second  threshold.   The  area  found  is  the
    contribution to the blue (or red) lobe of the flow.
2 IN$
        TASK\FILE "Input file name" IN$

    The input cube file name
2 OUT$
        TASK\FILE "Output file name" OUT$

    The output file name
2 RATIO$
        TASK\REAL "Ratio to peak value" RATIO$

    The ratio of the upper threshold used to define the  line  wing  to  the
    peak value of each spectrum.
2 THRESHOLD$
        TASK\REAL "Absolute detecttion threshold" THRESHOLD$

    The absolute level used to define the full line width.
1 END