cloudy trunk
Loading...
Searching...
No Matches
prt_lines_lv1_na_ar.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_na_ar place lines of elements sodium through argon into lines storage stack */
4#include "cddefines.h"
5#include "coolheavy.h"
6#include "sil.h"
7#include "phycon.h"
8#include "embesq.h"
9#include "taulines.h"
10#include "dense.h"
11#include "ionbal.h"
12#include "trace.h"
13#include "lines_service.h"
14#include "lines.h"
15
17{
18 long int ipnt;
19 double drec,
20 fac,
21 rec,
22 sum,
23 t4;
24
25 DEBUG_ENTRY( "lines_lv1_na_ar()" );
26
27 if( trace.lgTrace )
28 {
29 fprintf( ioQQQ, " lines_lv1_na_ar called\n" );
30 }
31
32 t4 = phycon.te/1e4;
33
35 " sodium sum of Na D lines");
36
38 " [NaIII] 7.319 micron ");
39
41 " [NaIV] 9.048 micron ");
42
44 " [NaIV] 21.29 micron ");
45
46 PntForLine(1365.,"Na 5",&ipnt);
47 lindst(CoolHeavy.c1365,1365,"Na 5",ipnt,'t',true ,
48 " [NaV] 1365, sum of 1365.1+1365.8; cs only guess ");
49
50 PntForLine(2067.,"Na 5",&ipnt);
51 lindst(CoolHeavy.c2067,2067,"Na 5",ipnt,'t',true ,
52 " [NaV] 2067, sum of 2066.9+2068.4; cs only guess ");
53
54 PntForLine(4017.,"Na 5",&ipnt);
55 lindst(CoolHeavy.c4017,4017,"Na 5",ipnt,'t',true,
56 " [NaV] 4017, sum of 4010.9+4016.7+4022.7; cs only guess " );
57
58 PntForLine(2569.,"Na 6",&ipnt);
59 lindst(CoolHeavy.c2569,2569,"Na 6",ipnt,'t',true,
60 " [Na VI] 2568.9 " );
61
62 PntForLine(1357.,"Na 6",&ipnt);
63 lindst(CoolHeavy.c1357,1357,"Na 6",ipnt,'t',true ,
64 " [Na VI] 1356.6 ");
65
66 PntForLine(2972.,"Na 6",&ipnt);
67 lindst(CoolHeavy.c2972/(1.+1./3.02),2972,"Na 6",ipnt,'t',true ,
68 " [Na VI] 2971.9 ");
69
70 PntForLine(2872.,"Na 6",&ipnt);
71 lindst(CoolHeavy.c2972/(1.+3.02),2872,"Na 6",ipnt,'t',true,
72 " [Na VI] 2872.7 ");
73
75 " [NaVI] 14.32 micron ");
76
78 " [NaVI] 8.62 micron ");
79
81 " [NaVII] 4.675 micron ");
82
84 "Magnesium I 4571, O I data for coll strength and trans prob ");
85
87 " Mg I 2853 ");
88
90 " Mg I 2026 ");
91
92 linadd(TauLines[ipT2796].Emis().xIntensity()+TauLines[ipT2804].Emis().xIntensity(),2798,
93 "TOTL",'i',"Mg II 2798 add both lines of multiplet together " );
94
95 /* sum of inward fracs of lines */
96 sum = TauLines[ipT2796].Emis().xIntensity()*TauLines[ipT2796].Emis().FracInwd() +
97 TauLines[ipT2804].Emis().xIntensity()*TauLines[ipT2804].Emis().FracInwd();
98
99 linadd(sum,2798,"Inwd",'i',
100 "inward part of Mg II 2798" );
101
103 "one member of Mg II multiplet");
104
106 "one member of Mg II multiplet");
107
109 " Mg IV 4.5 micron ");
110
112 " Mg V 13.5 micron emission");
113
115 " Mg V 5.6 micron emission ");
116
118 " [Mg 5] 2571, 2893 ");
119
121 " [Mg 5] 2417.5A, 3-2 in model atom" );
122
124 " [Mg 5] 1324.58A, 3-1 in model atom ");
125
126 PntForLine(1806.,"Mg 6",&ipnt);
127 lindst(CoolHeavy.c1806,1806,"Mg 6",ipnt,'t',true ,
128 " MG VI");
129
131 " [Mg 7] 2510, 2629, 2-1 transitions, together");
132
134 " [Mg 7] 3-2 transition, 2261.5 ");
135
137 " [Mg 7] 3-1 transition, 1189.82A ");
138
140 " Mg 7 IR line 5.50 microns ");
141
143 " Mg 7 IR line 9.03 microns ");
144
146 " [Mg 8] 3.03 micron ");
147
148 PutLine(TauLines[ipT705]," Mg 9 704.5 ");
149
150 linadd(TauLines[ipTMg610].Emis().xIntensity()+TauLines[ipTMg625].Emis().xIntensity(),615,"TOTL",'i',
151 " Mg 10 614.9 both of doublet, li seq 2s 2p" );
152
154 "");
156 "");
157
159 " part of Mg 10 destroyed by background opacity Mg 10 58 li seq 2s 3p ");
160 linadd(TauLines[ipTMg610].Emis().ots()*TauLines[ipTMg610].EnergyErg()+
161 TauLines[ipTMg625].Emis().ots()*TauLines[ipTMg625].EnergyErg(),615,"dest",'i',
162 "" );
163
165 " Aluminum Al I 3957 ");
166
168 " Al I 3090 ");
169
170 linadd(embesq.em2669+TauLines[ipT2670].Emis().xIntensity(),2665,"totl",'i',
171 "Al II 1671 total emission in Al II] 2669.7, 2660 doublet" );
173 " ");
174
175 linadd(embesq.em2669,2660,"Al 2",'i',
176 "emission in Al II] 2669 alone" );
177
178 linadd(TauLines[ipT1855].Emis().xIntensity()+TauLines[ipT1863].Emis().xIntensity(),1860,"TOTL",'i',
179 " Al III" );
180 sum = TauLines[ipT1855].Emis().xIntensity()*TauLines[ipT1855].Emis().FracInwd() +
181 TauLines[ipT1863].Emis().xIntensity()* TauLines[ipT1863].Emis().FracInwd();
182
183 linadd(sum,1860,"Inwd",'i',
184 " inward part of AlIII line" );
186 "");
188 "");
189
191 " [Al V] 2.905 micron ");
192
194 " [Al VI] 3.66 micron ");
195
197 " [Al VI] 9.12 micron");
198
199 PntForLine(2428.,"Al 6",&ipnt);
200 lindst(CoolHeavy.c2428/(1.+1./3.73),2428,"Al 6",ipnt,'t',true,
201 " [Al VI] 2428.4 " );
202
203 PntForLine(2601.,"Al 6",&ipnt);
204 lindst(CoolHeavy.c2428/(1.+3.73),2601,"Al 6",ipnt,'t',true,
205 " [Al VI] 2601.0");
206
207 PntForLine(1170.,"Al 6",&ipnt);
208 lindst(CoolHeavy.c1170,1170,"Al 6",ipnt,'t',true,
209 " [Al VI] 1169.86 ");
210
211 PntForLine(2125.,"Al 6",&ipnt);
212 lindst(CoolHeavy.c2125,2125,"Al 6",ipnt,'t',true ,
213 " [Al VI] 2124.95" );
214
216 " [Al VIII] 5.75 micron " );
217
219 " [Al VIII] 3.70 micron");
220
222 " [Al IX] 2.04 micron, no collision strength, A NIST ");
223
224 PutLine(TauLines[ipT639]," Al X ");
225
226 linadd(TauLines[ipTAl550].Emis().xIntensity()+TauLines[ipTAl568].Emis().xIntensity(),556,"TOTL",'i',
227 " Al 11, Li seq 2s2p" );
228
230 "");
231
233 "");
234
236 " Al 11, Li seq 2s3p ");
237
239 " silicon Silicon Si I 130m ");
240
242 " Si I 68m");
243
245 " Si I 2518A ");
246
248 " Si I 2215A ");
249
251 " Silicon II 35 micron ");
252
253 linadd(
254 TauLines[ipSi2_2334].Emis().xIntensity()+
255 TauLines[ipSi2_2329].Emis().xIntensity()+
256 TauLines[ipSi2_2350].Emis().xIntensity()+
257 TauLines[ipSi2_2344].Emis().xIntensity()+
258 TauLines[ipSi2_2336].Emis().xIntensity(),
259 2335,"TOTL",'i',
260 "total intensity of S IV] 1406, all lines in the multiplet" );
262 " ");
264 " ");
266 " ");
268 " ");
270 " ");
271
273 " SI II 1808, permitted resonance line, collisionally excited ");
274
276 " SI II 1527, permitted resonance line, collisionally excited ");
277
279 " SI II 1305, permitted resonance line, collisionally excited ");
280
282 " SI II 1260, permitted resonance line, collisionally excited ");
283
284 /* SI II 1260, rough guess of dielec contribution */
285 drec = dense.xIonDense[13][2]*dense.eden*7.6e-7/phycon.te32*
286 1.57e-11;
287
288 fac = emit_frac(TauLines[ipT1260]);
289 PntForLine(1260.,"Si 2",&ipnt);
290 lindst(drec*fac,1260,"diel",ipnt,'r',true,
291 " fac = (1.-TauLines[ipT1260].ColOvTot());" );
292
293
294 PntForLine(1909.,"Si 2",&ipnt);
295 lindst(drec*1260./1909.,1909,"diel",ipnt,'r',true,
296 " dielectronic recombination SiII 1909" );
297
298 rec = 1e-12*(0.1152/t4 - 0.3082 + 4.4734*t4 + 0.0207*t4*t4)/pow(t4,1.5)*
299 sexp(0.2981/t4);
300
301 /* >>chng 96 jul 8, added dielectronic recombination contribution to 1207 */
302 /*rec *= dense.xIonDense[13][3]*dense.eden*(1.-TauLines[ipT1207].ColOvTot())**/
304 1.65e-11;
305
306 PutExtra(MAX2(0.,rec));
307
309 " SI III 1207, collisional excitation and dielectronic recombination ");
310
311 linadd(MAX2(0.,rec),1207,"rec ",'i',
312 " Si III 1207, dielectronic recombination only" );
313
314 linadd(embesq.em1895+TauLines[ipT1895].Emis().xIntensity(),1888,"TOTL",'i',
315 " Si III] 1892+1883, total intensity of both lines" );
317 " ");
318
319 PntForLine(1883.,"Si 3",&ipnt);
320 lindst(embesq.em1895,1883,"Si 3",ipnt,'t',true ,
321 " Si III] 1883 by itself");
322
323 /*fac = (1.-TauLines[ipT1895].ColOvTot());*/
324 fac = emit_frac(TauLines[ipT1895]);
325
326 double p1895 = ionbal.PhotoRate_Shell[ipSILICON][1][2][0]*
327 dense.xIonDense[ipSILICON][1]*0.85;
328 linadd( p1895*1.05e-11*fac,1895,"PHOT",'i',
329 " photoproduction by inner shell removal" );
330
331 linadd(TauLines[ipT1403].Emis().xIntensity()+TauLines[ipT1394].Emis().xIntensity(),1397,"TOTL",'i',
332 " Si IV 1397, collisionally excited " );
333
334 sum = TauLines[ipT1403].Emis().xIntensity()*TauLines[ipT1403].Emis().FracInwd() +
335 TauLines[ipT1394].Emis().xIntensity()* TauLines[ipT1394].Emis().FracInwd();
336
337 linadd(sum,1397,"Inwd",'i',
338 " inward part of SiIV 1397" );
340 " ");
342 " ");
343
345 " SI VI 1.9641 micron ");
346
347 PntForLine(2148.,"Si 7",&ipnt);
348 lindst(sil.c2148,2148,"Si 7",ipnt,'t',true,
349 " SI VII, 2148, O III like, collisionally excited" );
350
352 " Si VII 2.48, 6.49 micron, collisionally excited ");
353
355 " Si VII 2.48, 6.49 micron, collisionally excited ");
356
357 PntForLine(1446.,"Si 8",&ipnt);
358 lindst(sil.c1446,1446,"Si 8",ipnt,'t',true,
359 " SI VIII 1446, OIII like, collisionally excited" );
360
361 PntForLine(1985.,"Si 9",&ipnt);
362 lindst(sil.c1985,1985,"Si 9",ipnt,'t',true,
363 " SI IX 1985, 2150, collisionally excited" );
364
365
366 PntForLine(949.,"Si 9",&ipnt);
367 lindst(sil.c949,949,"Si 9",ipnt,'t',true,
368 " collisionally excited" );
369
370 PntForLine(1815.,"Si 9",&ipnt);
371 lindst(sil.c1815,1815,"Si 9",ipnt,'t',true,
372 " collisionally excited " );
373
375 " SI 9, 3.86, 2.84 3P fine structure lines ");
376
378 " SI 9, 3.86, 2.84 3P fine structure lines ");
379
380 PntForLine(691.,"Si 9",&ipnt);
381 lindst(sil.c691,691,"Si 9",ipnt,'t',true ,
382 " both components of 5S-3P doublet");
383
385 " SI 10 606A, actually group of 4 intercombination lines ");
386
388 " [Si 10] 1.43 micron, collisionally excited ");
389
390 PntForLine(583.,"Si11",&ipnt);
391 lindst(sil.c583,581,"Si11",ipnt,'t',true,
392 " Si 11 582.9, collisionally excited >>chng 01 may 23, wavelength from 583 to 581" );
393
394 linadd(TauLines[ipTSi499].Emis().xIntensity()+TauLines[ipTSi521].Emis().xIntensity(),506,"TOTL",'i' ,
395 "emission total Si 12 506 + 499 ");
396
398 " Si 12 506 li seq 2s 2p ");
400 " ");
401
403 " Si 12 40.9A, li seq 2s 3p ");
404
406 " phosphorus [P II] 60.64 micron ");
407
409 " [P II] 32.87 micron ");
410
411 PntForLine(16400.,"P 2",&ipnt);
412 lindst(CoolHeavy.p2_32,16400,"P 2",ipnt,'t',true ,
413 " 3-2 1.64 micron ");
414
415 PntForLine(4669.,"P 2",&ipnt);
416 lindst(CoolHeavy.p2_31*0.75,4669,"P 2",ipnt,'t',true ,
417 " >>chng 01 may 15, add these lines P 2 3-1 4670, 4738 vac wl, 4669, 4737 air");
418
419 PntForLine(4737.,"P 2",&ipnt);
420 lindst(CoolHeavy.p2_31*0.25,4737,"P 2",ipnt,'t',true ,
421 " P 2 3-1 4670, 4738 vac wl, 4737, 4737 air");
422
423 PntForLine(11890.,"P 2",&ipnt);
424 lindst(CoolHeavy.p2_21*0.75,11890,"P 2",ipnt,'t',true,
425 " 2-1 1.147, 1.189 micron" );
426
427 PntForLine(11470.,"P 2",&ipnt);
428 lindst(CoolHeavy.p2_21*0.25,11470,"P 2",ipnt,'t',true ,
429 " [P II] 1.14 micron");
430
432 " [P III] 17.885 micron ");
433
435 " [P VII] 1.3745 micron ");
436
438 " sulphur S I 25m ");
439
441 " S I 56m ");
442
444 phycon.te10/phycon.te03,1807,"S 1R",'i',
445 " guesstimate of Sulphur I triplet excited state recombination rate.this is to check whether photoexcit of S II is ever important S I 1807 recombination " );
446
447 PntForLine(6731.,"S II",&ipnt);
448 lindst(CoolHeavy.c6731,6720,"S 2",ipnt,'t',false,
449 " S II 6731 + 6716 together " );
450
451 PntForLine(4070.,"S II",&ipnt);
453 "S 2",ipnt,'i',false,"S II 4070 +4078 together" );
454
455 PntForLine(10330.,"S 2",&ipnt);
456 lindst(CoolHeavy.c10330,10330,"S 2",ipnt,'t',true,
457 " S II N=3 lines, all four lines together " );
458
459 PntForLine(6731.,"S II",&ipnt);
460 lindst(CoolHeavy.S6733,6731,"S II",ipnt,'t',true,
461 " individual line from five level atom" );
462
463 PntForLine(6716.,"S II",&ipnt);
464 lindst(CoolHeavy.S6718,6716,"S II",ipnt,'t',true,
465 " individual line from five level atom" );
466
467 PntForLine(4070.,"S II",&ipnt);
468 lindst(CoolHeavy.S4070,4070,"S II",ipnt,'t',true,
469 " individual line from five level atom " );
470
471 PntForLine(4078.,"S II",&ipnt);
472 lindst(CoolHeavy.S4078,4078,"S II",ipnt,'t',true,
473 " individual line from five level atom" );
474
475 PntForLine(10330.,"S 2",&ipnt);
476 lindst(CoolHeavy.S10323,10323,"S II",ipnt,'t',false,
477 " individual line from five level atom " );
478
479 lindst(CoolHeavy.S10289,10289,"S II",ipnt,'t',false,
480 " individual line from five level atom" );
481
482 lindst(CoolHeavy.S10373,10373,"S II",ipnt,'t' ,false,
483 " individual line from five level atom ");
484
485 lindst(CoolHeavy.S10339,10339,"S II",ipnt,'t',false,
486 " individual line from five level atom " );
487
489 " resonance line near NV, collisionally excited ");
490
492 " S III fine structure 18.7 ");
493
495 " S III fine structure 34 ");
496
498 " S III] 1713.12, 1728.94 ");
499
500 PntForLine(9532.,"S 3",&ipnt);
501 lindst(CoolHeavy.c9532/(1.+1./2.48),9532,"S 3",ipnt,'t',true ,
502 " [S III] 9532 alone ");
503
504 PntForLine(9069.,"S 3",&ipnt);
505 lindst(CoolHeavy.c9532/(1.+2.48),9069,"S 3",ipnt,'t',true,
506 " [S III] 9069 alone" );
507
508 PntForLine(6312.,"S 3",&ipnt);
509 lindst(CoolHeavy.c6312,6312,"S 3",ipnt,'t',true ,
510 " [S III] 6312, trans-auroral temperature sensitive ");
511
512 PntForLine(3722.,"S 3",&ipnt);
513 lindst(CoolHeavy.c6312*0.59,3722,"S 3",ipnt,'t',true,
514 " [S III] 3722, same upper level as 6312" );
515
517 " WL, other data, from Ho + Henry Ap.J. 1984 ");
518
520 " S IV 10.5 micron, collisionally excited (label is 105) ");
521
522 linadd(
523 TauLines[ipS4_1405].Emis().xIntensity()+
524 TauLines[ipS4_1398].Emis().xIntensity()+
525 TauLines[ipS4_1424].Emis().xIntensity()+
526 TauLines[ipS4_1417].Emis().xIntensity()+
527 TauLines[ipS4_1407].Emis().xIntensity(),
528 1406,"TOTL",'i',
529 " total intensity of S IV] 1406, all lines in the multiplet " );
531 " ");
533 " ");
535 " ");
537 " ");
539 " ");
540
541 linadd(embesq.em1198+TauLines[ipT1198].Emis().xIntensity(),1198,"TOTL",'i',
542 " S V 1198] both lines together " );
543
545 " S V 1198] the stronger transition ");
546
547 linadd(embesq.em1198,1188,"S 5",'i',
548 " Be seq, weaker of the two transitions" );
549
551 " S V 786.5, collisionally excited ");
552
554 " chlorine lines [Cl I] 11 micron ");
555
557 " [Cl II] 33.281 micron ");
558
560 " [Cl II] 14.3678 micron ");
561
562 PntForLine(8578.7,"Cl 2",&ipnt);
563 lindst(CoolHeavy.c8579*0.791,8579,"Cl 2",ipnt,'t',true,
564 " Chlorine II 8578.7, 9123.6 doublet");
565
566 PntForLine(9123.6,"Cl 2",&ipnt);
567 lindst(CoolHeavy.c8579*0.209,9124,"Cl 2",ipnt,'t',true,
568 " Chlorine II 8578.7, 9123.6 doublet" );
569
570 PntForLine(6161.8,"Cl 2",&ipnt);
571 lindst(CoolHeavy.c6164,6162,"Cl 2",ipnt,'t',true,
572 " Chlorine II 6161.8 auroral line >>chng 03 feb 24, change wavelength from 6164 to correct 6161.8 " );
573
574 PntForLine(3677.9,"Cl 2",&ipnt);
575 lindst(CoolHeavy.c3679,3678,"Cl 2",ipnt,'t',true,
576 " Chlorine II 3677.9 auroral line >>chng 03 feb 24, to correct wavelength " );
577
578 linadd(CoolHeavy.c5525,5525,"TOTL",'i',
579 " Cl III 5519, 5539 doublet, both together " );
580
581 linadd(CoolHeavy.c3350,3350,"TOTL",'i',
582 " Cl III 3354, 3344 doublet, both together " );
583
584 linadd(CoolHeavy.c8494,8494,"TOTL",'i',
585 " Cl III 8504, 8436, 8552, 8483 multiplet, all together " );
586
587 PntForLine(5538.,"Cl 3",&ipnt);
588 lindst(CoolHeavy.Cl5539,5538,"Cl 3",ipnt,'t',true,
589 " Cl III 5538 " );
590
591 PntForLine(5518.,"Cl 3",&ipnt);
592 lindst(CoolHeavy.Cl5519,5518,"Cl 3",ipnt,'t',true,
593 " Cl III 5518" );
594
595 PntForLine(3354.,"Cl 3",&ipnt);
596 lindst(CoolHeavy.Cl3354,3354,"Cl 3",ipnt,'t',true,
597 " Cl III 3354 " );
598
599 PntForLine(3344.,"Cl 3",&ipnt);
600 lindst(CoolHeavy.Cl3344,3344,"Cl 3",ipnt,'t',true,
601 " Cl III 3344 " );
602
603 PntForLine(8504.,"Cl 3",&ipnt);
604 lindst(CoolHeavy.Cl8504,8504,"Cl 3",ipnt,'t',true ,
605 " Cl III 8504 ");
606
607 PntForLine(8436.,"Cl 3",&ipnt);
608 lindst(CoolHeavy.Cl8436,8436,"Cl 3",ipnt,'t',true,
609 " Cl III 8436" );
610
611 PntForLine(8552.,"Cl 3",&ipnt);
612 lindst(CoolHeavy.Cl8552,8552,"Cl 3",ipnt,'t',true,
613 " Cl III 8552 " );
614
615 PntForLine(8483.,"Cl 3",&ipnt);
616 lindst(CoolHeavy.Cl8483,8483,"Cl 3",ipnt,'t',true,
617 " Cl III 8483" );
618
620 " [Cl IV] fine structure line 20.354 microns");
621
623 " [Cl IV] fine structure line 11.741 microns ");
624
625 PntForLine(8047.,"Cl 4",&ipnt);
626 lindst(CoolHeavy.c8047*0.667,8047,"Cl 4",ipnt,'t',true,
627 " ClIV 8047" );
628
629 PntForLine(7532.,"Cl 4",&ipnt);
630 lindst(CoolHeavy.c8047*0.333,7532,"Cl 4",ipnt,'t',true,
631 " ClIV 7532" );
632
633 PntForLine(3119.,"Cl 4",&ipnt);
634 lindst(CoolHeavy.c3119,3119,"Cl 4",ipnt,'t',true,
635 " ClIV 3119" );
636
637 PntForLine(5324.,"Cl 4",&ipnt);
638 lindst(CoolHeavy.c5324,5324,"Cl 4",ipnt,'t',true,
639 " ClIV 5324" );
640
642 " Cl IX 7334A ");
643
645 " Argon II 7 micron ");
646
647 PntForLine(7135.,"Ar 3",&ipnt);
648 lindst(CoolHeavy.c7136/(1.+1./4.144),7135,"Ar 3",ipnt,'t',true ,
649 " Argon III 7135");
650
651 PntForLine(7751.,"Ar 3",&ipnt);
652 lindst(CoolHeavy.c7136/(1.+4.144),7751,"Ar 3",ipnt,'t',true,
653 " Argon III 7751" );
654
655 PntForLine(5192.,"Ar 3",&ipnt);
656 lindst(CoolHeavy.c5192,5192,"Ar 3",ipnt,'t',true,
657 " Argon III 5192" );
658
659 PntForLine(3109.,"Ar 3",&ipnt);
660 lindst(CoolHeavy.c3109*0.9894,3109,"Ar 3",ipnt,'t',true,
661 " Argon III 3109" );
662
663 PntForLine(3005.,"Ar 3",&ipnt);
664 lindst(CoolHeavy.c3109*(1.-0.9894),3005,"Ar 3",ipnt,'t',true ,
665 " Argon III 3005 ");
666
668 " Argon III 21.8, 9 micron lines");
669
671 " Argon III 21.8, 9 micron lines ");
672
673 linadd(CoolHeavy.Ar4740+CoolHeavy.Ar4711,4725,"TOTL",'i',
674 " Argon IV 4711 + 4740 together, 4740=90%" );
675
676 linadd(CoolHeavy.Ar2868+CoolHeavy.Ar2854,2860,"TOTL",'i',
677 " [AvIV] 2868, 2854 together " );
678
680 " [AvIV] auroral lines, 7237, 7331, 7171, 7263 " );
681
682 PntForLine(4740.,"Ar 4",&ipnt);
683 lindst(CoolHeavy.Ar4740,4740,"Ar 4",ipnt,'t',true,
684 " [Ar IV] 4740" );
685
686 PntForLine(4711.,"Ar 4",&ipnt);
687 lindst(CoolHeavy.Ar4711,4711,"Ar 4",ipnt,'t',true,
688 " [Ar IV] 4711" );
689
690 PntForLine(2868.,"Ar 4",&ipnt);
691 lindst(CoolHeavy.Ar2868,2868,"Ar 4",ipnt,'t',true,
692 " [Ar IV] 2868" );
693
694 PntForLine(2854.,"Ar 4",&ipnt);
695 lindst(CoolHeavy.Ar2854,2854,"Ar 4",ipnt,'t',true,
696 " [Ar IV] 2854" );
697
698
699 PntForLine(7263.,"Ar 4",&ipnt);
700 lindst(CoolHeavy.Ar7263,7263,"Ar 4",ipnt,'t',true,
701 " [Ar IV] 7263" );
702
703 PntForLine(7171.,"Ar 4",&ipnt);
704 lindst(CoolHeavy.Ar7171,7171,"Ar 4",ipnt,'t',true,
705 " [Ar IV] 7171" );
706
707 PntForLine(7331.,"Ar 4",&ipnt);
708 lindst(CoolHeavy.Ar7331,7331,"Ar 4",ipnt,'t',true,
709 " [Ar IV] 7331" );
710
711 PntForLine(7237.,"Ar 4",&ipnt);
712 lindst(CoolHeavy.Ar7237,7237,"Ar 4",ipnt,'t',true,
713 " [Ar IV] 7237" );
714
715 PntForLine(7005.,"Ar 5",&ipnt);
716 lindst(CoolHeavy.c7007/(1.+1./2.143),7005,"Ar 5",ipnt,'t',true,
717 " Argon V, 3P lines, 7005, collisionally excited" );
718
719 PntForLine(6435.,"Ar 5",&ipnt);
720 lindst(CoolHeavy.c7007/(1.+2.143),6435,"Ar 5",ipnt,'t',true ,
721 " Argon V, 3P lines, 6435, collisionally excited");
722
723 PntForLine(4626.,"Ar 5",&ipnt);
724 lindst( CoolHeavy.c4626 , 4626 ,"Ar 5",ipnt,'t',true,
725 " >>chng 01 mar 10, add following two lines Argon V, 4626" );
726
727 PntForLine(2691.,"Ar 5",&ipnt);
728 lindst( CoolHeavy.c2691 , 2691 ,"Ar 5",ipnt,'t',true ,
729 " Argon V, 2691");
730
732 " Argon V fine structure lines, 13.09, 7.903 micron line ");
733
735 " Argon V fine structure lines, 13.09, 7.903 micron line ");
736
738 " [Ar VI] 4.53 micron ");
739
740 if( trace.lgTrace )
741 {
742 fprintf( ioQQQ, " lines_lv1_na_ar returns\n" );
743 }
744 return;
745}
long ipxMg51325
long ipTAr13
long ipP0233
long ipTMg610
long ipCl1_11m
long ipTSi25
long ipTAl48
long ipAl6912
long ipTS34
long ipTAr22
long ipSi2_2334
long ipMgI2026
long ipfsMg790
long ipxNa0746
long ipSi10143
long ipAl8370
long ipAr06453
long ipAl8575
long ipTAl568
long ipAl09204
long ipSi2_2336
long ipCl973
long ipTSi65
long ipSi2_2329
long ipT1194
long ipAlI3957
long ipxMg52855
long ipT1394
long ipT1207
long ipSi1_68m
long ipTSi41
long ipCl04203
long ipT786
long ipT1305
long ipTS19
long ipP713
long ipAl6366
long ipT1863
long ipS1_25m
long ipTMg625
long ipTMg14
long ipTAl550
long ipT5895
long ipxNa6143
long ipSi2_2350
long ipAl529
long ipP0318
long ipfsNa373
long ipfsMg755
long ipT639
long ipxMg52417
long ipP0260
long ipTAr8
long ipxMg71190
long ipSii2215
long ipT1256
long ipT58
long ipT2804
long ipfsCl214
long ipxMg08303
long ipT4561
long ipT1527
long ipT1198
long ipMgI2853
long ipfsNa421
long ipxMg72261
long ipTAr7
long ipTS11
long ipfsCl233
long ipT705
long ipTMg4
long ipT1855
long ipS4_1407
long ipT1403
long ipCl04117
long ipTS1720
long ipTSi521
long ipTSi4
long ipTAr9
long ipT2796
long ipS1_56m
long ipT1260
long ipSii2518
long ipSi619
long ipSi2_2344
long ipAlI3090
long ipS4_1398
long ipT1895
long ipSi10_606
long ipT2670
long ipS4_1405
long ipTSi35
long ipTSi3
long ipxMg72569
long ipfsNa490
long ipS4_1424
long ipTMg6
long ipSi1_130m
long ipxNa6862
long ipTSi499
long ipS4_1417
long ipT1808
FILE * ioQQQ
Definition: cddefines.cpp:7
sys_float sexp(sys_float x)
Definition: service.cpp:914
const int ipSILICON
Definition: cddefines.h:318
#define MAX2
Definition: cddefines.h:782
const int ipSULPHUR
Definition: cddefines.h:320
#define DEBUG_ENTRY(funcname)
Definition: cddefines.h:684
EmissionList & Emis()
Definition: transition.h:329
double **** PhotoRate_Shell
Definition: ionbal.h:111
t_CoolHeavy CoolHeavy
Definition: coolheavy.cpp:5
t_dense dense
Definition: dense.cpp:24
t_embesq embesq
Definition: embesq.cpp:5
t_ionbal ionbal
Definition: ionbal.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 lines_lv1_na_ar(void)
t_sil sil
Definition: sil.cpp:5
double Cl5539
Definition: coolheavy.h:32
double p2_21
Definition: coolheavy.h:100
double Ar7263
Definition: coolheavy.h:22
double Ar7171
Definition: coolheavy.h:23
double Ar4711
Definition: coolheavy.h:19
double Ar7237
Definition: coolheavy.h:25
double Cl3344
Definition: coolheavy.h:35
double c10330
Definition: coolheavy.h:86
realnum S4070
Definition: coolheavy.h:105
double c5525
Definition: coolheavy.h:29
double c3109
Definition: coolheavy.h:17
double c4626
Definition: coolheavy.h:14
realnum S4078
Definition: coolheavy.h:106
double c6164
Definition: coolheavy.h:27
double c3679
Definition: coolheavy.h:28
realnum S10339
Definition: coolheavy.h:110
double c8047
Definition: coolheavy.h:40
double c2691
Definition: coolheavy.h:15
realnum S6733
Definition: coolheavy.h:103
double c1357
Definition: coolheavy.h:55
double c7136
Definition: coolheavy.h:12
double c6312
Definition: coolheavy.h:87
double c7007
Definition: coolheavy.h:13
double c3119
Definition: coolheavy.h:41
double Ar2868
Definition: coolheavy.h:20
double c2972
Definition: coolheavy.h:56
double Cl3354
Definition: coolheavy.h:34
double p2_31
Definition: coolheavy.h:102
realnum S10289
Definition: coolheavy.h:108
double c2428
Definition: coolheavy.h:10
double c1806
Definition: coolheavy.h:53
double Cl8552
Definition: coolheavy.h:38
double Cl8436
Definition: coolheavy.h:37
double c8579
Definition: coolheavy.h:26
double Cl8504
Definition: coolheavy.h:36
realnum S10323
Definition: coolheavy.h:107
double Ar2854
Definition: coolheavy.h:21
double c6731
Definition: coolheavy.h:85
realnum S10373
Definition: coolheavy.h:109
double Ar7331
Definition: coolheavy.h:24
double c2569
Definition: coolheavy.h:54
double Cl5519
Definition: coolheavy.h:33
double c5192
Definition: coolheavy.h:16
realnum S6718
Definition: coolheavy.h:104
double p2_32
Definition: coolheavy.h:101
double c3350
Definition: coolheavy.h:30
double c5324
Definition: coolheavy.h:42
double c2067
Definition: coolheavy.h:58
double c8494
Definition: coolheavy.h:31
double c9532
Definition: coolheavy.h:88
double c1170
Definition: coolheavy.h:9
double Ar4740
Definition: coolheavy.h:18
double Cl8483
Definition: coolheavy.h:39
double c1365
Definition: coolheavy.h:57
double c4017
Definition: coolheavy.h:59
double c2125
Definition: coolheavy.h:11
double eden
Definition: dense.h:190
double xIonDense[LIMELM][LIMELM+1]
Definition: dense.h:125
realnum em1198
Definition: embesq.h:16
realnum em1895
Definition: embesq.h:15
realnum em2669
Definition: embesq.h:14
double te03
Definition: phycon.h:59
double te
Definition: phycon.h:11
double sqrte
Definition: phycon.h:48
double te32
Definition: phycon.h:49
double te10
Definition: phycon.h:55
double c1985
Definition: sil.h:12
double c2148
Definition: sil.h:16
double c949
Definition: sil.h:10
double c1815
Definition: sil.h:11
double c1446
Definition: sil.h:15
double c583
Definition: sil.h:14
double c691
Definition: sil.h:13
bool lgTrace
Definition: trace.h:12
TransitionList TauLines("TauLines", &AnonStates)
t_trace trace
Definition: trace.cpp:5
void PutExtra(double Extra)
Definition: transition.cpp:468
void PutLine(const TransitionProxy &t, const char *chComment, const char *chLabelTemp)
Definition: transition.cpp:449
double emit_frac(const TransitionProxy &t)
Definition: transition.cpp:84