176 lines
8.0 KiB
Groff
176 lines
8.0 KiB
Groff
Reading input from stdin
|
|
|
|
================================
|
|
M O D E L A T M O S P H E R E
|
|
================================
|
|
|
|
TEFF = 35000.0
|
|
LOG G = 4.00
|
|
LTE = T
|
|
LTGRAY = T
|
|
|
|
FREQUENCIES:
|
|
NFREAD = 50
|
|
|
|
ATOMS: 8 elements configured
|
|
|
|
IONS:
|
|
H 1 (Z= 1, ion=0) - 9 levels, file: ./data/h1.dat
|
|
H 2 (Z= 1, ion=1) - 1 levels, file: (none)
|
|
He 1 (Z= 2, ion=0) - 14 levels, file: ./data/he1.dat
|
|
He 2 (Z= 2, ion=1) - 14 levels, file: ./data/he2.dat
|
|
He 3 (Z= 2, ion=2) - 1 levels, file: (none)
|
|
|
|
--- Reading atomic data files ---
|
|
Ion 1: H 1 <- ./data/h1.dat
|
|
Levels: 9, Continua: 8, Lines: 32
|
|
Level 1: G=0, NQUANT=0, IFWOP=0
|
|
Level 2: G=0, NQUANT=0, IFWOP=0
|
|
Level 3: G=0, NQUANT=0, IFWOP=0
|
|
... (6 more levels)
|
|
Ion 2: H 2 <- ground state only (fully ionized)
|
|
Ion 3: He 1 <- ./data/he1.dat
|
|
Levels: 14, Continua: 13, Lines: 72
|
|
Level 1: G=1, NQUANT=1, IFWOP=0
|
|
Level 2: G=3, NQUANT=2, IFWOP=0
|
|
Level 3: G=1, NQUANT=2, IFWOP=0
|
|
... (11 more levels)
|
|
Ion 4: He 2 <- ./data/he2.dat
|
|
Levels: 14, Continua: 13, Lines: 121
|
|
Level 1: G=0, NQUANT=1, IFWOP=0
|
|
Level 2: G=0, NQUANT=2, IFWOP=0
|
|
Level 3: G=0, NQUANT=3, IFWOP=0
|
|
... (11 more levels)
|
|
Ion 5: He 3 <- ground state only (fully ionized)
|
|
|
|
Total: 39 levels, 34 continua, 225 lines
|
|
|
|
--- Populating atomic data ---
|
|
Ion 1: nfirst=1, ntrans=40, ntranc=8
|
|
Ion 2: nfirst=10, ntrans=0, ntranc=0
|
|
Ion 3: nfirst=11, ntrans=85, ntranc=13
|
|
Ion 4: nfirst=25, ntrans=134, ntranc=13
|
|
Ion 5: nfirst=39, ntrans=0, ntranc=0
|
|
|
|
Total transitions: 259, continuum: 34
|
|
Total levels in atomic data: 39
|
|
|
|
--- Verification ---
|
|
Level 1 (H1 n=1): enion=8.7141e-11, g=2
|
|
Level 2 (H1 n=2): enion=2.1785e-11, g=8
|
|
Level 10 (He1 n=1): enion=2.1785e-11, g=1
|
|
|
|
--- Generating initial LTE grey atmosphere ---
|
|
Initial opacity estimate: 0.6654 cm²/g (computed from Teff=35000K)
|
|
DEBUG ELDENS iter 1: ane=3.010435e8, q=1.999638e-1, dqn=0.000000e0
|
|
DEBUG ELDENS coeffs iter 1: ae=5.627574e-12, gg=0.000000e0, e_val=0.000000e0, b_val=0.000000e0
|
|
DEBUG ELDENS params iter 1: d_val=1.480848e11, g2=1.480848e11, a_val=1.480848e11
|
|
DEBUG ELDENS matrix iter 1: R=[1.1000e0,0.0000e0,1.0000e0; -1.9996e-1,-1.4808e11,1.8334e0; -1.0000e0,1.4808e11,-8.3336e-1]
|
|
DEBUG ELDENS rhs iter 1: S=[2.5079e7,0.0000e0,0.0000e0]
|
|
DEBUG ELDENS lineqs: p=[1.090392e7, -1.090392e7, -1.090392e7]
|
|
DEBUG ELDENS update iter 1: ah=2.617811e8, anh=-1.090392e7, ane=2.901396e8, delne=-1.090392e7
|
|
DEBUG ELDENS iter 2: ane=2.901396e8, q=1.999648e-1, dqn=0.000000e0
|
|
DEBUG ELDENS coeffs iter 2: ae=-3.758164e-2, gg=-0.000000e0, e_val=-0.000000e0, b_val=-0.000000e0
|
|
DEBUG ELDENS params iter 2: d_val=1.536501e11, g2=1.536501e11, a_val=1.536501e11
|
|
DEBUG ELDENS matrix iter 2: R=[1.1000e0,0.0000e0,1.0000e0; -1.9996e-1,-1.5365e11,-5.7744e9; -1.0000e0,1.5365e11,5.7744e9]
|
|
DEBUG ELDENS rhs iter 2: S=[2.3988e7,-1.6754e18,1.6754e18]
|
|
DEBUG ELDENS lineqs: p=[9.235542e17, 8.466162e17, 8.466162e17]
|
|
DEBUG ELDENS update iter 2: ah=9.235542e17, anh=8.466162e17, ane=8.466162e17, delne=8.466162e17
|
|
DEBUG ELDENS iter 3: ane=8.466162e17, q=9.831719e-2, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[-1.131223e18, -1.839166e17, 6.362627e17]
|
|
DEBUG ELDENS iter 4: ane=1.482879e18, q=9.587030e-2, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[-5.527896e18, -5.094804e18, -4.337802e18]
|
|
DEBUG ELDENS iter 5: ane=-2.854923e18, q=0.000000e0, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[1.593445e19, 1.244316e19, 8.363850e18]
|
|
DEBUG ELDENS iter 6: ane=5.508926e18, q=7.231990e-2, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[-2.496604e18, 7.664910e18, 1.342711e19]
|
|
DEBUG ELDENS iter 7: ane=1.893604e19, q=2.035586e-2, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[-2.421280e19, -6.409854e18, 9.049457e17]
|
|
DEBUG ELDENS iter 8: ane=1.984098e19, q=1.868318e-2, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[-1.611619e19, -2.287177e19, -1.562108e19]
|
|
DEBUG ELDENS iter 9: ane=4.219901e18, q=8.040489e-2, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[-4.611725e20, -5.291496e20, -4.964130e20]
|
|
DEBUG ELDENS iter 10: ane=-4.921931e20, q=0.000000e0, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[1.657005e21, -1.928458e21, 7.873333e20]
|
|
DEBUG ELDENS after loop: ane=2.951402e20, an=6.020870e8
|
|
DEBUG ELDENS return: id=1, ane=2.951402e20, anp=2.951402e20, ahtot=1.163206e21, anerel=4.901952e11
|
|
DEBUG after ELDENS: ane=2.951402e20, an=6.020870e8, dens=-4.938669e-4
|
|
DEBUG ELDENS iter 1: ane=2.827072e8, q=1.999656e-1, dqn=0.000000e0
|
|
DEBUG ELDENS coeffs iter 1: ae=5.284795e-12, gg=0.000000e0, e_val=0.000000e0, b_val=0.000000e0
|
|
DEBUG ELDENS params iter 1: d_val=1.576896e11, g2=1.576896e11, a_val=1.576896e11
|
|
DEBUG ELDENS matrix iter 1: R=[1.1000e0,0.0000e0,1.0000e0; -1.9997e-1,-1.5769e11,1.8334e0; -1.0000e0,1.5769e11,-8.3336e-1]
|
|
DEBUG ELDENS rhs iter 1: S=[2.3552e7,0.0000e0,0.0000e0]
|
|
DEBUG ELDENS lineqs: p=[1.023994e7, -1.023994e7, -1.023994e7]
|
|
DEBUG ELDENS update iter 1: ah=2.458360e8, anh=-1.023994e7, ane=2.724673e8, delne=-1.023994e7
|
|
DEBUG ELDENS iter 2: ane=2.724673e8, q=1.999666e-1, dqn=0.000000e0
|
|
DEBUG ELDENS coeffs iter 2: ae=-3.758226e-2, gg=-0.000000e0, e_val=-0.000000e0, b_val=-0.000000e0
|
|
DEBUG ELDENS params iter 2: d_val=1.636159e11, g2=1.636159e11, a_val=1.636159e11
|
|
DEBUG ELDENS matrix iter 2: R=[1.1000e0,0.0000e0,1.0000e0; -1.9997e-1,-1.6362e11,-6.1491e9; -1.0000e0,1.6362e11,6.1491e9]
|
|
DEBUG ELDENS rhs iter 2: S=[2.2528e7,-1.6754e18,1.6754e18]
|
|
DEBUG ELDENS lineqs: p=[9.235706e17, 8.466299e17, 8.466299e17]
|
|
DEBUG ELDENS update iter 2: ah=9.235706e17, anh=8.466299e17, ane=8.466299e17, delne=8.466299e17
|
|
DEBUG ELDENS iter 3: ane=8.466299e17, q=9.831714e-2, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[-1.131242e18, -1.839184e17, 6.362738e17]
|
|
DEBUG ELDENS iter 4: ane=1.482904e18, q=9.587019e-2, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[-5.527914e18, -5.094814e18, -4.337800e18]
|
|
DEBUG ELDENS iter 5: ane=-2.854896e18, q=0.000000e0, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[1.593477e19, 1.244349e19, 8.364203e18]
|
|
DEBUG ELDENS iter 6: ane=5.509307e18, q=7.231754e-2, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[-2.497957e18, 7.663996e18, 1.342641e19]
|
|
DEBUG ELDENS iter 7: ane=1.893572e19, q=2.035647e-2, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[-2.421205e19, -6.410456e18, 9.043028e17]
|
|
DEBUG ELDENS iter 8: ane=1.984002e19, q=1.868489e-2, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[-1.611511e19, -2.287123e19, -1.562117e19]
|
|
DEBUG ELDENS iter 9: ane=4.218853e18, q=8.041147e-2, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[-4.610347e20, -5.289827e20, -4.962716e20]
|
|
DEBUG ELDENS iter 10: ane=-4.920528e20, q=0.000000e0, dqn=0.000000e0
|
|
DEBUG ELDENS lineqs: p=[1.656640e21, -1.928150e21, 7.872242e20]
|
|
DEBUG ELDENS after loop: ane=2.951715e20, an=5.654144e8
|
|
DEBUG ELDENS return: id=1, ane=2.951715e20, anp=2.951715e20, ahtot=1.162979e21, anerel=5.220445e11
|
|
DEBUG before pressure update: ane=2.951715e20, an=5.654144e8
|
|
Depth 1: T=28392K, ne=2.95e20, nH=8.68e20, rho=-4.94e-4
|
|
Quick estimate κ_R=4.0000e-1, Full LTE κ_R=4.0000e-1, κ_P=4.0000e-1
|
|
Components: κ_es=0.0000e0, κ_bf=0.0000e0, κ_ff=0.0000e0, κ_H-=0.0000e0
|
|
Depth 70: T=137404K, ne=1.29e236, nH=1.99e236, rho=2.17e212
|
|
Quick estimate κ_R=1.1870e8, Full LTE κ_R=1.1965e8, κ_P=3.0634e12
|
|
Components: κ_es=3.9751e-1, κ_bf=3.4214e-29, κ_ff=1.7235e-243, κ_H-=0.0000e0
|
|
Generated 70 depth points
|
|
Temperature range: 28392 K (surface) to 137404 K (bottom)
|
|
Electron density range: 2.95e20 to 1.29e236 cm^-3
|
|
|
|
--- Starting TLUSTY initialization ---
|
|
Initialization completed successfully
|
|
NN = 0
|
|
Success = true
|
|
|
|
--- Setting up frequency grid ---
|
|
Frequency range: 1.00e14 - 1.00e16 Hz
|
|
Number of frequency points: 50
|
|
|
|
--- Starting main iteration loop ---
|
|
|
|
=== Iteration 1 ===
|
|
Computing opacities...
|
|
Solving radiative transfer...
|
|
Updating populations...
|
|
Max flux error: 2.37e2
|
|
|
|
=== Iteration 2 ===
|
|
Computing opacities...
|
|
Solving radiative transfer...
|
|
Updating populations...
|
|
Max flux error: 2.99e276
|
|
|
|
=== Iteration 3 ===
|
|
Computing opacities...
|
|
Solving radiative transfer...
|
|
Updating populations...
|
|
Max flux error: 1.51e277
|
|
|
|
Main loop completed after 3 iterations
|
|
|
|
--- Writing model to fort.7 ---
|
|
Model written to fort.7
|
|
|
|
--- TLUSTY START completed ---
|