cloudy  trunk
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
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-2008 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 
20 void lines_lv1_k_zn(void)
21 {
22  long int i,
23  ipnt,
24  ilo,
25  ihi;
26 
27  double c10,
28  c14,
29  eff,
30  fela,
31  r14;
32 
33  DEBUG_ENTRY( "lines_lv1_k_zn()" );
34 
35  if( trace.lgTrace )
36  {
37  fprintf( ioQQQ, " lines_lv1_k_zn called\n" );
38  }
39 
41  " potassium K I 7745 ");
42 
44  " [K III] 4.62 microns ");
45 
47  " [KIV] 5.983 min ");
48 
50  " [KIV] 15.39 micron ");
51 
53  " [KVI] 8.823 micron ");
54 
56  " [KVI] 5.575 micron ");
57 
59  " [K VII] 3.189 microns ");
60 
62  " K 11 4249.99A ");
63 
65  " calcium Ca I 4228 ");
66 
67  linadd(ca.Cakh,3933,"Ca 2",'c',
68  " coll excited calcium k+h " );
69 
70  linadd(ca.Cair,8579,"Ca 2",'c' ,
71  " infrared triplet ");
72 
73  linadd(ca.c7306,7306,"Ca 2",'c',
74  " forbidden lines, 7291+7324 together " );
75 
76  linadd(ca.dCakh,3933,"Phot",'i' ,
77  " fraction H Ly-alpha destruction of excited levels ");
78 
79  linadd(ca.dCaf12,7306,"Phot",'i' ,
80  " fraction H Ly-alpha destruction of excited levels ");
81 
82  PntForLine(3934.,"Ca2K",&ipnt);
83  lindst(ca.Cak,3934,"Ca2K",ipnt,'i',true,
84  " individual lines from five level atom");
85 
86 
87  PntForLine(3969.,"Ca2H",&ipnt);
88  lindst(ca.Cah,3969,"Ca2H",ipnt,'i',true,
89  " individual lines from five level atom" );
90 
91 
92  PntForLine(8498.,"Ca2X",&ipnt);
93  lindst(ca.Cax,8498,"Ca2X",ipnt,'i',true,
94  " individual lines from five level atom " );
95 
96 
97  PntForLine(8542.,"Ca2Y",&ipnt);
98  lindst(ca.Cay,8542,"Ca2Y",ipnt,'i',true,
99  " individual lines from five level atom" );
100 
101 
102  PntForLine(8662.,"Ca2Z",&ipnt);
103  lindst(ca.Caz,8662,"Ca2Z",ipnt,'i',true,
104  " individual lines from five level atom" );
105 
106 
107  PntForLine(7291.,"CaF1",&ipnt);
108  lindst(ca.Caf1,7291,"CaF1",ipnt,'i',true,
109  " individual lines from five level atom" );
110 
111 
112  PntForLine(7324.,"CaF2",&ipnt);
113  lindst(ca.Caf2,7324,"CaF2",ipnt,'i',true,
114  " individual lines from five level atom" );
115 
116 
117  eff = dense.eden*dense.xIonDense[ipCALCIUM][2]*5.4e-21/(phycon.te/
119  linadd(eff,3933,"Rec ",'i',
120  " recombination contribution to CaII emission" );
121 
123  " Ca IV 3.2 micron ");
124 
126  " Ca V 4.16, 11.47 micron");
127 
129  " Ca V 4.16, 11.47 micron ");
130 
131 
132  PntForLine(6087.,"Ca 5",&ipnt);
133  lindst(ca.Ca6087,6087,"Ca 5",ipnt,'c',true ,
134  " Ca V optical and UV lines, collisional excitation, 3-level atom");
135 
136 
137  PntForLine(5311.,"Ca 5",&ipnt);
138  lindst(ca.c5311,5311,"Ca 5",ipnt,'c',true ,
139  " Ca V optical and UV lines, collisional excitation, 3-level atom");
140 
141 
142  PntForLine(2414.,"Ca 5",&ipnt);
143  lindst(ca.c2414,2414,"Ca 5",ipnt,'c',true ,
144  " Ca V optical and UV lines, collisional excitation, 3-level atom");
145 
146 
147  PntForLine(3997.,"Ca 5",&ipnt);
148  lindst(ca.c3997,3997,"Ca 5",ipnt,'c',true,
149  " Ca V optical and UV lines, collisional excitation, 3-level atom" );
150 
152  " [Ca VII] 4.09 microns" );
153 
155  " [Ca VII] 6.15 microns " );
156 
157 
158  PntForLine(5620.,"Ca 7",&ipnt);
159  lindst(ca.Ca5620,5620,"Ca 7",ipnt,'c',true,
160  " Ca VII optical and UV lines, collisional excitation, 3-level atom" );
161 
162 
163  PntForLine(4941.,"Ca 7",&ipnt);
164  lindst(ca.Ca4941,4941,"Ca 7",ipnt,'c',true,
165  " Ca VII optical and UV lines, collisional excitation, 3-level atom" );
166 
167 
168  PntForLine(2112.,"Ca 7",&ipnt);
169  lindst(ca.Ca2112,2112,"Ca 7",ipnt,'c',true ,
170  " Ca VII optical and UV lines, collisional excitation, 3-level atom");
171 
172 
173  PntForLine(3688.,"Ca 7",&ipnt);
174  lindst(ca.Ca3688,3688,"Ca 7",ipnt,'c',true,
175  " Ca VII optical and UV lines, collisional excitation, 3-level atom" );
176 
178  " [Ca VIII] 2.32 microns, A Saraph and Strey ");
179 
181  " [Ca 12] 3328.78A ");
182 
184  " Ca 18 Li seq 2s2p, 302, 345 separate ");
185 
187  " Ca 18 Li seq 2s2p, 302, 345 separate ");
188 
190  " Ca 18 Li seq 2s3p, ");
191 
192 
193  PntForLine(22.08e4,"Sc 2",&ipnt);
194  lindst(CoolHeavy.Sc22p08m,21,"Sc 2",ipnt,'c',true,
195  " Sc II 2.08 (1-3) " );
196 
197 
198  PntForLine(24.1e4,"Sc 2",&ipnt);
199  lindst(CoolHeavy.Sc24p1m,41,"Sc 2",ipnt,'c',true,
200  " Sc II 4.1 micron (1-2)" );
201 
202 
203  PntForLine(24.2e4,"Sc 2",&ipnt);
204  lindst(CoolHeavy.Sc24p2m,42,"Sc 2",ipnt,'c',true,
205  " Sc II 4.22 (2-3)" );
206 
207 
208  PntForLine(3933.,"Sc 3",&ipnt);
209  lindst(CoolHeavy.Sc33936,3933,"Sc 3",ipnt,'c',true,
210  " Sc III 3936" );
211 
213  " [Sc V] 1.46 microns ");
214 
215 
216  PntForLine(5054.,"Sc 6",&ipnt);
217  lindst(CoolHeavy.Sc45058,5054,"Sc 6",ipnt,'c',true ,
218  " Sc VI 5054 (1-2)");
219 
220 
221  PntForLine(3592.,"Sc 6",&ipnt);
222  lindst(CoolHeavy.Sc43595,3592,"Sc 6",ipnt,'c',true,
223  " Sc VI 3595 (2-3)" );
224 
225 
226  PntForLine(2100.,"Sc 6",&ipnt);
227  lindst(CoolHeavy.Sc42100,2100,"Sc 6",ipnt,'c',true,
228  " Sc VI 2100 (1-3)" );
229 
231  " [Sc 13] 2637.97A");
232 
233 
234  PntForLine(1.21e4,"Ti 3",&ipnt);
235  lindst(CoolHeavy.Ti3l21,12,"Ti 3",ipnt,'c',true,
236  " Ti III 1.21 micron, (actually multiplet) 2-1 transition from model atom " );
237 
238 
239  PntForLine(9594.,"Ti 3",&ipnt);
240  lindst(CoolHeavy.Ti3l31,9594,"Ti 3",ipnt,'c',true,
241  " Ti III 9594, 3-1 transition, (actually multiplet) from model atom" );
242 
243  PntForLine(4.57e4,"Ti 3",&ipnt);
244  lindst(CoolHeavy.Ti3l32,45,"Ti 3",ipnt,'c',true,
245  " Ti III 4.57 micron, 3-2 transition, (actually multiplet) from model atom" );
246 
248  " [Ti VI] 1.72 microns ");
249 
251  " [Ti XIV] 2117.79 ");
252 
253 
254  PntForLine(8823.,"V 3",&ipnt);
255  lindst(CoolHeavy.V38830,8823,"V 3",ipnt,'c',true ,
256  " V III 8823 ");
257 
258 
259  PntForLine(8507.,"V 3",&ipnt);
260  lindst(CoolHeavy.V38507,8507,"V 3",ipnt,'c',true,
261  " V III 8507" );
262 
263 
264  PntForLine(7735.,"V 4",&ipnt);
265  lindst(CoolHeavy.V47741,7735,"V 4",ipnt,'c',true,
266  " V IV 7741 1-3" );
267 
268 
269  PntForLine(9489.,"V 4",&ipnt);
270  lindst(CoolHeavy.V49496,9489,"V 4",ipnt,'c',true,
271  " V IV 9496 2-1 " );
272 
273 
274  PntForLine(4.19e4,"V 4",&ipnt);
275  lindst(CoolHeavy.V44p2m,42,"V 4",ipnt,'c',true,
276  " V IV 4.19 micron 3-2" );
277 
279  " [V VII] 1.304 microns ");
280 
282  " [V 15] 1721.38 ");
283 
284 
285  PntForLine(5828.,"Cr 3",&ipnt);
286  lindst(CoolHeavy.Cr3l21,5828,"Cr 3",ipnt,'c',true,
287  " [CrIII] multiplet blend at 5828A" );
288 
289  PntForLine(7267.,"Cr 4",&ipnt);
290  lindst(CoolHeavy.Cr4l21,7267,"Cr 4",ipnt,'c',true,
291  " [CrIV] 2 - 1 multiplet blend at 7272" );
292 
293 
294  PntForLine(6801.,"Cr 4",&ipnt);
295  lindst(CoolHeavy.Cr4l31,6801,"Cr 4",ipnt,'c',true,
296  " [CrIV] 3 - 1 multiplet blend at 6806" );
297 
298 
299  PntForLine(7979.,"Cr 5",&ipnt);
300  lindst(CoolHeavy.Cr5l21,7979,"Cr 5",ipnt,'c',true,
301  " [CrV] 2 - 1 multiplet blend at 7985" );
302 
303  PntForLine(6577.,"Cr 5",&ipnt);
304  lindst(CoolHeavy.Cr5l31,6577,"Cr 5",ipnt,'c',true,
305  " [CrV] 3 - 1 multiplet blend at 6582" );
306 
307 
308  PntForLine(3.75e4,"Cr 5",&ipnt);
309  lindst(CoolHeavy.Cr5l32,37,"Cr 5",ipnt,'c',true,
310  " [CrV] 3 - 2 multiplet blend at 3.75 microns " );
311 
313  " [Cr VIII] 1.01 microns ");
314 
316  " [Cr 16] 1410.60 ");
317 
319  " [Mn IX] 7968.5 A ");
320 
322  " [Mn 17] 1169.59 ");
323 
324  /* bob Rubin's UV line
325  * f2 = dense.xIonDense(26,4)*sexp(50 764./te)*0.45*cdsqte/6.*7.01e-12
326  * call linadd( f2 , 2837 , 'BobR' , 'i')
327  * f2 = dense.xIonDense(26,4)*sexp(55 989./te)*0.384*cdsqte/6.*7.74e-12
328  * call linadd( f2 , 2568 , 'BobR' , 'i') */
329 
330  /* iron */
331 
333  " Fe 1 24m ");
334 
336  " Fe 1 35m ");
337 
339  " Fe 1 54m ");
340 
342  " Fe 1 111m ");
343 
345  " Fe 1 3884 ");
346 
348  " Fe 1 3729 ");
349 
351  " Fe 1 3457 ");
352 
354  " Fe 1 3021 ");
355 
357  " Fe 1 2966 ");
358 
360  "total of all Fe 2, both simple UV and large atom together ");
361 
362  linadd(MAX2(0.,-FeII.Fe2_large_cool-FeII.Fe2_UVsimp_cool),0,"Fe2h",'h' ,
363  "total of all Fe 2 heating, both simple UV and large atom together ");
364 
365  linadd(FeII.for7,4300,"Fe 2",'i' ,
366  " Fe 2 forbidden 2-1 transition from Netzer's atom ");
367 
369  " 2400 ");
371  " 6200 ");
373  " 1080 ");
375  " 1500 ");
377  " 11500 ");
379  " 2500 ");
381  " 2300 ");
383  " 8900 ");
384 
385  /* >>chng 03 nov 02, changed wl in lindst calls to be wl in microns,
386  * this will change wl in printout */
387  /* results of FeIILevelPops Fe 2 atom,
388  * not all are printed, first is */
389  /* total cooling due to 16 level atom
390  * >>chng 05 dec 04, rm this since no longer meaningful
391  linadd(FeII.Fe2_16levl_cool,0,"Fe 2",'c' );*/
392 
393 
394  PntForLine(16640.,"Fe 2",&ipnt);
395  lindst(FeII.fe21308,16640,"Fe 2",ipnt,'i',true,
396  " Fe 2 1.664 microns 8-13 " );
397 
398 
399  PntForLine(15990.,"Fe 2",&ipnt);
400  lindst(FeII.fe21207,15990,"Fe 2",ipnt,'i',true,
401  " Fe 2 1.599 microns 7-12" );
402 
403 
404  PntForLine(15340.,"Fe 2",&ipnt);
405  lindst(FeII.fe21106,15340,"Fe 2",ipnt,'i',true ,
406  " Fe 2 1.534 microns 6-11 ");
407 
408 
409  PntForLine(16440.,"Fe 2",&ipnt);
410  lindst(FeII.fe21006,16440,"Fe 2",ipnt,'i',true,
411  " Fe 2 1.644 microns 6-10 " );
412 
413 
414  PntForLine(12790.,"Fe 2",&ipnt);
415  lindst(FeII.fe21204,12790,"Fe 2",ipnt,'i',true,
416  " Fe 2 1.279 microns 12-4" );
417 
418 
419  PntForLine(12950.,"Fe 2",&ipnt);
420  lindst(FeII.fe21103,12950,"Fe 2",ipnt,'i',true,
421  " Fe 2 1.295 microns 11-3 " );
422 
423 
424  PntForLine(13280.,"Fe 2",&ipnt);
425  lindst(FeII.fe21104,13280,"Fe 2",ipnt,'i',true,
426  " Fe 2 1.328 microns 11-4" );
427 
428 
429  PntForLine(12570.,"Fe 2",&ipnt);
430  lindst(FeII.fe21001,12570,"Fe 2",ipnt,'i',true,
431  " Fe 2 1.257 microns 10-1 " );
432 
433 
434  PntForLine(13210.,"Fe 2",&ipnt);
435  lindst(FeII.fe21002,13210,"Fe 2",ipnt,'i',true,
436  " Fe 2 1.321 microns 10-2 " );
437 
438 
439  PntForLine(259880.,"Fe 2",&ipnt);
440  lindst(FeII.fe20201,259880,"Fe 2",ipnt,'i',true,
441  " Fe 2 25.988 microns 2-1" );
442 
443  PntForLine(353480.,"Fe 2",&ipnt);
444  lindst(FeII.fe20302,353480,"Fe 2",ipnt,'i',true,
445  " Fe 2 35.348 microns 3-2" );
446 
447 
448  PntForLine(179360.,"Fe 2",&ipnt);
449  lindst(FeII.fe20706,179360,"Fe 2",ipnt,'i',true,
450  " Fe 2 17.936 microns 7-6, label is 178 to be unique " );
451 
452 
453  PntForLine(245180.,"Fe 2",&ipnt);
454  lindst(FeII.fe20807,245180,"Fe 2",ipnt,'i',true,
455  " Fe 2 24.518 microns 8-7" );
456 
457  PntForLine(357760.,"Fe 2",&ipnt);
458  lindst(FeII.fe20908,357760,"Fe 2",ipnt,'i',true,
459  " Fe 2 35.776 microns 9-8" );
460 
461  PntForLine(18100.,"Fe 2",&ipnt);
462  lindst(FeII.fe21007,18100,"Fe 2",ipnt,'i',true,
463  " Fe 2 1.810 microns 10-7" );
464 
465 
466  PntForLine(16770.,"Fe 2",&ipnt);
467  lindst(FeII.fe21107,16770,"Fe 2",ipnt,'i',true,
468  " Fe 2 1.677 microns 11-7 " );
469 
470 
471  PntForLine(18000.,"Fe 2",&ipnt);
472  lindst(FeII.fe21108,18000,"Fe 2",ipnt,'i',true,
473  " Fe 2 1.800 microns 11-8" );
474 
475 
476  PntForLine(17120.,"Fe 2",&ipnt);
477  lindst(FeII.fe21208,17120,"Fe 2",ipnt,'i',true,
478  " Fe 2 1.712 microns 12-8 " );
479 
480 
481  PntForLine(17980.,"Fe 2",&ipnt);
482  lindst(FeII.fe21209,17980,"Fe 2",ipnt,'i',true,
483  " Fe 2 1.798 microns 12-9" );
484 
485 
486  PntForLine(229020.,"Fe 2",&ipnt);
487  lindst(FeII.fe21110,229020,"Fe 2",ipnt,'i',true,
488  " Fe 2 22.902 microns 11-10 " );
489 
490 
491  PntForLine(346600.,"Fe 2",&ipnt);
492  lindst(FeII.fe21211,346600,"Fe 2",ipnt,'i',true,
493  " Fe 2 34.660 microns 12-11" );
494 
495  /* >>chng 97 mar 14, added following 4 lines after consultation with Verner
496  * these are typically within 5-9% of their large model atom. These are
497  * predicted with both large and small atom */
498 
499 
500  PntForLine(8617.,"Fe 2",&ipnt);
501  lindst(FeII.fe21406,8617,"Fe 2",ipnt,'i',true,
502  " Fe 2 8617A 14-06," );
503 
504 
505  PntForLine(8892.,"Fe 2",&ipnt);
506  lindst(FeII.fe21507,8892,"Fe 2",ipnt,'i',true,
507  " Fe 2 8892A 15-07" );
508 
509 
510  PntForLine(9227.,"Fe 2",&ipnt);
511  lindst(FeII.fe21508,9227,"Fe 2",ipnt,'i',true,
512  " Fe 2 9227A 15-08 " );
513 
514 
515  PntForLine(9268.,"Fe 2",&ipnt);
516  lindst(FeII.fe21609,9268,"Fe 2",ipnt,'i',true,
517  " Fe 2 9268A 16-09" );
518 
519  /* following lines are only non-zero when large atom is enabled. When
520  * running with small FeII atoms these are set zero in routine */
521 
522 
523  PntForLine(5159.,"Fe 2",&ipnt);
524  lindst(FeII.fe25to6,5159,"Fe 2",ipnt,'i',true,
525  " Fe 2 5158.81 25-6" );
526 
527 
528  PntForLine(5262.,"Fe 2",&ipnt);
529  lindst(FeII.fe27to7,5262,"Fe 2",ipnt,'i',true,
530  " Fe 2 5261.61 27-7" );
531 
532  PntForLine(5334.,"Fe 2",&ipnt);
533  lindst(FeII.fe28to8,5334,"Fe 2",ipnt,'i',true,
534  " Fe 2 5333.65 28-8 " );
535 
536 
537  PntForLine(5376.,"Fe 2",&ipnt);
538  lindst(FeII.fe29to9,5376,"Fe 2",ipnt,'i',true,
539  " Fe 2 5376.47 29-9 " );
540 
541  PntForLine(4815.,"Fe 2",&ipnt);
542  lindst(FeII.fe32to6,4815,"Fe 2",ipnt,'i',true,
543  " Fe 2 4814.55 32-6 " );
544 
545  PntForLine(4905.,"Fe 2",&ipnt);
546  lindst(FeII.fe33to7,4905,"Fe 2",ipnt,'i',true,
547  " Fe 2 4905.35 33-7" );
548 
549 
550  PntForLine(4347.,"Fe 2",&ipnt);
551  lindst(FeII.fe37to7,4347,"Fe 2",ipnt,'i',true,
552  " Fe 2 4346.85 37-7 " );
553 
554 
555  PntForLine(4353.,"Fe 2",&ipnt);
556  lindst(FeII.fe39to8,4353,"Fe 2",ipnt,'i',true,
557  " Fe 2 4352.78 39-8 " );
558 
559 
560  PntForLine(4372.,"Fe 2",&ipnt);
561  lindst(FeII.fe40to9,4372,"Fe 2",ipnt,'i',true,
562  " Fe 2 4372.43 40-9 " );
563 
564 
565  PntForLine(4244.,"Fe 2",&ipnt);
566  lindst(FeII.fe37to6,4244,"Fe 2",ipnt,'i',true,
567  " Fe 2 4243.98 37-6" );
568 
569 
570  PntForLine(4277.,"Fe 2",&ipnt);
571  lindst(FeII.fe39to7,4277,"Fe 2",ipnt,'i',true,
572  " Fe 2 4276.83 39-7" );
573 
574 
575  PntForLine(4320.,"Fe 2",&ipnt);
576  lindst(FeII.fe40to8,4320,"Fe 2",ipnt,'i',true,
577  " Fe 2 4319.62 40-8" );
578 
579 
580  PntForLine(4358.,"Fe 2",&ipnt);
581  lindst(FeII.fe41to9,4358,"Fe 2",ipnt,'i',true,
582  " Fe 2 4358.37 41-9 " );
583 
584 
585  PntForLine(4177.,"Fe 2",&ipnt);
586  lindst(FeII.fe39to6,4177,"Fe 2",ipnt,'i',true,
587  " Fe 2 4177.21 39-6 " );
588 
589 
590  PntForLine(4245.,"Fe 2",&ipnt);
591  lindst(FeII.fe40to7,4245,"Fe 2",ipnt,'i',true,
592  " Fe 2 4244.81 40-7 " );
593 
594 
595  PntForLine(4306.,"Fe 2",&ipnt);
596  lindst(FeII.fe41to8,4306,"Fe 2",ipnt,'i',true,
597  " Fe 2 4305.90 41-8" );
598 
599  /* >>chng 99 jun 18 new lines added as per Katya email, so all orion lines are in */
600 
601  PntForLine(4114.,"Fe 2",&ipnt);
602  lindst(FeII.fe42to6,4114,"Fe 2",ipnt,'i',true,
603  " Fe 2 4114.48 42-6 " );
604 
605 
606  PntForLine(4179.,"Fe 2",&ipnt);
607  lindst(FeII.fe43to7,4179,"Fe 2",ipnt,'i',true,
608  " Fe 2 4178.95 7- 43" );
609 
610 
611  PntForLine(4211.,"Fe 2",&ipnt);
612  lindst(FeII.fe42to7,4211,"Fe 2",ipnt,'i',true,
613  " Fe 2 4211.10 7- 42 " );
614 
615 
616  PntForLine(4314.,"Fe 2",&ipnt);
617  lindst(FeII.fe80to28,4314,"Fe 2",ipnt,'i',true,
618  " Fe 2 4314.29 28- 80" );
619 
620  PntForLine(4359.,"Fe 2",&ipnt);
621  lindst(FeII.fe36to2,4359,"Fe 2",ipnt,'i',true,
622  " Fe 2 4359.34 2- 36" );
623 
624 
625  PntForLine(4414.,"Fe 2",&ipnt);
626  lindst(FeII.fe36to3,4414,"Fe 2",ipnt,'i',true,
627  " Fe 2 4413.78 3- 36 " );
628 
629 
630  PntForLine(4416.,"Fe 2",&ipnt);
631  lindst(FeII.fe32to1,4416,"Fe 2",ipnt,'i',true ,
632  " Fe 2 4416.27 1- 32");
633 
634 
635  PntForLine(4458.,"Fe 2",&ipnt);
636  lindst(FeII.fe33to2,4458,"Fe 2",ipnt,'i',true,
637  " Fe 2 4457.95 2- 33 " );
638 
639 
640  PntForLine(4475.,"Fe 2",&ipnt);
641  lindst(FeII.fe36to5,4475,"Fe 2",ipnt,'i',true,
642  " Fe 2 4474.91 5- 36 " );
643 
644 
645  PntForLine(4493.,"Fe 2",&ipnt);
646  lindst(FeII.fe32to2,4493,"Fe 2",ipnt,'i',true,
647  " Fe 2 4492.64 2- 32" );
648 
649 
650  PntForLine(4515.,"Fe 2",&ipnt);
651  lindst(FeII.fe33to3,4515,"Fe 2",ipnt,'i',true,
652  " Fe 2 4514.90 3- 33" );
653 
654 
655  PntForLine(4728.,"Fe 2",&ipnt);
656  lindst(FeII.fe30to3,4728,"Fe 2",ipnt,'i',true,
657  " Fe 2 4728.07 3- 30" );
658 
659 
660  PntForLine(4775.,"Fe 2",&ipnt);
661  lindst(FeII.fe33to6,4775,"Fe 2",ipnt,'i',true,
662  " Fe 2 4774.74 6- 33" );
663 
664  PntForLine(4890.,"Fe 2",&ipnt);
665  lindst(FeII.fe24to2,4890,"Fe 2",ipnt,'i',true,
666  " Fe 2 4889.63 2- 24" );
667 
668 
669  PntForLine(4947.,"Fe 2",&ipnt);
670  lindst(FeII.fe32to7,4947,"Fe 2",ipnt,'i',true,
671  " Fe 2 4947.37 7- 32" );
672 
673 
674  PntForLine(4951.,"Fe 2",&ipnt);
675  lindst(FeII.fe35to8,4951,"Fe 2",ipnt,'i',true,
676  " Fe 2 4950.74 8- 35" );
677 
678 
679  PntForLine(4973.,"Fe 2",&ipnt);
680  lindst(FeII.fe34to8,4973,"Fe 2",ipnt,'i',true,
681  " Fe 2 4973.39 8- 34" );
682 
683 
684  PntForLine(5112.,"Fe 2",&ipnt);
685  lindst(FeII.fe27to6,5112,"Fe 2",ipnt,'i',true,
686  " Fe 2 5111.63 6- 27" );
687 
688 
689  PntForLine(5220.,"Fe 2",&ipnt);
690  lindst(FeII.fe28to7,5220,"Fe 2",ipnt,'i',true,
691  " Fe 2 5220.06 7- 28" );
692 
693 
694  PntForLine(5269.,"Fe 2",&ipnt);
695  lindst(FeII.fe30to8,5269,"Fe 2",ipnt,'i',true,
696  " Fe 2 5268.880 8- 30 " );
697 
698 
699  PntForLine(5273.,"Fe 2",&ipnt);
700  lindst(FeII.fe24to6,5273,"Fe 2",ipnt,'i',true,
701  " Fe 2 5273.35 6- 24" );
702 
703  PntForLine(5297.,"Fe 2",&ipnt);
704  lindst(FeII.fe29to8,5297,"Fe 2",ipnt,'i',true,
705  " Fe 2 5296.840 8- 29 " );
706 
707 
708  PntForLine(5433.,"Fe 2",&ipnt);
709  lindst(FeII.fe24to7,5433,"Fe 2",ipnt,'i',true,
710  " Fe 2 5433.13 7- 24" );
711 
712 
713  PntForLine(5527.,"Fe 2",&ipnt);
714  lindst(FeII.fe22to7,5527,"Fe 2",ipnt,'i',true,
715  " Fe 2 5527.340 7- 22" );
716 
717 
718  PntForLine(5747.,"Fe 2",&ipnt);
719  lindst(FeII.fe38to11,5747,"Fe 2",ipnt,'i',true,
720  " Fe 2 5746.96 11- 38" );
721 
722 
723  PntForLine(6440.,"Fe 2",&ipnt);
724  lindst(FeII.fe19to8,6440,"Fe 2",ipnt,'i',true,
725  " Fe 2 6440.400 8- 19" );
726 
727  PntForLine(7155.,"Fe 2",&ipnt);
728  lindst(FeII.fe17to6,7155,"Fe 2",ipnt,'i',true,
729  " Fe 2 7155.140 6- 17 " );
730 
731 
732  PntForLine(7172.,"Fe 2",&ipnt);
733  lindst(FeII.fe18to7,7172,"Fe 2",ipnt,'i',true,
734  " Fe 2 7171.980 7- 18" );
735 
736 
737  PntForLine(7388.,"Fe 2",&ipnt);
738  lindst(FeII.fe18to8,7388,"Fe 2",ipnt,'i',true,
739  " Fe 2 7388.160 8- 18" );
740 
741 
742  PntForLine(7453.,"Fe 2",&ipnt);
743  lindst(FeII.fe17to7,7453,"Fe 2",ipnt,'i',true,
744  " Fe 2 7452.540 7- 17" );
745 
746  /* option to save all intensities predicted by large FeII atom,
747  * code is in FeIILevelPops */
748  FeIIAddLines();
749 
750  /* some bands from the large model Fe II atom
751  * only done when full large atom is turned on */
752  if( FeII.lgFeIILargeOn )
753  {
754  for( i=0; i < nFeIIBands; i++ )
755  {
756  /* [i][0] is center wavelength, [i][1] and [i][2] are upper and
757  * lower bounds in Angstroms. These are set in FeIIZero
758  * units are erg s-1 cm-3 */
759  eff = FeIISumBand(FeII_Bands[i][1],FeII_Bands[i][2]);
760 
761  linadd(eff,FeII_Bands[i][0],"Fe2b",'i' ,
762  " total Fe II emission in Fe II bands, as defined in bands_Fe2.dat ");
763  }
765  if( LineSave.ipass > 0 )
766  {
767  for( i=0; i < nFeIIConBins; i++ )
768  {
769  /* [i][0] is total intensity in cell, [i][1] and [i][2] are lower and
770  * upper bounds in Angstroms. these are set in FeIIZero *
771  * find total emission from large FeII atom, integrated over band */
772  FeII_Cont[i][0] +=
774  /*fprintf(ioQQQ,"DEBUG feii\t%li\t%.2e\n", i, FeII_Cont[i][0]);*/
775  }
776  }
777  else
778  {
779  for( i=0; i < nFeIIConBins; i++ )
780  {
781  /* initialize arrays */
782  FeII_Cont[i][0] = 0.;
783  }
784  }
785  }
786  else
787  {
788  /* must still add the bands in case of grid where FeII is turned off and on */
789  for( i=0; i < nFeIIBands; i++ )
790  {
791 
792  linadd(0., FeII_Bands[i][0] ,"Fe2b",'i' ,
793  " total Fe II emission in Fe II bands, as defined in bands_Fe2.dat ");
794  }
795  }
796 
798  " anomalous Fe 2 transition at 1787, RMT 191");
799 
800  linadd(fe.Fe3CoolTot,0,"Fe3c",'c' ,
801  " chng 05 dec 16, FeIII code created by Kevin Blagrave Fe3c 0 - total cooling due to 14-level Fe 3 atom ");
802 
803  /* Fe 3 14-level atom
804  * following from print statements within loop */
805  /* Fe 3 22.92m from Blagrave 14-level atom */
806  /* Fe 3 13.53m from Blagrave 14-level atom */
807  /* Fe 3 33.03m from Blagrave 14-level atom */
808  /* Fe 3 10.72m from Blagrave 14-level atom */
809  /* Fe 3 20.15m from Blagrave 14-level atom */
810  /* Fe 3 51.67m from Blagrave 14-level atom */
811  /* Fe 3 9.732m from Blagrave 14-level atom */
812  /* Fe 3 16.91m from Blagrave 14-level atom */
813  /* Fe 3 34.66m from Blagrave 14-level atom */
814  /* Fe 3 105.3m from Blagrave 14-level atom */
815  /* Fe 3 5152A from Blagrave 14-level atom */
816  /* Fe 3 5271A from Blagrave 14-level atom */
817  /* Fe 3 5356A from Blagrave 14-level atom */
818  /* Fe 3 5412A from Blagrave 14-level atom */
819  /* Fe 3 5440A from Blagrave 14-level atom */
820  /* Fe 3 4986A from Blagrave 14-level atom */
821  /* Fe 3 5097A from Blagrave 14-level atom */
822  /* Fe 3 5177A from Blagrave 14-level atom */
823  /* Fe 3 5230A from Blagrave 14-level atom */
824  /* Fe 3 5256A from Blagrave 14-level atom */
825  /* Fe 3 15.47m from Blagrave 14-level atom */
826  /* Fe 3 4925A from Blagrave 14-level atom */
827  /* Fe 3 5033A from Blagrave 14-level atom */
828  /* Fe 3 5111A from Blagrave 14-level atom */
829  /* Fe 3 5162A from Blagrave 14-level atom */
830  /* Fe 3 5188A from Blagrave 14-level atom */
831  /* Fe 3 11.16m from Blagrave 14-level atom */
832  /* Fe 3 40.04m from Blagrave 14-level atom */
833  /* Fe 3 4881A from Blagrave 14-level atom */
834  /* Fe 3 4988A from Blagrave 14-level atom */
835  /* Fe 3 5064A from Blagrave 14-level atom */
836  /* Fe 3 5114A from Blagrave 14-level atom */
837  /* Fe 3 5139A from Blagrave 14-level atom */
838  /* Fe 3 9.282m from Blagrave 14-level atom */
839  /* Fe 3 23.21m from Blagrave 14-level atom */
840  /* Fe 3 55.20m from Blagrave 14-level atom */
841  /* Fe 3 4833A from Blagrave 14-level atom */
842  /* Fe 3 4937A from Blagrave 14-level atom */
843  /* Fe 3 5012A from Blagrave 14-level atom */
844  /* Fe 3 5061A from Blagrave 14-level atom */
845  /* Fe 3 5085A from Blagrave 14-level atom */
846  /* Fe 3 7.789m from Blagrave 14-level atom */
847  /* Fe 3 15.69m from Blagrave 14-level atom */
848  /* Fe 3 25.79m from Blagrave 14-level atom */
849  /* Fe 3 48.41m from Blagrave 14-level atom */
850  /* Fe 3 4714A from Blagrave 14-level atom */
851  /* Fe 3 4813A from Blagrave 14-level atom */
852  /* Fe 3 4884A from Blagrave 14-level atom */
853  /* Fe 3 4931A from Blagrave 14-level atom */
854  /* Fe 3 4954A from Blagrave 14-level atom */
855  /* Fe 3 5.543m from Blagrave 14-level atom */
856  /* Fe 3 8.638m from Blagrave 14-level atom */
857  /* Fe 3 11.01m from Blagrave 14-level atom */
858  /* Fe 3 13.76m from Blagrave 14-level atom */
859  /* Fe 3 19.22m from Blagrave 14-level atom */
860  /* Fe 3 4659A from Blagrave 14-level atom */
861  /* Fe 3 4755A from Blagrave 14-level atom */
862  /* Fe 3 4825A from Blagrave 14-level atom */
863  /* Fe 3 4870A from Blagrave 14-level atom */
864  /* Fe 3 4893A from Blagrave 14-level atom */
865  /* Fe 3 4.859m from Blagrave 14-level atom */
866  /* Fe 3 7.085m from Blagrave 14-level atom */
867  /* Fe 3 8.608m from Blagrave 14-level atom */
868  /* Fe 3 10.20m from Blagrave 14-level atom */
869  /* Fe 3 12.92m from Blagrave 14-level atom */
870  /* Fe 3 39.41m from Blagrave 14-level atom */
871  /* Fe 3 4608A from Blagrave 14-level atom */
872  /* Fe 3 4702A from Blagrave 14-level atom */
873  /* Fe 3 4770A from Blagrave 14-level atom */
874  /* Fe 3 4814A from Blagrave 14-level atom */
875  /* Fe 3 4836A from Blagrave 14-level atom */
876  /* Fe 3 4.356m from Blagrave 14-level atom */
877  /* Fe 3 6.063m from Blagrave 14-level atom */
878  /* Fe 3 7.146m from Blagrave 14-level atom */
879  /* Fe 3 8.208m from Blagrave 14-level atom */
880  /* Fe 3 9.884m from Blagrave 14-level atom */
881  /* Fe 3 20.34m from Blagrave 14-level atom */
882  /* Fe 3 42.06m from Blagrave 14-level atom */
883  /* Fe 3 4574A from Blagrave 14-level atom */
884  /* Fe 3 4668A from Blagrave 14-level atom */
885  /* Fe 3 4734A from Blagrave 14-level atom */
886  /* Fe 3 4778A from Blagrave 14-level atom */
887  /* Fe 3 4800A from Blagrave 14-level atom */
888  /* Fe 3 4.077m from Blagrave 14-level atom */
889  /* Fe 3 5.535m from Blagrave 14-level atom */
890  /* Fe 3 6.423m from Blagrave 14-level atom */
891  /* Fe 3 7.269m from Blagrave 14-level atom */
892  /* Fe 3 8.554m from Blagrave 14-level atom */
893  /* Fe 3 15.41m from Blagrave 14-level atom */
894  /* Fe 3 25.31m from Blagrave 14-level atom */
895  /* Fe 3 63.56m from Blagrave 14-level atom */
896  for( ihi=1; ihi<NLFE3; ++ihi )
897  {
898  for( ilo=0; ilo<ihi; ++ilo )
899  {
900  /* emission in these lines */
901  PntForLine(fe.Fe3_wl[ihi][ilo],"Fe 3",&ipnt);
902 # if 0
903  fprintf( ioQQQ,"\t/* FeIII ");
904  prt_wl( ioQQQ , (realnum)(fe.Fe3_wl[ihi][ilo]+0.5) );
905  fprintf( ioQQQ," from Blagrave 14-level atom */\n ");
906 # endif
907  lindst( fe.Fe3_emiss[ihi][ilo] , (realnum)(fe.Fe3_wl[ihi][ilo]+0.5) , "Fe 3",ipnt,'c',true,
908  " " );
909  }
910  }
911 
912  /*>>chng 05 dec 18, following are now in the above */
913  /* sum of 3p and 3g states together */
914  /* linadd(CoolHeavy.c5270,0,"Fe 3",'c' ); */
915 
916  /* Fe 3 5270, predictions from Garstang et al 78
917  PntForLine(5270.,"Fe 3",&ipnt);
918  lindst(CoolHeavy.c5270*0.2090,5270,"Fe 3",ipnt,'c',true );*/
919 
920  /* Fe 3 5270, predictions from Garstang et al 78
921  PntForLine(4658.,"Fe 3",&ipnt);
922  lindst(CoolHeavy.c5270*0.3667,4658,"Fe 3",ipnt,'c',true ); */
923 
924  PutLine(&TauLines[ipT1122]," Fe 3 1122 entire multiplet");
925 
926  linadd(fe.Fe4CoolTot,0,"Fe4c",'c',
927  " Fe4c 0 - total cooling due to 12-level Fe 4 atom " );
928 
929 
930  PntForLine(3096.,"Fe 4",&ipnt);
931  lindst(fe.fe40401,3096,"Fe 4",ipnt,'i',true,
932  " Fe 4 3096.A, 4-1 and 5-1 transitions together" );
933 
934 
935  PntForLine(2836.,"Fe 4",&ipnt);
936  lindst(fe.fe42836,2836,"Fe 4",ipnt,'i',true,
937  " Fe 4 2835.7A, 6-1 transition, 4P5/2 - 6S5/2 " );
938 
939 
940  PntForLine(2829.,"Fe 4",&ipnt);
941  lindst(fe.fe42829,2829,"Fe 4",ipnt,'i',true,
942  " Fe 4 2829.4A, 7-1 transition, 4P3/2 - 6S5/2" );
943 
944 
945  PntForLine(2567.,"Fe 4",&ipnt);
946  lindst(fe.fe42567,2567,"Fe 4",ipnt,'i',true,
947  " Fe 4 2567.6+ 2567.4. 11-1 and 12-1 transitions" );
948 
949 
950  PntForLine(2.774e4,"Fe 4",&ipnt);
951  lindst(fe.fe41207,277,"Fe 4",ipnt,'i',true,
952  " Fe 4 2.774 microns 12-7 transition " );
953 
954 
955  PntForLine(2.714e4,"Fe 4",&ipnt);
956  lindst(fe.fe41206,271,"Fe 4",ipnt,'i',true,
957  " Fe 4 2.714 microns 12-6 transition " );
958 
959 
960  PntForLine(2.716e4,"Fe 4",&ipnt);
961  lindst(fe.fe41106,272,"Fe 4",ipnt,'i',true,
962  " Fe 4 2.716 microns 11-6 transition" );
963 
964 
965  PntForLine(2.806e4,"Fe 4",&ipnt);
966  lindst(fe.fe41007,281,"Fe 4",ipnt,'i',true,
967  " Fe 4 2.806 microns 10-7 transition " );
968 
969 
970  PntForLine(2.865e4,"Fe 4",&ipnt);
971  lindst(fe.fe41008,287,"Fe 4",ipnt,'i',true ,
972  " Fe 4 2.865 microns 10-8 transition");
973 
974 
975  PntForLine(2.836e4,"Fe 4",&ipnt);
976  lindst(fe.fe40906,284,"Fe 4",ipnt,'i',true,
977  " Fe 4 2.836 microns 9-6 transition" );
978 
979 
980  PntForLine(3892.,"Fe 5",&ipnt);
981  lindst(CoolHeavy.c3892,3892,"Fe 5",ipnt,'c',true,
982  " Fe 5 3892+3839" );
983 
984  linadd(CoolHeavy.c5177,0,"Fe 6",'c' ,
985  " all of 2G lines together first ");
986 
987 
988  PntForLine(5177.,"Fe 6",&ipnt);
989  lindst(CoolHeavy.c5177*0.354,5177,"Fe 6",ipnt,'c',true,
990  " Fe 6 5177, approximate correct " );
991 
992  linadd(fe.Fe7CoolTot,0,"Fe7c",'c' ,
993  " Fe7c 0 - total cooling due to n-level Fe 7 atom ");
994 
995  /* >>chng 04 nov 04, move to multi-level system */
996  for( ilo=0; ilo<NLFE7-1; ++ilo )
997  {
998  /* must not do 1-0 or 2-1, which are transferred lines */
999  for( ihi=MAX2(3,ilo+1); ihi<NLFE7; ++ihi )
1000  {
1001 
1002  PntForLine(fe.Fe7_wl[ihi][ilo],"Fe 7",&ipnt);
1003  lindst( fe.Fe7_emiss[ihi][ilo] , (realnum)(fe.Fe7_wl[ihi][ilo]+0.5) , "Fe 7",ipnt,'c',true,
1004  " emission in these lines" );
1005  }
1006  }
1007 # if 0
1008  PntForLine(5721.,"Fe 7",&ipnt);
1009  lindst( fe.Fe7_5721 , 5721 , "Fe 7",ipnt,'c',true,
1010  " " );
1011 
1012  PntForLine(6601.,"Fe 7",&ipnt);
1013  lindst( fe.Fe7_6601 , 6601 , "Fe 7",ipnt,'c',true,
1014  " " );
1015 
1016  PntForLine(3760.,"Fe 7",&ipnt);
1017  lindst( fe.Fe7_3760 , 3760 , "Fe 7",ipnt,'c',true,
1018  " " );
1019 
1020  PntForLine(3588.,"Fe 7",&ipnt);
1021  lindst( fe.Fe7_3588 , 3588 , "Fe 7",ipnt,'c',true,
1022  " " );
1023 # endif
1024 
1026  " [Fe 7] 9.51 micron ");
1027 
1029  " [Fe 7] 7.81 micron ");
1030 
1031  /* [Fe 7] 6087
1032  PntForLine(6087.,"Fe 7",&ipnt);
1033  lindst(CoolHeavy.c6087,6087,"Fe 7",ipnt,'c',true );*/
1034 
1035  /* [Fe 7] 5722
1036  PntForLine(5722.,"Fe 7",&ipnt);
1037  lindst(CoolHeavy.Fe5722,5722,"Fe 7",ipnt,'c',true );*/
1038 
1040  " Be-seq lines ");
1041 
1042  PntForLine(242.,"Fe 7",&ipnt);
1043  lindst(CoolHeavy.c242,242,"Fe 7",ipnt,'c',true,
1044  " Fe 9 242 j=1 slower decay");
1045 
1047  " the E1 transition that can pump [Fe X] ");
1048 
1049  /* optically thin Fe X pumping */
1050  eff = 69.4/(69.4 + 0.27*dense.cdsqte);
1051 
1052  /* coll excitation of 352 which decays to excited state of 6374
1053  * assumes 17/56 (ratio of A's) go to excited state */
1054  c10 = TauLines[ipT352].Emis->phots*eff*.01676*TauLines[ipT352].EnergyErg*352/6374.;
1055 
1056  /* Fe 10 and Fe 14 from Mason 75
1057  * total (coll, pumped)
1058  * call linadd( C6374+C10 , 6374 , 'Fe10','i') */
1059  PutExtra( c10 );
1060 
1061 
1063  " [Fe 10] 6375, collisions with pumping too ");
1064 
1065  /* collisional contribution
1066  * call linadd( C6374 , 6374 , 'Coll','c')
1067  * collisions of E1 line, plus pumped by continuum fluorescence
1068  * call linadd( C10 , 6374 , ' 352','c')
1069  * Fe XI 7892, 6.08 micron
1070  PutLine(&TauLines[ipTFe07]);*/
1071 
1083  /* Fe 11 2649 collisional excitation
1084  PntForLine(2649.,"Fe11",&ipnt);
1085  lindst(CoolHeavy.c2649,2649,"Fe11",ipnt,'c',true ); */
1086 
1087  /* Fe 11 1467 collisional excitation
1088  PntForLine(1467.,"Fe11",&ipnt);
1089  lindst(CoolHeavy.c1467,1467,"Fe11",ipnt,'c',true );*/
1090 
1091  linadd(fe.Fe11CoolTot,0,"Fe11",'c' ,
1092  " >>chng 05 dec 18, add Fe 11 Fe11 0 - total cooling due to 5-level Fe 11 atom ");
1093 
1094  /* Fe 11 5-level atom */
1095  for( ihi=1; ihi<NLFE11; ++ihi )
1096  {
1097  for( ilo=0; ilo<ihi; ++ilo )
1098  {
1099  PntForLine(fe.Fe11_wl[ihi][ilo],"Fe11",&ipnt);
1100  lindst( fe.Fe11_emiss[ihi][ilo] , (realnum)(fe.Fe11_wl[ihi][ilo]+0.5) , "Fe11",ipnt,'c',true,
1101  " emission in these lines" );
1102  }
1103  }
1104 
1105  PntForLine(1242.,"Fe12",&ipnt);
1106  lindst(CoolHeavy.c1242,1242,"Fe12",ipnt,'c',true ,
1107  " Fe 12, 1242, 1349 together, collisional excitation");
1108 
1109  PntForLine(2170.,"Fe12",&ipnt);
1110  lindst(CoolHeavy.c2170,2170,"Fe12",ipnt,'c',true ,
1111  " Fe 12, 2170, 2406 together, collisional excitation");
1112 
1113 
1114  PntForLine(2568.,"Fe12",&ipnt);
1115  lindst(CoolHeavy.c2568,2568,"Fe12",ipnt,'c',true,
1116  " Fe12 2904, 2567, 3567, 3073 together, collisional excitation" );
1117 
1118  /* >>chng 05 dec 18, add Fe 13 */
1119  /* Fe13 0 - total cooling due to 5-level Fe 13 atom */
1120  linadd(fe.Fe13CoolTot,0,"Fe13",'c' ,
1121  "total cooling due to Fe 13 model atom ");
1122 
1123  /* Fe 13 5-level atom */
1124  for( ihi=1; ihi<NLFE13; ++ihi )
1125  {
1126  for( ilo=0; ilo<ihi; ++ilo )
1127  {
1128 
1129  PntForLine(fe.Fe13_wl[ihi][ilo],"Fe13",&ipnt);
1130  lindst( fe.Fe13_emiss[ihi][ilo] , (realnum)(fe.Fe13_wl[ihi][ilo]+0.5) , "Fe13",ipnt,'c',true ,
1131  " Fe 13 emission");
1132  }
1133  }
1134 
1135  /* Fe 14 optically thin in line 344 */
1136  eff = 60.3/(60.3 + 0.23*dense.cdsqte/4.);
1137  r14 = dense.xIonDense[ipIRON][14-1]*fe.pfe14*rfield.flux[fe.ipfe10-1]*
1138  3.75e-12*eff/(1. + 24./.63*
1139  TauLines[ipT347].Emis->Pesc)*esc_PRD_1side(TauLines[ipT347].Emis->TauIn,1e-4);
1140  PntForLine(5303.,"Fe14",&ipnt);
1141 
1142  lindst(CoolHeavy.c5303+r14,5303,"Fe14",ipnt,'i',true ," total emission in Fe 14 5304");
1143 
1144  linadd(CoolHeavy.c5303,5303,"Coll",'c' ,
1145  " Fe 14 5304 contribution from collisional excitation ");
1146 
1147  linadd(r14,5303,"Pump",'r' ,
1148  " Fe 14 5304 continuum fluorescense ");
1149 
1152  c14 = 0.;
1153  linadd(c14,5303," 347",'c' ,
1154  " collisional excitation of E1 line ");
1155 
1157  " Fe 18 974.86A ");
1158 
1159 
1160  PntForLine(7047.,"Fe19",&ipnt);
1161  lindst(CoolHeavy.c7082,7047,"Fe19",ipnt,'c',true,
1162  " O-like Fe19, 3P ground term, 7046.72A vacuum wl, 1328.90A >>chng 01 aug 10, updated wavelengths Fe 19 7047 '85 " );
1163 
1164 
1165  PntForLine(1328.,"Fe19",&ipnt);
1166  lindst(CoolHeavy.c1328,1329,"Fe19",ipnt,'c',true,
1167  " Fe 19 1329" );
1168 
1169 
1170  PntForLine(592.,"Fe19",&ipnt);
1171  lindst(CoolHeavy.c592,592,"Fe19",ipnt,'c',true,
1172  " Fe 19 from loulergue et al '85" );
1173 
1174 
1175  PntForLine(1118.,"Fe19",&ipnt);
1176  lindst(CoolHeavy.c1118,1118,"Fe19",ipnt,'c',true ,
1177  " Fe 19 from loulergue et al '85");
1178 
1180  " next two 3p ground state lines, collisional excitation ");
1181 
1183  " collisional excitation ");
1184 
1186  " Fe20 721.40A, 578");
1188  " ");
1189 
1190  linadd(
1191  TauLines[ipFe22_247].Emis->xIntensity+
1196  260,"TOTL",'i',
1197  " Fe 22 845.6A total intensity of Fe22, all lines in the multiplet " );
1199  "Fe 22 247");
1201  "Fe 22 217");
1203  "Fe 22 348");
1205  "Fe 22 292");
1206  PutLine(&TauLines[ipFe22_253],
1207  "Fe 22 253");
1208 
1209  /* Fe 23 1909-like 262.6 */
1210  PntForLine(263.,"Fe23",&ipnt);
1211  lindst(CoolHeavy.c263,263,"Fe23",ipnt,'c',true,
1212  "Fe 23 1909-like 262.6" );
1213 
1214 
1216  " Fe 24 only 192 of 255, 192 Li seq doublet, collisional excitation ");
1217 
1219  " Fe 24 255 of 255, 192 Li seq doublet, collisional excitation ");
1220 
1222  " Fe 24 Li seq 2s3p collisional excitation ");
1223 
1224  /* recombination Ka */
1225  if( dense.lgElmtOn[ipIRON] )
1226  {
1227  /* these lines added to outlin in metdif - following must be false
1228  * fela = xLyaHeavy(nelem,nelem)*dense.xIonDense(nelem,nelem+1) */
1229  fela = Transitions[ipH_LIKE][ipIRON][ipH2p][ipH1s].Emis->xIntensity;
1230  }
1231  else
1232  {
1233  fela = 0.;
1234  }
1235 
1236  /* >>chng 02 jan 14, add grain fe to this sum */
1237  /* total intensity of K-alpha line */
1238  /*linadd((fe.fekcld+fe.fegrain)*1.03e-8+(fe.fekhot+fela)*1.11e-8,2,"FeKa",'i' );*/
1239  if( dense.lgElmtOn[ipIRON] )
1240  {
1241  lindst((fe.fekcld+fe.fegrain)*1.03e-8+(fe.fekhot+fela)*1.11e-8,2,"FeKa",
1242  Transitions[ipH_LIKE][ipIRON][ipH2p][ipH1s].ipCont,'i',false,
1243  "total intensity of K-alpha line" );
1244  }
1245 
1246  linadd(fela*1.11e-8,2,"FeLr",'i' ,
1247  " recombination from fully stripped ion ");
1248 
1249  /* >>chng 03 aug 14, label changed from TotH to AugH to be like rest total hot iron Ka; */
1250  linadd((fe.fekhot+fela)*1.11e-8,2,"AugH",'i' ,
1251  " Auger hot iron, assumes case b for H and He-like ");
1252 
1253  linadd(fe.fekcld*1.03e-8,2,"AugC",'i',
1254  " Auger production of cold iron, less than or 17 times ionized " );
1255 
1256  linadd(fe.fegrain*1.03e-8,2,"AugG",'i' ,
1257  " grain production of cold iron ");
1258 
1260  " [Co XI] 5168. A ");
1261 
1263  " nickel [Ni I] 7m ");
1264 
1265  /* nickel*/
1266 
1267 
1269  " [Ni I] 11m ");
1270 
1272  " [Ni XII] 4230.8 A ");
1273 
1274  /* copper */
1275 
1276  /* zinc */
1278  "zinc iv 3.625 microns, cs and A just made up ");
1279 
1280  if( trace.lgTrace )
1281  {
1282  fprintf( ioQQQ, " lines_lv1_k_zn returns\n" );
1283  }
1284  return;
1285 }

Generated for cloudy by doxygen 1.8.4