root/trunk/tsuite/auto/blr_n12_p19.in

Revision 2539, 8.5 kB (checked in by peter, 4 weeks ago)

Merge the newsolvers branch.

This is the first merge from this branch implementing the new electron temperature and
density solvers, as well as a new method for tracking the convergence of the line optical
depths through different iterations.

In iter_track.h and iter_track.cpp the new classes iter_track and iter_track_basic are
implemented. The iter_track class iterates for the zero-point of a function based on the
van Wijngaarden-Dekker-Brent method, and is used for heavy-duty tracking of convergence
towards the root, complete with a full history of function evaluations. The
iter_track_basic class is a light-weight version that is based on a bisection search for
the zero-point and it does not keep track of the function evaluations.

The new electron density and temperature solvers are based on the iter_track class and as
a result now always bracket the solution before converging onto the solution. This should
improve the quality of the solution. The temperature solver now checks the accuracy of the
electron temperature in addition to the quality of the heating-cooling match. The old
solvers have been removed.

This merge also implements a new method based on the iter_track_basic class to track
convergence of the line optical depths through different iterations. This should improve
convergence when oscillations in the optical depth occur.

In cddefines.h new versions of the safe_div function are added that do not crash on 0/0,
but return a user-supplied argument instead.

The command PUNCH TEGRID has been removed since it is largely redundant. Use the command
PUNCH TEMPERATURE HISTORY instead. The commands SET TEMP/EDEN SOLVER have been removed
since they no longer serve any purpose (only one solver remains).

  • Property svn:eol-style set to native
