root/trunk/tsuite/auto/blr_n09_p18_Z20.in

Revision 2614, 6.2 kB (checked in by gary, 8 hours ago)

small changes in iter/zn, iterations, and zones, to account for small changes that have built up over last month

  • Property svn:eol-style set to native
Line 
1title BLR model, density 1e09 cm-3, flux of H-ion phots 1e18 cm2 s-1, Z=20
2c
3c commands controlling continuum =========
4table agn
5phi(h) 18
6c
7c commands for density & abundances =========
8hden 9
9abundances starburst 20
10c
11c commands controlling geometry  =========
12stop column density 23
13c
14c other commands for details     =========
15iterate to convergence
16c these are to try to speed things up
17init "c84.ini"
18no level2 lines
19c
20c commands controlling output    =========
21normalize to "totl" 1216 100
22print h-like hydrogen departure coefficients
23print line faint 1  // get rid of lots of faint lines
24punch asserts "blr_n09_p18_Z20.asr" last
25punch overview "blr_n09_p18_Z20.ovr"
26punch dr "blr_n09_p18_Z20.dr"
27punch convergence reason "blr_n09_p18_Z20.cvr"
28punch convergence error "blr_n09_p18_Z20.cve"
29c
30c commands giving the asserts    =========
31// strong blr lines
32// >>chng 04 mar 04, from 3.696 to 3.806, chng convergence criteria
33// >>chng 04 sep 06, from 3.806 to 3.792, bug nej Auger
34// >>chng 04 oct 02, from 3.792 to 3.778, drift down
35assert line luminosity "totl"  1035 3.778 //total of sum of both O VI comp
36//
37//total Hydrogen Lyman-alpha 1215.67
38// >>chng 04 jun 01, from 6.796 to 6.900, unified calc of Lya line width
39// >>chng 06 nov 29, from 6.900 to 6.944, Badnell DR turned on
40// >>chng 07 jul 25, from 6.944 to 6.987, bugfix pvh FeII tau r1290
41assert line luminosity "totl"  1216 6.987
42//
43//total of sum of both N V comp
44// >>chng 06 nov 29, from 5.621 to 5.566, Badnell DR turned on
45assert line luminosity "totl"  1240 5.566
46//
47//total of sum of both Si IV comp
48// >>chng 04 jul 28, from 5.731 to 5.753, drift
49// >>chng 06 aug 09, from 5.753 to 5.777 Meudon 06, Bowen change
50// >>chng 06 nov 29, from 5.777 to 5.817, Badnell DR turned on
51assert line luminosity "totl"  1397 5.817
52//
53// >>chng 06 nov 29, from 5.163 to 5.214, Badnell DR turned on
54assert line luminosity "totl"  1402 5.214 //total O IV] 1402
55//
56// >>chng 07 jan 20, from 5.980 to 5.955, e2 bug
57assert line luminosity "totl"  1549 5.955 //total of sum of both C IV comp
58//
59// >>chng 06 nov 29, from 5.394 to 5.415, Badnell DR turned on
60// >>chng 07 oct 22, from 5.415 to 5.388, resolve l-levels of h-like sequence
61// >>chng 08 feb 13, from 5.388 to 5.449, r1844 - mostly iso PopOpc and rt_line_one_tauinc changes
62assert line luminosity "he 2"  1640 5.449 //total He II Balmer-alpha 1640
63//
64// >>chng 06 jul 25, from 5.464 to 5.382 HN update OIII collision data
65// >>chng 06 nov 29, from 5.382 to 5.418, Badnell DR turned on
66assert line luminosity "o  3"  1666 5.418 //total O III] 1666
67//
68// >>chng 04 sep 23, from 6.062 to 6.041, slow ddrift down
69// >>chng 06 nov 29, from 6.041 to 6.113, Badnell DR turned on
70assert line luminosity "totl"  1909 6.113 //total of sum of CIII] 1908.73+
71//
72// >>chng 04 jun 01, from 6.378 to 6.343, unified calc of Lya line width
73// >>chng 07 jul 25, from 6.343 to 6.321, bugfix pvh FeII tau r1290
74// >>chng 07 oct 22, from 6.321 to 6.278, resolve l-levels of h-like sequence
75// >>chng 07 oct 26, from 6.278 to 6.323, due to changes in r1472
76// >>chng 07 nov 17, from 6.323 to 6.282, changes in r1512
77// >>chng 07 nov 23, from 6.282 to 6.316, result with 32 bit vs03
78assert line luminosity "totl"  1888 6.316 //total of sum of Si III] 1892+
79//
80// >>chng 04 jun 01, from 5.432 to 5.372, unified calc of Lya line width
81// >>chng 07 jul 25, from 5.372 to 5.352, bugfix pvh FeII tau r1290
82assert line luminosity "totl"  1860 5.352 //total sum of both Al III comp
83//
84//total C II] 2324.69 + 2328.12
85// >>chng 04 may 19, from 5.715 to 5.686 RP update helike sequence
86// >>chng 04 jun 01, from 5.686 to 5.634, unified calc of Lya line width
87// >>chng 04 jul 02, from 5.634 to 5.914, O H ct
88// >>chng 04 jun 03, from 5.914 to 5.622, NA update CO, H chem net, incr err
89// >>chng 04 jul 07, from 5.62 to 6.12, oh ct integ
90// >>chng 04 jul 16, from 6.12 to 5.95, chng order of eval of heav elem, inc er
91// >>chng 04 jul 29, from 5.96 to 6.22, impr conv h chem, incr error
92// >>chng 04 sep 25, from 6.22 to 6.21, incr temp toler 0.02 to 0.005
93// >>chng 05 jul 17, had been missing, reset
94// >>chng 07 nov 17, from 5.926 to 6.006, changes in r1512
95// >>chng 07 nov 23, from 6.006 to 5.917, result with 32 bit vs03
96// >>chng 08 apr 30, from 5.92 to 5.77, H, He topoff change
97// >>chng 08 nov 13, from 5.77 to 5.683, fix Boltzmann factor in CT rates (r2455)
98assert line luminosity "totl"  2326 5.683 error 0.2
99//
100//total of sum of both Mg II comp
101// >>chng 04 jun 01, from 6.629 to 6.578, unified calc of Lya line width
102// >>chng 04 jun 12, from 6.578 to 6.616 chng cont shield fcn
103// >>chng 04 jul 02, from 6.615 to 6.745, O H ct
104// >>chng 04 jun 03, from 6.745 to 6.614, NA update CO, H chem net, incr err
105// >>chng 04 jul 07, from 6.61 to 6.88, oh ct integ
106// >>chng 04 jul 16, from 6.88 to 6.76, chng order of eval of heav elem, inc er
107// >>chng 04 jul 29, from 6.76 to 6.94, impr conv h chem, incr error
108// >>chng 04 sep 25, from 6.94 to 6.79, incr temp toler 0.02 to 0.005
109// >>chng 04 dec 08, from 6.79 to 6.61 co-ion feedback, check converg
110// >>chng 05 mar 25, from 6.61 to 6.74 renorm H pops in level solver
111// >>chng 08 apr 30, from 6.74 to 6.63, H, He topoff change
112assert line luminosity "totl"  2798 6.63 error 0.2
113//
114// >>chng 04 sep 25, from 6.6 to 8.5 incr temp toler 0.02 to 0.005
115// >>chng 08 dec 08, from 8.5 to 6.686, upper limit -> equals (r2515)
116// >>chng 08 dec 13, from 6.686 to 18.5, merging newsolvers branch
117assert itrzn 18.5 error 0.02
118// >>chng 07 dec 06, from 5 to 7, changes in h-like collisions
119// >>chng 07 dec 19, from 7 to 9, changes in r1698
120// >>chng 08 dec 08, from 9 to 3, upper limit -> equals (r2515)
121assert niter 3 error 0.01
122//
123// >>chng 03 nov 10, from 390 to 450, improve he0 zoning
124// >>chng 03 dec 07, from 450 to 550, dr improve metals resolution
125// >>chng 05 dec 22, from 550 to 420 (371), rationalize zone asserts
126// >>chng 08 apr 30, from 420 to 460 (444) H, He topoff change
127// >>chng 08 dec 08, from 460 to 455, upper limit -> equals (r2515)
128// >>chng 08 dec 13, from 455 to 400, merging newsolvers branch
129// >>chng 09 jan 08, from 400 to 407, minor changes in iter/zn, iterations, and zone that crept in over last month
130assert nzone 407 error 0.01
131//
132c blr_n09_p18_Z20.in
133c class blr
134c ========================================
135
136This is one of the 5 models that sample the LOC plane.
Note: See TracBrowser for help on using the browser.