Input of the model atmosphere

The input of the model atmosphere is driven by the parameter IGRAY in the file FT05.DAT. If the value of the parameter IGRAY is less then zero then there is no input of the model atmosphere and a plane-parallel LTE grey model is calculated as the first model. The values of the parameter IGRAY larger or equal to 2 enable the input of models calculated by other programs. The value IGRAY=1 means that the input model is the one calculated previously by this program. The output of the model atmosphere file is called MODELO.DAT.

XLUM0, RADIU0, XMSTA0, TAUMIN, TAUMAX, LNLINP, ITER0, LSPHIN
ND
ID, DM(ID), DR(ID), TEMP(ID), ELEC(ID), DENS(ID), TAUR(ID)
NL, LBG0, LNLLEV

following lines only for LBG0=LT

ID, BFAC(ID,IL)
NL

ID, BFAC(ID,IL)

XLUM0

stellar kuminosity (in ergs) - does not affect calculation, since the value from FT05.DAT is used

RADIU0

stellar radius (in cm) - does not affect calculation, the value from FT05.DAT is used

XMSTA0

stellar mass (in g) - does not affect calculation, the value from FT05.DAT is used

TAUMIN

minimum optical depth - does not affect calculation

TAUMAX

maximum optical depth - does not affect calculation

LNLINP

T - input of the NLTE model atmosphere

F - input of the LTE model atmosphere, b-factors are set to 1

ITER0

number of iterations needed for the calculation of the model - does not affect calculation

LSPHIN

T - input of the spherically symmetric model atmosphere

F - input of the plane-parallel model atmosphere (all radii are 0); if spherically symmetric model is calculated from a plane-parallel one, the first estimate of radius is calculated by program during initial phase

number of depth points of the model

ND lines, for each depth point 1:

index of the depth point - does not affect calculation

DM(ID)

column mass depth

DR(ID)

radius (in cm); for a plane-parallel atmosphere 0

TEMP(ID)

temperature

ELEC(ID)

electron density

DENS(ID)

density - does not affect calculation, density is recalculated by program

TAUR(ID)

Rosseland optical depth - does not affect calculation, Rosseland optical depth is recalculated by program

the rest of the input only for LNLINP=.TRUE.

for LBG0=.FALSE. number of NLTE levels, for LBG0=.TRUE. number of background NLTE levels; must be consistent with the value determined from input of atomic data, otherwise the program complains and stops

LBG0

T - background levels considered (background levels do not change its population numbers during calculation)

F - no background levels

LNLLEV

only for LBG0=.TRUE.

T - input first of background NLTE levels; then of the NLTE levels to be solved

F - input of background NLTE levels only

only for LBG0=.TRUE.; NL/5+1 times ND lines

index of the depth point - does not affect calculation

BFAC(ID,IL)

departure coefficients (b-factors) for each background NLTE energy level and depth point

only for LBG0=.TRUE.

number of NLTE levels; must be consistent with the value determined from input of atomic data, otherwise the program complains and stops

index of the depth point - does not affect calculation

BFAC(ID,IL)

departure coefficients (b-factors) for each NLTE energy level considered

Back to the main page of the code ATA
Last update: 27.8.2001