title cool HII region model from Lexington Meeting
punch overview last 11 "coolhii.ovr"
punch results last 12 "coolhii.rlt"
punch asserts last "coolhii.asr"
normalize to "Ca B" 4861
init file="ism.ini"
c Hbeta 4.93E36, L(total)4.30xHbeta
iterate   ; must iterate since fine structure lines are optically thick
dielec kludge 0 ; turn off guesses of dielectronic recombination
sphere
black body, T=20,000K 
q(h) 49
hden = 2
radius = 18.477121
abund He-1 C-3.6576 N-4.39794 O-3.481146 ne-4.30103 mg-8
continue si-8 s-5.04576 ar-8 fe-8
plot continuum range .1
print faint .01
print last iteration
print line sum :this is the sum of lines in table 2 of the Lexington meeting
N  2 6584
N  2 6548
TOTL 3727
Ne 2  128
S  2 6720
S  3  187
S  3  334
S  3 9532
S  3 9069
end of lines
assert line "TOTL" 4861 1
assert line "Q(H)" 4861 > 1.
assert line "Q(H)" 4861 < 1.05
assert line "N  2" 6584 0.6765
assert line "TOTL" 3727 1.16
assert line "Ne 2"  128 0.266
assert line "S  2" 6720 0.455
assert line "S  3"  187 0.458
assert line "S  3"  334 0.655
assert line "S  3" 9532 0.451
assert line "Stoy"    0 4.54
assert ionization fraction helium 2 volume 0.019
assert nzone < 140
assert itrzn < 6
c coolhii.in