flare

 FLARE N1 N2 OUT CYCLES WMOD HEAD -- Models spectra as the sum of a constant
                  spectrum plus variable amount of a flare spectrum. 
                  Iteratively estimates constant and low state spectra
                  by minimising chi**2 with restriction that sum of flare 
                  multipliers = 0 and sum over flare spectrum = 1

 Parameters

       N1, N2  -- Spectra to be used.
       OUT     -- First output slot. Will store constant then final flare
                  spectrum.
       CYCLES  -- Number of cycles to go through (10 or so should do)
       WMOD    -- If the model is accurate and you have low photon
                  numbers you may want to modify the weights on the
                  points. This is done by estimating a background count
                  by taking lowest value in whole of spectrum (suited to
                  HST), and then estimting variances in counts according to
                  V(J) = B + R(J)*ABS(L(J)+M(I)*F(J)) where R(J) is the
                  count/flux ratio, L(J) is the low state spectrum, F(J)
                  the flare spectrum and M(I) the multipliers. ABS
                  prevents silly errors. This option only applies for
                  CYCLES > 3 to get some stability. NB THIS IS A PERMANENT
                  CHANGE TO THE ERRORS ON THE SELECTED SPECTRA.
       BACK    -- Estimate background rate in counts/pixel to provide 
                  softening of weights at low count rates. Only if you
                  answer yes to WMOD. Note that the FOS on HST has a fairly
                  constant rate per pixel.
       HEAD    -- Name of Real type header parameter to store multipliers
                  in spectra used to form flare.  to ignore. The
                  parameter will be stored as HEAD_D, and HEAD_E.

This command belongs to the class: time_series