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
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
38 " potassium K I 7745 ");
39
41 " [K III] 4.62 microns ");
42
44 " [KIV] 5.983 min ");
45
47 " [KIV] 15.39 micron ");
48
50 " [K VII] 3.189 microns ");
51
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
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
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
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
198 " Fe 1 24m ");
199
201 " Fe 1 35m ");
202
204 " Fe 1 54m ");
205
207 " Fe 1 111m ");
208
210 " Fe 1 3884 ");
211
213 " Fe 1 3729 ");
214
216 " Fe 1 3457 ");
217
219 " Fe 1 3021 ");
220
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
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
532 " [Co XI] 5168. A ");
533
535 " nickel [Ni I] 7m ");
536
537 /* nickel*/
538
539
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}
long ipCo11527
long ipTFe34
long ipxK03462
long ipNi1_7m
long ipTFe26
long ipTCa3
long ipFe1_24m
long ipFe1_35m
long ipNi1_11m
long ipCaI4228
long ipT1122
long ipT191
long ipxK04598
long ipKI7745
long ipSc05231
long ipFe1_111m
long ipTuv3
long ipFe1_54m
long ipSc13264
long ipTFe46
long ipTFe56
long ipFeI2966
long ipxK07319
long ipFeI3021
long ipFeI3729
long ipxK04154
long ipTFe35
long ipFeI3884
long ipFeI3457
long ipTr48
long ipTFe16
long int nFeIIBands
realnum ** FeII_Bands
long int nFeIIConBins
void FeIIAddLines(void)
double FeIISumBand(realnum wl1, realnum wl2, double *SumBandInward)
realnum ** FeII_Cont
t_FeII FeII
Definition atomfeii.cpp:5
t_ca ca
Definition ca.cpp:5
long int nzone
Definition cddefines.cpp:14
FILE * ioQQQ
Definition cddefines.cpp:7
const int ipIRON
Definition cddefines.h:330
float realnum
Definition cddefines.h:103
#define MAX2
Definition cddefines.h:782
const int ipCALCIUM
Definition cddefines.h:324
#define DEBUG_ENTRY(funcname)
Definition cddefines.h:684
long & ipCont() const
Definition transition.h:450
t_CoolHeavy CoolHeavy
Definition coolheavy.cpp:5
t_dense dense
Definition dense.cpp:24
t_fe fe
Definition fe.cpp:5
#define NLFE3
Definition fe.h:10
t_iso_sp iso_sp[NISO][LIMELM]
Definition iso.cpp:8
const int ipH1s
Definition iso.h:27
const int ipH2p
Definition iso.h:29
const int ipH_LIKE
Definition iso.h:62
t_LineSave LineSave
Definition lines.cpp:5
void linadd(double xInten, realnum wavelength, const char *chLab, char chInfo, const char *chComment)
void PntForLine(double wavelength, const char *chLabel, long int *ipnt)
void lindst(double xInten, realnum wavelength, const char *chLab, long int ipnt, char chInfo, bool lgOutToo, const char *chComment)
t_phycon phycon
Definition phycon.cpp:6
void prt_wl(FILE *ioOUT, realnum wl)
Definition prt.cpp:13
void lines_lv1_k_zn(void)
t_radius radius
Definition radius.cpp:5
TransitionList Fe2LevN("Fe2LevN", &Fe2LevNStates)
multi_arr< int, 2 > ipFe2LevN
Definition taulines.cpp:34
TransitionList TauLines("TauLines", &AnonStates)
t_trace trace
Definition trace.cpp:5
void PutLine(const TransitionProxy &t, const char *chComment, const char *chLabelTemp)