cloudy  trunk
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 
21 void 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 }
t_ca::c7306
double c7306
Definition: ca.h:23
t_ca::Caf1
realnum Caf1
Definition: ca.h:19
ipTFe16
long ipTFe16
Definition: atmdat_readin.cpp:87
lines.h
ipNi1_11m
long ipNi1_11m
Definition: atmdat_readin.cpp:97
PntForLine
void PntForLine(double wavelength, const char *chLabel, long int *ipnt)
Definition: lines_service.cpp:583
t_dense::eden
double eden
Definition: dense.h:190
ipTFe34
long ipTFe34
Definition: atmdat_readin.cpp:87
ipTr48
long ipTr48
Definition: atmdat_readin.cpp:87
ipTuv3
long ipTuv3
Definition: atmdat_readin.cpp:86
dense
t_dense dense
Definition: dense.cpp:24
ipTCa3
long ipTCa3
Definition: atmdat_readin.cpp:83
ipSc13264
long ipSc13264
Definition: atmdat_readin.cpp:84
t_CoolHeavy::Cr4l21
double Cr4l21
Definition: coolheavy.h:44
ioQQQ
FILE * ioQQQ
Definition: cddefines.cpp:7
EmissionProxy::xIntensity
double & xIntensity() const
Definition: emission.h:483
ipTFe46
long ipTFe46
Definition: atmdat_readin.cpp:88
t_ca::Cax
realnum Cax
Definition: ca.h:16
realnum
float realnum
Definition: cddefines.h:103
t_radius::dVeffAper
double dVeffAper
Definition: radius.h:87
rfield.h
t_ca::dCakh
realnum dCakh
Definition: ca.h:28
t_fe::fe42567
double fe42567
Definition: fe.h:20
t_fe::fekcld
realnum fekcld
Definition: fe.h:35
t_CoolHeavy::Sc22p08m
double Sc22p08m
Definition: coolheavy.h:89
t_dense::lgElmtOn
bool lgElmtOn[LIMELM]
Definition: dense.h:146
ipT1122
long ipT1122
Definition: atmdat_readin.cpp:88
ipFeI3884
long ipFeI3884
Definition: atmdat_readin.cpp:85
ipIRON
const int ipIRON
Definition: cddefines.h:330
phycon
t_phycon phycon
Definition: phycon.cpp:6
trace.h
t_fe::fekhot
realnum fekhot
Definition: fe.h:34
t_ca::Cak
realnum Cak
Definition: ca.h:14
ca
t_ca ca
Definition: ca.cpp:5
lines_service.h
TransitionProxy::ipCont
long & ipCont() const
Definition: transition.h:450
ipNi1_7m
long ipNi1_7m
Definition: atmdat_readin.cpp:97
ipFeI3729
long ipFeI3729
Definition: atmdat_readin.cpp:85
t_ca::Cakh
double Cakh
Definition: ca.h:24
iso.h
t_ca::Cah
realnum Cah
Definition: ca.h:15
fe
t_fe fe
Definition: fe.cpp:5
t_FeII::for7
double for7
Definition: atomfeii.h:259
TransitionProxy
Definition: transition.h:23
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
cddrive.h
nzone
long int nzone
Definition: cddefines.cpp:14
LineSave
t_LineSave LineSave
Definition: lines.cpp:5
FeII_Cont
realnum ** FeII_Cont
Definition: atom_feii.cpp:113
t_CoolHeavy::V38830
double V38830
Definition: coolheavy.h:96
FeIISumBand
double FeIISumBand(realnum wl1, realnum wl2, double *SumBandInward)
Definition: atom_feii.cpp:1355
radius
t_radius radius
Definition: radius.cpp:5
FeIIAddLines
void FeIIAddLines(void)
Definition: atom_feii.cpp:1574
t_FeII::nFeIILevel_malloc
long int nFeIILevel_malloc
Definition: atomfeii.h:193
TransitionProxy::Emis
EmissionList::reference Emis() const
Definition: transition.h:408
ipxK03462
long ipxK03462
Definition: atmdat_readin.cpp:80
linadd
void linadd(double xInten, realnum wavelength, const char *chLab, char chInfo, const char *chComment)
Definition: lines_service.cpp:316
dense.h
t_fe::fe40401
double fe40401
Definition: fe.h:21
coolheavy.h
t_FeII::Fe2_UVsimp_cool
double Fe2_UVsimp_cool
Definition: atomfeii.h:258
trace
t_trace trace
Definition: trace.cpp:5
ipSc05231
long ipSc05231
Definition: atmdat_readin.cpp:84
Fe2LevN
TransitionList Fe2LevN("Fe2LevN", &Fe2LevNStates)
cddefines.h
fe.h
t_CoolHeavy::Cr3l21
double Cr3l21
Definition: coolheavy.h:48
ipCALCIUM
const int ipCALCIUM
Definition: cddefines.h:324
t_LineSave::ipass
long int ipass
Definition: lines.h:75
t_fe::fegrain
realnum fegrain
Definition: fe.h:38
t_fe::fe41007
double fe41007
Definition: fe.h:25
radius.h
ipTFe56
long ipTFe56
Definition: atmdat_readin.cpp:88
ipCaI4228
long ipCaI4228
Definition: atmdat_readin.cpp:81
FeII
t_FeII FeII
Definition: atomfeii.cpp:5
ipFe2LevN
multi_arr< int, 2 > ipFe2LevN
Definition: taulines.cpp:34
ipFeI3457
long ipFeI3457
Definition: atmdat_readin.cpp:85
t_ca::Caz
realnum Caz
Definition: ca.h:18
t_fe::Fe3_wl
double ** Fe3_wl
Definition: fe.h:31
MAX2
#define MAX2
Definition: cddefines.h:782
t_CoolHeavy::Sc33936
double Sc33936
Definition: coolheavy.h:92
t_CoolHeavy::Cr4l31
double Cr4l31
Definition: coolheavy.h:43
ipFe1_24m
long ipFe1_24m
Definition: atmdat_readin.cpp:96
ipKI7745
long ipKI7745
Definition: atmdat_readin.cpp:80
t_dense::xIonDense
double xIonDense[LIMELM][LIMELM+1]
Definition: dense.h:125
TauLines
TransitionList TauLines("TauLines", &AnonStates)
ipTFe26
long ipTFe26
Definition: atmdat_readin.cpp:87
nFeIIBands
long int nFeIIBands
Definition: atom_feii.cpp:116
t_fe::fe42829
double fe42829
Definition: fe.h:19
ipxK04598
long ipxK04598
Definition: atmdat_readin.cpp:80
t_fe::Fe4CoolTot
double Fe4CoolTot
Definition: fe.h:17
ipH2p
const int ipH2p
Definition: iso.h:29
t_fe::Fe3CoolTot
double Fe3CoolTot
Definition: fe.h:30
ipTFe35
long ipTFe35
Definition: atmdat_readin.cpp:88
rt.h
PutLine
void PutLine(const TransitionProxy &t, const char *chComment, const char *chLabelTemp)
Definition: transition.cpp:449
ipFe1_54m
long ipFe1_54m
Definition: atmdat_readin.cpp:96
t_fe::fe41008
double fe41008
Definition: fe.h:26
t_CoolHeavy::Sc24p2m
double Sc24p2m
Definition: coolheavy.h:91
t_phycon::te10
double te10
Definition: phycon.h:55
t_CoolHeavy::Cr5l31
double Cr5l31
Definition: coolheavy.h:45
t_CoolHeavy::Cr5l21
double Cr5l21
Definition: coolheavy.h:46
t_ca::dCaf12
realnum dCaf12
Definition: ca.h:29
t_ca::Caf2
realnum Caf2
Definition: ca.h:20
t_FeII::Fe2_large_cool
double Fe2_large_cool
Definition: atomfeii.h:250
ipFe1_111m
long ipFe1_111m
Definition: atmdat_readin.cpp:96
t_fe::fe41106
double fe41106
Definition: fe.h:24
t_fe::fe41206
double fe41206
Definition: fe.h:23
t_CoolHeavy::c3892
double c3892
Definition: coolheavy.h:49
prt_wl
void prt_wl(FILE *ioOUT, realnum wl)
Definition: prt.cpp:13
t_fe::fe41207
double fe41207
Definition: fe.h:22
t_fe::fe40906
double fe40906
Definition: fe.h:27
ipxK04154
long ipxK04154
Definition: atmdat_readin.cpp:81
ipFeI3021
long ipFeI3021
Definition: atmdat_readin.cpp:86
ipT191
long ipT191
Definition: atmdat_readin.cpp:89
t_iso_sp::trans
TransitionProxy trans(const long ipHi, const long ipLo)
Definition: iso.h:444
taulines.h
CoolHeavy
t_CoolHeavy CoolHeavy
Definition: coolheavy.cpp:5
ipCo11527
long ipCo11527
Definition: atmdat_readin.cpp:89
FeII_Bands
realnum ** FeII_Bands
Definition: atom_feii.cpp:106
phycon.h
t_ca::Cay
realnum Cay
Definition: ca.h:17
t_CoolHeavy::V38507
double V38507
Definition: coolheavy.h:97
ipH1s
const int ipH1s
Definition: iso.h:27
lines_lv1_k_zn
void lines_lv1_k_zn(void)
Definition: prt_lines_lv1_k_zn.cpp:21
t_CoolHeavy::Sc24p1m
double Sc24p1m
Definition: coolheavy.h:90
iso_sp
t_iso_sp iso_sp[NISO][LIMELM]
Definition: iso.cpp:8
t_CoolHeavy::Sc42100
double Sc42100
Definition: coolheavy.h:93
NLFE3
#define NLFE3
Definition: fe.h:10
t_CoolHeavy::Sc43595
double Sc43595
Definition: coolheavy.h:95
nFeIIConBins
long int nFeIIConBins
Definition: atom_feii.cpp:119
t_phycon::te
double te
Definition: phycon.h:11
ipxK07319
long ipxK07319
Definition: atmdat_readin.cpp:81
ipFeI2966
long ipFeI2966
Definition: atmdat_readin.cpp:86
t_CoolHeavy::Cr5l32
double Cr5l32
Definition: coolheavy.h:47
t_CoolHeavy::Sc45058
double Sc45058
Definition: coolheavy.h:94
ipFe1_35m
long ipFe1_35m
Definition: atmdat_readin.cpp:96
t_fe::Fe3_emiss
double ** Fe3_emiss
Definition: fe.h:31
t_ca::Cair
double Cair
Definition: ca.h:22
atomfeii.h
DEBUG_ENTRY
#define DEBUG_ENTRY(funcname)
Definition: cddefines.h:684
ipH_LIKE
const int ipH_LIKE
Definition: iso.h:62
t_fe::fe42836
double fe42836
Definition: fe.h:18
t_trace::lgTrace
bool lgTrace
Definition: trace.h:12
ca.h
Generated by   doxygen 1.8.17