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())**/
303 rec *= dense.xIonDense[13][3]*dense.eden*emit_frac(TauLines[ipT1207])*
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
443 linadd(1.75e-22*dense.eden*dense.xIonDense[ipSULPHUR][1]/phycon.sqrte/
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);
452 lindst(CoolHeavy.S4070+CoolHeavy.S4078,4074,
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
679 linadd(CoolHeavy.Ar7237+CoolHeavy.Ar7331+CoolHeavy.Ar7171+CoolHeavy.Ar7263,7250,"TOTL",'i',
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
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
TransitionList TauLines("TauLines", &AnonStates)
t_trace trace
Definition trace.cpp:5
void PutExtra(double Extra)
void PutLine(const TransitionProxy &t, const char *chComment, const char *chLabelTemp)
double emit_frac(const TransitionProxy &t)