root/trunk/tsuite/auto/blr_rnfa.in

title table 1 of Rees et al. ApJ 347, 648
// >>refer      blr     model   Rees, M. J., Netzer, H., &
// >>refercon   Ferland, G. J. 1989, ApJ, 347, 640
c
c commands controlling continuum =========
table agn
ionization parameter -2
c
c commands for density & abundances =========
hden 10
init file="c84.ini"
abundances he -1 c -3.328 n -4.0088 o -3.0809 ne -4 na -20
continue mg -4.3768 al -5.5686 si -4.36653 s -4.76955 
continue a -5.4202 ca -5.6383 fe -4.4815 ni -20
c
c commands controlling geometry  =========
stop column density 23
c
c other commands for details     =========
iterate to convergence
c
c commands controlling output    =========
normalise "H  1" 1216 = 100
print line faint -1
punch overview last "blr_rnfa.ovr"
punch continuum last "blr_rnfa.con" units microns 
punch transmitted continuum last "blr_rnfa.trn" units rydbergs 
punch element nitrogen last "blr_rnfa.nit"
punch dr last "blr_rnfa.dr"
punch results last "blr_rnfa.rlt"
punch asserts last "blr_rnfa.asr"
c
c commands giving the asserts    =========
// >>chng 06 aug 09, from 7.645 to 7.667 Meudon 06, Bowen change
// >>chng 08 apr 30, from 7.667 to 7.645, H, He topoff change
assert line luminosity "TOTL" 1216  7.645
//
// >>chng 00 oct 03, from 3.18 to 3.411 rebin uv continuum H- opacity bug
// >>chng 01 mar 09, from 3.41 to 2.95, change in H redis fcn
// >>chng 01 dec 28, from 3.08 to 3.41, auger yields
// >>chng 02 may 09, from 3.41 to 3.191, e scat escp prob
// >>chng 03 apr 11, from 3.19 to 3.64, PopOpc bug
// >>chng 03 jun 09, from 3.64 to 3.82, larger H atom, finer zones
// >>chng 03 aug 22, from 3.82 to 3.48, line overlap
// >>chng 04 jun 12, from 3.48 to 3.67, chng cont shield fcn
// >>chng 04 jul 02, from 3.67 to 4.02 O H ct
// >>chng 04 jul 03, from 4.02 to 3.66, NA update H, CO chem net, inc err
// >>chng 04 jul 07, from 3.66 to 4.07, oh ct integ
// >>chng 04 dec 08, from 4.07 to 3.63 co-ion feedback, check converg
// >>chng 06 jul 18, from 3.63 to 3.96 RP major upgrade VS coll strengths
// >>chng 07 jul 27, from 3.96 to 3.54, repair common error in several collision routines
// >>chng 07 oct 22, from 3.54 to 3.16, resolve l-levels of h-like sequence
// >>chng 08 mar 17, from 3.16 to 2.36, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert line "TOTL" 4861      2.36 error 0.1
//
// >>chng 03 apr 11, from 12.6 to 13.7 PopOpc bug
// >>chng 04 jul 02, from 13.7 to 15.5 O H ct
// >>chng 04 jul 03, from 15.5 to 13.9 NA update H, CO chem net
// >>chng 04 jul 07, from 13.9 to 15.9 oh ct integ
// >>chng 04 dec 08, from 15.9 to 13.7 co-ion feedback, check converg
// >>chng 07 oct 22, from 13.7 to 15.6, resolve l-levels of h-like sequence
// >>chng 08 apr 20, incr error to more realistic 10%
// >>chng 08 dec 13, from 15.6 to 14.1, merging newsolvers branch
assert line "H  1" 6563      14.1 error 0.1
//
// >>chng 00 oct 03, from 1.43 to 1.322 rebin uv continuum H- opacity bug
// >>chng 01 mar 09, from 1.32 to 1.44 change in H redis fcn
// >>chng 04 jul 02, from 1.445 to 1.584 O H ct
// >>chng 04 jul 03, from 1.58 to 1.48 NA update H, CO chem net
// >>chng 04 jul 07, from 1.48 to 1.62 oh ct integ
// >>chng 04 dec 08, from 1.62 to 1.46 co-ion feedback, check converg
// >>chng 05 oct 01, from 1.46 to 1.52, drift down
// >>chng 07 jul 27, from 1.52 to 1.44, repair common error in several collision routines
// >>chng 07 oct 22, from 1.44 to 1.89, resolve l-levels of h-like sequence
// >>chng 08 mar 17, from 1.89 to 3.12, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 apr 03, from 3.12 to 3.41 include 10830 in fine opacity array
// >>chng 08 apr 10, from 3.41 to 3.21, change default h-like levels, other minor changes
// >>chng 08 aug 04, from 3.21 to 3.06, speed hack bug fix
assert line "H  1" 1.875m        3.06
//
// >>chng 00 oct 03, from 12.1 to 12.86 rebin uv continuum H- opacity bug
// >>chng 01 mar 09, from 12.9 to 14.2 change in H redis fcn
// >>chng 01 dec 28, from 14.8 to 15.6 auger yields
// >>chng 03 apr 11, from 15.6 to 13.5 PopOpc bug
// >>chng 04 jul 02, from 13.5 to 14.8 O H ct
// >>chng 04 jul 03, from 14.8 to 13.2 NA update H, CO chem net
// >>chng 04 jul 07, from 13.2 to 15.2 oh ct integ
// >>chng 04 dec 08, from 15.2 to 13.1 co-ion feedback, check converg
// >>chng 05 mar 25, from 13.1 to 14.7 renorm H pops in level solver
// >>chng 06 jul 18, from 14.7 to 13.5 RP major upgrade VS coll strengths
// >>chng 07 oct 25, from 13.5 to 12.8, bugfix stark broadening incorrect
// >>chng 08 mar 17, from 12.8 to 14.1, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 apr 10, from 14.1 to 15.9, change default h-like levels, other minor changes
// >>chng 08 apr 20, incr error to more realistic 10%
assert line "Ba C"    0      15.9 error 0.1
//
// >>chng 00 oct 03, from 4.54 to 4.86 rebin uv continuum H- opacity bug
// >>chng 01 mar 09, from 4.86 to 5.40 change in H redis fcn
// >>chng 01 dec 28, from 5.63 to 5.98 auger yields
// >>chng 03 apr 11, from 5.98 to 5.15 PopOpc bug
// >>chng 04 jul 02, from 5.15 to 5.70 O H ct
// >>chng 04 jul 03, from 5.70 to 5.04 NA update H, CO chem net
// >>chng 04 jul 07, from 5.04 to 5.87 oh ct integ
// >>chng 04 dec 08, from 5.87 to 4.98 co-ion feedback, check converg
// >>chng 05 mar 25, from 4.98 to 5.67 renorm H pops in level solver
// >>chng 06 jul 18, from 5.67 to 5.20 RP major upgrade VS coll strengths
// >>chng 07 oct 25, from 5.20 to 4.93, bugfix stark broadening incorrect
// >>chng 08 mar 17, from 4.93 to 5.43, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 apr 10, from 5.43 to 6.14, change default h-like levels, other minor changes
// >>chng 08 apr 20, incr error to more realistic 10%
assert line "Pa C"    0      6.14 error 0.1
//
// >>chng 01 mar 09, from 8.6 to 9.3 change in H redis fcn
// >>chng 02 may 27, from 9.45 10.08 better gaunts
// >>chng 02 jun 21, from 10.08 to 9.55, better gaunts
// >>chng 03 apr 11, from 9.55 to 8.50 PopOpc bug
// >>chng 04 jul 07, from 8.50 to 9.05 oh ct integ
// >>chng 04 dec 08, from 9.05 to 8.25 co-ion feedback, check converg
// >>chng 05 mar 25, from 8.25 to 8.87 renorm H pops in level solver
// >>chng 06 jul 18, from 8.87 to 8.23 RP major upgrade VS coll strengths
// >>chng 08 mar 20, from 8.23 to 8.63, value drifted but did not throw assert
// >>chng 08 apr 10, from 8.63 to 9.50, change default h-like levels, other minor changes
assert line "H FF"    0      9.50
//
// >>chng 00 oct 03, from 40.58 to 0.609 rebin uv continuum H- opacity bug
// >>chng 01 mar 09, from 0.609 to 0.689 change in H redis fcn
// >>chng 01 dec 28, from 0.690 to 0.734 auger yields
// >>chng 02 dec 11, from 0.734 to 1.235, move to new hei atom
// >>chng 03 apr 11, from 1.23 to 1.105 PopOpc bug
// >>chng 04 jul 07, from 1.105 to 1.174 oh ct integ
// >>chng 04 dec 08, from 1.174 to 1.098 co-ion feedback, check converg
// >>chng 05 mar 12, from 1.098 to 0.855, RP merge He-like code
// >>chng 05 oct 01, from 0.855 to 1.015 RP update levels
// >>chng 06 feb 03, from 1.015 to 1.063, guess dr to Badnell mean
// >>chng 06 jul 18, from 1.063 to 0.897 RP major upgrade VS coll strengths
// >>chng 07 jul 27, from 0.897 to 1.009, repair common error in several collision routines
// >>chng 07 oct 22, from 1.009 to 0.945, resolve l-levels of h-like sequence
// >>chng 08 mar 22, from 0.945 to 0.890, bugfix r1879
// >>chng 08 apr 10, from 0.890 to 1.056, change default h-like levels, other minor changes
assert line "He 1" 5876      1.056 error 0.1
//
// >>chng 00 oct 03, from 3.1 to 3.987 rebin uv continuum H- opacity bug
// >>chng 01 mar 09, from 3.99 to 4.19 change in H redis fcn
// >>chng 02 dec 11, from 4.194 to 6.45 move to new hei atom
// >>chng 03 apr 11, from 6.45 to 5.82 PopOpc bug
// >>chng 04 may 19, from 5.82 to 5.31, RP update helike sequence
// >>chng 05 mar 12, from 5.31 to 5.97 RP merge He-like code
// >>chng 05 dec 20, from 5.97 to 6.35, RP change default helium levels
// >>chng 06 jul 18, from 6.35 to 5.97 RP major upgrade VS coll strengths
// >>chng 07 oct 26, from 5.97 to 5.69, very close to botch after recent changes
// >>chng 08 mar 20, from 5.69 to 5.96, value drifted but did not throw assert
// >>chng 08 mar 22, from 5.96 to 6.26 bugfix r1879
// >>chng 08 apr 03, from 6.26 to 5.44 include 10830 in fine opacity array
assert line "TOTL" 10830     5.44 error 0.1
//
// >>chng 01 mar 09, from 0.39 to 0.42 change in H redis fcn
// >>chng 06 jul 18, from 0.420 to 0.397 RP major upgrade VS coll strengths
// >>chng 07 oct 22, from 0.397 to 0.356, resolve l-levels of h-like sequence
// >>chng 07 nov 29, from 0.356 to 0.375, due to changes in r1579
// >>chng 08 mar 17, from 0.375 to 0.499, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert line "He 2" 4686      0.499
//
// >>chng 02 nov 20, from 3.4 to 3.59, new he ots diff
// >>chng 06 aug 09, from 3.59 to 3.36 Meudon 06, Bowen change
// >>chng 08 mar 17, from 3.36 to 3.89, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
assert line "He 2" 1640      3.89
//
// >>chng 00 aug 18 from 1.01 to 0.829 result of large bo-seq atom
// >>chng 00 dec 20 from 0.829 to 0.881 result of code cleanup
// >>chng 01 mar 09, from 0.88 to 0.95 change in H redis fcn
// >>chng 01 may 19, from 0.881 to 1.002 CT network upgrade
// >>chng 02 nov 20, from 1.002 to 0.897 new he ots diff
// >>chng 04 jul 02, from 0.897 to 0.988 O H ct
// >>chng 08 apr 30, from 0.988 to 1.099 H, He topoff change
assert line "c  2" 1335      1.099 error 0.1
//
// >>chng 00 aug 18 from 0.96 to 1.048 result of large bo-seq atom
// >>chng 00 oct 03, from 1.048 to 1.582 rebin uv continuum H- opacity bug
// >>chng 00 nov 30, from 1.582 to 1.455, changed sec ioniz efficiency
// >>chng 02 jul 23, from 1.71 to 1.81, improved brems opacity
// >>chng 03 oct 16, from 1.81 to 1.90, feii cooling, ots cleanup
// >>chng 04 may 13, from 1.90 to 1.80 dest albedo chng, incr err
// >>chng 04 jun 12, from 1.80 to 2.04, chng cont shield fcn
// >>chng 04 jul 02, from 2.04 to 2.56 O H ct
// >>chng 04 jul 03, from 2.56 to 2.09 NA update H, CO chem net
// >>chng 04 jul 07, from 2.09 to 2.77 oh ct integ
// >>chng 04 dec 08, from 2.77 to 2.03 co-ion feedback, check converg
// >>chng 05 mar 25, from 2.03 to 2.33 renorm H pops in level solver
// >>chng 08 jul 08, from 2.33 to 2.10 damping wing pump rates
assert line "totl" 2326      2.10 error 0.1 
//
// >>chng 02 apr 26, from 2.52 to 2.366, he iso-seq cooling enabled
// >>chng 04 jun 12, from 2.37 to 2.52, chng cont shield fcn
// >>chng 06 nov 29, from 2.52 to 3.09 Badnell DR turned on
assert line "c  3" 977A      3.09
//
// >>chng 01 mar 09, from 9.66 to 10.2 change in H redis fcn
// >>chng 02 may 27, from 10.2 9.66 improved brems gaunts
// >>chng 06 nov 29, from 9.66 to 11.1 Badnell DR turned on
// >>chng 08 oct 18, from 11.1 to 11.61, include Badnell UTA data
assert line "TOTL" 1909      11.61
//
// >>chng 02 jul 23, from 42.9 to 40.7 improved brems opacity
// >>chng 06 nov 29, from 40.7 to 37.3 Badnell DR turned on
assert line "TOTL" 1549      37.3
//
// >>chng 00 aug 18 from 1.26 to 1.04 result of large bo-seq atom
// >>chng 06 jul 18, from 1.04 to 0.987 RP major upgrade VS coll strengths
// >>chng 06 nov 29, from 0.987 to 1.043 Badnell DR turned on
assert line "totl" 1750      1.043
//
// >>chng 06 nov 29, from 1.14 to 1.24 Badnell DR turned on
assert line "TOTL" 1486      1.24
//
// >>chng 04 apr 16, from 0.222 to 0.234, drift up
// >>chng 06 jul 18, from 0.234 to 0.221 RP major upgrade VS coll strengths
// >>chng 06 aug 09, from 0.221 to 0.234 Meudon 06, Bowen change
// >>chng 06 nov 29, from 0.234 to 0.212 Badnell DR turned on
// >>chng 08 nov 08, from 0.212 to 0.2217, fix VS80 elec impact coll. str. (r2444)
assert line "TOTL" 1240      0.2217
//
// >>chng 01 mar 09, from 0.065 to 0.077 change in H redis fcn
// >>chng 01 jul 04, from 0.080 to 0.076, chng diffuse fields
// >>chng 01 dec 28, from 0.076 to 0.084 auger yields
// >>chng 02 may 12, from 0.084 to 0.079, elec scat escape
// >>chng 03 feb 28, from 0.079 to 0.083, better hei treatment
// >>chng 03 apr 11, from 0.083 to 0.0677 PopOpc bug
// >>chng 03 jun 09, from 0.0677 to 0.0717, larger H atom, finer zones
// >>chng 03 aug 22, from 0.0717 to 0.0674 line overlap
// >>chng 04 apr 24, from 0.0674 to 0.0710, slow drift, incr error
// >>chng 04 nov 15, from 0.0710 to 0.372, incorr upper level for line
// >>chng 04 dec 08, from 0.377 to 0.278 co-ion feedback, check converg
// >>chng 05 mar 25, from 0.278 to 0.315 renorm H pops in level solver
// >>chng 07 oct 22, from 0.315 to 0.238, resolve l-levels of h-like sequence
// >>chng 08 mar 17, from 0.238 to 0.395, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 apr 30, from 0.395 to 0.349 H, He topoff change
// >>chng 08 aug 04, from 0.349 to 0.317 speed hack bug fix
assert line "6lev" 8446      0.317 error 0.1
//
// amount of Lya lost to Fe 2
// >>chng 03 oct 16, from 7.76 to 8.27, feii cooling, ots cleanup
// >>chng 04 apr 28, from 8.27 to 7.56, He+ - H ct rate update
// >>chng 04 jun 01, from 7.56 to 8.98, unified calc of Lya line width
// >>chng 04 jun 12, from 8.98 to 10.86, chng cont shield fcn
// >>chng 04 jul 02, from 10.86 to 12.97 O H ct
// >>chng 04 jul 03, from 12.97 to 11.16 NA update H, CO chem net
// >>chng 04 jul 07, from 11.16 to 14.05 oh ct integ
// >>chng 04 jul 29, from 14.05 to 14.99 oh ct integ
// this is not an observed line, but an intermediate result
// very sensitive - do not assert i
// assert line "Fe 2" 1216 14.99 
//
// >>chng 03 oct 16, from 4.58 to 2.83, feii cooling, ots cleanup
// >>chng 04 may 24, from 2.83 to 2.69, slow drift
// >>chng 04 jun 01, from 2.69 to 2.93, unified calc of Lya line width
// >>chng 04 jun 12, from 2.93 to 3.23, chng cont shield fcn
// >>chng 04 jul 02, from 3.23 to 3.55 O H ct
// >>chng 04 jul 03, from 3.55 to 3.26 NA update H, CO chem net
// >>chng 04 jul 07, from 3.26 to 3.70 oh ct integ
// >>chng 04 dec 08, from 3.70 to 3.19 co-ion feedback, check converg
// >>chng 05 mar 25, from 3.19 to 3.38 renorm H pops in level solver
// >>chng 06 nov 29, from 3.38 to 3.22 Badnell DR turned on
// >>chng 07 dec 11, from 3.22 to 3.30 drift, take ave
// >>chng 08 mar 17, from 3.30 to 3.48, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 jul 08, from 3.48 to 3.26 damping wing pump rates
assert line "Fe 2" 6200 3.26 
//
// >>chng 03 aug 15, from 1.61 to 1.74, add guess low T dr for fe
// >>chng 03 oct 16, from 1.74 to 2.00, feii cooling, ots cleanup
// >>chng 04 jun 01, from 2.00 to 2.14, unified calc of Lya line width
// >>chng 04 jun 12, from 2.13 to 2.28, chng cont shield fcn
// >>chng 05 mar 25, from 2.28 to 2.14 renorm H pops in level solver
// >>chng 07 oct 03, from 2.14 to 2.04, drift down
// >>chng 07 oct 22, from 2.04 to 2.15, resolve l-levels of h-like sequence
// >>chng 08 apr 10, from 2.15 to 1.98, change default h-like levels, other minor changes
// >>chng 08 sep 25, from 1.98 to 1.888, remove damper on line pump (r2367)
assert line "Fe 2" 2400 1.888
//
// >>chng 00 oct 03, from 0.915 to 1.21 rebin uv continuum H- opacity bug
// >>chng 01 mar 09, from 1.224 to 1.47 change in H redis fcn
// >>chng 01 jul 04, from 1.65 to 1.544 chng diffuse fields
// >>chng 01 dec 28, from 1.544 to 1.945 auger yields
// >>chng 03 feb 14, from 1.945 to 2.164, e + H3+ dr rate, overlap block
// >>chng 03 apr 11, from 2.164 to 1.727 PopOpc bug
// >>chng 03 oct 16, from 1.727 to 1.633 feii cooling, ots cleanup
// >>chng 03 nov 08, from 1.633 to 1.534, diff reset bug
// >>chng 04 mar 04, from 1.534 to 1.633, improv conr, incr error
// >>chng 04 jun 01, from 1.633 to 1.212 unified calc of Lya line width
// >>chng 04 jun 12, from 1.21 to 1.54, chng cont shield fcn
// >>chng 04 jul 02, from 1.54 to 1.81 O H ct
// >>chng 04 jul 03, from 1.81 to 1.61 NA update H, CO chem net
// >>chng 04 jul 07, from 1.61 to 1.89 oh ct integ
// >>chng 04 dec 08, from 1.89 to 1.57 co-ion feedback, check converg
// >>chng 06 jul 18, from 1.57 to 1.73 RP major upgrade VS coll strengths
// >>chng 07 oct 22, from 1.73 to 1.26, resolve l-levels of h-like sequence
// >>chng 08 mar 17, from 1.26 to 2.09, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
// >>chng 08 apr 30, from 2.09 to 1.88 H, He topoff change
// >>chng 08 nov 13, from 1.88 to 1.705, fix Boltzmann factor in CT rates (r2455)
assert line "6lev" 1304      1.705 error 0.1
//
// >>chng 06 jul 18, from 5.90 to 5.61 RP major upgrade VS coll strengths
// >>chng 06 jul 25, from 5.61 to 4.98 HN update OIII collision data
// >>chng 07 jan 23, from 4.98 to 5.21 update Gu UTA data file
assert line "TOTL" 1665      5.21
//
// >>chng 00 aug 18 from 1.4 to 1.169 result of large bo-seq atom
// >>chng 06 nov 29, from 1.17 to 1.45 Badnell DR turned on
// >>chng 08 apr 20, incr error to more realistic 10%
assert line "totl" 1402      1.45 error 0.1
//
// >>chng 00 oct 03, from 18.4 to 21.4 rebin uv continuum H- opacity bug
// >>chng 04 jul 02, from 22.5 to 25.6 O H ct
// >>chng 04 jul 03, from 25.6 to 23.5 NA update H, CO chem net
// >>chng 04 jul 07, from 23.5 to 26.4 oh ct integ
// >>chng 04 dec 08, from 26.4 to 23.1 co-ion feedback, check converg
// >>chng 05 mar 25, from 23.1 to 24.6 renorm H pops in level solver
// >>chng 06 jul 18, from 24.6 to 23.3 RP major upgrade VS coll strengths
// >>chng 06 nov 24, from 23.3 to 21.4, turn on Badnell RR
// >>chng 08 apr 03, from 21.4 to 22.4 include 10830 in fine opacity array
// >>chng 08 oct 18, from 22.4 to 21.31, include Badnell UTA data
assert line "TOTL" 2798      21.31
//
// >>chng 01 jul 04, from 1.05 to 0.987 chng diffuse fields
// >>chng 02 may 27, from 0.987 0.933 improved brems gaunts
// >>chng 04 jun 12, from 0.933 to 0.989, chng cont shield fcn
// >>chng 06 nov 29, from 0.989 to 1.04 Badnell DR turned on
// >>chng 08 apr 30, from 1.04 to 0.99 H, He topoff change
assert line "si 3" 1207      0.99 error 0.1
//
// >>chng 00 nov 02, from 8.5 to 9.14 due to updated transit prob
// >>chng 02 nov 20, from 9.14 to 8.63 new he ots diff
// >>chng 06 aug 09, from 8.63 to 8.06 Meudon 06, Bowen change
// >>chng 07 jan 23, from 8.06 to 8.43 update Gu UTA data file
assert line "TOTL" 1888      8.43
//
// >>chng 02 apr 26, from 4.098 to 3.877 he iso-seq cooling enabled
// >>chng 04 jul 03, from 3.88 to 4.07 NA update H, CO chem net
// >>chng 06 aug 09, from 4.07 to 4.29 Meudon 06, Bowen change
// >>chng 06 nov 29, from 4.29 to 4.49 Badnell DR turned on
assert line "TOTL" 1397      4.49
//
// >>chng 00 oct 03, from 2.77 to 3.45 rebin uv continuum H- opacity bug
// >>chng 01 jul 04, from 3.16 to 2.96 chng diffuse fields
// >>chng 01 dec 28, from 2.96 to 2.04 auger yields
// >>chng 02 jul 23, from 2.04 to 1.92 improved brems opacity
// >>chng 02 nov 08, from 1.92 to 2.064, fix recom cooling, increase error
// >>chng 03 apr 11, from 2.064 to 2.295 PopOpc bug
// >>chng 03 oct 16, from 2.295 to 1.963 feii cooling, ots cleanup
// >>chng 04 jun 12, from 1.963 to 1.67, chng cont shield fcn
// >>chng 05 apr 17, from 1.67 to 1.60, slow drift
// >>chng 08 mar 22, from 1.60 to 1.43 bugfix r1879
assert line "Fe2c"    0     1.43 error 0.1
//
// >>chng 06 jul 18, from 0.102 to 0.0971 RP major upgrade VS coll strengths
assert line "FeKa"    1.78      0.0971
//
assert ionization Hydrogen 1 -0.019
//
// >>chng 00 oct 03, from -1.53 to -1.418 rebin uv continuum H- opacity bug
// >>chng 01 dec 28, from -1.36 to -1.319 auger yields
// >>chng 02 nov 20, from -1.319 to -1.341 new he ots diff
// >>chng 03 dec 07, from -1.341 to -1.365 improve ots conver, finer zone
// >>chng 04 jun 01, from -1.365 to -1.400 unified calc of Lya line width
// >>chng 04 jul 02, from -1.400 to -1.262 O H ct
// >>chng 04 jul 03, from -1.262 to -1.386 NA update H, CO chem net
// >>chng 04 jul 07, from -1.386 to -1.243 oh ct integ
// >>chng 04 jul 29, from -1.243 to -1.196 chng conv chem net, inc error
// >>chng 04 dec 08, from -1.196 to -1.389 co-ion feedback, check converg
// >>chng 05 mar 25, from -1.389 to -1.27 renorm H pops in level solver
// >>chng 08 aug 05, from -1.27 to -1.31 drift down
assert ionization Hydrogen 2 -1.31 error 0.1
//
assert ionization Helium   1 -0.011
//
// >>chng 01 may 06, from -1.571 to -1.690, chng in He+ H ct rates
// >>chng 01 jul 04, from -1.698 to -1.676 chng diffuse fields
// >>chng 02 nov 20, from -1.676 to -1.708 new he ots diff
// >>chng 03 apr 11, from -1.708 to -1.745 PopOpc bug
// >>chng 04 apr 16, from -1.745 to -1.652, drift up
// >>chng 04 apr 28, from -1.652 to -1.616, He+ - H ct rate update
// >>chng 04 jun 01, from -1.616 to -1.593 unified calc of Lya line width
// >>chng 04 jul 02, from -1.593 to -1.673 O H ct
// >>chng 04 jul 03, from -1.673 to -1.380 NA update H, CO chem net, inc err
// >>chng 08 mar 22, from -1.380 to -1.78 bugfix r1879
// >>chng 08 apr 30, from -1.78 to -1.28 H, He topoff change
assert ionization Helium   2 -1.72 err 0.1
//
// >>chng 06 nov 29, from -2.672 to -2.651 Badnell DR turned on
// >>chng 08 apr 30, from -2.651 to -2.631 H, He topoff change
assert ionization Helium   3 -2.631
//
assert ionization Carbon   2 -0.006
//
// >>chng 01 mar 09, from -2.05 to -1.999 change in H redis fcn
// >>chng 01 may 19, from -1.999 to -1.973 CT network upgrade
// >>chng 01 dec 28, from -1.973 to -1.946 auger yields
// >>chng 02 nov 20, from -1.946 to -1.989 new he ots diff
// >>chng 03 apr 11, from -1.989 to -2.027 PopOpc bug
// >>chng 06 jul 18, from -02.027to -2.064 RP major upgrade VS coll strengths
// >>chng 06 nov 29, from -2.672 to -2.002 Badnell DR turned on
// >>chng 08 apr 10, from -2.002 to -1.955, change default h-like levels, other minor changes
assert ionization Carbon   3 -1.955
//
// >>chng 01 may 19, from -2.336 to -2.313 CT network upgrade
// >>chng 03 apr 11, from -2.313 to -2.342 PopOpc bug
// >>chng 06 nov 29, from -2.342 to -2.422 Badnell DR turned on
// >>chng 08 apr 30, from -2.42 to -2.40 H, He topoff change
assert ionization Carbon   4 -2.40 error 0.1
//
// >>chng 04 sep 06, from -3.252 to -3.300 bug nej Auger
// >>chng 07 jan 23, from -3.30 to -3.34 update Gu UTA data file
assert ionization Carbon   5 -3.34 error 0.1
//
// >>chng 00 oct 03, from -0.213 to -0.190 rebin uv continuum H- opacity bug
// >>chng 04 jul 02, from -0.190 to -0.167 O H ct
// >>chng 08 mar 20, from -0.167 to -0.147 value drift but did not throw assert
assert ionization Nitrogen 1 -0.147
//
// >>chng 00 oct 03, from -0.427 to -0.467 rebin uv continuum H- opacity bug
// >>chng 03 oct 16, from -0.467 to -0.489 feii cooling, ots cleanup
// >>chng 04 jul 02, from -0.489 to -0.544 O H ct
// >>chng 08 apr 30, from -0.54 to -0.57 H, He topoff change
assert ionization Nitrogen 2 -0.57 error 0.1
//
// >>chng 01 mar 09, from -2.08 to -2.02 change in H redis fcn
// >>chng 02 dec 10, from -1.966 to -2.0153, HeI 23P A
// >>chng 03 aug 22, from -2.02 to -2.06, line overlap
// >>chng 08 apr 10, from -2.06 to -2.01, change default h-like levels, other minor changes
assert ionization Nitrogen 3 -2.01 error 0.1
//
// >>chng 01 may 19, from -2.318 to -2.286 CT network upgrade
// >>chng 03 apr 11, from -2.286 to -2.315 PopOpc bug
// >>chng 08 apr 10, from -2.315 to -2.286, change default h-like levels, other minor changes
assert ionization Nitrogen 4 -2.286
//
// >>chng 03 jun 02, from -3.554 to -3.577, larger size default H 
// >>chng 06 nov 29, from -3.577 to -3.637 Badnell DR turned on
assert ionization Nitrogen 5 -3.637
//
assert ionization Oxygen   1 -0.019
//
// >>chng 00 oct 03, from -1.813 to -1.635 rebin uv continuum H- opacity bug
// >>chng 03 nov 08, from -1.529 to -1.552, diff reset bug
// >>chng 04 apr 28, from -1.552 to -1.578, He+ - H ct rate update
// >>chng 04 jun 01, from -1.578 to -1.605 unified calc of Lya line width
// >>chng 04 jun 12, from -1.605 to -1.580 chng cont shield fcn
// >>chng 04 jul 02, from -1.580 to -1.417 O H ct
// >>chng 04 jul 03, from -1.417 to -1.646 NA update H, CO chem net, inc err
// >>chng 04 jul 07, from -1.646 to -1.449 oh ct integ
// >>chng 04 jul 29, from -1.449 to -1.390 chng conv chem net, inc error
// >>chng 04 dec 08, from -1.39 to -1.65 co-ion feedback, check converg
// >>chng 05 mar 25, from -1.65 to -1.49, renorm H pops in level solver
// >>chng 08 jul 08, from -1.49 to -1.56 damping wing pump rates
assert ionization Oxygen   2 -1.56 error 0.15
//
// >>chng 01 mar 09, from -1.93 to -1.89 change in H redis fcn
// >>chng 02 dec 10, from -1.836 to -1.8828, HeI 23P A
// >>chng 03 oct 29, from -1.883 to -1.931, drift down
// >>chng 08 apr 10, from -1.931 to -1.883, change default h-like levels, other minor changes
assert ionization Oxygen   3 -1.883 error 0.1
//
// >>chng 06 nov 29, from -2.818 to -2.746 Badnell DR turned on
assert ionization Oxygen   4 -2.746
//
// >>chng 03 dec 10, from -3.576 to -3.598, dr chng, slow drift
assert ionization Oxygen   5 -3.598
//
// >>chng 03 dec 12, from -4.875 to -4.898, slow drift, change dr 
// >>chng 04 sep 06, from -4.898 to -4.927 bug nej Auger
// >>chng 06 aug 09, from -4.927 to -4.902 Meudon 06, Bowen change
// >>chng 06 nov 29, from -4.902 to -4.948 Badnell DR turned on
assert ionization Oxygen   6 -4.948
//
assert ionization Iron     2 -0.007
//
// >>chng 00 oct 03, from -2.627 to -2.503 rebin uv continuum H- opacity bug
// >>chng 00 nov 30, from -2.503 to -2.564 changed sec ioniz effic
// >>chng 03 apr 12, from -2.412 to -2.437 drift down 
// >>chng 03 aug 15, from -2.437 to -2.384 add guess low T dr for fe
// >>chng 04 apr 28, from -2.384 to -2.425, He+ - H ct rate update
// >>chng 04 jul 02, from -2.425 to -2.294 O H ct
// >>chng 04 jul 03, from -2.294 to -2.421 NA update H, CO chem net
// >>chng 04 jul 07, from -2.421 to -2.275 oh ct integ
// >>chng 04 jul 28, from -2.275 to -2.297 imp conv mole
// >>chng 04 jul 29, from -2.297 to -2.230 chng conv chem net, inc error
// >>chng 04 dec 08, from -2.23 to -2.42 co-ion feedback, check converg
// >>chng 05 mar 25, from -2.42 to -2.30, renorm H pops in level solver
// >>chng 08 jul 08, from -2.30 to -2.35 damping wing pump rates
// >>chng 08 nov 13, from -2.35 to -2.649, fix Boltzmann factor in CT rates (r2455)
assert ionization Iron     3 -2.649 error 0.1
//
// >>chng 01 mar 09, from -1.94 to -1.95 change in H redis fcn
// >>chng 02 nov 20, from -1.888 to -1.935 new he ots diff
// >>chng 03 apr 11, from -1.935 to -1.982 PopOpc bug
// >>chng 08 apr 10, from -1.982 to -1.917, change default h-like levels, other minor changes
assert ionization Iron     4 -1.917 error 0.1
//
// >>chng 03 jun 02, from -2.891 to -2.919, larger size default H 
// >>chng 03 aug 15, from -2.919 to -3.108 add guess low T dr for fe
// >>chng 03 nov 23, from -3.108 to -2.998, Kraemer et al. guess of dr
// >>chng 06 feb 03, from -2.998 to -2.965, guess dr to Badnell mean
// >>chng 06 jul 18, from -2.965 to -2.943 RP major upgrade VS coll strengths
// >>chng 06 nov 26, from -2.943 to -2.921, use Gu et al. UTA data
// >>chng 07 apr 30, from -2.921 to -2.952, log DR average, use mix S DR
// >>chng 08 apr 30, from -2.95 to -2.92 H, He topoff change
// >>chng 08 sep 22, from -2.92 to -2.868, remove Ali et al. DR guess
assert ionization Iron     5 -2.868 error 0.1
//
// >>chng 03 jun 02, from -2.856 to -2.882, larger size default H 
// >>chng 03 aug 15, from -2.882 to -2.906 add guess low T dr for fe
// >>chng 03 nov 23, from -2.906 to -3.065, Kraemer et al. guess of dr
// >>chng 03 dec 07, from -3.065 to -3.090 improve ots conver, finer zone
// >>chng 04 feb 09, from -3.090 to -3.065, chng ion_solver auger, incr error
// >>chng 04 sep 06, from -3.065 to -3.117 bug nej Auger
// >>chng 06 nov 29, from -3.117 to -3.071 Badnell DR turned on
assert ionization Iron     6 -3.071 error 0.1
//
// >>chng 03 oct 29, from -3.362 to -3.408, drift down
// >>chng 03 nov 23, from -3.408 to -3.585, Kraemer et al. guess of dr
// >>chng 06 feb 03, from -3.585 to -3.846, guess dr to Badnell mean
assert ionization Iron     7 -3.846 error 0.1
//
// >>chng 03 nov 10, from -4.338 to -4.385, chng zoning
// >>chng 03 nov 23, from -4.385 to -4.543, Kraemer et al. guess of dr
// >>chng 04 sep 06, from -4.543 to -4.633 bug nej Auger
// >>chng 06 feb 03, from -4.633 to -5.105, guess dr to Badnell mean
// >>chng 06 nov 26, from -5.105 to -5.063, use Gu et al. UTA data
// >>chng 06 nov 26, do not assert this very minor species 
// assert ionization Iron     8 -5.063 error 0.1
//
// >>chng 03 oct 16, from 3.817 to 3.840 feii cooling, ots cleanup
assert temperature Hydrogen 1  3.840
//
// >>chng 00 oct 03, from 4.019 to 3.992 rebin uv continuum H- opacity bug
// >>chng 04 jul 07, from 3.980 to 3.956 oh ct integ
// >>chng 04 dec 08, from 3.956 to 3.988 co-ion feedback, check converg
// >>chng 05 mar 25, from 3.988 to 3.957, renorm H pops in level solver
// >>chng 08 aug 04, from 3.957 to 3.976 speed hack bug fix
assert temperature Hydrogen 2  3.976
//
// >>chng 03 oct 16, from 3.817 to 3.840 feii cooling, ots cleanup
assert temperature Helium   1  3.840
//
// >>chng 01 may 06, from 4.009 to 4.058 chng in He+ H ct rates
// >>chng 04 apr 16, from 4.063 to 4.034, drift up
// >>chng 04 jul 03, from 4.034 to 3.960 NA update H, CO chem net
// >>chng 08 mar 22, from 3.960 to 4.089 bugfix r1879
assert temperature Helium   2  4.089
//
assert temperature Helium   3  4.210
//
// >>chng 03 oct 16, from 3.818 to 3.841 feii cooling, ots cleanup
assert temperature Carbon   1  3.841
//
// >>chng 03 oct 16, from 3.818 to 3.841 feii cooling, ots cleanup
assert temperature Carbon   2  3.841
//
assert temperature Carbon   3  4.137
//
assert temperature Carbon   4  4.171
//
assert temperature Carbon   5  4.211
//
assert temperature Nitrogen 1  3.829
//
// >>chng 03 sep 03, from 3.825 to 3.856, heavy brems cooling
assert temperature Nitrogen 2  3.856
//
assert temperature Nitrogen 3  4.153
//
assert temperature Nitrogen 4  4.172
//
assert temperature Nitrogen 5  4.211
//
assert temperature Nitrogen 6  4.205
//
// >>chng 03 oct 16, from 3.817 to 3.840 feii cooling, ots cleanup
assert temperature Oxygen   1  3.840
//
assert temperature Oxygen   2  3.854
//
assert temperature Oxygen   3  4.154 
//
assert temperature Oxygen   4  4.205
//
assert temperature Oxygen   5  4.214
//
assert temperature Oxygen   6  4.215
//
// >>chng 03 oct 16, from 3.818 to 3.841 feii cooling, ots cleanup
assert temperature Iron     2  3.841
//
// >>chng 03 aug 15, from 3.938 to 3.959 add guess low T dr for fe
// >>chng 04 jul 29, from 3.959 to 3.937 chng conv chem net
// >>chng 04 dec 08, from 3.937 to 3.974 co-ion feedback, check converg
// >>chng 05 mar 25, from 3.974 to 3.946, renorm H pops in level solver
// >>chng 08 nov 13, from 3.946 to 4.012, fix Boltzmann factor in CT rates (r2455)
assert temperature Iron     3  4.012
//
assert temperature Iron     4  4.150 
//
// >>chng 07 jan 23, from 4.196 to 4.218, update Gu UTA data file
assert temperature Iron     5  4.218
//
assert temperature Iron     6  4.208
//
assert temperature Iron     7  4.214
//
assert temperature Iron     8  4.214
//
//
// >>chng 03 jun 09, from 7 to 6 larger H atom, finer zones
// >>chng 04 sep 25, from 6 to 8, incr temp toler 0.02 to 0.005
// >>chng 08 dec 08, from 8 to 7.452, upper limit -> equals (r2515)
// >>chng 08 dec 13, from 7.452 to 15.39, merging newsolvers branch
// >>chng 09 jan 08, from 15.39 to 15.76, minor changes in iter/zn, iterations, and zone that crept in over last month
assert itrzn 15.76 error 0.02
//
// >>chng 03 jun 09, from 250 to 280 larger H atom, finer zones
// >>chng 03 dec 07, from 280 to 320 improve ots conver, finer zone
// >>chng 08 dec 08, from 320 to 303, upper limit -> equals (r2515)
// >>chng 09 jan 08, from 303 to 300, minor changes in iter/zn, iterations, and zone that crept in over last month
assert nzone 300 error 0.01
//
c blr_rnfa.in
c class blr
c ========================================
c 

This is the lower density cloud computed in Rees et al. (1989). 
Table 1 of that paper lists the predictions, which were a mean of those
of Hagai Netzer's ION and roughly version 76 of CLOUDY.  The lines
are generally still in good agreement with the predictions of that paper. 
In particular the changes in the line fluxes shown in Figure 1
of that paper are reproduced quite well.