______________________________________________________________________________
SF
SADFIT is a module which can fit a function to a row of a SAD file.
The functions presently supported are, a power law, a polynomial, and a
series of modeled absorption and/or emission lines with continuum. Nonlinear
weighted least squares is used for the solutions.
______________________________________________________________________________
Enter an option: 21-APR-1992 11:45:28.47
OPTION=SF
Enter the filename: (<CR> for default)
FILE=SPECTRUM2
There is one map in the file
Reading from the file
1 rows of 60 pixels in R4 format
Enter the type of fit:
I = Power law plus constant
P = Polynomial
S = Line profiles and continuum for spectral data (default)
CHOICE=
S
Enter the weight scheme:
E = Equal weights
R = Weights contained in a row
S = Weights based on signal-to-noise ratio (default)
WEIGHT_SCHEME=
S
Enter the readout noise in electrons and the gain in electrons per ADU:
(defaults = 10.0 10.0)
NOISE_PARAMETERS=
10.0000
10.0000
After the plot is drawn, press C to turn on the cursor. Then position the
cursor at the left end of the continuum to be fit and press the <CR>. For
each line to be fit, position the cursor at the left side of the line
(imagined to be unblended) at half amplitude, the peak or bottom of the
line for an emission or absorption line, respectively, and the imagined
continuum level in the center of the line. Press the <CR> at each
location. Next, position the cursor at the right end of the
continuum to be fit and press the <CR>. Mistakes can be voided by pressing
the - key. Press Q to turn off the cursor and <CR> to proceede.
Press the <CR> when through reading these instructions
THROUGH=
T
Enter the plot devices:
G = Graphics terminal (VT100/Retrographics and Visual 550, default)
G1 = Graphics terminal (GraphOn)
G2 = Graphics terminal (Codonics)
G3 = Graphics terminal (Micro-term, Versaterm)
F = File
T = Tektronix 4662 flat bed plotter (1200 baud)
DEVICES=
G
4001.51 920.798 These are the coordinates of each position
4005.47 812.880 at which the <CR> was pressed.
4006.73 715.519
4006.73 913.173
4010.19 844.551
4011.19 765.959
4011.19 924.317
4012.45 780.622
4013.30 616.399
4013.30 924.317
4017.61 917.865
3 lines were chosen to be fit between the limits of
4000.45 and 4016.05 .
53 data points will be used in the nonlinear least squares fit
Enter the degree of the polynomial which will be fit to the continuum:
(default = 3)
DEGREE=
3
Please wait for the solution
Mean error for an observation of average weight =
9.22081 +- 1.03092
The polynomial solution for the continuum is
DEGREE COEFFICIENT MEAN ERROR
0 0.41083E+07 0.95548E+08
1 -3169.4 71508.
2 0.81437 17.840
3 -0.69684E-04 0.14836E-02
Line Number 1
Optical Depth = 0.243063 +- 0.775325E-02
Sigma = 0.811531 +- 0.329358E-01
Position = 4006.56 +- 0.253076E-01
Line Number 2
Optical Depth = 0.179165 +- 0.811902E-02
Sigma = 0.747693 +- 0.579495E-01
Position = 4011.16 +- 0.572024E-01
Line Number 3
Optical Depth = 0.379496 +- 0.877719E-02
Sigma = 0.779997 +- 0.303247E-01
Position = 4013.45 +- 0.289311E-01
Do you wish to plot the residuals? (default = Y)
PLOT_RESIDUALS=
T
Enter the plot devices:
G = Graphics terminal (VT100/Retrographics and Visual 550, default)
G1 = Graphics terminal (GraphOn)
G2 = Graphics terminal (Codonics)
G3 = Graphics terminal (Micro-term, Versaterm)
F = File
T = Tektronix 4662 flat bed plotter (1200 baud)
DEVICES=
G
Do you wish to plot the data and function? (default = Y)
PLOT=
T
Enter the plot devices:
G = Graphics terminal (VT100/Retrographics and Visual 550, default)
G1 = Graphics terminal (GraphOn)
G2 = Graphics terminal (Codonics)
G3 = Graphics terminal (Micro-term, Versaterm)
F = File
T = Tektronix 4662 flat bed plotter (1200 baud)
DEVICES=
G
Do you want to write the function to a SAD file? (default = Y)
WRITE=N
______________________________________________________________________________