Context Navigation

agn_blr_albedo.in

Revision 2578, 1.9 kB (checked in by gary, 3 weeks ago)

series of bug fixes after email from Michel Fioc - wavelengths of lines entered into lindst did not agree with energy pointed to by ipCont. added assert that the specified line wavelength and the continuum energy of ipCont are in rough agreement. this exposed a number of bugs which are fixed in this rev.

the wavelength in the argument to lindst is the one printed in the output. this can't change due to the code's history or spectroscopic convention. The continuum bands are even worse since the "ID" for the band was set by observers longs ago and often is nowhere near the center of the band. in this case a negative wavelength says not to do the check at all.

Really the printed wavelength and the continuum index should be derived from the same data. this goes back to the problem of how lines are identified in the output that we discussed at Imperial early this year.

data files -
Fe Ka wavelength changed from 2A to 1.78A - more accurate

cont_createpointers.cpp - ipCont for satellite lines was not allowed to go above ionization potential of species.
do not allow wavelength label for continuum bands to lie outside the band

lines_service.cpp - assert that wavelength of line is in rough agreement with line energies dervied from continuum index

prt_lines_continuum.cpp negative wavelength says not to check that label wavelength for continuum bands is in agreement with continuum index - wavelegth label is set by observational papers and not in agreement with center of band

Fe 22, update level energies and line wavelengths

Property svn:eol-style set to native

Line
1	title measure rayleigh scattering of Lya
2	c model from Korista, K., & Ferland, G. 1998, ApJ, 495, 672
3	c
4	c commands controlling continuum =========
5	agn 6.683 -1.20 -1.20 -0.90
6	ionization parameter 1.0
7	c
8	c commands for density & abundances =========
9	init file = "ism.ini"
10	abundances old solar 84
11	hden 11.0
12	c
13	c commands controlling geometry =========
14	stop total column density = 23.75
15	c
16	c other commands for details =========
17	iterate
18	c
19	c commands controlling output =========
20	normalize to "FeKa" 1.78
21	print line faint -1
22	print diffuse continua
23	print lines inward
24	punch continuum last "agn_blr_albedo.con" units kev
25	punch dr "agn_blr_albedo.dr"
26	punch reflected continuum last "agn_blr_albedo.ref"
27	punch emitted continuum last "agn_blr_albedo.emt"
28	punch asserts last "agn_blr_albedo.asr"
29	c
30	c commands giving the asserts =========
31	c must assert the luminosity since change in this
32	c line would appear to change the rest of the spectrum
33	assert line luminosity "FeKa" 1.78 10.070
34	assert line "nFnu" 1000 26.2
35	assert line "nInu" 1000 64.663
36	assert line "nFnu" 925.1 26.0
37	assert line "nInu" 925.1 64.05
38	assert line "nFnu" 897.8 24.8
39	//
40	// >>chng 03 dec 08, from 57.7 to 60.7, chng zoning, ots
41	assert line "nInu" 897.8 60.7
42	//
43	// >>chng 04 sep 30, from 5 to 7, conv_base conv heat cool
44	// >>chng 08 dec 08, from 7 to 4.95, upper limit -> equals (r2515)
45	// >>chng 08 dec 13, from 4.95 to 18.33, merging newsolvers branch
46	assert itrzn 18.33 error 0.02
47	// >>chng 08 dec 08, from 130 to 99, upper limit -> equals (r2515)
48	assert nzone 99 error 0.01
49	//
50	c agn_blr_albedo.in
51	c class blr
52	c ========================================
53
54	This model computes the albedo of a fairly
55	standard BLR cloud. This is the type of
56	model that was presented in the BLR albedo
57	paper by Korista & Ferland, 1998, ApJ 495, 672.
58
59	The print diffuse continua command enters continuum
60	fluxs into the emission line stack. The asserts then
61	check that these continua have the expected brightness.
62

Note: See TracBrowser for help on using the browser.

Download in other formats:

Original Format