root/trunk/tsuite/auto/blr_n09_p22.in

Revision 2539, 2.8 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 1e20 cm2 s-1
2c
3c commands controlling continuum =========
4table agn
5phi(h) 22
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_p22.asr" last
23punch overview "blr_n09_p22.ovr" last
24punch dr "blr_n09_p22.dr" last
25punch convergence reason "blr_n09_p22.cvr"
26c
27c commands giving the asserts    =========
28// strong blr lines
29//
30//total of sum of both O VI comp
31// >>chng 06 feb 03, from 3.579 to 4.515, guess dr to Badnell mean
32// >>chng 06 apr 22, from 4.515 to 4.675, RP use Badnell DR for He-like
33// >>chng 06 nov 29, from 4.675 to 4.714, Badnell DR turned on
34// >>chng 07 jan 20, from 4.714 to 4.754, e2 bug
35// >>chng 07 may 21, from 4.754 to 4.725, fix bug in He-like DR extrapolation
36assert line luminosity "totl"  1035 4.725
37//
38//total Hydrogen Lyman-alpha 1215.67
39// >>chng 06 feb 03, from 6.259 to 6.764, guess dr to Badnell mean
40// >>chng 06 nov 29, from 6.764 to 6.811, Badnell DR turned on
41assert line luminosity "totl"  1216 6.811
42//
43//total of sum of both N V comp
44// >>chng 06 feb 03, from 1.613 to 2.647, guess dr to Badnell mean
45// >>chng 06 apr 22, from 2.647 to 2.851, RP use Badnell DR for He-like
46// >>chng 06 nov 29, from 2.851 to 2.894, Badnell DR turned on
47// >>chng 07 jan 20, from 2.894 to 2.921, e2 bug
48// >>chng 07 may 21, from 2.921 to 2.884, fix bug in He-like DR extrapolation
49assert line luminosity "totl"  1240 2.884
50//
51//total He II Balmer-alpha 1640
52// >>chng 06 feb 03, from 5.513 to 5.974, guess dr to Badnell mean
53// >>chng 06 nov 29, from 5.974 to 6.017, Badnell DR turned on
54// >>chng 07 oct 22, from 6.017 to 5.606, resolve l-levels of h-like sequence
55// >>chng 08 dec 13, from 5.606 to 5.584, merging newsolvers branch
56assert line luminosity "he 2"  1640 5.584
57//
58// >>chng 08 dec 08, from 6 to 4.533, upper limit -> equals (r2515)
59// >>chng 08 dec 13, from 4.533 to 13.07, merging newsolvers branch
60assert itrzn 13.07 error 0.02
61// >>chng 08 dec 08, from 4 to 3, upper limit -> equals (r2515)
62assert niter 3 error 0.01
63//
64// >>chng 06 feb 03, from 7.020 to 6.628, guess dr to Badnell mean
65// >>chng 06 nov 29, from 6.628 to 6.591, Badnell DR turned on
66assert hydrogen 2 temperature 6.591
67//
68// >>chng 03 dec 07, from 280 to 330, dr resolve metal changes
69// >>chng 05 dec 22, from 330 to 20 (13), very hot, why did it ever need 330?
70// >>chng 08 dec 08, from 20 to 15, upper limit -> equals (r2515)
71assert nzone 15 error 0.01
72//
73c blr_n09_p22.in
74c class blr
75c ========================================
76
77This is one of the models that sample the LOC plane.
Note: See TracBrowser for help on using the browser.