Input for LINFIT

LINE.DAT

In the following, green quantities in italics have to be read in, blue quantities in boldface keep their values until they are changed.
For blends, line data for each component must be in consecutive lines.

Part 1

Start of file
IJMA Number of ions to analyze
NAMIJ (1...IJMA) Names of all ions to analyze, e.g., Fe I: 2600, Fe II: 2601
Hydrogen (100) has always to be included and listed first.
IFATM Listing of atmosphere: 0: no listing, 1: listing
IFDEP Same as IFATM for departure coefficients

Control flags (Part 2)

The following lines consist of a control flag and the corresponding data: the lines can be in any sequence that makes astrophysical sense.

Control
flag
DataDefaultMeaning
70 IFSPHA 0 0: plan parallel,
1: spherical approximation

60 IFSCAT 0 0: treat scattering as true absorption,
1: correct treatment

80 NLTE 0 0: LTE,
1: non-LTE (requires file DEPART.DAT)

50 Delta log C6 0 global scaling factor for van der Waals broadening
(may appear within blend)

51 Delta log C4 0 global scaling factor for Stark broadening
(may appear within blend)

52 Delta log gammarad 0 global scaling factor for radiative damping
(may appear within blend)

53 NAMK ABU ATM.DAT
value
new abundance ABU for element/ion NAMK;
effectively changes log gf of all lines of element NAMK

54 FIT_ITEM... Items to vary (fit) for spectrum fits (IA > 8); is in general the control flag of item to fit---several can be specified on one line.
Implemented so far: 13 ximacro
14 v sin i
55 ximicro
>=98 abundance;
same as NAMK
Effect of several 54 lines is cumulative.

55 XI_MICRO ATM.DAT
value
new ximicro (always depth independent);
changes vDop only---no pressure change

40 cos theta 1 cos theta; if -1 calculate flux Flambda

30 KMA 1 number of blended lines,
for synthesis: number of all lines to synthesize

if -1: linfor counts blended lines until it finds another 30 -1 line.

31 ---.--- do not iterate log gf for this line in the blend; must precede line data immediately

21 WLIM 0 don't print lines with Wlambda calc < WLIM (mÅ)

20 IA 6
1:compute Wlambda from log gf,
and output of each iteration
2: iteration of log gf until Wlambda calc = Wlambda obs
or synthetic spectrum matches observation (requires file OBS.OUT)
4: Coarse analysis, get coarse correction to Wlambda obs from curve of growth
8: make spectrum synthesis
numbers can be added, e.g., 6 to start from curve of growth or 10 to match synthetic and observed spectrum

10 LAMMIN LAMMAX DELLAM ITV [YPMIN YPMAX [XPMIN XPMAX [PID] ] ]
none
for synthesis: lambdamin, lambdamax, Delta lambda; all [Å]
ymin, ymax, xmin, xmax for graphics plot
if LAMMIN < 0 (might not work): wavenumbers [cm-1]
ITV = 0: only printout
1: also crude ASCII screen plot
2: graphics plot on screen
3: like 2 but on paper
For graphics plot:
ymin = ymax: scale ordinate that observation fits in plot (default)
xmin = xmax: scale abscissa that both observation and synthesis fit in plot (default)
PID = 1: mark significant lines in plot (default)
!= 1: no line id's in plot
ymin etc. may be omitted as indicated by the square brackets

Environment variables LINLTE_SCREEN and LINLTE_PAPER allow control over plot devices and LINLTE_PRINT over the host print command for paper plots (as a Fortran FORMAT statement). LINLTE_SCALE sets a scale factor for paper plots.

11 IFABS 0 for synthesis:
0: relative flux/intensity in PLOT.OUT
(continuum normalized to 10000)
1: absolute flux/intensity F/Ilambda [erg/s/cm^2/cm]
= 1/pi Flambda = 4 Hlambda
2: absolute flux/intensity Fnu/Inu
3: relative flux/intensity (continuum normalized to 1)
Numbers can be combined to get more than one of the above, e.g., 31 to get Flambda and relative flux (unity). Restriction: leading zero as in 02 does not work. Use 20 instead.
12 ISYN 0 depth dependent output for wavelength point ISYN during synthesis if requested by line data, i.e., multiplet number <0

13 XIMAC FALT.DAT macro turbulence [km s-1] for synthesis

14 VSINI BETACL 0 0 projected rotational velocity v sin i and limb darkening parameter betacl for synthesis

15 IRV [RV] | [XRMIN XRMAX] radial velocity correction of observation
0 0
IRV = 0correct by RV [km/s] specified here
1: determine RV automatically from cross-correlation over interval XRMIN,XRMAX.

If XRMIN=XRMAX or omitted use whole range of overlap between observation and synthesis.

RV or XRMIN,XRMAX may be omitted as indicated by the square brackets

