root/trunk/tsuite/slow/feii_blr_n11_p20.in

Revision 2539, 5.4 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
6hden 11
7phi(h) 20
8stop column density 23
9iterate convergence
10//
11// affects appearance of printout
12normalize to "totl" 1216 100
13print line faint 1  // get rid of lots of faint lines
14//
15// these are to try to speed things up
16init "c84.ini"
17//
18punch asserts "feii_blr_n11_p20.asr" last
19punch feii lines "feii_blr_n11_p20.lin" last
20punch feii column densities "feii_blr_n11_p20.fe2col" last
21punch feii populations all "feii_blr_n11_p20.lev" last
22punch overview "feii_blr_n11_p20.ovr" last
23punch dr "feii_blr_n11_p20.dr" last
24//
25// strong blr lines
26// >>chng 06 aug 20 from 6.921 to 6.949, bowen meudon trickledown
27// >>chng 06 nov 29, from 6.949 to 6.879, Badnell DR
28assert line luminosity "totl"  1035 6.879 //total of sum of both O VI comp
29//
30// >>chng 06 nov 29, from 8.960 to 8.981, Badnell DR
31// >>chng 08 apr 30, from 8.981 to 8.955 H, He topoff change
32assert line luminosity "totl"  1216 8.955 //total Hydrogen Lyman-alpha 1215.67
33//
34// >>chng 06 aug 20 from 7.479 to 7.519, bowen meudon trickledown
35// >>chng 06 nov 29, from 7.519 to 7.439, Badnell DR
36assert line luminosity "totl"  1240 7.439 //total of sum of both N V comp
37//
38//total of sum of both Si IV comp
39// >>chng 03 nov 23, from 7.747 to 7.776, Kraemer et al. guess of dr
40// >>chng 05 jul 17, from 7.78 to 7.74, first run in 6 mo
41// >>chng 06 jul 24, from 7.74 to 7.81, RP update H-like collision rates
42// >>chng 07 dec 09, from 7.81 to 7.76 cumulative RP changes in iso seq
43// >>chng 08 apr 30, from 7.76 to 7.71 H, He topoff change
44assert line luminosity "totl"  1397 7.71 error 0.1
45//
46// >>chng 06 aug 20 from 7.309 to 7.350, bowen meudon trickledown
47// >>chng 06 nov 29, from 7.350 to 7.411, Badnell DR
48// >>chng 08 mar 17, from 7.411 to 7.380 RP update line overlap
49assert line luminosity "totl"  1402 7.380 //total O IV] 1402
50//
51// >>chng 06 aug 20 from 8.607 to 8.641, bowen meudon trickledown
52// >>chng 07 aug 09, from 8.641 to 8.611, R1305 update H coll
53// >>chng 08 mar 17, from 8.611 to 8.583 RP update line overlap
54assert line luminosity "totl"  1549 8.583 //total of sum of both C IV comp
55//
56//total He II Balmer-alpha 1640
57// >>chng 05 oct 20, from 7.806 to 7.823, drift
58assert line luminosity "he 2"  1640 7.823
59//
60//total O III] 1666
61// >>chng 05 jul 17, from 7.20 to 7.14, first run in 6 mo
62// >>chng 06 jul 24, from 7.14 to 7.24, RP update H-like collision rates
63// >>chng 07 aug 09, from 7.24 to 7.19 R1305 update H coll
64// >>chng 07 dec 09, from 7.19 to 7.14 cumulative RP changes in iso seq
65// >>chng 08 mar 17, from 7.14 to 7.07 RP update line overlap
66assert line luminosity "o  3"  1666 7.07 error 0.1
67//
68//total of sum of CIII] 1908.73+
69// >>chng 03 oct 04, from 7.044 to 7.066, atom2 ots
70// >>chng 05 jul 17, from 7.04 to 6.98, first run in 6 mo
71// >>chng 06 jul 24, from 6.98 to 7.10, RP update H-like collision rates
72// >>chng 06 nov 29, from 7.10 to 7.20, Badnell DR
73// >>chng 07 aug 09, from 7.20 to 7.15 R1305 update H coll
74// >>chng 07 dec 09, from 7.15 to 7.08 cumulative RP changes in iso seq
75// >>chng 08 mar 17, from 7.08 to 6.99 RP update line overlap
76assert line luminosity "totl"  1909 6.99 error 0.1
77//
78//total of sum of Si III] 1892+
79// >>chng 05 jul 17, from 7.55 to 7.48, first run in 6 mo
80// >>chng 06 jul 24, from 7.48 to 7.60, RP update H-like collision rates
81// >>chng 07 aug 09, from 7.60 to 7.56 R1305 update H coll
82// >>chng 07 dec 09, from 7.56 to 7.49 cumulative RP changes in iso seq
83// >>chng 08 mar 17, from 7.49 to 7.42 RP update line overlap
84// >>chng 08 apr 30, from 7.42 to 7.38 H, He topoff change
85assert line luminosity "totl"  1888 7.38 error 0.1
86//
87//total sum of both Al III comp
88// >>chng 05 jul 17, from 6.91 to 6.85, first run in 6 mo
89// >>chng 06 jul 24, from 6.85 to 6.97, RP update H-like collision rates
90// >>chng 07 aug 09, from 6.97 to 6.90 R1305 update H coll
91// >>chng 07 dec 09, from 6.90 to 6.83 cumulative RP changes in iso seq
92// >>chng 08 mar 17, from 6.83 to 6.76 RP update line overlap
93// >>chng 08 apr 30, from 6.76 to 6.72 H, He topoff change
94assert line luminosity "totl"  1860 6.72 error 0.1
95//
96//total C II] 2324.69 + 2328.12
97// >>chng 05 jul 17, from 6.23 to 6.17, first run in 6 mo
98// >>chng 06 jul 24, from 6.17 to 6.31, RP update H-like collision rates
99// >>chng 07 aug 09, from 6.31 to 6.26 R1305 update H coll
100// >>chng 07 dec 09, from 6.26 to 6.15 cumulative RP changes in iso seq
101// >>chng 08 mar 17, from 6.15 to 5.98 RP update line overlap
102// >>chng 08 apr 30, from 5.98 to 5.81 H, He topoff change
103// >>chng 08 jul 08, from 5.81 to 5.92 damping wing pump rates
104assert line luminosity "totl"  2326 5.92 error 0.1
105//
106//total of sum of both Mg II comp
107// >>chng 05 oct 20, from 7.650 to 7.627, drift
108// >>chng 06 jul 24, from 7.63 to 7.69, RP update H-like collision rates
109// >>chng 07 aug 09, from 7.69 to 7.65 R1305 update H coll
110// >>chng 07 dec 09, from 7.65 to 7.58 cumulative RP changes in iso seq
111// >>chng 08 mar 17, from 7.58 to 7.52 RP update line overlap
112// >>chng 08 apr 30, from 7.52 to 7.48 H, He topoff change
113assert line luminosity "totl"  2798 7.48 error 0.1
114//
115// >>chng 08 dec 08, from 7 to 6.297, upper limit -> equals (r2515)
116// >>chng 08 dec 13, from 6.297 to 14.12, merging newsolvers branch
117assert itrzn 14.12 error 0.02
118assert niter 4 error 0.01
119// >>chng 03 nov 10, from 360 to 400, improve he0 zoning
120// >>chng 08 apr 30, from < 400 to 347 H, He topoff change
121assert nzone 347 error 0.01
122
123This is one of the 5 models that sample the LOC plane.
Note: See TracBrowser for help on using the browser.