title test ots, inward fractions for pure hydrogen, open geo, filling factor
c must iterate to get inward line fractions, that are asserted below
iterate
diffuse ots
filling 0.01
print diffuse continuum
print last iteration
print line inward
ionization -2
blackbody 40,000
stop efrac 0.05 %stop when 5 percent ionized 
hydrogen levels large
hden 5
constant temper 7500
init file "honly.ini"  
punch overview  "honlyotsopen.ovr" last
punch asserts "honlyotsopen.asr" last
* not iterating does not affect the solution
no level2
no induced processes (OK) 
normalize to "Ca B" 4861 = 1 % normalize to exact result 
c this is an open geometry, so do not expect H atom to exactly be case B
c due to Lyman line leakage.  we iterated, so have a good estimate of
c inward and outward fractions, which we test 
assert line "Q(H)" 4861 > 1.0,  %>1 due to lyman line leaking
assert line "Q(H)" 4861 < 1.1,  %>but not a lot less
* 99jun26, error had been just under 2%, increased above it when optical
* depths of 2s and 2p were brought together, so increased error
assert line "H  1" 4861 = 1, error 0.03 %no surprise if not unity, since open
assert line "Inwd" 4861 =0.5 error 0.03 %no surprise if not unity, since open
assert line "nFnu" 3644 = 19.1 error 0.02 %total Balmer continuum rel to h-beta
assert line "InwT" 3644 = 9.55 error 0.02 %inward Balmer continuum rel to h-beta
assert line "BaC " 3646 = 18.26 error 0.02%residual Balmer continuum rel to h-b
assert line "cout" 3646 = 9.07 error 0.02 %Balmer continuum rel to h-beta
assert line "cref" 3646 = 9.07 error 0.02 %Balmer continuum rel to h-beta
c honlyotsopen.in