root/trunk/tsuite/auto/blr_f92.in

title standard blr cloud in Ferland et al. 1992
c
c commands controlling continuum =========
table agn
ionization parameter -1
c
c commands for density & abundances =========
hden 11
abundances old solar 84 
c
c commands controlling geometry  =========
stop column density 25
c
c other commands for details     =========
atom h-like element hydrogen levels resolved 10
atom h-like element helium levels resolved 10
* stop zone 20
* element iron ionization -2 0 -2 -2 -2 - 2 
iterate to convergence
c
c commands controlling output    =========
print h-like departure hydrogen 
print h-like populations hydrogen 
print h-like departure helium
print h-like populations helium
plot continuum
// >>chng 02 nov 20, Lya now 1.049 fainter, make this smaller
// normalize to "H  1" 1216 = 111.12
normalize to "H  1" 1216 = 105.90
print line faint -2
punch overview last "blr_f92.ovr"
punch convergence reason "blr_f92.cvr"
punch element nitrogen last "blr_f92.nit"
punch element calcium last "blr_f92.ca"
punch asserts last "blr_f92.asr"
punch dr "blr_f92.dr"
punch transmitted continuum "blr_f92.trn" last
c
c commands giving the asserts    =========
// >>chng 01 jul 04, from 9.327 to 9.371, diffuse fields
// >>chng 02 oct 29, from 9.371 to 9.348, slow drift
// >>chng 02 nov 20, from 9.348 to 9.327, new he ots diff
// >>chng 05 feb 18, from 9.327 to 9.349 incld outwrd cont in line pump
// >>chng 07 jul 25, from 9.349 to 9.369, bugfix pvh FeII tau r1290
// >>chng 07 oct 22, from 9.369 to 9.305, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 9.305 to 9.367, iteration had been disabled, turn back on
// >>chng 08 mar 17, from 9.367 to 9.341, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert line luminosity "TOTL" 1216  9.341
//
// >>chng 01 mar 09 from 1.00 to 0.84 change in H redis function
// >>chng 01 mar 16 from 0.84 to 0.894 change in atomic data
// >>chng 01 jul 04, from 0.895 to 0.829, diffuse fields
// >>chng 02 may 09, from 0.829 to 1.007, e scat esc prob
// >>chng 03 apr 11, from 1.007 to 1.499, PopOpc bug
// >>chng 03 sep 09, from 1.499 to 1.421, diffuse fields */
// >>chng 04 jun 12, from 1.421 to 1.326 chng cont shield fcn
// >>chng 05 mar 12, from 1.326 to 1.121, RP merge He-like code
// >>chng 06 jul 18, from 1.121 to 1.798, RP major upgrade VS coll strengths
// >>chng 06 jul 18, from 1.798 to 1.417, repair common error in several collision routines
// >>chng 07 oct 22, from 1.417 to 2.002, resolve l-levels of h-like sequence
// >>chng 07 oct 25, from 2.002 to 2.400, bugfix stark broadening incorrect
// >>chng 07 dec 19, from 2.400 to 2.147, Stark changes in r1698
// >>chng 08 feb 13, from 2.147 to 1.351, iteration had been disabled, turn back on
// >>chng 08 mar 17, from 1.351 to 0.880, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 mar 22, from 0.880 to 0.938, cutoff between volatile/static opacities was broken
// >>chng 08 apr 20, incr error to more realistic 10%
assert line "H  1" 4861      0.938 error 0.1
//
// >>chng 02 nov 08, add next assert
// >>chng 02 dec 11, from 8.457 to 8.469, new hei
// >>chng 04 apr 28, from 8.469 to 8.435, He+ - H ct rate update
// >>chng 04 jul 03, from 8.435 to 8.472, NA update CO, H chem net
// >>chng 07 jan 20, from 8.472 to 8.512, e2 bug
// >>chng 07 oct 22, from 8.512 to 8.455, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 8.455 to 8.488, iteration had been disabled, turn back on
// >>chng 08 mar 22, from 08.488to 8.467 bugfix r1879
assert line luminosity "H- H" 0  8.467
//
// >>chng 02 nov 08, add next assert
// >>chng 02 dec 11, from 8.503 to 8.528, new hei
// >>chng 04 apr 28, from 8.528 to 8.490, He+ - H ct rate update
// >>chng 04 jul 03, from 8.490 to 8.552, NA update CO, H chem net
// >>chng 06 jul 18, from 8.585 to 8.580, RP major upgrade VS coll strengths
// >>chng 07 oct 22, from 8.580 to 8.489, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 8.489 to 8.534, iteration had been disabled, turn back on
// >>chng 08 mar 22, from 8.534 to 8.510, cutoff between volatile/static opacities was broken
// >>chng 08 jul 08, from 8.510 to 8.48, damping wing pump rates
assert line luminosity "H-FB" 0  8.48 error 0.1
//
// >>chng 02 nov 08, add next assert
// >>chng 06 nov 29, from 9.214 to 9.235, Badnell DR turned on
// >>chng 07 oct 22, from 9.235 to 9.206, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 9.206 to 9.233, iteration had been disabled, turn back on
assert line luminosity "TotM" 0  9.233
//
// >>chng 01 mar 09 from 2.086 to 4.695 change in H redis function
// >>chng 01 mar 16 from 4.695 to 5.01 change in atomic data
// >>chng 01 may 19 from 5.01 to 6.75 ct upgrade, more ionization in piz
// >>chng 01 may 21 bad opacity in HeI 7065, from 6.75 to 6.17
// >>chng 01 may 30, from 6.17 to 6.84, change in line redistribution fcn
// >>chng 01 jun 15, from 6.84 to 4.95, O ct
// >>chng 01 jul 04, from 4.95 to 4.585 diffuse fields
// >>chng 02 may 09, from 4.255 to 4.998 e scat esc prob
// >>chng 03 feb 13, from 4.998 to 4.743, new e + H3+ dr rate
// >>chng 03 apr 11, from 4.74 to 6.54 PopOpc bug
// >>chng 03 Aug 22, from 6.54 to 6.98, line overlap, incr error
// >>chng 05 mar 12, from 6.98 to 6.16 RP merge He-like code
// >>chng 06 jul 18, from 6.16 to 7.23, RP major upgrade VS coll strengths
// >>chng 07 oct 22, from 7.23 to 9.18, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 9.18 to 6.42, iteration had been disabled, turn back on
// >>chng 08 mar 17, from 6.42 to 4.92, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 mar 22, from 4.92 5.44 bugfix r1879
assert line "H  1" 6563      5.44 error 0.1
//
// >>chng 01 mar 09 from 0.119 to 0.271 change in H redis function
// >>chng 01 may 19 from 0.271 to 0.359 ct upgrade, more ionization in piz
// >>chng 01 may 21 bad opacity in HeI 7065, from 0.359 to 0.324
// >>chng 01 may 30, from 0.324 to 0.362 change in line redistribution fcn
// >>chng 01 jun 15, from 0.362 to 0.276 O ct
// >>chng 01 jul 04, from 0.276 to 0.254 diffuse fields
// >>chng 02 may 09, from 0.254 to 0.285 e scat esc prob
// >>chng 04 may 19, from 0.285 to 0.270, RP helike update
// >>chng 05 mar 12, from 0.270 to 0.240 RP merge He-like code
// >>chng 06 jul 18, from 0.240 to 0.354, RP major upgrade VS coll strengths
// >>chng 06 nov 29, from 0.354 to 0.336, Badnell DR turned on
// >>chng 07 jan 20, from 0.336 to 0.356, e2 bug
// >>chng 06 jul 18, from 0.356 to 0.296, repair common error in several collision routines
// >>chng 07 oct 22, from 0.296 to 0.394, resolve l-levels of h-like sequence
// >>chng 07 oct 25, from 0.394 to 0.435, bugfix stark broadening incorrect
// >>chng 08 feb 13, from 0.435 to 0.273, iteration had been disabled, turn back on
// >>chng 08 mar 17, from 0.273 to 0.493, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 mar 22, from 0.493 to 0.522, cutoff between volatile/static opacities was broken
// >>chng 08 apr 20, incr error to more realistic 10%
assert line "H  1" 1.875m       0.522 0.1
//
// >>chng 00 oct 03, from 102.9 to 111, cont rebin, H- opac bug
// >>chng 01 dec 28 to 110 auger yields
// >>chng 04 may 19, from 110 to 104, RP update he-like seq
// >>chng 08 mar 17, from 104 to 113, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert line "Ba C"    0      113
//
// >>chng 00 oct 03, from 54.3 to 93.1, cont rebin, H- opac bug
// >>chng 01 dec 28 to 85.01, auger yields
// >>chng 02 apr 02, from 85.09 to 80.8, better bound compton scattering
// >>chng 02 may 31, from 80.8 to 85.0, better convergence
// >>chng 04 apr 28, from 85.0 to 77.11 He+ - H ct rate update
// >>chng 04 jun 12, from 77.1 to 81.3 chng cont shield fcn
// >>chng 04 jul 03, from 81.3 to 86.5 NA update CO, H chem net
// >>chng 06 aug 09, from 86.5 to 82.6 Meudon 06, Bowen change
// >>chng 07 jul 25, from 82.6 to 78.6 bugfix pvh FeII tau r1290
// >>chng 07 oct 22, from 78.6 to 83.7, resolve l-levels of h-like sequence
// >>chng 08 mar 17, from 83.7 to 94.8, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 dec 13, from 94.8 to 90.09, merging newsolvers branch
assert line "Pa C"    0      90.09
//
// >>chng 00 oct 03, from 106 to 143.8 cont rebin, H- opac bug
// >>chng 01 may 19, from 143.8 to 292.8 ct upgrade, more ionization in piz
// >>chng 01 may 21, from 292.8 to 266.2 bad opacity in HeI 7065,
// >>chng 01 may 30, from 266.2 to 291 change in line redistribution fcn
// >>chng 01 jun 15, from 291 to 138, O ct
// >>chng 02 jan 21, from 138 to 146, chng in he trans
// >>chng 02 apr 02, from 145 to 132, better bound compton scattering,
//                   increase error
// >>chng 02 jul 23, from 132 to 116, better brems opacity
// >>chng 03 feb 28, from 116 to 128, hei coll
// >>chng 04 apr 28, from 128 to 111 He+ - H ct rate update
// >>chng 04 jul 03, from 111 to 127 NA update CO, H chem net
// >>chng 07 oct 22, from 127 to 107, resolve l-levels of h-like sequence
assert line "H FF"    0      107 error 0.1
//
// >>chng 02 nov 20, add this line
// >>chng 03 feb 28, from 0.924 to 1.549, hei coll
// >>chng 03 apr 25, from 1.549 to 1.382, bring degenerate lines together
// >>chng 03 sep 26, from 1.382 to 1.050, update HeI cs
// >>chng 04 may 19, from 1.050 to 0.735, RP update he-like seq
// >>chng 04 jul 03, from 0.735 to 0.902 NA update CO, H chem net
// >>chng 05 mar 12, from 0.902 to 0.422 RP merge He-like code
// >>chng 05 oct 01, from 0.422 to 0.448, RA chng levels
// >>chng 06 jul 18, from 0.448 to 0.897, RP major upgrade VS coll strengths
// >>chng 06 jul 18, from 0.897 to 0.745, repair common error in several collision routines
// >>chng 07 oct 22, from 0.745 to 1.080, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 1.080 to 0.805, iteration had been disabled, turn back on
// >>chng 08 mar 17, from 0.805 to 0.667, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 mar 22, from 0.667 to 0.589 bugfix r1879
// >>chng 08 apr 30, from 0.589 to 0.562 H, He topoff change
// >>chng 08 dec 13, from 0.562 to 0.5014, merging newsolvers branch
assert line "He 1" 5876      0.5014 error 0.1
//
// >>chng 00 oct 03, from 0.512 to 0.555 cont rebin, H- opac bug
// >>chng 01 mar 09 from 0.55 to 0.672 change in H redis function
// >>chng 02 dec 11, rm line new hei
// assert line "TOTL" 10830     0.677
// >>chng 02 dec 11, add this line
// >>chng 03 feb 28, from 1.249 to 1.413, hei coll
// >>chng 03 sep 26, from 1.413 to 1.310, update HeI cs
// >>chng 04 may 19, from 1.310 to 1.155, RP update he-like seq
// >>chng 04 jun 12, from 1.155 to 1.214 chng cont shield fcn
// >>chng 04 jul 03, from 1.214 to 1.529 NA update CO, H chem net
// >>chng 06 jul 18, from 1.529 to 1.731, RP major upgrade VS coll strengths
// >>chng 06 jul 18, from 1.731 to 1.561, repair common error in several collision routines
// >>chng 07 oct 22, from 1.561 to 2.311, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 2.311 to 1.577, iteration had been disabled, turn back on
// >>chng 08 mar 22, from 1.577 to 1.195 bugfix r1879
// >>chng 08 apr 03, from 1.195 to 0.937, include 10830 in fine opacity array
// >>chng 08 dec 13, from 0.937 to 0.8904, merging newsolvers branch
assert line "TOTL" 10830     0.8904
//
// >>chng 01 mar 09 from 1.23 to 1.35 change in H redis function
// chng 02 jan 21, from 1.35 to 1.437, chng in he trans
// >>chng 03 feb 28, from 1.437 to 1.520, hei coll
// >>chng 03 apr 11, from 1.52 to 1.36 PopOpc bug
// >>chng 04 jun 12, from 1.36 to 1.46 chng cont shield fcn
// >>chng 06 jul 18, from 1.46 to 1.37, RP major upgrade VS coll strengths
// >>chng 06 nov 29, from 1.37 to 1.30, Badnell DR turned on
// >>chng 07 oct 22, from 1.30 to 1.19, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 1.19 to 1.28, iteration had been disabled, turn back on
// >>chng 08 mar 17, from 1.28 to 1.44, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 nov 08, from 1.44 to 1.364, fix VS80 elec impact coll. str. (r2444)
assert line "He 2" 4686      1.364
//
// >>chng 01 mar 16 from 8.67 to 9.16 change in H atomic data
// >>chng 02 jun 13, from 9.16 to 9.71, slow drift
// >>chng 06 nov 29, from 9.71 to 8.98 Badnell DR turned on
// >>chng 08 feb 13, from 8.98 to 9.12, iteration had been disabled, turn back on
assert line "He 2" 1640      9.12
//
// >>chng 01 jul 04, from 3.58 to 4.227 diffuse fields
// >>chng 02 jun 15, from 4.227 to 4.014, slow drift
// >>chng 02 jul 01, from 4.014 to 2.162, better corr for dest of line phot
// >>chng 02 nov 20, from 2.162 to 1.794, new he diff ots, increase error
// >>chng 07 oct 22, from 1.794 to 2.242, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 2.242 to 1.863, iteration had been disabled, turn back on
assert line "c  2" 1335      1.863 error 0.1
//
// >>chng 99 aug 18, from 0.262 to 0.6555 result of bo-sequence improvement
// >>chng 00 oct 03, from 0.656 to 0.505 cont rebin, H- opac bug
// >>chng 01 mar 09 from 0.505 to 0.54 change in H redis function
// >>chng 01 jul 04, from 0.518 to 0.452 diffuse fields
// >>chng 01 dec 28 to 0.424 auger yields
// >>chng 02 jan 23, from 0.424 to 0.452, better he-like c
// >>chng 02 mar 04, from 0.452 to 0.428, increase errors, hei atomic data, 
// >>chng 02 apr 19, from 0.428 to 0.396, slow drive down, increase errors
// >>chng 02 may 09, from 0.396 to 0.496 e scat esc prob
// >>chng 02 jul 23, from 0.498 to 0.427, better brems opacity
// >>chng 02 nov 08, from 0.427 to 0.477, rec cool increase error
// >>chng 03 jan 28, from 0.477 to 0.547, drift up
// >>chng 03 may 23, from 0.547 to 0.449, H excit state escp prob
// >>chng 04 may 13, from 0.449 to 0.514 dest albedo chng
// >>chng 06 jul 18, from 0.514 to 0.653, RP major upgrade VS coll strengths
// >>chng 06 jul 18, from 0.653 to 0.500, repair common error in several collision routines
// >>chng 07 oct 26, decrease error from 0.15 to 0.10
// >>chng 07 oct 26, from 0.500 to 0.569, very close to botch after recent changes
// >>chng 08 feb 13, from 0.569 to 0.407, iteration had been disabled, turn back on
// >>chng 08 mar 22, from 0.407 to 0.367 bugfix r1879
assert line "totl" 2326      0.367 error 0.1
//
// >>chng 01 jul 04, from 2.26 to 1.99 diffuse fields
// >>chng 02 jan 26, from 1.99 to 2.09, slow drift up 
// >>chng 02 mar 07, from 2.09 to 1.96, 3-bod rec he-like, increase error
// >>chng 02 jul 19, from 1.96 to 1.76, better two photon
// >>chng 02 oct 04, from 1.76 to 1.58, slow drift down
// >>chng 02 nov 08, from 1.58 to 1.797, rec cool increase error
// >>chng 06 jul 18, from 1.79 to 2.13, RP major upgrade VS coll strengths
// >>chng 06 nov 29, from 2.13 to 2.44 Badnell DR turned on
// >>chng 08 nov 04, from 2.44 to 2.104, line RT bug fixes (r2429-35)
assert line "c  3"  977A     2.104 error 0.15
//
// >>chng 01 jun 15, from 2.79 to 2.505, O ct
// >>chng 02 jan 23, from 2.505 to 2.640, better he-like c
// >>chng 02 mar 04, from 2.64 to 2.492 increase errors, hei atomic data, 
// >>chng 02 may 12, from 2.492 to 2.704, elec scat escape
// >>chng 02 jul 23, from 2.704 to 2.28 better brems opacity
// >>chng 02 nov 08, from 2.28 to 2.79 rec cool increase error
// >>chng 04 oct 02, from 2.79 to 2.40, drift down
// >>chng 05 mar 12, from 2.40 to 2.01 RP merge He-like code
// >>chng 06 jul 18, from 2.01 to 3.01, RP major upgrade VS coll strengths
// >>chng 07 oct 03, from 3.01 to 2.59, drift down
// >>chng 08 feb 13, from 2.59 to 1.96, iteration had been disabled, turn back on
// >>chng 08 mar 17, from 1.96 to 1.59, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert line "TOTL" 1909      1.59 error 0.15
//
// >>chng 01 jul 04, from 51.2 to 45.4 diffuse fields
// >>chng 02 nov 08, from 45.4 to 47.93 rec cool
// >>chng 05 mar 12, from 47.93 to 45.25 RP merge He-like code
// >>chng 06 jul 18, from 45.25 to 47.33, RP major upgrade VS coll strengths
// >>chng 07 oct 12, from 47.33 to 45.17, drift down, no out tau first iter
// >>chng 07 oct 22, from 45.17 to 55.09, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 55.09 to 43.30, iteration had been disabled, turn back on
// >>chng 08 mar 17, from 43.34 to 46.30, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert line "TOTL" 1549      46.30
//
// >>chng 99 aug 18, from 1.07 to 0.892 result of bo-sequence improvement
// >>chng 01 jun 15, from 0.892 to 0.791, O ct rate
// >>chng 01 jul 04, from 0.791 to 0.856 diffuse fields
// >>chng 02 jan 25, from 0.856 to 0.903, slow drift up
// >>chng 02 mar 04, from 0.903 to 0.844 increase errors, hei atomic data, 
// >>chng 02 apr 26, from 0.844 to 0.861, he-iso cooling
// >>chng 02 jul 23, from 0.861 to 0.73 better brems opacity
// >>chng 02 nov 08, from 0.73 to 0.92 rec cool increase error
// >>chng 03 may 16, from 0.92 to 0.78, h line rt cleanup
// >>chng 05 mar 12, from 0.78 to 0.64 RP merge He-like code
// >>chng 06 jul 18, from 0.64 to 0.96 RP major upgrade VS coll strengths
// >>chng 06 jul 18, from 0.96 to 0.81, repair common error in several collision routines
// >>chng 07 oct 22, from 0.81 to 0.56, resolve l-levels of h-like sequence
// >>chng 08 mar 17, from 0.56 to 0.38, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert line "totl" 1750      0.38 error 0.15
//
// >>chng 01 jun 15, from 1.84 to 1.747 O ct rate
// >>chng 01 jul 04, from 1.747 to 1.636 diffuse fields
// >>chng 02 may 31, from 1.636 to 1.728 better convergence
// >>chng 02 jul 23, from 1.728 to 1.62 better brems opacity
// >>chng 02 nov 08, from 1.62 to 1.75 rec cool increase error
// >>chng 05 mar 12, from 1.75 to 1.56 RP merge He-like code
// >>chng 06 jul 18, from 1.56 to 1.74 RP major upgrade VS coll strengths
// >>chng 08 mar 17, from 1.74 to 1.16, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert line "TOTL" 1486      1.16 error 0.1
//
// >>chng 02 mar 07, from 7.84 to 7.43, 3-bod rec he-like
// >>chng 02 may 12, from 7.43 to 7.843, elec scat escape
// >>chng 03 may 02, from 7.843 to 8.25, small changes
// >>chng 04 jul 17, from 8.25 to 7.84, drift down
// >>chng 06 nov 29, from 7.84 to 7.05 Badnell DR turned on
// >>chng 08 mar 17, from 7.05 to 7.59, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert line "TOTL" 1240      7.59
//
// >>chng 00 oct 03, from 0.022 to 0.027 cont rebin, H- opac bug
// >>chng 01 mar 16 from 0.027 to 0.029 change in H atomic data
// >>chng 01 apr 01, from 0.29 to 0.28, H redistribution
// >>chng 01 may 01, from 0.028 to 0.026, reorder evaluated of rec rates
// >>chng 01 jul 21, from 0.027 to 0.024, very sensitive
// >>chng 01 dec 28 to 0.260 auger yields
// >>chng 04 jul 03, from 0.260 to 0.203 NA update CO, H chem net
// >>chng 04 nov 15, from 0.203 to 0.294 incorrect upper level for line
assert line "6lev" 8446      0.294 error 0.2
//
// >>chng 00 oct 03, from 0.021 to 0.019 cont rebin, H- opac bug
// >>chng 01 jul 04, from 0.020 to 0.014 diffuse fields
// >>chng 01 aug 04, from 0.014 to 0.015, very sensitive, error 0.1
// >>chng 01 dec 28 to 01.833 auger yields
// >>chng 03 may 02, from 1.833 to 1.56, small changes
// >>chng 03 aug 15, from 1.56 to 1.78, change in fe ionization
// >>chng 04 may 19, from 1.78 to 1.53, RP update he-like seq
// >>chng 04 jul 03, from 1.53 to 1.38 NA update CO, H chem net
// >>chng 06 jul 18, from 1.38 to 1.52 RP major upgrade VS coll strengths
// >>chng 06 nov 29, from 1.52 to 1.33 Badnell DR turned on
// >>chng 07 oct 22, from 1.33 to 1.49, resolve l-levels of h-like sequence
// >>chng 07 oct 26, from 1.49 to 1.33, due to changes in r1472
// >>chng 07 dec 19, from 1.33 to 1.47, changes in r1698
// >>chng 08 feb 13, from 1.47 to 1.18, iteration had been disabled, turn back on
assert line "6lev" 1304      1.18 error 0.1
//
// >>chng 01 jun 15, from 5.405 to 4.92, O ct
// >>chng 01 jul 04, from 4.92 to 4.67 diffuse fields
// >>chng 02 jan 25, from 4.67 to 4.91, slow drift up
// >>chng 02 mar 04, from 4.91 to 4.647 increase errors, hei atomic data, 
// >>chng 02 may 29, from 4.65 to 4.99, brems opacity
// >>chng 02 jul 23, from 4.99 to 4.36 better brems opacity
// >>chng 02 nov 08, from 4.36 to 5.16 rec cool increase error
// >>chng 03 may 15, from 5.16 to 4.59, h line rt cleanup
// >>chng 05 mar 12, from 4.59 to 4.02 RP merge He-like code
// >>chng 06 jul 18, from 4.02 to 5.45 RP major upgrade VS coll strengths
// >>chng 06 nov 29, from 5.45 to 4.94 Badnell DR turned on
// >>chng 06 jul 18, from 4.94 to 4.24, repair common error in several collision routines
// >>chng 08 feb 13, from 4.24 to 3.40, iteration had been disabled, turn back on
// >>chng 08 mar 17, from 3.40 to 2.67, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert line "TOTL" 1665      2.67 error 0.10
//
// >>chng 99 aug 18, from 4.49 to 4.34 result of bo-sequence improvement
// >>chng 01 jul 04, from 4.34 to 3.93 diffuse fields
// >>chng 02 jun 13, from 3.93 to 4.14, slow drift
// >>chng 03 jun 02, from 4.14 to 3.93, larger default H atom
// >>chng 03 Aug 22, from 3.93 to 4.16, line overlap
// >>chng 05 feb 18, from 4.16 to 3.95, incld outwrd cont in line pump
// >>chng 06 aug 09, from 3.95 to 4.16 Meudon 06, Bowen change
// >>chng 06 nov 29, from 4.16 to 4.72 Badnell DR turned on
// >>chng 08 mar 17, from 4.72 to 3.75, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert line "totl" 1402      3.75
//
// >>chng 00 oct 03, from 3.66 to 3.407 cont rebin, H- opac bug
// >>chng 01 jul 04, from 3.41 to 3.02 diffuse fields
// >>chng 02 apr 26, from 3.02 to 2.898 he-iso cooling
// >>chng 02 may 12, from 2.898 to 3.205, elec scat escape
// >>chng 02 jul 23, from 3.205 to 2.89 better brems opacity
// >>chng 02 nov 08, from 2.89 to 3.19 rec cool increase error
// >>chng 03 may 07, from 3.19 to 3.52, Lya Pesc bug
// >>chng 04 may 13, from 3.52 to 3.09, dest albedo chng
// >>chng 07 oct 22, from 3.09 to 4.48, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 4.48 to 2.23, iteration had been disabled, turn back on
assert line "TOTL" 2798      2.23 error 0.15
//
// >>chng 01 jun 15, from 0.613 to 0.556 O ct
// >>chng 01 dec 28 to 0.519 auger yields
// >>chng 02 jan 21, from 0.519 to 0.555, chng in he trans
// >>chng 02 mar 04, from 0.555 to 0.527 increase errors, hei atomic data, 
// >>chng 02 may 29, from 0.527 to 0.564, brems opacity
// >>chng 02 jun 15, from 0.564 to 0.514, slow drift
// >>chng 02 jul 23, from 0.514 to 0.403 better brems opacity
// >>chng 02 nov 08, from 0.403 to 0.5104 rec cool increase error
// >>chng 06 jul 18, from 0.510 to 0.663 RP major upgrade VS coll strengths
// >>chng 06 jul 18, from 0.663 to 0.531, repair common error in several collision routines
// >>chng 08 feb 13, from 0.531 to 0.432, iteration had been disabled, turn back on
assert line "si 3" 1207      0.432 error 0.15
//
// >>chng 00 oct 03, from 2.45 to 2.32 cont rebin, H- opac bug
// >>chng 00 nov 02, from 2.32 to 2.937, updated A for transition
// >>chng 01 jun 15, from 2.937 to 2.746 O ct
// >>chng 01 jul 04, from 2.746 to 2.578 diffuse fields
// >>chng 02 apr 26, from 2.578 to 2.394 he-iso cooling
// >>chng 02 may 12, from 2.394 to 2.713, elec scat escape
// >>chng 02 may 27, from 2.713 to 2.554, improved brems cooling 
// >>chng 02 may 29, from 2.554 to 2.734, improved brems cooling 
// >>chng 02 jun 13, from 2.73 to 2.54, slow drift
// >>chng 02 jul 23, from 2.54 to 2.23 better brems opacity
// >>chng 02 nov 08, from 2.23 to 2.73 rec cool increase error
// >>chng 03 may 15, from 2.73 to 2.41, h line rt cleanup
// >>chng 05 mar 12, from 2.41 to 1.93 RP merge He-like code
// >>chng 06 jul 18, from 1.93 to 2.76 RP major upgrade VS coll strengths
// >>chng 07 jul 25, from 2.76 to 2.48 bugfix pvh FeII tau r1290
// >>chng 06 jul 18, from 2.48 to 2.12, repair common error in several collision routines
// >>chng 07 oct 22, from 2.12 to 2.48, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 2.48 to 1.62, iteration had been disabled, turn back on
// >>chng 08 nov 04, from 1.62 to 1.446, line RT bug fixes (r2429-35)
assert line "TOTL" 1888      1.446 error 0.1
//
// >>chng 01 jun 15, from 6.05 to 5.623 O ct
// >>chng 01 dec 28 to 5.287 auger yields
// >>chng 02 jan 21, from 5.287 to 5.57 chng in he trans
// >>chng 02 mar 07, from 5.57 to 5.20, 3-bod rec he-like, increase error
// >>chng 06 jul 18, from 5.50 to 6.61 RP major upgrade VS coll strengths
// >>chng 06 nov 29, from 6.61 to 5.98 Badnell DR turned on
// >>chng 06 jul 18, from 5.98 to 5.34, repair common error in several collision routines
// >>chng 08 feb 13, from 5.34 to 4.55, iteration had been disabled, turn back on
assert line "TOTL" 1397      4.55 error 0.1
//
// >>chng 00 oct 03, from 9.34 to 1.16 cont rebin, H- opac bug
// >>chng 01 may 19 from 1.16 to 0.582 ct upgrade, more ionization in piz
// >>chng 01 may 21 bad opacity in HeI 7065, from 0.582 to 0.653
// >>chng 01 may 30, from 0.653 to 0.578 change in line redistribution fcn
// >>chng 01 jun 15, from 0.578 to 1.130 O ct
// >>chng 01 jul 04, from 1.130 to 0.985 diffuse fields
// >>chng 01 dec 28 to 0.823 auger yields
// >>chng 02 jan 29, from 0.823 to 0.870, slow drift 
// >>chng 02 mar 04, from 0.870 to 0.826 increase errors, hei atomic data, 
// >>chng 02 apr 02, from 0.826 to 0.738, better bound compton scattering,
//                   increase error
// >>chng 02 may 09, from 0.738 to 0.902 e scat esc prob
// >>chng 02 jul 23, from 0.902 to 0.732 better brems opacity
// >>chng 02 nov 08, from 0.732 to 0.8101 rec cool increase error
// >>chng 03 feb 28, from 0.810 to 0.939, hei coll
// >>chng 03 oct 04, from 0.939 to 0.782, atom2 ots
// >>chng 04 jun 01, from 0.78 to 0.213, unified calc of Lya line width
// >>chng 04 jun 12, from 0.213 to 0.0077 chng cont shield fcn
// >>chng 05 mar 12, from 0.0077 to 0.0052 RP merge He-like code
// >>chng 06 jul 18, from 0.0052 to 0.0112 RP major upgrade VS coll strengths
// >>chng 07 jul 25, from 0.0112 to 0.107 bugfix pvh FeII tau r1290
// >>chng 07 jul 27, from 0.107 to 0.0743, repair common error in several collision routines
// >>chng 07 oct 22, from 0.0743 to 0.0960, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 0.0960 to 0.0365, iteration had been disabled, turn back on
// >>chng 08 mar 17, from 0.0365 to 0.0271, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 nov 08, from 0.0271 to 0.0216, fix VS80 elec impact coll. str. (r2444)
assert line "Fe2c"    0      0.0216 error 0.15
//
// the amount of Lya lost to FeII
// >>chng 03 aug 15, from 3.225 to 5.19, add low-T dr guess for fe
// >>chng 04 apr 28, from 5.195 to 4.35 He+ - H ct rate update
// >>chng 04 may 19, from 4.35 to 4.06, RP update he-like seq
// >>chng 04 jun 01, from 4.06 to 12.9, unified calc of Lya line width
// >>chng 04 jul 03, from 12.9 to 15.0 NA update CO, H chem net
// >>chng 05 mar 27, from 15 to 13.2, incr err, should remove?
// >>chng 06 jul 18, from 13.2 to 16.3 RP major upgrade VS coll strengths
// >>chng 07 jul 25, from 16.3 to 12.8 bugfix pvh FeII tau r1290
// >>chng 08 feb 13, from 12.8 to 9.57, iteration had been disabled, turn back on
// >>chng 08 oct 18, from 9.57 to 8.233, include Badnell UTA data
assert line "Fe 2" 1216 8.233 error 0.15
//
// >>chng 03 aug 15, from 0.338 to 0.701, add low-T dr guess for fe
// >>chng 04 apr 28, from 0.701 to 0.660 He+ - H ct rate update
// >>chng 04 may 13, from 0.660 to 0.760 dest albedo chng, incr error
// >>chng 04 jun 01, from 0.76 to 0.94 unified calc of Lya line width
// >>chng 06 nov 29, from 0.94 to 0.84 Badnell DR turned on
// >>chng 07 jul 25, from 0.84 to 0.111 bugfix pvh FeII tau r1290
// >>chng 08 feb 13, from 0.111 to 0.124, iteration had been disabled, turn back on
// >>chng 08 mar 17, from 0.124 to 0.180, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 sep 22, from 0.18 to 0.1624, remove Ali et al. DR guess
// >>chng 08 nov 13, from 0.1624 to 0.1781, fix Boltzmann factor in CT rates (r2455)
assert line "Fe 2" 2500 0.1781 error 0.1
//
// >>chng 03 aug 15, from 0.194, to 0.487 low-T dr guess for fe
// >>chng 04 apr 28, from 0.487 to 0.443 He+ - H ct rate update
// >>chng 04 may 13, from 0.443 to 0.488 dest albedo chng, incr error
// >>chng 04 jun 01, from 0.49 to 1.07 unified calc of Lya line width
// >>chng 07 jul 25, from 1.07 to 0.578 bugfix pvh FeII tau r1290
// >>chng 07 oct 22, from 0.578 to 0.431, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 0.431 to 0.447, iteration had been disabled, turn back on
// >>chng 08 mar 17, from 0.447 to 0.643, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 mar 22, from 0.64 to 0.76 bugfix r1879
assert line "Fe 2" 2300 0.76 error 0.1 
//
// >>chng 01 mar 09 from 1.45 to 1.54, change in H redis function
// >>chng 01 jul 04, from 1.54 to 1.355 diffuse fields
// >>chng 03 feb 28, from 1.355 to 1.429, hei coll
// >>chng 03 nov 08, from 1.429 to 1.277, diff reset bug
// >>chng 06 mar 11, from 1.277 to 1.215, small drift down
// >>chng 07 oct 22, from 1.215 to 1.385, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from 1.201 to 1.274, iteration had been disabled, turn back on
assert line "FeKa"    1.78      1.274
//
// >>chng 01 may 19 from -0.029 to -0.092 ct upgrade, more ionization in piz
// >>chng 01 jun 15, from -0.092 to -0.045 O ct
assert ionization Hydrogen 1 -0.045
//
// >>chng 00 oct 03, from -1.19 to -0.989 cont rebin, H- opac bug
// >>chng 01 may 19 from -0.989 to -0.720 ct upgrade, more ionization in piz
// >>chng 01 may 21 bad opacity in HeI 7065, from -0.720 to -0.758
// >>chng 01 may 30, from -0.758 to -0.720 change in line redistribution fcn
// >>chng 01 jun 15, from -0.720 to -1.007 O ct
// >>chng 02 apr 02, from -1.007 to -1.031 better bound compton scattering,
//                   increase error 
// >>chng 07 oct 22, from -1.031 to -1.094, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -1.094 to -1.033, iteration had been disabled, turn back on
assert ionization Hydrogen 2 -1.033 error 0.1
//
assert ionization Helium   1 -0.002
//
// >>chng 00 oct 03, from -2.27 to -2.317 cont rebin, H- opac bug
// >>chng 00 oct 03, from -2.317 to 2.285 update sec ionizaton rates
// >>chng 01 may 19 from -2.285 to -2.343 ct upgrade, more ionization in piz
// >>chng 01 jun 15, from -2.343 to -2.302 O ct
// >>chng 01 jul 04, from -2.302 to -2.227 diffuse fields
// >>chng 02 jun 13, from -2.27 to -2.257, slow drift
// >>chng 02 nov 20, from -2.257 to -2.287 new he diff ots, increase error
// >>chng 04 may 19, from -2.29 to -2.24, RP update he-like seq
// >>chng 04 jun 12, from -2.24 to -2.20 chng cont shield fcn
// >>chng 04 jul 03, from -2.20 to -1.59 NA update CO, H chem net
// >>chng 08 feb 13, from -1.59 to -1.64, iteration had been disabled, turn back on
// >>chng 08 mar 22, from -1.64 to -2.24 bugfix r1879
assert ionization Helium   2 -2.24 error 0.1
//
// >>chng 01 jul 04, from -3.642 to -3.593 diffuse fields
// >>chng 03 jun 02, from -3.593 to -3.633 larger default H atom
// >>chng 06 jul 18, from -3.633 to -3.604 RP major upgrade VS coll strengths
// >>chng 07 oct 22, from -3.604 to -3.650, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -3.650 to -3.601, iteration had been disabled, turn back on
// >>chng 08 apr 30, from -3.601 to -3.573 H, He topoff change
assert ionization Helium   3 -3.573 error 0.1
//
// >>chng 00 oct 03, from -0.221 to -0.136 cont rebin, H- opac bug
// >>chng 03 oct 16, from -0.134 to -0.109 feii cooling, ots converg
// >>chng 08 feb 13, from -0.109 to -0.141, iteration had been disabled, turn back on
// >>chng 08 jul 08, from -0.141 to -0.121 damping wing pump rates
assert ionization Carbon   1 -0.121
//
// >>chng 00 oct 03, from -0.404 to -0.579 cont rebin, H- opac bug
// >>chng 02 apr 02, from -0.583 to -0.607 better bound compton scattering,
// >>chng 02 may 27, from -0.607 to -0.634 improved brems cooling 
// >>chng 02 dec 11, from -0.634 to -0.604, new hei
// >>chng 03 oct 16, from -0.604 to -0.661 feii cooling, ots converg
// >>chng 04 apr 28, from -0.661 to -0.699 He+ - H ct rate update
// >>chng 04 jul 03, from -0.699 to -0.657 NA update CO, H chem net
// >>chng 05 aug 03, from -0.657 to -0.690, NA update CO, H chem net, inc err
// >>chng 06 jul 18, from -0.690 to -0.627 RP major upgrade VS coll strengths
// >>chng 07 oct 03, from -0.627 to -0.670, drift down
// >>chng 07 oct 22, from -0.670 to -0.724, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -0.724 to -0.566, iteration had been disabled, turn back on
// >>chng 08 mar 17, from -0.566 to -0.618, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert ionization Carbon   2 -0.618 error 0.1
//
// >>chng 01 jul 04, from -2.471 to -2.367 diffuse fields
// >>chng 02 may 12, from -2.367 to -2.390
// >>chng 02 jul 23, from -2.390 to -2.421 better brems opacity
// >>chng 02 nov 20, from -2.421 to -2.459 new he diff ots,
// >>chng 04 jun 12, from -2.459 to -2.432 chng cont shield fcn
// >>chng 06 nov 29, from -2.432 to -2.373 Badnell DR turned on
// >>chng 07 oct 22, from -2.373 to -2.527, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -2.527 to -2.393, iteration had been disabled, turn back on
assert ionization Carbon   3 -2.393
//
// >>chng 01 may 30, from -2.984 to -3.006 change in line redistribution fcn
// >>chng 05 feb 18, from -3.006 to -2.989 incld outwrd cont in line pump
// >>chng 06 nov 29, from -2.989 to -3.084 Badnell DR turned on
assert ionization Carbon   4 -3.084
//
// >>chng 01 jul 04, from -3.752 to -3.705 diffuse fields
// >>chng 03 jun 02, from -3.705 to -3.741 larger default H atom
// >>chng 06 jul 18, from -3.741 to -3.716 RP major upgrade VS coll strengths
// >>chng 07 oct 22, from -3.716 to -3.772, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -3.772 to -3.729, iteration had been disabled, turn back on
// >>chng 08 apr 30, from -3.729 to -3.703 H, He topoff change
assert ionization Carbon   5 -3.703 error 0.1
//
assert ionization Nitrogen 1 -0.039
//
// >>chng 00 oct 03, from -1.236 to -1.072 cont rebin, H- opac bug
// >>chng 01 may 19 from -1.072 to -0.848 ct upgrade, more ionization in piz
// >>chng 01 may 21 bad opacity in HeI 7065, from -0.848 to -0.880
// >>chng 01 may 30, from -0.880 to -0.848 change in line redistribution fcn
// >>chng 01 jun 15, from -0.848 to -1.085, ct O rate
// >>chng 02 apr 02, from -1.085 to -1.118 better bound compton scattering,
//                   increase error
// >>chng 07 oct 22, from -1.120 to -1.184, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -1.184 to -1.141, iteration had been disabled, turn back on
assert ionization Nitrogen 2 -1.141 error 0.1
//
// >>chng 01 jul 04, from -2.503 to -2.393 diffuse fields
// >>chng 02 jun 13, from -2.393 to -2.434, slow drift
// >>chng 02 nov 20, from -2.434 to -2.479 new he diff ots,
// >>chng 04 jun 12, from -2.479 to -2.449 chng cont shield fcn
// >>chng 06 nov 29, from -2.449 to -2.401 Badnell DR turned on
// >>chng 07 oct 22, from -2.401 to -2.559, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -2.559 to -2.420, iteration had been disabled, turn back on
assert ionization Nitrogen 3 -2.420
//
// >>chng 09 dec 10, from -2.933 to -2.953, drift down
// >>chng 07 oct 22, from -2.953 to -2.977, resolve l-levels of h-like sequence
assert ionization Nitrogen 4 -2.977
//
// >>chng 01 jul 04, from -3.996 to -3.951 diffuse fields
// >>chng 03 jun 02, from -3.951 to -3.982 larger default H atom
// >>chng 06 jul 18, from -3.982 to -3.962 RP major upgrade VS coll strengths
// >>chng 06 nov 29, from -3.962 to -4.033 Badnell DR turned on
// >>chng 07 oct 22, from -4.033 to -4.079, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -4.079 to -4.034, iteration had been disabled, turn back on
// >>chng 08 feb 13, from -4.034 to -4.009, iteration had been disabled, turn back on
assert ionization Nitrogen 5 -4.009 error 0.1
//
// >>chng 01 jul 04, from -4.249 to -4.215 diffuse fields
// >>chng 01 dec 28 to -4.193, auger yields
// >>chng 03 dec 08, from -4.193 to -4.216, chng zoning, ots
// >>chng 06 nov 29, from -4.216 to -4.281 Badnell DR turned on
// >>chng 07 oct 22, from -4.281 to -4.327, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -4.327 to -4.283, iteration had been disabled, turn back on
assert ionization Nitrogen 6 -4.283
//
// >>chng 01 may 19 from -0.027 to -0.085 ct upgrade, more ionization in piz
// >>chng 01 jun 15, from -0.085 to -0.042 ct O rate
assert ionization Oxygen   1 -0.042
//
// >>chng 00 oct 03, from -1.255 to -1.041 cont rebin, H- opac bug
// >>chng 01 may 19 from -1.041 to -0.760 ct upgrade, more ionization in piz
// >>chng 01 may 21 bad opacity in HeI 7065, from -0.760 to -0.800
// >>chng 01 may 30, from -0.800 to -0.761 change in line redistribution fcn
// >>chng 01 jun 15, from -0.761 to -1.060 ct O rate
// >>chng 02 apr 02, from -1.060 to -1.089 better bound compton scattering,
//                   increase error
// >>chng 04 dec 08, from -1.089 to -1.134, co-ion feedback, check converg
// >>chng 07 oct 22, from -1.134 to -1.203, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -1.203 to -1.144, iteration had been disabled, turn back on
assert ionization Oxygen   2 -1.144 error 0.1
//
// >>chng 01 jul 04, from -2.397 to -2.316 diffuse fields
// >>chng 02 jun 13, from -2.316 to -2.352, slow drift
// >>chng 02 nov 20, from -2.352 to -2.378 new he diff ots,
// >>chng 04 jun 14, from -2.378 to -2.357, slow drift
// >>chng 06 nov 29, from -2.357 to -2.326 Badnell DR turned on
// >>chng 07 oct 22, from -2.326 to -2.451, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -2.451 to -2.347, iteration had been disabled, turn back on
assert ionization Oxygen   3 -2.347
//
// >>chng 00 oct 03, from -3.626 to -3.657 cont rebin, H- opac bug
// >>chng 01 dec 28 to -3.631, auger yields
// >>chng 06 nov 29, from -3.631 to -3.572 Badnell DR turned on
// >>chng 07 oct 22, from -3.572 to -3.600, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -3.600 to -3.576, iteration had been disabled, turn back on
assert ionization Oxygen   4 -3.576
//
// >>chng 01 jul 04, from -3.970 to -3.921 diffuse fields
// >>chng 03 jun 02, from -3.921 to -3.947 larger default H atom
// >>chng 06 aug 09, from -3.947 to -3.924 Meudon 06, Bowen change
// >>chng 07 mar 20, from -3.924 to -3.946, due to chng in r987
// >>chng 07 oct 22, from -3.946 to -3.991, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -3.991 to -3.942, iteration had been disabled, turn back on
// >>chng 08 feb 13, from -3.942 to -3.917, iteration had been disabled, turn back on
assert ionization Oxygen   5 -3.917
//
// >>chng 01 jul 04, from -4.449 to -4.415 diffuse fields
// >>chng 06 feb 03, from -4.415 to -4.438, guess dr to Badnell mean
// >>chng 06 nov 29, from -4.438 to -4.515 Badnell DR turned on
// >>chng 07 oct 22, from -4.515 to -4.556, resolve l-levels of h-like sequence
// >>chng 08 feb 13, from -4.556 to -4.515, iteration had been disabled, turn back on
assert ionization Oxygen   6 -4.515
//
// >>chng 02 apr 24, from 3.836 to 3.813, slow drift down
// >>chng 05 feb 18, from 3.813 to 3.832 incld outwrd cont in line pump
// >>chng 08 mar 22, from 3.832 to 3.810 bugfix r1879
assert temperature Hydrogen 1   3.810
//
// >>chng 00 oct 03, from 3.871 to 3.847 cont rebin, H- opac bug
assert temperature Hydrogen 2   3.847
//
// >>chng 02 apr 24, from 3.836 to 3.813, slow drift down
// >>chng 05 feb 18, from 3.813 to 3.832 incld outwrd cont in line pump
// >>chng 08 mar 22, from 3.832 to 3.810 bugfix r1879
assert temperature Helium  1   3.810
//
// >>chng 04 jun 12, from 4.079 to 4.056 chng cont shield fcn
// >>chng 04 jul 03, from 4.056 to 3.901 NA update CO, H chem net
// >>chng 08 mar 22, from 3.901 to 4.044 bugfix r1879
assert temperature Helium  2   4.044
//
assert temperature Helium  3   4.336
//
// >>chng 02 apr 24, from 3.835 to 3.812, slow drift down
// >>chng 05 feb 18, from 3.812 to 3.835 incld outwrd cont in line pump
// >>chng 08 mar 22, from 3.835 to 3.809 bugfix r1879
assert temperature Carbon  1   3.809
//
// >>chng 02 may 27, from 3.838 to 3.816 improved brems cooling 
// >>chng 05 feb 18, from 3.816 to 3.835 incld outwrd cont in line pump
// >>chng 08 mar 22, from 3.835 to 3.811 bugfix r1879
assert temperature Carbon  2   3.811 
//
// >>chng 05 mar 12, from 4.098 to 4.074 RP merge He-like code
// >>chng 06 jul 18, from 4.074 to 4.106 RP major upgrade VS coll strengths
// >>chng 07 oct 22, from 4.106 to 4.077, resolve l-levels of h-like sequence
// >>chng 08 apr 30, from 4.077 to 4.054 H, He topoff change
assert temperature Carbon  3   4.054 
//
assert temperature Carbon  4   4.184 
//
assert temperature Carbon  5   4.334 
//
assert temperature Carbon  6   4.339
//
// >>chng 02 apr 24, from 3.836 to 3.813, slow drift down
// >>chng 05 feb 18, from 3.813 to 3.831 incld outwrd cont in line pump
// >>chng 08 mar 22, from 3.831 to 3.810 bugfix r1879
assert temperature Nitrogen 1   3.810 
//
// >>chng 01 may 19 from 3.845 to 3.817 ct upgrade, more ionization in piz
// >>chng 05 feb 18, from 3.817 to 3.837 incld outwrd cont in line pump
// >>chng 08 mar 22, from 3.837 to 3.810 bugfix r1879
assert temperature Nitrogen 2   3.810
//
// >>chng 02 oct 04, from 4.102 to 4.079 slow drift down
// >>chng 03 feb 28, from 4.079 to 4.102, hei coll
// >>chng 05 mar 12, from 4.102 to 4.073 RP merge He-like code
// >>chng 06 jul 18, from 4.073 to 4.106 RP major upgrade VS coll strengths
// >>chng 07 oct 22, from 4.106 to 4.074, resolve l-levels of h-like sequence
// >>chng 08 mar 22, from 4.074 to 4.054, cutoff between volatile/static opacities was broken
assert temperature Nitrogen 3   4.054
//
assert temperature Nitrogen 4   4.182 
//
assert temperature Nitrogen 5   4.333 
//
assert temperature Nitrogen 6   4.335
//
// >>chng 02 apr 24, from 3.836 to 3.813, slow drift down
// >>chng 05 feb 18, from 3.813 to 3.832 incld outwrd cont in line pump
// >>chng 08 mar 22, from 3.832 to 3.810 bugfix r1879
assert temperature Oxygen  1   3.810
//
// >>chng 01 may 19 from 3.844 to 3.817 ct upgrade, more ionization in piz
// >>chng 05 feb 18, from 3.817 to 3.836 incld outwrd cont in line pump
// >>chng 08 mar 22, from 3.836 to 3.810 bugfix r1879
assert temperature Oxygen  2   3.810
//
// >>chng 05 mar 12, from 4.119 to 4.095 RP merge He-like code
// >>chng 06 jul 18, from 4.095 to 4.122 RP major upgrade VS coll strengths
// >>chng 07 oct 22, from 4.122 to 4.093, resolve l-levels of h-like sequence
// >>chng 08 apr 30, from 4.093 to 4.073 H, He topoff change
assert temperature Oxygen  3   4.073 error 0.1
//
assert temperature Oxygen  4   4.247
//
assert temperature Oxygen  5   4.336
//
assert temperature Oxygen  6   4.340
//
// >>chng 04 jun 13, from < 4 to < 5, chng cont shield fcn
// >>chng 08 dec 08, from 5 to 3, upper limit -> equals (r2515)
assert niter 3 error 0.01
//
//
// >>chng 02 nov 08, from 7.5 to 8.0, had been around 6.5
// >>chng 02 nov 24, from 8 to 9.5, caused by using iso he rec effic
// >>chng 02 dec 13, from 9.5 to 7.5, better conv
// >>chng 03 may 22, from 7.5 to 6.5, actually 5.2
// >>chng 03 dec 08, from 7.5 to 5 chng zoning, ots
// >>chng 04 feb 15, from 5 to 6, converge ionization loops auger
// >>chng 04 dec 08, from 6 to 7.5, co-ion feedback, check converg
// >>chng 08 dec 08, from 7.5 to 7.035, upper limit -> equals (r2515)
// >>chng 08 dec 13, from 7.035 to 15.32, merging newsolvers branch
assert itrzn 15.32 error 0.02
//
//
// >>chng 02 jun 13, inc from 500 to 520, finer zoning
// >>chng 02 oct 28, inc from 520 to 540, ? ?
// >>chng 02 dec 13, from 540 to 500 better conv
// >>chng 03 jun 08, from 500 to 550, N0 zoning
// >>chng 03 dec 08, from 550 to 850 chng zoning, ots
// >>chng 05 mar 27, from 850 to 600 (actually 502)
// >>chng 08 dec 08, from 600 to 551, upper limit -> equals (r2515)
assert nzone 551 error 0.01
//
//
c blr_f92.ini
c class blr 
c ========================================

This is similar to one of the BLR models presented in Ferland et al. (1992)
for the well-studied Seyfert galaxy NGC 5548.  It has a very large column
density and is marginally optically thick to electron scattering.  The spectrum
is given relative to Lya, and the intensity of this line is reset to 
produce a spectrum that is on the same intensity scale as that paper.