title parispn.in Meudon Planetary nebula
init file="ism.ini"
punch overview "parispn.ovr"
punch physical conditions "parispn.phy"
punch asserts"parispn.asr"
sphere
c standard" PN model of the Pequignot Meudon Conferance
dielec kludge 0
black body, T=150,000K radius = 10
hden = 3.4771213
radius = 17
normalize to "Ca b" 4861
abund he -1 C-3.523 N-4. O-3.222 ne-3.824 mg-4.523 
continue si-4.523 s-4.824 ar-10 fe-10
print faint .01
normalize to "Ca b" 4861
// 99jun15, conv lyman lines change following uv line, error increased as relt
assert line "Q(H)" 4861 > 1.
assert line "Q(H)" 4861 < 1.1
assert line "BaC " 3646 14.4
assert line "H  1" 6563 2.898 error 0.02
assert line "H  1" 4861 1. error 0.02
assert line "H  1"   18 0.32
assert line "H  1"   12 0.158
assert line "TOTL" 5876 0.11
assert line "TOTL" 10830 0.916
assert line "He 2" 1640 2.35
assert line "He 2" 4686 0.34
assert line "C  2" 2326 0.30
assert line "TOTL" 1909 1.90
assert line "TOTL" 1549 2.60
assert line "N  2" 6584 1.03
assert line "N  3" 1750 0.12 error 0.1
assert line "N  3" 570  0.126
assert line "TOTL" 1486 0.24 error 0.1
assert line "TOTL" 1240 0.18 error 0.1
assert line "O  1" 6300 0.121
assert line "TOTL" 3727 2.33
assert line "O  3" 5007 16.8
assert line "TOTL" 4363 0.172
assert line "O  3" 518  1.305
assert line "O  4" 258  3.49
assert line "O  4" 1402  0.25 error 0.1
assert line "TOTL" 1218  0.20 error 0.1
assert line "Ne 2" 128  0.0277
assert line "Ne 3" 156  1.88
assert line "Ne 3" 3869  2.15
assert line "Ne 4" 2424  0.80
// chng 99 nov 18, small change in code chanced this line
assert line "Ne 5" 3426  0.59 error 0.1
assert line "Ne 5" 243  1.05
assert line "TOTL" 2798  2.33
assert line "Mg 4" 44  0.122
assert line "Si 2" 348  0.161
assert line "Si 2" 2335  0.16
assert line "TOTL" 1888  0.50 error 0.1
assert line "TOTL" 1397  0.23
assert line "S  2" 6720  0.36
assert line "S  3" 187  0.48
assert line "S  3" 9532  1.36
assert line "S  4" 105  2.08
assert ionization fraction volume helium 2 0.64
assert nzone < 130
assert itrzn < 6
assert temperature hydrogen 2 volume 12,300 0.02
assert temperature helium 3 volume 14,500 0.02
c parispn.in