title "New" Paris meeting HII region
init file="ism.ini"
c "standard" HII region model of the Pequignot Meudon Conferance
sphere
black body, T=40,000K radius = 12.113943
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
iterate   ; must iterate since fine structure lines are opticall thick
*print last
print faint 0.01
normalize to "Ca B" 4861
punch overview last "parishii.ovr"
punch results last "parishii.rlt"
punch asserts last "parishii.asr"
print line sum
TOTL 5876
C  2 2326
TOTL 1909
N  2  121
N  2 6584
N  2 6548
N  3  570
TOTL 3727
O  3 5007
O  3 4959
O  3  518
O  3  883
Ne 2  128
Ne 3  156
Ne 3 3869
Ne 3 3968
S  2 6720
S  3  187
S  3  334
S  3 9532
S  3 9069
S  4  105
end of lines
dielectronic kludge 0
assert line "TOTL" 4861 1
assert line "Q(H)" 4861 > 1.
assert line "Q(H)" 4861 < 1.05
assert line "totl" 5876 0.1175
assert line "C  2" 2326 0.15
assert line "TOTL" 1909 0.0693
assert line "N  2"  121 0.0278
assert line "N  2" 6584 0.578
assert line "N  3"  570 0.303
assert line "TOTL" 3727 2.06
assert line "O  3"  518 1.263
assert line "O  3"  883 1.156
assert line "O  3" 5007 1.753
assert line "Ne 2"  128 0.176
assert line "Ne 3"  156 0.294
assert line "Ne 3" 3869 0.0671
assert line "S  2" 6720 0.122
assert line "S  3"  187 0.617
assert line "S  3"  334 0.87
assert line "S  3" 9532 0.914
assert line "S  4"  105 0.548
assert line "Stoy"    0 12.28
assert ionization fraction helium 2 volume 0.932
assert nzone < 150
assert itrzn < 6
assert hydrogen 2 temperature volume 7651 0.01
c parishii.in