root/trunk/tsuite/auto/blr_n09_p18.in

Revision 2539, 4.6 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 1e09 cm-3, flux of H-ion phots 1e18 cm2 s-1
2c
3c commands controlling continuum =========
4table agn
5phi(h) 18
6c
7c commands for density & abundances =========
8hden 9
9c
10c commands controlling geometry  =========
11stop column density 23
12c
13c other commands for details     =========
14init "c84.ini"
15iterate convergence
16c these are to try to speed things up
17no level2 lines
18c
19c commands controlling output    =========
20normalize to "totl" 1216 100
21print line faint 1  // get rid of lots of faint lines
22punch asserts "blr_n09_p18.asr" last
23punch overview "blr_n09_p18.ovr" last
24punch dr "blr_n09_p18.dr" last
25punch convergence reason "blr_n09_p18.cvr"
26c
27c commands giving the asserts    =========
28// strong blr lines
29// >>chng 06 aug 09, from 4.913 to 4.943 Meudon 06, Bowen change
30// >>chng 06 nov 29, from 4.943 to 4.857, Badnell DR turned on
31assert line luminosity "totl"  1035 4.857 //total of sum of both O VI comp
32//
33// >>chng 06 aug 09, from 7.195 to 7.215 Meudon 06, Bowen change
34assert line luminosity "totl"  1216 7.215 //total Hydrogen Lyman-alpha 1215.67
35//
36// >>chng 06 aug 09, from 5.420 to 5.449 Meudon 06, Bowen change
37// >>chng 06 nov 29, from 5.449 to 5.370, Badnell DR turned on
38assert line luminosity "totl"  1240 5.370 //total of sum of both N V comp
39//
40// >>chng 03 nov 23, from 5.813 to 5.836, Kraemer et al. guess of dr
41// >>chng 06 jul 25, from 5.836 to 5.856 HN update OIII collision data
42// >>chng 07 jan 20, from 5.856 to 5.860, e2 bug
43assert line luminosity "totl"  1397 5.860 //total of sum of both Si IV comp
44//
45// >>chng 06 aug 09, from 5.498 to 5.518 Meudon 06, Bowen change
46// >>chng 06 nov 29, from 5.518 to 5.573, Badnell DR turned on
47assert line luminosity "totl"  1402 5.573 //total O IV] 1402
48//
49// >>chng 06 nov 29, from 6.908 to 6.883, Badnell DR turned on
50assert line luminosity "totl"  1549 6.883 //total of sum of both C IV comp
51//
52// >>chng 08 mar 17, from 5.738 to 5.777, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
53assert line luminosity "he 2"  1640 5.777 //total He II Balmer-alpha 1640
54//
55// >>chng 06 jul 25, from 5.756 to 5.698 HN update OIII collision data
56assert line luminosity "o  3"  1666 5.698 //total O III] 1666
57//
58// >>chng 06 nov 29, from 6.178 to 6.314, Badnell DR turned on
59assert line luminosity "totl"  1909 6.314 //total of sum of CIII] 1908.73+
60//
61// >>chng 06 nov 29, from 5.734 to 5.766, Badnell DR turned on
62assert line luminosity "totl"  1888 5.766 //total of sum of Si III] 1892+
63//
64//total sum of both Al III comp
65// >>chng 06 nov 24, from 4.619 to 4.586, turn on Badnell RR
66// >>chng 06 nov 29, from 4.586 to 4.544, Badnell DR turned on
67assert line luminosity "totl"  1860 4.544
68//
69//total C II] 2324.69 + 2328.12
70// >>chng 04 jun 01, from 5.126 to 5.095, unified calc of Lya line width
71// >>chng 04 jun 12, from 5.096 to 5.164 chng cont shield fcn
72// >>chng 04 jul 02, from 5.164 to 5.295, O H ct
73// >>chng 04 jul 03, from 5.295 to 5.171, O H ct, incr error
74// >>chng 04 jul 07, from 5.17 to 5.31, oh ct integt
75// >>chng 04 dec 08, from 5.31 to 5.17, co-ion feedback, check converg
76// >>chng 07 oct 22, from 5.17 to 5.21, resolve l-levels of h-like sequence
77// >>chng 08 jul 08, from 5.21 to 5.16 damping wing pump rates
78assert line luminosity "totl"  2326 5.16 error 0.1
79//
80//total of sum of both Mg II comp
81// >>chng 04 jun 01, from 6.215 to 6.185, unified calc of Lya line width
82// >>chng 04 jun 12, from 6.185 to 6.237 chng cont shield fcn
83// >>chng 04 jul 02, from 6.237 to 6.334, O H ct
84// >>chng 04 jul 03, from 6.334 to 6.241, O H ct, incr error
85// >>chng 04 jul 07, from 6.24 to 6.35, oh ct integt
86// >>chng 04 dec 08, from 6.35 to 6.24, co-ion feedback, check converg
87// >>chng 05 mar 25, from 6.24 to 6.27 renorm H pops in level solver
88// >>chng 08 jul 08, from 6.27 to 6.22 damping wing pump rates
89// >>chng 08 nov 13, from 6.22 to 6.175, fix Boltzmann factor in CT rates (r2455)
90assert line luminosity "totl"  2798 6.175 error 0.1
91//
92// >>chng 04 sep 25, from 6 to 8 incr temp toler 0.02 to 0.005
93// >>chng 04 sep 30, from 8 to 10, conv_base conv heat cool
94// >>chng 08 dec 08, from 10 to 7.338, upper limit -> equals (r2515)
95// >>chng 08 dec 13, from 7.338 to 15.3, merging newsolvers branch
96assert itrzn 15.3 error 0.02
97// >>chng 08 dec 08, from 4 to 3, upper limit -> equals (r2515)
98assert niter 3 error 0.01
99//
100// >>chng 03 dec 07, from 280 to 330, dr resolve metal changes
101// >>chng 08 dec 08, from 330 to 320, upper limit -> equals (r2515)
102// >>chng 08 dec 13, from 320 to 314, merging newsolvers branch
103assert nzone 314 error 0.01
104//
105c blr_n09_p18.in
106c class blr
107c ========================================
108
109This is one of the 5 models that sample the LOC plane.
Note: See TracBrowser for help on using the browser.