root/trunk/tsuite/auto/blr_n13_p22.in

Revision 2539, 6.9 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 1e13 cm-3, flux of H-ion phots 1e22 cm2 s-1
2c
3c commands controlling continuum =========
4table agn
5phi(h) 22
6c
7c commands for density & abundances =========
8hden 13
9c
10c commands controlling geometry  =========
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    =========
20print h-like departure hydrogen
21print h-like departure helium
22normalize to "totl" 1216 100
23print line faint 1  // get rid of lots of faint lines
24punch asserts "blr_n13_p22.asr" last
25punch overview "blr_n13_p22.ovr" last
26punch dr "blr_n13_p22.dr" last
27c
28c commands giving the asserts    =========
29// strong blr lines
30// >>chng 06 jul 18, from 9.690 to 9.530, RP major upgrade VS coll strengths
31// >>chng 06 aug 09, from 9.690 to 9.559 Meudon 06, Bowen change
32// >>chng 06 nov 29, from 9.559 to 9.452, Badnell DR turned on
33// >>chng 08 feb 20, from 9.452 to 9.474, grain charge transfer changes in ion_solver and iso*
34// >>chng 08 mar 17, from 9.474 to 9.400, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
35// >>chng 08 mar 20, from 9.400 to 9.411, dbl precision gets diff answer, take average, widen
36assert line luminosity "totl"  1035 9.411 0.10 //total of sum of both O VI comp
37//
38// >>chng 07 oct 02, increase error, slightly different results with all doubles
39assert line luminosity "totl"  1216 9.690 0.10 //total Hydrogen Lyman-alpha 1215.67
40//
41// >>chng 05 jul 08, from 9.642 to 9.659 drift down
42// >>chng 06 jul 18, from 9.659 to 9.777, RP major upgrade VS coll strengths
43// >>chng 06 aug 09, from 9.777 to 9.806 Meudon 06, Bowen change
44// >>chng 06 nov 29, from 9.806 to 9.734, Badnell DR turned on
45// >>chng 08 mar 17, from 9.734 to 9.687, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
46// >>chng 08 dec 13, from 9.687 to 9.664, merging newsolvers branch
47assert line luminosity "totl"  1240 9.664 //total of sum of both N V comp
48//
49// >>chng 06 jul 18, from 9.554 to 9.624, RP major upgrade VS coll strengths
50// >>chng 08 mar 17, from 9.624 to 9.590, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
51assert line luminosity "totl"  1397 9.590 //total of sum of both Si IV comp
52//
53// >>chng 06 jul 18, from 8.711 to 8.758, RP major upgrade VS coll strengths
54// >>chng 06 aug 09, from 8.758 to 8.779 Meudon 06, Bowen change
55// >>chng 06 nov 29, from 8.779 to 8.804, Badnell DR turned on
56// >>chng 07 oct 22, from 8.804 to 8.780, resolve l-levels of h-like sequence
57// >>chng 08 mar 17, from 8.780 to 8.736, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
58// >>chng 08 mar 20, from 8.736 to 8.747, dbl precision gets diff answer, take average, widen
59assert line luminosity "totl"  1402 8.747 0.10  //total O IV] 1402
60//
61// >>chng 06 jul 18, from 9.856 to 9.912, RP major upgrade VS coll strengths
62// >>chng 07 oct 03, from 9.912 to 9.893 drift down
63// >>chng 08 mar 17, from 9.893 to 9.870, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
64assert line luminosity "totl"  1549 9.870  //total of sum of both C IV comp
65//
66// >>chng 04 jun 12, from 9.290 to 9.255 chng cont shield fcn
67// >>chng 06 jul 18, from 9.255 to 9.512, RP major upgrade VS coll strengths
68// >>chng 07 oct 22, from 9.512 to 9.676, resolve l-levels of h-like sequence
69// >>chng 07 dec 19, from 9.676 to 9.727, changes in r1698
70// >>chng 07 dec 23, from 9.727 to 9.681, small correction to h-like collision strengths
71// >>chng 08 mar 17, from 9.681 to 9.431, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
72// >>chng 08 nov 08, from 9.431 to 9.405, fix VS80 elec impact coll. str. (r2444)
73assert line luminosity "he 2"  1640 9.405  //total He II Balmer-alpha 1640
74//
75// >>chng 06 jul 18, from 8.160 to 8.239, RP major upgrade VS coll strengths
76// >>chng 06 nov 29, from 8.239 to 8.219, Badnell DR turned on
77// >>chng 07 oct 22, from 8.219 to 8.194, resolve l-levels of h-like sequence
78// >>chng 08 mar 17, from 8.194 to 8.157, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
79assert line luminosity "o  3"  1666 8.157  //total O III] 1666
80//
81// >>chng 04 may 19, from 7.752 to 7.724 RP update helike sequence
82// >>chng 04 jun 12, from 7.724 to 7.750 chng cont shield fcn
83// >>chng 06 jul 18, from 7.750 to 7.804, RP major upgrade VS coll strengths
84// >>chng 06 nov 29, from 7.804 to 7.933, Badnell DR turned on
85// >>chng 07 oct 22, from 7.933 to 7.896, resolve l-levels of h-like sequence
86// >>chng 08 mar 17, from 7.896 to 7.869, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
87assert line luminosity "totl"  1909 7.869  //total of sum of CIII] 1908.73+
88//
89// >>chng 06 jul 18, from 8.571 to 8.627, RP major upgrade VS coll strengths
90// >>chng 06 nov 29, from 8.627 to 8.745, Badnell DR turned on
91// >>chng 07 nov 29, from 8.745 to 8.721, due to changes in r1579
92// >>chng 08 mar 20, from 8.721 to 8.710, dbl precision gets diff answer, take average, widen
93assert line luminosity "totl"  1888 8.710 0.10  //total of sum of Si III] 1892+
94//
95// >>chng 06 jul 18, from 8.968 to 9.018, RP major upgrade VS coll strengths
96// >>chng 06 nov 29, from 9.018 to 8.975, Badnell DR turned on
97// >>chng 07 nov 29, from 8.975 to 8.952, due to changes in r1579
98// >>chng 08 dec 13, from 8.952 to 8.932, merging newsolvers branch
99assert line luminosity "totl"  1860 8.932  //total sum of both Al III comp
100//
101// >>chng 04 apr 28, from 7.044 to 7.020, He+ - H ct rate update
102// >>chng 04 may 10, from 7.020 to 7.047, drift up
103// >>chng 04 set 25, from 7.047 to 7.026, chng te conv from 0.02 to 0.005
104// >>chng 06 jul 18, from 7.026 to 7.078, RP major upgrade VS coll strengths
105// >>chng 07 jan 20, from 7.078 to 7.027, e2 bug
106// >>chng 07 jul 27, from 7.027 to 6.995, repair common error in several collision routines
107// >>chng 07 oct 22, from 6.995 to 6.930, resolve l-levels of h-like sequence
108// >>chng 08 mar 20, from 6.930 to 6.920, dbl precision gets diff answer, take average, widen
109// >>chng 08 apr 30, from 6.920 to 6.877 H, He topoff change
110assert line luminosity "totl"  2326 6.877 0.10 //total C II] 2324.69 + 2328.12
111//
112// >>chng 06 jul 18, from 8.985 to 9.017, RP major upgrade VS coll strengths
113// >>chng 06 nov 29, from 9.017 to 8.982, Badnell DR turned on
114// >>chng 08 oct 18, from 8.982 to 8.955, include Badnell UTA data
115assert line luminosity "totl"  2798 8.955  //total of sum of both Mg II comp
116//
117// >>chng 04 sep 30, from 6 to 8, conv_base conv heat cool
118// >>chng 08 dec 08, from 8 to 7.394, upper limit -> equals (r2515)
119// >>chng 08 dec 13, from 7.394 to 14.38, merging newsolvers branch
120assert itrzn 14.38 error 0.02
121// >>chng 08 dec 08, from 4 to 3, upper limit -> equals (r2515)
122assert niter 3 error 0.01
123//
124// >>chng 03 nov 10, from 420 to 450, improve he0 zoning
125// >>chng 03 dec 07, from 450 to 530, dr improve metals resolution
126// >>chng 07 oct 29, from 530 to 395, change from "<" to "="
127assert nzone 395 error 0.01
128//
129c blr_n13_p22.in
130c class blr
131c ========================================
132
133This is one of the 5 models that sample the LOC plane.
Note: See TracBrowser for help on using the browser.