cloudy trunk
Loading...
Searching...
No Matches
prt_lines_lv1_k_zn.cpp
Go to the documentation of this file.
1/* This file is part of Cloudy and is copyright (C)1978-2013 by Gary J. Ferland and
2 * others. For conditions of distribution and use see copyright notice in license.txt */
3/*lines_lv1_k_zn place lines of elements potassium and heavier into lines storage stack */
4#include "cddefines.h"
5#include "cddrive.h"
6#include "coolheavy.h"
7#include "ca.h"
8#include "fe.h"
9#include "rfield.h"
10#include "dense.h"
11#include "phycon.h"
12#include "radius.h"
13#include "taulines.h"
14#include "trace.h"
15#include "lines_service.h"
16#include "rt.h"
17#include "atomfeii.h"
18#include "lines.h"
19#include "iso.h"
20
21void lines_lv1_k_zn(void)
22{
23 long int i,
24 ipnt,
25 ilo,
26 ihi;
27
28 double eff, fela;
29
30 DEBUG_ENTRY( "lines_lv1_k_zn()" );
31
32 if( trace.lgTrace )
33 {
34 fprintf( ioQQQ, " lines_lv1_k_zn called\n" );
35 }
36
37 PutLine(TauLines[ipKI7745],
38 " potassium K I 7745 ");
39
40 PutLine(TauLines[ipxK03462],
41 " [K III] 4.62 microns ");
42
43 PutLine(TauLines[ipxK04598],
44 " [KIV] 5.983 min ");
45
46 PutLine(TauLines[ipxK04154],
47 " [KIV] 15.39 micron ");
48
49 PutLine(TauLines[ipxK07319],
50 " [K VII] 3.189 microns ");
51
52 PutLine(TauLines[ipCaI4228],
53 " calcium Ca I 4228 ");
54
55 linadd(ca.Cakh,3933,"Ca 2",'i',
56 " coll excited calcium k+h " );
57
58 linadd(ca.Cair,8579,"Ca 2",'i' ,
59 " infrared triplet ");
60
61 linadd(ca.c7306,7306,"Ca 2",'i',
62 " forbidden lines, 7291+7324 together " );
63
64 linadd(ca.dCakh,3933,"Phot",'i' ,
65 " fraction H Ly-alpha destruction of excited levels ");
66
67 linadd(ca.dCaf12,7306,"Phot",'i' ,
68 " fraction H Ly-alpha destruction of excited levels ");
69
70 PntForLine(3934.,"Ca2K",&ipnt);
71 lindst(ca.Cak,3934,"Ca2K",ipnt,'t',true,
72 " individual lines from five level atom");
73
74
75 PntForLine(3969.,"Ca2H",&ipnt);
76 lindst(ca.Cah,3969,"Ca2H",ipnt,'t',true,
77 " individual lines from five level atom" );
78
79
80 PntForLine(8498.,"Ca2X",&ipnt);
81 lindst(ca.Cax,8498,"Ca2X",ipnt,'t',true,
82 " individual lines from five level atom " );
83
84
85 PntForLine(8542.,"Ca2Y",&ipnt);
86 lindst(ca.Cay,8542,"Ca2Y",ipnt,'t',true,
87 " individual lines from five level atom" );
88
89
90 PntForLine(8662.,"Ca2Z",&ipnt);
91 lindst(ca.Caz,8662,"Ca2Z",ipnt,'t',true,
92 " individual lines from five level atom" );
93
94
95 PntForLine(7291.,"CaF1",&ipnt);
96 lindst(ca.Caf1,7291,"CaF1",ipnt,'t',true,
97 " individual lines from five level atom" );
98
99
100 PntForLine(7324.,"CaF2",&ipnt);
101 lindst(ca.Caf2,7324,"CaF2",ipnt,'t',true,
102 " individual lines from five level atom" );
103
104 eff = dense.eden*dense.xIonDense[ipCALCIUM][2]*5.4e-21/(phycon.te/
105 phycon.te10/phycon.te10);
106 linadd(eff,3933,"Rec ",'i',
107 " recombination contribution to CaII emission" );
108
109 PutLine(TauLines[ipTCa3],
110 " Ca IV 3.2 micron ");
111
112
113 PntForLine(22.08e4,"Sc 2",&ipnt);
114 lindst(CoolHeavy.Sc22p08m,22.08e4,"Sc 2",ipnt,'t',true,
115 " Sc II 2.08 (1-3) " );
116
117
118 PntForLine(24.1e4,"Sc 2",&ipnt);
119 lindst(CoolHeavy.Sc24p1m,24.1e4,"Sc 2",ipnt,'t',true,
120 " Sc II 4.1 micron (1-2)" );
121
122
123 PntForLine(24.2e4,"Sc 2",&ipnt);
124 lindst(CoolHeavy.Sc24p2m,24.2e4,"Sc 2",ipnt,'t',true,
125 " Sc II 4.22 (2-3)" );
126
127
128 PntForLine(3933.,"Sc 3",&ipnt);
129 lindst(CoolHeavy.Sc33936,3933,"Sc 3",ipnt,'t',true,
130 " Sc III 3936" );
131
132 PutLine(TauLines[ipSc05231],
133 " [Sc V] 1.46 microns ");
134
135
136 PntForLine(5054.,"Sc 6",&ipnt);
137 lindst(CoolHeavy.Sc45058,5054,"Sc 6",ipnt,'t',true ,
138 " Sc VI 5054 (1-2)");
139
140
141 PntForLine(3592.,"Sc 6",&ipnt);
142 lindst(CoolHeavy.Sc43595,3592,"Sc 6",ipnt,'t',true,
143 " Sc VI 3595 (2-3)" );
144
145
146 PntForLine(2100.,"Sc 6",&ipnt);
147 lindst(CoolHeavy.Sc42100,2100,"Sc 6",ipnt,'t',true,
148 " Sc VI 2100 (1-3)" );
149
150 PutLine(TauLines[ipSc13264],
151 " [Sc 13] 2637.97A");
152
153 PntForLine(8823.,"V 3",&ipnt);
154 lindst(CoolHeavy.V38830,8823,"V 3",ipnt,'t',true ,
155 " V III 8823 ");
156
157
158 PntForLine(8507.,"V 3",&ipnt);
159 lindst(CoolHeavy.V38507,8507,"V 3",ipnt,'t',true,
160 " V III 8507" );
161
162 PntForLine(5828.,"Cr 3",&ipnt);
163 lindst(CoolHeavy.Cr3l21,5828,"Cr 3",ipnt,'t',true,
164 " [CrIII] multiplet blend at 5828A" );
165
166 PntForLine(7267.,"Cr 4",&ipnt);
167 lindst(CoolHeavy.Cr4l21,7267,"Cr 4",ipnt,'t',true,
168 " [CrIV] 2 - 1 multiplet blend at 7272" );
169
170
171 PntForLine(6801.,"Cr 4",&ipnt);
172 lindst(CoolHeavy.Cr4l31,6801,"Cr 4",ipnt,'t',true,
173 " [CrIV] 3 - 1 multiplet blend at 6806" );
174
175
176 PntForLine(7979.,"Cr 5",&ipnt);
177 lindst(CoolHeavy.Cr5l21,7979,"Cr 5",ipnt,'t',true,
178 " [CrV] 2 - 1 multiplet blend at 7985" );
179
180 PntForLine(6577.,"Cr 5",&ipnt);
181 lindst(CoolHeavy.Cr5l31,6577,"Cr 5",ipnt,'t',true,
182 " [CrV] 3 - 1 multiplet blend at 6582" );
183
184
185 PntForLine(3.75e4,"Cr 5",&ipnt);
186 lindst(CoolHeavy.Cr5l32,3.75e4,"Cr 5",ipnt,'t',true,
187 " [CrV] 3 - 2 multiplet blend at 3.75 microns " );
188
189 /* bob Rubin's UV line
190 * f2 = dense.xIonDense(26,4)*sexp(50 764./te)*0.45*cdsqte/6.*7.01e-12
191 * call linadd( f2 , 2837 , 'BobR' , 'i')
192 * f2 = dense.xIonDense(26,4)*sexp(55 989./te)*0.384*cdsqte/6.*7.74e-12
193 * call linadd( f2 , 2568 , 'BobR' , 'i') */
194
195 /* iron */
196
197 PutLine(TauLines[ipFe1_24m],
198 " Fe 1 24m ");
199
200 PutLine(TauLines[ipFe1_35m],
201 " Fe 1 35m ");
202
203 PutLine(TauLines[ipFe1_54m],
204 " Fe 1 54m ");
205
206 PutLine(TauLines[ipFe1_111m],
207 " Fe 1 111m ");
208
209 PutLine(TauLines[ipFeI3884],
210 " Fe 1 3884 ");
211
212 PutLine(TauLines[ipFeI3729],
213 " Fe 1 3729 ");
214
215 PutLine(TauLines[ipFeI3457],
216 " Fe 1 3457 ");
217
218 PutLine(TauLines[ipFeI3021],
219 " Fe 1 3021 ");
220
221 PutLine(TauLines[ipFeI2966],
222 " Fe 1 2966 ");
223
224 linadd(MAX2(0.,FeII.Fe2_large_cool+FeII.Fe2_UVsimp_cool),0,"Fe2c",'c' ,
225 "total of all Fe 2 cooling, both simple UV and large atom together ");
226
227 linadd(MAX2(0.,-FeII.Fe2_large_cool-FeII.Fe2_UVsimp_cool),0,"Fe2h",'h' ,
228 "total of all Fe 2 heating, both simple UV and large atom together ");
229
230 linadd(FeII.for7,4300,"Fe 2",'i' ,
231 " Fe 2 forbidden 2-1 transition from Netzer's atom ");
232
233 linadd(TauLines[ipTuv3].Emis().xIntensity(), TauLines[ipTuv3].WLAng(), "Fe 2", 'i',
234 " 2400 in simple Wills, Netzer, Wills FeII");
235 linadd(TauLines[ipTr48].Emis().xIntensity(), TauLines[ipTr48].WLAng(), "Fe 2", 'i',
236 " 6200 in simple Wills, Netzer, Wills FeII");
237 linadd(TauLines[ipTFe16].Emis().xIntensity(), TauLines[ipTFe16].WLAng(), "Fe 2", 'i',
238 " 1080 in simple Wills, Netzer, Wills FeII");
239 linadd(TauLines[ipTFe26].Emis().xIntensity(), TauLines[ipTFe26].WLAng(), "Fe 2", 'i',
240 " 1500 in simple Wills, Netzer, Wills FeII");
241 linadd(TauLines[ipTFe34].Emis().xIntensity(), TauLines[ipTFe34].WLAng(), "Fe 2", 'i',
242 " 11500 in simple Wills, Netzer, Wills FeII");
243 linadd(TauLines[ipTFe35].Emis().xIntensity(), TauLines[ipTFe35].WLAng(), "Fe 2", 'i',
244 " 2500 in simple Wills, Netzer, Wills FeII");
245 linadd(TauLines[ipTFe46].Emis().xIntensity(), TauLines[ipTFe46].WLAng(), "Fe 2", 'i',
246 " 2300 in simple Wills, Netzer, Wills FeII");
247 linadd(TauLines[ipTFe56].Emis().xIntensity(), TauLines[ipTFe56].WLAng(), "Fe 2", 'i',
248 " 8900 in simple Wills, Netzer, Wills FeII");
249
250 /* option to save all intensities predicted by large FeII atom,
251 * code is in FeIILevelPops */
252 FeIIAddLines();
253 /* we were called by lines, and we want to zero out Fe2SavN */
254 for( long ipLo=0; ipLo < (FeII.nFeIILevel_malloc - 1); ipLo++ )
255 {
256 for( long ipHi=ipLo + 1; ipHi < FeII.nFeIILevel_malloc; ipHi++ )
257 {
258 /* only evaluate real transitions */
259 TransitionProxy tr = Fe2LevN[ipFe2LevN[ipHi][ipLo]];
260 if( tr.ipCont() > 0 )
261 PutLine( tr ," Fe II emission" );
262 }
263 }
264
265 /* emission bands from the model Fe II atom */
266 for( i=0; i < nFeIIBands; i++ )
267 {
268 double SumBandInward;
269 /* [i][0] is center wavelength, [i][1] and [i][2] are upper and
270 * lower bounds in Angstroms. These are set in FeIIZero
271 * units are erg s-1 cm-3 */
272 eff = FeIISumBand(FeII_Bands[i][1],FeII_Bands[i][2],
273 &SumBandInward);
274
275 linadd(eff,FeII_Bands[i][0],"Fe2b",'i' ,
276 " total Fe II emission in Fe II bands, as defined in FeII_bands.ini ");
277 linadd(SumBandInward,FeII_Bands[i][0],"Inwd",'i' ,
278 " inward Fe II emission in Fe II bands, as defined in FeII_bands.ini ");
279 }
280
281 // integrate the pseudo continuum of FeII emission
282 if( LineSave.ipass > 0 )
283 {
284 // initialize
285 if( nzone == 1 )
286 {
287 for( i=0; i < nFeIIConBins; i++ )
288 {
289 /* initialize arrays */
290 FeII_Cont[i][1] = 0.;
291 FeII_Cont[i][2] = 0.;
292 }
293 }
294
295 // integrate
296 for( i=0; i < nFeIIConBins; i++ )
297 {
298 double SumBandInward;
299 /* [i][0] is total intensity in cell, [i][1] and [i][2] are lower and
300 * upper bounds in Angstroms. these are set in FeIIZero *
301 * find total emission from large FeII atom, integrated over band */
302 double TotalFeII = FeIISumBand(FeII_Cont[i][0],FeII_Cont[i+1][0],
303 &SumBandInward);
304 FeII_Cont[i][1] += (realnum)(SumBandInward*radius.dVeffAper);
305 FeII_Cont[i][2] += (realnum)(MAX2(0.,TotalFeII-SumBandInward)*radius.dVeffAper);
306 /*fprintf(ioQQQ,"DEBUG feii\t%li\t%.2e\n", i, FeII_Cont[i][0]);*/
307 }
308 }
309
310 PutLine(TauLines[ipT191],
311 " anomalous Fe 2 transition at 1787, RMT 191");
312
313 linadd(fe.Fe3CoolTot,0,"Fe3c",'c' ,
314 " chng 05 dec 16, FeIII code created by Kevin Blagrave Fe3c 0 - total cooling due to 14-level Fe 3 atom ");
315
316 /* Fe 3 14-level atom
317 * following from print statements within loop */
318 /* Fe 3 22.92m from Blagrave 14-level atom */
319 /* Fe 3 13.53m from Blagrave 14-level atom */
320 /* Fe 3 33.03m from Blagrave 14-level atom */
321 /* Fe 3 10.72m from Blagrave 14-level atom */
322 /* Fe 3 20.15m from Blagrave 14-level atom */
323 /* Fe 3 51.67m from Blagrave 14-level atom */
324 /* Fe 3 9.732m from Blagrave 14-level atom */
325 /* Fe 3 16.91m from Blagrave 14-level atom */
326 /* Fe 3 34.66m from Blagrave 14-level atom */
327 /* Fe 3 105.3m from Blagrave 14-level atom */
328 /* Fe 3 5152A from Blagrave 14-level atom */
329 /* Fe 3 5271A from Blagrave 14-level atom */
330 /* Fe 3 5356A from Blagrave 14-level atom */
331 /* Fe 3 5412A from Blagrave 14-level atom */
332 /* Fe 3 5440A from Blagrave 14-level atom */
333 /* Fe 3 4986A from Blagrave 14-level atom */
334 /* Fe 3 5097A from Blagrave 14-level atom */
335 /* Fe 3 5177A from Blagrave 14-level atom */
336 /* Fe 3 5230A from Blagrave 14-level atom */
337 /* Fe 3 5256A from Blagrave 14-level atom */
338 /* Fe 3 15.47m from Blagrave 14-level atom */
339 /* Fe 3 4925A from Blagrave 14-level atom */
340 /* Fe 3 5033A from Blagrave 14-level atom */
341 /* Fe 3 5111A from Blagrave 14-level atom */
342 /* Fe 3 5162A from Blagrave 14-level atom */
343 /* Fe 3 5188A from Blagrave 14-level atom */
344 /* Fe 3 11.16m from Blagrave 14-level atom */
345 /* Fe 3 40.04m from Blagrave 14-level atom */
346 /* Fe 3 4881A from Blagrave 14-level atom */
347 /* Fe 3 4988A from Blagrave 14-level atom */
348 /* Fe 3 5064A from Blagrave 14-level atom */
349 /* Fe 3 5114A from Blagrave 14-level atom */
350 /* Fe 3 5139A from Blagrave 14-level atom */
351 /* Fe 3 9.282m from Blagrave 14-level atom */
352 /* Fe 3 23.21m from Blagrave 14-level atom */
353 /* Fe 3 55.20m from Blagrave 14-level atom */
354 /* Fe 3 4833A from Blagrave 14-level atom */
355 /* Fe 3 4937A from Blagrave 14-level atom */
356 /* Fe 3 5012A from Blagrave 14-level atom */
357 /* Fe 3 5061A from Blagrave 14-level atom */
358 /* Fe 3 5085A from Blagrave 14-level atom */
359 /* Fe 3 7.789m from Blagrave 14-level atom */
360 /* Fe 3 15.69m from Blagrave 14-level atom */
361 /* Fe 3 25.79m from Blagrave 14-level atom */
362 /* Fe 3 48.41m from Blagrave 14-level atom */
363 /* Fe 3 4714A from Blagrave 14-level atom */
364 /* Fe 3 4813A from Blagrave 14-level atom */
365 /* Fe 3 4884A from Blagrave 14-level atom */
366 /* Fe 3 4931A from Blagrave 14-level atom */
367 /* Fe 3 4954A from Blagrave 14-level atom */
368 /* Fe 3 5.543m from Blagrave 14-level atom */
369 /* Fe 3 8.638m from Blagrave 14-level atom */
370 /* Fe 3 11.01m from Blagrave 14-level atom */
371 /* Fe 3 13.76m from Blagrave 14-level atom */
372 /* Fe 3 19.22m from Blagrave 14-level atom */
373 /* Fe 3 4659A from Blagrave 14-level atom */
374 /* Fe 3 4755A from Blagrave 14-level atom */
375 /* Fe 3 4825A from Blagrave 14-level atom */
376 /* Fe 3 4870A from Blagrave 14-level atom */
377 /* Fe 3 4893A from Blagrave 14-level atom */
378 /* Fe 3 4.859m from Blagrave 14-level atom */
379 /* Fe 3 7.085m from Blagrave 14-level atom */
380 /* Fe 3 8.608m from Blagrave 14-level atom */
381 /* Fe 3 10.20m from Blagrave 14-level atom */
382 /* Fe 3 12.92m from Blagrave 14-level atom */
383 /* Fe 3 39.41m from Blagrave 14-level atom */
384 /* Fe 3 4608A from Blagrave 14-level atom */
385 /* Fe 3 4702A from Blagrave 14-level atom */
386 /* Fe 3 4770A from Blagrave 14-level atom */
387 /* Fe 3 4814A from Blagrave 14-level atom */
388 /* Fe 3 4836A from Blagrave 14-level atom */
389 /* Fe 3 4.356m from Blagrave 14-level atom */
390 /* Fe 3 6.063m from Blagrave 14-level atom */
391 /* Fe 3 7.146m from Blagrave 14-level atom */
392 /* Fe 3 8.208m from Blagrave 14-level atom */
393 /* Fe 3 9.884m from Blagrave 14-level atom */
394 /* Fe 3 20.34m from Blagrave 14-level atom */
395 /* Fe 3 42.06m from Blagrave 14-level atom */
396 /* Fe 3 4574A from Blagrave 14-level atom */
397 /* Fe 3 4668A from Blagrave 14-level atom */
398 /* Fe 3 4734A from Blagrave 14-level atom */
399 /* Fe 3 4778A from Blagrave 14-level atom */
400 /* Fe 3 4800A from Blagrave 14-level atom */
401 /* Fe 3 4.077m from Blagrave 14-level atom */
402 /* Fe 3 5.535m from Blagrave 14-level atom */
403 /* Fe 3 6.423m from Blagrave 14-level atom */
404 /* Fe 3 7.269m from Blagrave 14-level atom */
405 /* Fe 3 8.554m from Blagrave 14-level atom */
406 /* Fe 3 15.41m from Blagrave 14-level atom */
407 /* Fe 3 25.31m from Blagrave 14-level atom */
408 /* Fe 3 63.56m from Blagrave 14-level atom */
409 for( ihi=1; ihi<NLFE3; ++ihi )
410 {
411 for( ilo=0; ilo<ihi; ++ilo )
412 {
413 /* emission in these lines */
414 PntForLine(fe.Fe3_wl[ihi][ilo],"Fe 3",&ipnt);
415# if 0
416 fprintf( ioQQQ,"\t/* FeIII ");
417 prt_wl( ioQQQ , (realnum)(fe.Fe3_wl[ihi][ilo]+0.5) );
418 fprintf( ioQQQ," from Blagrave 14-level atom */\n" );
419# endif
420 lindst( fe.Fe3_emiss[ihi][ilo] , (realnum)(fe.Fe3_wl[ihi][ilo]+0.5) , "Fe 3",ipnt,'t',true,
421 " " );
422 }
423 }
424
425 /*>>chng 05 dec 18, following are now in the above */
426 /* sum of 3p and 3g states together */
427 /* linadd(CoolHeavy.c5270,0,"Fe 3",'c' ); */
428
429 /* Fe 3 5270, predictions from Garstang et al 78
430 PntForLine(5270.,"Fe 3",&ipnt);
431 lindst(CoolHeavy.c5270*0.2090,5270,"Fe 3",ipnt,'c',true );*/
432
433 /* Fe 3 5270, predictions from Garstang et al 78
434 PntForLine(4658.,"Fe 3",&ipnt);
435 lindst(CoolHeavy.c5270*0.3667,4658,"Fe 3",ipnt,'c',true ); */
436
437 PutLine(TauLines[ipT1122]," Fe 3 1122 entire multiplet");
438
439 linadd(fe.Fe4CoolTot,0,"Fe4c",'i',
440 " Fe4c 0 - total cooling due to 12-level Fe 4 atom " );
441
442
443 PntForLine(3096.,"Fe 4",&ipnt);
444 lindst(fe.fe40401,3096,"Fe 4",ipnt,'t',true,
445 " Fe 4 3096.A, 4-1 and 5-1 transitions together" );
446
447
448 PntForLine(2836.,"Fe 4",&ipnt);
449 lindst(fe.fe42836,2836,"Fe 4",ipnt,'t',true,
450 " Fe 4 2835.7A, 6-1 transition, 4P5/2 - 6S5/2 " );
451
452
453 PntForLine(2829.,"Fe 4",&ipnt);
454 lindst(fe.fe42829,2829,"Fe 4",ipnt,'t',true,
455 " Fe 4 2829.4A, 7-1 transition, 4P3/2 - 6S5/2" );
456
457
458 PntForLine(2567.,"Fe 4",&ipnt);
459 lindst(fe.fe42567,2567,"Fe 4",ipnt,'t',true,
460 " Fe 4 2567.6+ 2567.4. 11-1 and 12-1 transitions" );
461
462
463 PntForLine(2.774e4,"Fe 4",&ipnt);
464 lindst(fe.fe41207,2.774e4,"Fe 4",ipnt,'t',true,
465 " Fe 4 2.774 microns 12-7 transition " );
466
467
468 PntForLine(2.714e4,"Fe 4",&ipnt);
469 lindst(fe.fe41206,2.714e4,"Fe 4",ipnt,'t',true,
470 " Fe 4 2.714 microns 12-6 transition " );
471
472
473 PntForLine(2.716e4,"Fe 4",&ipnt);
474 lindst(fe.fe41106,2.716e4,"Fe 4",ipnt,'t',true,
475 " Fe 4 2.716 microns 11-6 transition" );
476
477
478 PntForLine(2.806e4,"Fe 4",&ipnt);
479 lindst(fe.fe41007,2.806e4,"Fe 4",ipnt,'t',true,
480 " Fe 4 2.806 microns 10-7 transition " );
481
482
483 PntForLine(2.865e4,"Fe 4",&ipnt);
484 lindst(fe.fe41008,2.865e4,"Fe 4",ipnt,'t',true ,
485 " Fe 4 2.865 microns 10-8 transition");
486
487
488 PntForLine(2.836e4,"Fe 4",&ipnt);
489 lindst(fe.fe40906,2.836e4,"Fe 4",ipnt,'t',true,
490 " Fe 4 2.836 microns 9-6 transition" );
491
492
493 PntForLine(3892.,"Fe 5",&ipnt);
494 lindst(CoolHeavy.c3892,3892,"Fe 5",ipnt,'t',true,
495 " Fe 5 3892+3839" );
496 /* recombination Ka */
497 if( dense.lgElmtOn[ipIRON] )
498 {
499 /* these lines added to outlin in metdif - following must be false
500 * fela = xLyaHeavy(nelem,nelem)*dense.xIonDense(nelem,nelem+1) */
501 fela = iso_sp[ipH_LIKE][ipIRON].trans(ipH2p,ipH1s).Emis().xIntensity();
502 }
503 else
504 {
505 fela = 0.;
506 }
507
508 /* >>chng 02 jan 14, add grain fe to this sum */
509 /* total intensity of K-alpha line */
510 /*linadd((fe.fekcld+fe.fegrain)*1.03e-8+(fe.fekhot+fela)*1.11e-8,2,"FeKa",'i' );*/
511 if( dense.lgElmtOn[ipIRON] )
512 {
513 lindst((fe.fekcld+fe.fegrain)*1.03e-8+(fe.fekhot+fela)*1.11e-8,1.78f,"FeKa",
514 iso_sp[ipH_LIKE][ipIRON].trans(ipH2p,ipH1s).ipCont(),'t',false,
515 "total intensity of K-alpha line" );
516 }
517
518 linadd(fela*1.11e-8,2,"FeLr",'i' ,
519 " recombination from fully stripped ion ");
520
521 /* >>chng 03 aug 14, label changed from TotH to AugH to be like rest total hot iron Ka; */
522 linadd((fe.fekhot+fela)*1.11e-8,2,"AugH",'i' ,
523 " Auger hot iron, assumes case b for H and He-like ");
524
525 linadd(fe.fekcld*1.03e-8,2,"AugC",'i',
526 " Auger production of cold iron, less than or 17 times ionized " );
527
528 linadd(fe.fegrain*1.03e-8,2,"AugG",'i' ,
529 " grain production of cold iron ");
530
531 PutLine(TauLines[ipCo11527],
532 " [Co XI] 5168. A ");
533
534 PutLine(TauLines[ipNi1_7m],
535 " nickel [Ni I] 7m ");
536
537 /* nickel*/
538
539
540 PutLine(TauLines[ipNi1_11m],
541 " [Ni I] 11m ");
542
543 /* copper */
544
545 if( trace.lgTrace )
546 {
547 fprintf( ioQQQ, " lines_lv1_k_zn returns\n" );
548 }
549 return;
550}
ipCo11527
long ipCo11527
Definition: atmdat_readin.cpp:89
ipTFe34
long ipTFe34
Definition: atmdat_readin.cpp:87
ipxK03462
long ipxK03462
Definition: atmdat_readin.cpp:80
ipNi1_7m
long ipNi1_7m
Definition: atmdat_readin.cpp:97
ipTFe26
long ipTFe26
Definition: atmdat_readin.cpp:87
ipTCa3
long ipTCa3
Definition: atmdat_readin.cpp:83
ipFe1_24m
long ipFe1_24m
Definition: atmdat_readin.cpp:96
ipFe1_35m
long ipFe1_35m
Definition: atmdat_readin.cpp:96
ipNi1_11m
long ipNi1_11m
Definition: atmdat_readin.cpp:97
ipCaI4228
long ipCaI4228
Definition: atmdat_readin.cpp:81
ipT1122
long ipT1122
Definition: atmdat_readin.cpp:88
ipT191
long ipT191
Definition: atmdat_readin.cpp:89
ipxK04598
long ipxK04598
Definition: atmdat_readin.cpp:80
ipKI7745
long ipKI7745
Definition: atmdat_readin.cpp:80
ipSc05231
long ipSc05231
Definition: atmdat_readin.cpp:84
ipFe1_111m
long ipFe1_111m
Definition: atmdat_readin.cpp:96
ipTuv3
long ipTuv3
Definition: atmdat_readin.cpp:86
ipFe1_54m
long ipFe1_54m
Definition: atmdat_readin.cpp:96
ipSc13264
long ipSc13264
Definition: atmdat_readin.cpp:84
ipTFe46
long ipTFe46
Definition: atmdat_readin.cpp:88
ipTFe56
long ipTFe56
Definition: atmdat_readin.cpp:88
ipFeI2966
long ipFeI2966
Definition: atmdat_readin.cpp:86
ipxK07319
long ipxK07319
Definition: atmdat_readin.cpp:81
ipFeI3021
long ipFeI3021
Definition: atmdat_readin.cpp:86
ipFeI3729
long ipFeI3729
Definition: atmdat_readin.cpp:85
ipxK04154
long ipxK04154
Definition: atmdat_readin.cpp:81
ipTFe35
long ipTFe35
Definition: atmdat_readin.cpp:88
ipFeI3884
long ipFeI3884
Definition: atmdat_readin.cpp:85
ipFeI3457
long ipFeI3457
Definition: atmdat_readin.cpp:85
ipTr48
long ipTr48
Definition: atmdat_readin.cpp:87
ipTFe16
long ipTFe16
Definition: atmdat_readin.cpp:87
nFeIIBands
long int nFeIIBands
Definition: atom_feii.cpp:116
FeII_Bands
realnum ** FeII_Bands
Definition: atom_feii.cpp:106
nFeIIConBins
long int nFeIIConBins
Definition: atom_feii.cpp:119
FeIIAddLines
void FeIIAddLines(void)
Definition: atom_feii.cpp:1574
FeIISumBand
double FeIISumBand(realnum wl1, realnum wl2, double *SumBandInward)
Definition: atom_feii.cpp:1355
FeII_Cont
realnum ** FeII_Cont
Definition: atom_feii.cpp:113
FeII
t_FeII FeII
Definition: atomfeii.cpp:5
ca
t_ca ca
Definition: ca.cpp:5
nzone
long int nzone
Definition: cddefines.cpp:14
ioQQQ
FILE * ioQQQ
Definition: cddefines.cpp:7
ipIRON
const int ipIRON
Definition: cddefines.h:330
realnum
float realnum
Definition: cddefines.h:103
MAX2
#define MAX2
Definition: cddefines.h:782
ipCALCIUM
const int ipCALCIUM
Definition: cddefines.h:324
DEBUG_ENTRY
#define DEBUG_ENTRY(funcname)
Definition: cddefines.h:684
EmissionProxy::xIntensity
double & xIntensity() const
Definition: emission.h:483
TransitionProxy::ipCont
long & ipCont() const
Definition: transition.h:450
TransitionProxy::Emis
EmissionList::reference Emis() const
Definition: transition.h:408
t_iso_sp::trans
TransitionProxy trans(const long ipHi, const long ipLo)
Definition: iso.h:444
CoolHeavy
t_CoolHeavy CoolHeavy
Definition: coolheavy.cpp:5
dense
t_dense dense
Definition: dense.cpp:24
fe
t_fe fe
Definition: fe.cpp:5
NLFE3
#define NLFE3
Definition: fe.h:10
iso_sp
t_iso_sp iso_sp[NISO][LIMELM]
Definition: iso.cpp:8
ipH1s
const int ipH1s
Definition: iso.h:27
ipH2p
const int ipH2p
Definition: iso.h:29
ipH_LIKE
const int ipH_LIKE
Definition: iso.h:62
LineSave
t_LineSave LineSave
Definition: lines.cpp:5
linadd
void linadd(double xInten, realnum wavelength, const char *chLab, char chInfo, const char *chComment)
Definition: lines_service.cpp:316
PntForLine
void PntForLine(double wavelength, const char *chLabel, long int *ipnt)
Definition: lines_service.cpp:583
lindst
void lindst(double xInten, realnum wavelength, const char *chLab, long int ipnt, char chInfo, bool lgOutToo, const char *chComment)
Definition: lines_service.cpp:468
lines_service.h
phycon
t_phycon phycon
Definition: phycon.cpp:6
prt_wl
void prt_wl(FILE *ioOUT, realnum wl)
Definition: prt.cpp:13
lines_lv1_k_zn
void lines_lv1_k_zn(void)
Definition: prt_lines_lv1_k_zn.cpp:21
radius
t_radius radius
Definition: radius.cpp:5
t_CoolHeavy::Cr5l32
double Cr5l32
Definition: coolheavy.h:47
t_CoolHeavy::Cr3l21
double Cr3l21
Definition: coolheavy.h:48
t_CoolHeavy::Sc42100
double Sc42100
Definition: coolheavy.h:93
t_CoolHeavy::Cr4l21
double Cr4l21
Definition: coolheavy.h:44
t_CoolHeavy::Cr5l21
double Cr5l21
Definition: coolheavy.h:46
t_CoolHeavy::Cr4l31
double Cr4l31
Definition: coolheavy.h:43
t_CoolHeavy::Sc24p1m
double Sc24p1m
Definition: coolheavy.h:90
t_CoolHeavy::Sc45058
double Sc45058
Definition: coolheavy.h:94
t_CoolHeavy::Cr5l31
double Cr5l31
Definition: coolheavy.h:45
t_CoolHeavy::Sc22p08m
double Sc22p08m
Definition: coolheavy.h:89
t_CoolHeavy::V38507
double V38507
Definition: coolheavy.h:97
t_CoolHeavy::c3892
double c3892
Definition: coolheavy.h:49
t_CoolHeavy::Sc33936
double Sc33936
Definition: coolheavy.h:92
t_CoolHeavy::Sc24p2m
double Sc24p2m
Definition: coolheavy.h:91
t_CoolHeavy::V38830
double V38830
Definition: coolheavy.h:96
t_CoolHeavy::Sc43595
double Sc43595
Definition: coolheavy.h:95
t_FeII::for7
double for7
Definition: atomfeii.h:259
t_FeII::nFeIILevel_malloc
long int nFeIILevel_malloc
Definition: atomfeii.h:193
t_FeII::Fe2_large_cool
double Fe2_large_cool
Definition: atomfeii.h:250
t_FeII::Fe2_UVsimp_cool
double Fe2_UVsimp_cool
Definition: atomfeii.h:258
t_LineSave::ipass
long int ipass
Definition: lines.h:75
t_ca::Cax
realnum Cax
Definition: ca.h:16
t_ca::Cair
double Cair
Definition: ca.h:22
t_ca::dCaf12
realnum dCaf12
Definition: ca.h:29
t_ca::Cakh
double Cakh
Definition: ca.h:24
t_ca::Cak
realnum Cak
Definition: ca.h:14
t_ca::Caf2
realnum Caf2
Definition: ca.h:20
t_ca::c7306
double c7306
Definition: ca.h:23
t_ca::Cay
realnum Cay
Definition: ca.h:17
t_ca::Cah
realnum Cah
Definition: ca.h:15
t_ca::dCakh
realnum dCakh
Definition: ca.h:28
t_ca::Caf1
realnum Caf1
Definition: ca.h:19
t_ca::Caz
realnum Caz
Definition: ca.h:18
t_dense::eden
double eden
Definition: dense.h:190
t_dense::lgElmtOn
bool lgElmtOn[LIMELM]
Definition: dense.h:146
t_dense::xIonDense
double xIonDense[LIMELM][LIMELM+1]
Definition: dense.h:125
t_fe::fe40906
double fe40906
Definition: fe.h:27
t_fe::Fe3CoolTot
double Fe3CoolTot
Definition: fe.h:30
t_fe::fe41207
double fe41207
Definition: fe.h:22
t_fe::fe41008
double fe41008
Definition: fe.h:26
t_fe::fe41206
double fe41206
Definition: fe.h:23
t_fe::Fe4CoolTot
double Fe4CoolTot
Definition: fe.h:17
t_fe::Fe3_emiss
double ** Fe3_emiss
Definition: fe.h:31
t_fe::Fe3_wl
double ** Fe3_wl
Definition: fe.h:31
t_fe::fe42567
double fe42567
Definition: fe.h:20
t_fe::fe40401
double fe40401
Definition: fe.h:21
t_fe::fe41106
double fe41106
Definition: fe.h:24
t_fe::fegrain
realnum fegrain
Definition: fe.h:38
t_fe::fe42829
double fe42829
Definition: fe.h:19
t_fe::fekcld
realnum fekcld
Definition: fe.h:35
t_fe::fe42836
double fe42836
Definition: fe.h:18
t_fe::fekhot
realnum fekhot
Definition: fe.h:34
t_fe::fe41007
double fe41007
Definition: fe.h:25
t_phycon::te
double te
Definition: phycon.h:11
t_phycon::te10
double te10
Definition: phycon.h:55
t_radius::dVeffAper
double dVeffAper
Definition: radius.h:87
t_trace::lgTrace
bool lgTrace
Definition: trace.h:12
Fe2LevN
TransitionList Fe2LevN("Fe2LevN", &Fe2LevNStates)
ipFe2LevN
multi_arr< int, 2 > ipFe2LevN
Definition: taulines.cpp:34
TauLines
TransitionList TauLines("TauLines", &AnonStates)
trace
t_trace trace
Definition: trace.cpp:5
PutLine
void PutLine(const TransitionProxy &t, const char *chComment, const char *chLabelTemp)
Definition: transition.cpp:449
Generated by doxygen 1.9.5