first line
|
Teff |
logg |
RSU |
lambda ref |
|
F8.2 |
F6.2 |
F8.2 |
F9.2 |
|
Effective Temperature |
Gravity
log g1
|
Radius in solar radii
For calculation of pressure stratification:
- If RSU <= 0:
-
plane parallel stratification with
log g=const
(Parameter is not really used in this case)
- RSU >0:
- spherical approximation with
log g= log g1 +
2 log (R( tau1)/R(tau) )
where subscript 1 stands for the outermost layer.
|
reference wavelength [Å] for optical depth scale |
2nd line
|
I | number of depth points
(I3)
- I=0
- use 30 point standard optical depth scale
- I>0
- read monochromatic optical depth scale from next lines,
5 numbers per line (5E10.3)
|
(2+I)th line
|
M | control flag
(I3)
M = 0 |
scaled solar HM model (built in):
T(tau) = TSun(tau) ·Teff/5780
|
M = 1
|
theta (5040/T) values for each depth point follow, I many, 5 per line
(5F7.4)
|
M = -1
|
T values for each depth point follow, I many, 5 per line (5F7.0)
|
|
(2+I+M)th line
|
TUR | microturbulence
(F6.2)
TUR < 0 |
depth dependent microturbulence [km/s]; values follow for every depth point, I
many, 5 per line (5F7.2)
|
TUR > 0
|
micro=TUR [km/s] for all depth points
|
|
next line
|
MODE, IONST, IPRINT, KPRINT, LPR1, LPR2, LPRD, KPID, PLMIN, PLMAX, KPKT
|
|
|
7 to 9 control flags (free format)
- a) MODE
- =1: IONDIS computes particle concentrations ZETAJ only for the usually
important elements
- =2: IONDIS computes particle concentrations for all elements it knows
about
- (computation of pressure stratification (pgas and
DENSNC) is always done with important
elements only)
- b) IONST
- ignored (not used anymore)
- c) IPRINT
- =1: no output from IONDIS
- =2: IONDIS prints ionization equilibria and similar for depth points
specified below
- d) KPRINT
- =1: no output from OPALAM
- =2: OPALAM prints details for opacities for all
wavelengths
- =4: make a file suitable for plotting with program OPAPLOT
- e-g) LPR1,LPR2,LPRD
- first depth point, last depth point, and increment for printing of details
- h-k) KPID,PLMIN,PLMAX,KPKT (optional
if KPRINT!=4)
- KPID: # of depth point for plot
- PLMIN,PLMAX: wavelength interval
- KPKT: # of wavelength points for plot
|
(3+I+M)th line
|
NLAKS | # of
wavelengths for computation of continuous opacities
(free format)
NLAKS=0: |
use the following 12 standard wavelengths (good for optical and IR only):
3648, 4000, 5000, 6000, 7000,
8000, 10000, 12500, 15000, 16420,
20000, 25000 [Å]
|
NLAKS>0: |
list the desired NLAKS wavelengths on the following lines, 5 per line
(5F10.2)
|
|
last line
|
LCONTR | # of depth
points again for consistency check. Must be 30 if I=0. (I3)
|