title case A
c Seaton, M.J, 1959 MN 119, 90, 
c 4pi j(beta) 5.56E-26
c for this model total H-beta=4.745
c b(2)=3.73E-3 (3)=3.69E-2 (4)=0.091 (5)=0.145 (6)=0.193
c n.b. very different results if not l-mixed
c
// this should make the predicted totl Hbeta near unity for case a
normalize to "CaBo" 4861 = 2.23
no induced processes (OK) ;N.B. case A is a fiction; no incuded transitions
hydrogen collisions off ; must turn off all collisions 
stop lyman optical depth -6
constant temperature = 10,000K
stop zone 1
black body, T=50,000K
ionization parameter -2
hden = 13;
set dr 0
init file="honly.ini"
iterate
hydrogen levels large
print h-like levels
print departure coefficients
print last iteration
punch asserts "casea.asr" last
assert line "H  1" 6563 1.91
assert line "H  1" 4861 1
assert line "H  1" 4340 0.589
assert line "H  1" 4102 0.378
assert hydrogen 1 ionization -3.225
c casea.in