90 IFHOP 0
0: no Hydrogen lines,
1: include Hydrogen lines (Paschen to Lyman)--- different theory than explicit specification in line section (don't use both!)

This works properly only for syntheses.

9 IFLIMITS IFAUTO 0 0
IFLIMITS0: fit whole range
1: select range with cursor;
sets plot device [ITV] to 2 [screen] unless ITV=3 [paper] was set (see flag 10)
2: read range from file LIMITS.DAT
3: 1 and 2 combined: cursor with initial range read from file LIMITS.DAT
IFAUTO0: interactive termination of fit
>0: number of iterations

0 MLAM COSTHE Continuum calculation at wavelength MLAM
MLAM < 0: depth dependent output
cos theta < 0: Flambda

-1 Comment none ignore rest of line

-4 file name include line data from different file -- can be nested. (File name must be in single quotes.) Switches back to previous file on end of included file.

May be used within blends.

-5 ---.--- end of calculations; exit program

>=98 line data---see next page

Line data

line data, i.e., IZ > 98, all data must be in a single line

NAMK element/ion: z*100 + ion.stage, e.g.,
Fe I:2600
Fe II: 2601
MgH: 1200100

MULT multiplet number, if < 0: depth dependent output

LAMK lambda [Å], if < 0: [cm-1] (might not work)

CHIK chi [eV], excitation potential of lower level

GFLG log gf

DRRCA |Delta r2/a02| [Bohr radii] difference of mean square radii r2 of upper and lower level for van der Waals broadening; computes -log C6 = -log Delta r2/a02 + 32.3867

if <0: use Unsöld formula with
LU DIU LO DIO
llow Delta Elow lup Delta Eup
orbital quantum number of lower state excitation energy of parent term of lower state same for upper state

C4LG -log C4 for Stark broadening;
if -1: no broadening
if 0: approximation according to Griem, 1968, Phys.Rev. 165, 258 (ions) and Cowley, 1971, Obs 91, 139 (neutrals)
if -2: He I lines:
approximation according to Griem, 1962, Phys.Rev. 125, 177, for phase shift, and Bennett & Griem, Univ. Maryland Report 71, for Stark parameter. Available only for 2829, 2945, 3888, 3964, 4121, 4388, 4437, 4713, 4922, 5015, 5048, 5875, 6678, 7065Å
DRRCA and RAD ignored.

Other lines: not permitted

if -3: He I lines:
use broadening tables by Barnard, Cooper, Smith (1974, JQSRT 14, 1025) for 4471Å, and Shamey (1969, PhD thesis) for 4026, 4387, 4922Å. (requires files BCSS_4471 [4471Å] and BCSS_HE1)
DRRCA and RAD have no meaning.

Other lines: not permitted

if -4: He I lines:
use broadening tables by Beauchamp (priv. comm., following Beauchamp, Wesemael, and Bergeron 1997, ApJS, 108, 559) (requires file BWB)
DRRCA and RAD have no meaning.

Other lines: not permitted

anything else:not permitted

RAD gammarad [108 s-1]; if -1 gammarad class, i.e., gammarad = 2.22·1015/lambda2

WNOTE Wlambda obs [mÅ]; ignored for synthesis and all but last component of a blend

ILOW number or label of lower level in non-LTE model atom; can be omitted if LTE

IUP number or label of upper level in non-LTE model atom; can be omitted if LTE

Hydrogen and Helium II lines:

slightly changed meaning of line data, two theories available:

Simpler Stark broadening (Traving, Griem...):
(Hydrogen only)

GFLG log K + 17, see Traving 1962, ApJ 135, 439

DRRCA nlow principal quantum number

C4LG nup

RAD Ce/CH, resonance broadening, see Cayrel, Traving 1960, Z.Astrophys 50, 239 if 0, no resonance broadening

DLOGC6, DLOGC4, DLOGGR functionless

Unified Theory (Vidal, Cooper, Smith for Hydrogen; Schöning & Butler for He II)---requires file VCS.DAT (Hydrogen) or VCSSB.DAT (Helium):
LAMK lambda [Å]; if 0 and both nlow!=0 (LU) and nup!=0 (DIU): program computes wavelength;      if nlow=0 or nup=0: taken as approximate wavelength to fix nlow and/or nup and to compute exact wavelength---see also LU and DIU below.

CHIK chi [eV]; if 0: program computes excitation potential

GFLG log gf; if 0: program computes f-value

DRRCA -1 to select VCS theory

LU nlow; if 0 and LAMK != 0: program computes lower principal quantum number from approximate wavelength and sets LAMK to exact wavelength

DIU nup; if 0 and LAMK != 0: program computes upper principal quantum number from approximate wavelength and sets LAMK to exact wavelength

LO Hydrogen only
-1: do fudge convolution (i.e., add) of Stark profile with Lorentz profile for resonance broadening
0: do correct convolution, very time consuming, most likely worse than fudge

DIO no meaning

C4LG Hydrogen only
-log C3 [108 sec-1]; if 0: program computes resonance broadening,
C3 = e2/4 pi me ( f1i / nu1i + f1j / nu1j )
if -1: no broadening

RAD Hydrogen only
natural damping [108 sec-1]; if -1: program computes natural damping,
gammarad = \sum_{k=1}^{i-1} Aik + \sum_{k=1}^{j-1} Ajk

DLOGC6 functionless
DLOGC4 and DLOGGR usual meaning but for log C3 and gammarad