root/trunk/tsuite/slow/feii_blr_n13_p22_Z20.in

Revision 2539, 3.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 FeII model
2table agn
3print lines column
4print lines sort intensity
5atom feii
6abundances starburst 20
7hden 13
8phi(h) 22
9stop column density 23
10iterate convergence
11//
12// affects appearance of printout
13normalize to "totl" 1216 100
14print line faint 1  // get rid of lots of faint lines
15//
16// these are to try to speed things up
17init "c84.ini"
18//
19punch asserts "feii_blr_n13_p22_Z20.asr" last
20punch feii lines "feii_blr_n13_p22_Z20.lin" last
21punch feii column densities "feii_blr_n13_p22_Z20.fe2col" last
22punch feii populations all "feii_blr_n13_p22_Z20.lev" last
23punch overview "feii_blr_n13_p22_Z20.ovr" last
24punch dr "feii_blr_n13_p22_Z20.dr"
25punch heating "feii_blr_n13_p22_Z20.het"
26punch cooling "feii_blr_n13_p22_Z20.col"
27//
28// strong blr lines
29// >>chng 06 aug 09, from 9.122 to 9.148 Meudon 06, Bowen change
30// >>chng 06 nov 29, from 9.148 to 9.048, Badnell DR
31assert line luminosity "totl"  1035 9.048  //total of sum of both O VI comp
32//
33// >>chng 06 aug 09, from 9.296 to 9.330 Meudon 06, Bowen change
34// >>chng 07 dec 09, from 9.330 to 9.264, cumulative RP changes in iso seq
35// >>chng 08 mar 17, from 9.264 to 9.242, RP update line overlap
36// >>chng 08 aug 10, from 9.242 to 9.223, drift down
37assert line luminosity "totl"  1216 9.223  //total Hydrogen Lyman-alpha 1215.67
38//
39// >>chng 06 nov 29, from 9.717 to 9.693, Badnell DR
40// >>chng 08 dec 13, from 9.693 to 9.67, merging newsolvers branch
41assert line luminosity "totl"  1240 9.67  //total of sum of both N V comp
42//
43// >>chng 06 nov 29, from 9.877 to 9.850, Badnell DR
44assert line luminosity "totl"  1397 9.850  //total of sum of both Si IV comp
45//
46// >>chng 06 nov 29, from 8.662 to 8.714, Badnell DR
47// >>chng 08 mar 17, from 8.714 to 8.656, RP update line overlap
48assert line luminosity "totl"  1402 8.656  //total O IV] 1402
49//
50// >>chng 07 jan 20, from 9.763 to 9.740, e2 bug
51// >>chng 08 mar 17, from 9.740 to 9.704, RP update line overlap
52assert line luminosity "totl"  1549 9.704  //total of sum of both C IV comp
53//
54//total He II Balmer-alpha 1640
55// >>chng 06 jul 24, from 9.155 to 9.275, RP update H-like collision rates
56// >>chng 07 dec 09, from 9.275 to 9.387, cumulative RP changes in iso seq
57// >>chng 07 dec 31, from 9.387 to 9.361 RP changes??
58// >>chng 08 mar 17, from 9.361 to 9.248, RP update line overlap
59// >>chng 08 oct 18, from 9.248 to 9.226, include Badnell UTA data
60assert line luminosity "he 2"  1640 9.226 
61//
62// >>chng 07 dec 09, from 8.903 to 8.878, cumulative RP changes in iso seq
63// >>chng 08 mar 17, from 8.878 to 8.707, RP update line overlap
64assert line luminosity "o  3"  1666 8.707  //total O III] 1666
65//
66// >>chng 08 mar 17, from 8.493 to 8.210, RP update line overlap
67assert line luminosity "totl"  1909 8.210  //total of sum of CIII] 1908.73+
68//
69// >>chng 06 nov 29, from 9.207 to 9.246, Badnell DR
70assert line luminosity "totl"  1888 9.246  //total of sum of Si III] 1892+
71//
72// >>chng 07 dec 09, from 9.586 to 9.563, cumulative RP changes in iso seq
73assert line luminosity "totl"  1860 9.563  //total sum of both Al III comp
74//
75//total C II] 2324.69 + 2328.12
76// >>chng 05 jul 17, from 8.33 to 8.36, first run in 6 mo
77// >>chng 07 dec 09, from 8.36 to 8.31 cumulative RP changes in iso seq
78// >>chng 08 mar 17, from 8.31 to 8.24 RP update line overlap
79// >>chng 08 mar 22, from 8.24 to 8.28 bugfix r1879
80assert line luminosity "totl"  2326 8.28 error 0.1 
81//
82assert line luminosity "totl"  2798 9.392  //total of sum of both Mg II comp
83//
84// >>chng 08 dec 08, from 9 to 6.583, upper limit -> equals (r2515)
85// >>chng 08 dec 13, from 6.583 to 16.86, merging newsolvers branch
86assert itrzn 16.86 error 0.02
87// >>chng 08 dec 08, from 4 to 3, upper limit -> equals (r2515)
88assert niter 3 error 0.01
89//
90// >>chng 05 oct 07, from 500 to 550, ? ?
91// >>chng 08 dec 08, from 550 to 516, upper limit -> equals (r2515)
92// >>chng 08 dec 13, from 516 to 502, merging newsolvers branch
93assert nzone 502 error 0.01
94
95This is one of the 5 models that sample the LOC plane.
Note: See TracBrowser for help on using the browser.