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 ---