Navigation

Changeset 1996

Timestamp:

04/30/08 09:48:07 (2 weeks ago)

Author:

gary

Message:

further rationalization

Files:

r1994	r1996
32	32	print departure coefficients
33	33	print last iteration
34		punch asserts "~~h_casea~~.asr" last
	34	punch asserts "limit_casea_h_den13.asr" last
35	35	c
36	36	c commands giving the asserts =========
…	…
51	51	// >>chng 07 oct 22, from -3.295 to -3.341, resolve l-levels of h-like sequence
52	52	assert hydrogen 1 ionization -3.341
53		c ~~h_casea~~.in
	53	c limit_casea_h_den13.in
54	54	c class limit
55	55	c ========================================

r1994	r1996
14	14	print populaions h-like
15	15	print departure h-like
	16	punch asserts "limit_caseb_h_hs87.asr"
16	17	c
17	18	assert Case B H-like hydrogen error 0.15 range 1,000 30,000A
…	…
55	56	assert line "H 1" 1.588m 0.0033 error 0.02
56	57	assert line "H 1" 1.570m 0.0027 error 0.02
	58	c
	59	c limit_caseb_h_hs87.in
	60	c class limit
	61

r1994	r1996
29	29	atom h-like levels resolved 20
30	30	c
31		punch results "~~h_caseb~~_lot.rlt" last
32		punch asserts "~~h_caseb~~_lot.asr" last
	31	punch results "limit_caseb_h_lot.rlt" last
	32	punch asserts "limit_caseb_h_lot.asr" last
33	33	c
34	34	assert Case B H-like element hydrogen 0.05 range 1,000 40,000A
35	35
36		c ~~h_caseb~~_lot.in
	36	c limit_caseb_h_lot.in
37	37	c class limit
38	38	c ========================================

r1994	r1996
14	14	c
15	15	c other commands for details =========
	16	c this will not affect the main model, but will test CR heating rate
	17	cosmic ray background
16	18	case b hummer no photoionization no Pdest
17	19	no scattering escape
…	…
22	24	c commands controlling output =========
23	25	normalize to "Ca B" 4861
24		punch results "~~h_caseb~~_n8.rlt"
25		punch ionizing continuum "~~h_caseb~~_n8.cion"
26		punch asserts "~~h_caseb~~_n8.asr"
	26	punch results "limit_caseb_h_n8.rlt"
	27	punch ionizing continuum "limit_caseb_h_n8.cion"
	28	punch asserts "limit_caseb_h_n8.asr"
27	29	c
28	30	assert case b h-like element hydrogen 0.10 range 1,000A to 40,000A
	31	assert csupra -18 error 0.1
	32	c
	33	c limit_caseb_h_n8.in
29	34	c class limit
30	35	c ========================================

r1994	r1996
23	23	normalize to "Ca B" 4861 = 1 % normalize to exact result
24	24	//
25		punch overview "~~h_lym_thin~~.ovr" last
26		punch fine optical depths last "~~h_lym_thin~~.fin" range 0.9 1.02 every 1
27		punch optical depths last "~~h_lym_thin~~.opt"
28		punch continuum "~~h_lym_thin~~.con" last units microns
29		punch asserts "~~h_lym_thin~~.asr" last
30		punch dr "~~h_lym_thin~~.dr" last
	25	punch overview "limit_casec_h-den2.ovr" last
	26	punch fine optical depths last "limit_casec_h-den2.fin" range 0.9 1.02 every 1
	27	punch optical depths last "limit_casec_h-den2.opt"
	28	punch continuum "limit_casec_h-den2.con" last units microns
	29	punch asserts "limit_casec_h-den2.asr" last
	30	punch dr "limit_casec_h-den2.dr" last
31	31	c
32	32	c commands giving the asserts =========
…	…
139	139	assert itrzn < 6.5
140	140	//
141		c ~~honlyoutsp~~.in
	141	c limit_casec_h-den2.in
142	142	c class limit
143	143	c ========================================

r1994	r1996
23	23	print h-like populations
24	24	print departure coefficients
25		punch asserts "~~h_casec~~.asr" last
	25	punch asserts "limit_casec_h_den5.asr" last
26	26	c
27	27	c commands giving the asserts =========
…	…
42	42	assert ionization helium 3 -0.942
43	43	//
44		c ~~h_casec~~.in
	44	c limit_casec_h_den5.in
45	45	c class limit
46	46	c ========================================

r1994	r1996
24	24	print line faint -1
25	25	print last
26		punch asserts "h_induc.asr" last
	26	punch asserts "limit_h_induc.asr" last
27	27	c
28	28	c commands giving the asserts =========
29	29	assert h-like departure coef hydrogen, average 1
30	30	assert h-like departure coef helium, average 1
31		c h_induc.in
	31	c limit_h_induc.in
32	32	c class limit
33	33	c ========================================