Line 
1title BLR model, density 1e12 cm-3, flux of H-ion phots 1e19 cm2 s-1
2c
3c commands for density & abundances =========
4c
5c commands controlling continuum =========
6table agn
7phi(h) 19
8c
9c commands controlling geometry  =========
10hden 12
11stop column density 23
12c
13c other commands for details     =========
14iterate convergence
15c these are to try to speed things up
16init "c84.ini"
17no level2 lines
18c
19c commands controlling output    =========
20normalize to "totl" 1216 100
21print line faint 1  // get rid of lots of faint lines
22//
23punch asserts "blr_n12_p19.asr" last
24punch overview "blr_n12_p19.ovr" last
25punch dr "blr_n12_p19.dr" last
26c
27c commands giving the asserts    =========
28//
29// strong blr lines
30// >>chng 06 jul 18, from 7.757 to 7.780, RP major upgrade VS coll strengths
31// >>chng 06 aug 09, from 7.780 to 7.808 Meudon 06, Bowen change
32// >>chng 07 oct 22, from 7.808 to 7.785, resolve l-levels of h-like sequence
33// >>chng 08 mar 17, from 7.785 to 7.756, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
34assert line luminosity "totl"  1216 7.756 //total Hydrogen Lyman-alpha 1215.67
35//
36// >>chng 04 sep 25, from 4.989 to 4.961, incr temp toler 0.02 to 0.005
37// >>chng 06 aug 09, from 4.961 to 5.112 Meudon 06, Bowen change
38// >>chng 06 nov 29, from 5.112 to 5.156, Badnell DR turned on
39// >>chng 07 jul 25, from 5.156 to 5.130, bugfix pvh FeII tau r1290
40// >>chng 07 oct 22, from 5.130 to 5.105, resolve l-levels of h-like sequence
41// >>chng 07 dec 19, from 5.105 to 5.080, changes in r1698
42// >>chng 08 mar 17, from 5.080 to 5.036, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
43// >>chng 08 oct 18, from 5.036 to 5.098, include Badnell UTA data
44// >>chng 08 nov 08, from 5.098 to 5.077, fix VS80 elec impact coll. str. (r2444)
45assert line luminosity "totl"  1397 5.077 //total of sum of both Si IV comp
46//
47// >>chng 06 aug 09, from 3.145 to 3.301 Meudon 06, Bowen change
48// >>chng 06 nov 29, from 3.301 to 3.347, Badnell DR turned on
49// >>chng 07 jul 25, from 3.347 to 3.315, bugfix pvh FeII tau r1290
50// >>chng 07 jul 27, from 3.315 to 3.337, repair common error in several collision routines
51// >>chng 07 nov 29, from 3.337 to 3.308, due to changes in r1579
52// >>chng 08 mar 17, from 3.308 to 3.246, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
53// >>chng 08 aug 05 from 3.246 to 3.264, drift down
54// >>chng 08 oct 18, from 3.264 to 3.308, include Badnell UTA data
55// >>chng 08 nov 08, from 3.308 to 3.285, fix VS80 elec impact coll. str. (r2444)
56assert line luminosity "totl"  1402 3.285 //total O IV] 1402
57//
58// >>chng 06 aug 09, from 5.392 to 5.515 Meudon 06, Bowen change
59// >>chng 06 nov 29, from 5.515 to 5.485, Badnell DR turned on
60// >>chng 07 jan 20, from 5.485 to 5.452, e2 bug
61// >>chng 07 oct 22, from 5.452 to 5.420, resolve l-levels of h-like sequence
62// >>chng 08 mar 17, from 5.420 to 5.365, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
63// >>chng 08 oct 18, from 5.365 to 5.409, include Badnell UTA data
64// >>chng 08 nov 08, from 5.409 to 5.388, fix VS80 elec impact coll. str. (r2444)
65assert line luminosity "totl"  1549 5.388 //total of sum of both C IV comp
66//
67// >>chng 08 mar 17, from 6.732 to 6.799, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
68// >>chng 08 jun 25, from 6.799 to 6.823 (12.1) add source/sink to ioniz ladder
69assert line luminosity "he 2"  1640 6.823 //total He II Balmer-alpha 1640
70//
71// >>chng 05 mar 12, from 4.513 to 4.482, RP merge He-like code
72// >>chng 06 jul 18, from 4.482 to 4.534, RP major upgrade VS coll strengths
73// >>chng 06 aug 09, from 4.534 to 4.635 Meudon 06, Bowen change
74// >>chng 07 sep 15, from 4.635 to 4.613, drift down
75// >>chng 07 oct 22, from 4.613 to 4.574, resolve l-levels of h-like sequence
76// >>chng 08 mar 17, from 4.574 to 4.495, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
77// >>chng 08 oct 18, from 4.495 to 4.525, include Badnell UTA data
78// >>chng 08 nov 08, from 4.525 to 4.505, fix VS80 elec impact coll. str. (r2444)
79assert line luminosity "o  3"  1666 4.505 //total O III] 1666
80//
81// >>chng 05 mar 12, from 4.511 to 4.482, RP merge He-like code
82// >>chng 06 jul 18, from 4.482 to 4.531, RP major upgrade VS coll strengths
83// >>chng 06 aug 09, from 4.531 to 4.604 Meudon 06, Bowen change
84// >>chng 06 nov 29, from 4.604 to 4.632, Badnell DR turned on
85// >>chng 07 jul 27, from 4.632 to 4.604, repair common error in several collision routines
86// >>chng 07 oct 22, from 4.604 to 4.567, resolve l-levels of h-like sequence
87// >>chng 08 mar 17, from 4.567 to 4.497, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
88// >>chng 08 oct 18, from 4.497 to 4.521, include Badnell UTA data
89assert line luminosity "totl"  1909 4.521 //total of sum of CIII] 1908.73+
90//
91// >>chng 05 mar 12, from 6.384 to 6.333, RP merge He-like code
92// >>chng 06 jul 18, from 6.333 to 6.429, RP major upgrade VS coll strengths
93// >>chng 06 aug 09, from 6.429 to 6.469 Meudon 06, Bowen change
94// >>chng 07 jul 25, from 6.469 to 6.444, bugfix pvh FeII tau r1290
95// >>chng 07 sep 15, from 6.444 to 6.423, drift down
96// >>chng 07 oct 22, from 6.423 to 6.379, resolve l-levels of h-like sequence
97// >>chng 07 dec 19, from 6.379 to 6.352, changes in r1698
98// >>chng 08 mar 17, from 6.352 to 6.283, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
99// >>chng 08 nov 04, from 6.283 to 6.303, line RT bug fixes (r2429-35)
100// >>chng 08 nov 08, from 6.303 to 6.277, fix VS80 elec impact coll. str. (r2444)
101assert line luminosity "totl"  1888 6.277 //total of sum of Si III] 1892+
102//
103// >>chng 05 mar 12, from 5.590 to 5.562, RP merge He-like code
104// >>chng 06 jul 18, from 5.562 to 5.634, RP major upgrade VS coll strengths
105// >>chng 06 aug 09, from 5.634 to 5.696 Meudon 06, Bowen change
106// >>chng 06 nov 24, from 5.696 to 5.674, turn on Badnell RR
107// >>chng 06 nov 29, from 5.674 to 5.520, Badnell DR turned on
108// >>chng 07 jul 27, from 5.520 to 5.487, repair common error in several collision routines
109// >>chng 07 oct 22, from 5.487 to 5.453, resolve l-levels of h-like sequence
110// >>chng 08 mar 17, from 5.453 to 5.387, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
111// >>chng 08 oct 18, from 5.387 to 5.423, include Badnell UTA data
112assert line luminosity "totl"  1860 5.423 //total sum of both Al III comp
113//
114//total C II] 2324.69+2328.12
115// >>chng 03 sep 07, from 5.686 to 5.709, ? ?
116// >>chng 03 sep 10, from 5.709 to 5.771, changes in H2* population
117// >>chng 03 sep 12, from 5.771 to 5.725, changes in H2* population, incr err
118// >>chng 04 jun 01, from 5.725 to 5.816, unified calc of Lya line width
119// >>chng 04 jul 03, from 5.816 to 5.891, NA update CO, H chem net
120// >>chng 05 mar 12, from 5.891 to 5.853, RP merge He-like code
121// >>chng 06 jul 18, from 5.853 to 5.970, RP major upgrade VS coll strengths
122// >>chng 07 jul 27, from 5.970 to 5.913, repair common error in several collision routines
123// >>chng 08 mar 19, from 5.913 to 5.86, drift down with gcc/cygwin and icl
124// >>chng 08 mar 22, from 5.86 to 5.78 bugfix r1879
125// >>chng 08 apr 10, from 5.78 to 5.67, change default h-like levels, other minor changes
126assert line luminosity "totl"  2326 5.67 err 0.1
127//
128// >>chng 04 jul 20, from 7.533 to 7.552, drift up
129// >>chng 05 mar 12, from 7.552 to 7.531, RP merge He-like code
130// >>chng 06 jul 18, from 7.531 to 7.583, RP major upgrade VS coll strengths
131// >>chng 07 jul 27, from 7.583 to 7.559, repair common error in several collision routines
132// >>chng 07 oct 22, from 7.559 to 7.532, resolve l-levels of h-like sequence
133// >>chng 08 mar 17, from 7.532 to 7.484, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
134// >>chng 08 apr 10, from 7.484 to 7.462, change default h-like levels, other minor changes
135// >>chng 08 nov 13, from 7.462 to 7.442, fix Boltzmann factor in CT rates (r2455)
136assert line luminosity "totl"  2798 7.442 //total of sum of both Mg II comp
137//
138// >>chng 04 feb 27, from < 6 to < 5, dr on chng eden
139// >>chng 04 jun 04, from 5 to 6.5,
140// >>chng 04 sep 25, from 6.5 to 7.5, incr temp toler 0.02 to 0.005
141// >>chng 04 sep 30, from 7.5 to 9, conv_base conv heat cool
142// >>chng 08 dec 08, from 9 to 8.132, upper limit -> equals (r2515)
143// >>chng 08 dec 13, from 8.132 to 15.66, merging newsolvers branch
144assert itrzn 15.66 error 0.02
145//
146// >>chng 04 sep 08, from <4 to <5, small changes in logic
147// >>chng 08 dec 08, from 5 to 3, upper limit -> equals (r2515)
148assert niter 3 error 0.01
149//
150// >>chng 04 feb 27, from 330 to 410, dr on chng eden
151// >>chng 07 oct 29, from 410 to 304, change from "<" to "="
152// >>chng 08 dec 13, from 304 to 301, merging newsolvers branch
153assert nzone 301 error 0.01
154c
155c blr_n12_p19.in
156c class blr
157c ========================================
158
159This is one of the 5 models that sample the LOC plane.
Note: See TracBrowser for help on using the browser.