cloudy  trunk
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 
16 void lines_lv1_na_ar(void)
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 }
t_CoolHeavy::p2_32
double p2_32
Definition: coolheavy.h:101
t_CoolHeavy::c10330
double c10330
Definition: coolheavy.h:86
t_CoolHeavy::c1170
double c1170
Definition: coolheavy.h:9
ipS4_1405
long ipS4_1405
Definition: atmdat_readin.cpp:90
ipfsNa490
long ipfsNa490
Definition: atmdat_readin.cpp:57
ipTS19
long ipTS19
Definition: atmdat_readin.cpp:77
ipP713
long ipP713
Definition: atmdat_readin.cpp:75
emit_frac
double emit_frac(const TransitionProxy &t)
Definition: transition.cpp:84
t_CoolHeavy::c1365
double c1365
Definition: coolheavy.h:57
t_CoolHeavy::Ar4711
double Ar4711
Definition: coolheavy.h:19
ipAl529
long ipAl529
Definition: atmdat_readin.cpp:66
lines.h
PntForLine
void PntForLine(double wavelength, const char *chLabel, long int *ipnt)
Definition: lines_service.cpp:583
t_CoolHeavy::Cl8552
double Cl8552
Definition: coolheavy.h:38
t_ionbal::PhotoRate_Shell
double **** PhotoRate_Shell
Definition: ionbal.h:111
ipAl09204
long ipAl09204
Definition: atmdat_readin.cpp:67
ipT1808
long ipT1808
Definition: atmdat_readin.cpp:69
t_dense::eden
double eden
Definition: dense.h:190
ipSi2_2350
long ipSi2_2350
Definition: atmdat_readin.cpp:94
t_CoolHeavy::Cl3354
double Cl3354
Definition: coolheavy.h:34
t_CoolHeavy::c6731
double c6731
Definition: coolheavy.h:85
t_CoolHeavy::c5192
double c5192
Definition: coolheavy.h:16
ipS4_1407
long ipS4_1407
Definition: atmdat_readin.cpp:90
t_sil::c583
double c583
Definition: sil.h:14
t_CoolHeavy::c3679
double c3679
Definition: coolheavy.h:28
t_CoolHeavy::c2125
double c2125
Definition: coolheavy.h:11
dense
t_dense dense
Definition: dense.cpp:24
ipTS11
long ipTS11
Definition: atmdat_readin.cpp:77
ipxMg52417
long ipxMg52417
Definition: atmdat_readin.cpp:60
t_CoolHeavy::Cl5539
double Cl5539
Definition: coolheavy.h:32
t_CoolHeavy::Ar2854
double Ar2854
Definition: coolheavy.h:21
PutExtra
void PutExtra(double Extra)
Definition: transition.cpp:468
ipTSi25
long ipTSi25
Definition: atmdat_readin.cpp:73
ioQQQ
FILE * ioQQQ
Definition: cddefines.cpp:7
t_CoolHeavy::c6312
double c6312
Definition: coolheavy.h:87
ipAlI3090
long ipAlI3090
Definition: atmdat_readin.cpp:64
ipCl973
long ipCl973
Definition: atmdat_readin.cpp:78
ipSii2215
long ipSii2215
Definition: atmdat_readin.cpp:69
t_CoolHeavy::c8494
double c8494
Definition: coolheavy.h:31
t_embesq::em1895
realnum em1895
Definition: embesq.h:15
ipSii2518
long ipSii2518
Definition: atmdat_readin.cpp:68
ipTSi41
long ipTSi41
Definition: atmdat_readin.cpp:72
ipTSi499
long ipTSi499
Definition: atmdat_readin.cpp:72
ipT2804
long ipT2804
Definition: atmdat_readin.cpp:59
ipfsMg755
long ipfsMg755
Definition: atmdat_readin.cpp:64
ipS4_1398
long ipS4_1398
Definition: atmdat_readin.cpp:90
t_phycon::te03
double te03
Definition: phycon.h:59
ipSi2_2336
long ipSi2_2336
Definition: atmdat_readin.cpp:94
ipAl6912
long ipAl6912
Definition: atmdat_readin.cpp:66
t_CoolHeavy::p2_21
double p2_21
Definition: coolheavy.h:100
t_CoolHeavy::S10373
realnum S10373
Definition: coolheavy.h:109
phycon
t_phycon phycon
Definition: phycon.cpp:6
trace.h
ipS1_56m
long ipS1_56m
Definition: atmdat_readin.cpp:95
sil
t_sil sil
Definition: sil.cpp:5
ipTAr7
long ipTAr7
Definition: atmdat_readin.cpp:79
ipfsNa373
long ipfsNa373
Definition: atmdat_readin.cpp:57
t_CoolHeavy::Cl8483
double Cl8483
Definition: coolheavy.h:39
embesq.h
lines_service.h
ipT1198
long ipT1198
Definition: atmdat_readin.cpp:76
ipTMg6
long ipTMg6
Definition: atmdat_readin.cpp:63
ipTMg14
long ipTMg14
Definition: atmdat_readin.cpp:63
ipTAr22
long ipTAr22
Definition: atmdat_readin.cpp:79
ipAlI3957
long ipAlI3957
Definition: atmdat_readin.cpp:64
ipT58
long ipT58
Definition: atmdat_readin.cpp:63
t_CoolHeavy::c5324
double c5324
Definition: coolheavy.h:42
ipTAl48
long ipTAl48
Definition: atmdat_readin.cpp:68
ipP0260
long ipP0260
Definition: atmdat_readin.cpp:74
ipxMg72569
long ipxMg72569
Definition: atmdat_readin.cpp:62
t_CoolHeavy::c8047
double c8047
Definition: coolheavy.h:40
t_sil::c2148
double c2148
Definition: sil.h:16
t_CoolHeavy::Ar4740
double Ar4740
Definition: coolheavy.h:18
ipS1_25m
long ipS1_25m
Definition: atmdat_readin.cpp:95
ipCl1_11m
long ipCl1_11m
Definition: atmdat_readin.cpp:95
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
t_CoolHeavy::c3350
double c3350
Definition: coolheavy.h:30
t_CoolHeavy::S4070
realnum S4070
Definition: coolheavy.h:105
ipfsNa421
long ipfsNa421
Definition: atmdat_readin.cpp:57
ipT1895
long ipT1895
Definition: atmdat_readin.cpp:70
t_CoolHeavy::c1357
double c1357
Definition: coolheavy.h:55
ipT1305
long ipT1305
Definition: atmdat_readin.cpp:71
t_CoolHeavy::Ar7331
double Ar7331
Definition: coolheavy.h:24
ipT1207
long ipT1207
Definition: atmdat_readin.cpp:70
sexp
sys_float sexp(sys_float x)
Definition: service.cpp:914
ipSi10143
long ipSi10143
Definition: atmdat_readin.cpp:72
t_sil::c1985
double c1985
Definition: sil.h:12
ipT4561
long ipT4561
Definition: atmdat_readin.cpp:60
t_CoolHeavy::c3119
double c3119
Definition: coolheavy.h:41
t_CoolHeavy::p2_31
double p2_31
Definition: coolheavy.h:102
ipMgI2853
long ipMgI2853
Definition: atmdat_readin.cpp:58
ipfsMg790
long ipfsMg790
Definition: atmdat_readin.cpp:64
t_CoolHeavy::c5525
double c5525
Definition: coolheavy.h:29
ipT1260
long ipT1260
Definition: atmdat_readin.cpp:71
linadd
void linadd(double xInten, realnum wavelength, const char *chLab, char chInfo, const char *chComment)
Definition: lines_service.cpp:316
dense.h
ipP0233
long ipP0233
Definition: atmdat_readin.cpp:74
ipTAl568
long ipTAl568
Definition: atmdat_readin.cpp:68
t_CoolHeavy::Cl8436
double Cl8436
Definition: coolheavy.h:37
ipT1194
long ipT1194
Definition: atmdat_readin.cpp:75
t_CoolHeavy::Ar7171
double Ar7171
Definition: coolheavy.h:23
ipSILICON
const int ipSILICON
Definition: cddefines.h:318
t_CoolHeavy::c4017
double c4017
Definition: coolheavy.h:59
coolheavy.h
ipxNa6143
long ipxNa6143
Definition: atmdat_readin.cpp:58
ipT2670
long ipT2670
Definition: atmdat_readin.cpp:65
trace
t_trace trace
Definition: trace.cpp:5
ipxMg08303
long ipxMg08303
Definition: atmdat_readin.cpp:62
cddefines.h
ipSULPHUR
const int ipSULPHUR
Definition: cddefines.h:320
ipSi2_2334
long ipSi2_2334
Definition: atmdat_readin.cpp:94
t_CoolHeavy::c2067
double c2067
Definition: coolheavy.h:58
t_CoolHeavy::c4626
double c4626
Definition: coolheavy.h:14
ipSi10_606
long ipSi10_606
Definition: atmdat_readin.cpp:49
ipCl04117
long ipCl04117
Definition: atmdat_readin.cpp:78
t_CoolHeavy::c9532
double c9532
Definition: coolheavy.h:88
ipT705
long ipT705
Definition: atmdat_readin.cpp:60
ipTAr13
long ipTAr13
Definition: atmdat_readin.cpp:79
ipTSi65
long ipTSi65
Definition: atmdat_readin.cpp:73
t_CoolHeavy::c2691
double c2691
Definition: coolheavy.h:15
ipMgI2026
long ipMgI2026
Definition: atmdat_readin.cpp:59
t_CoolHeavy::c7007
double c7007
Definition: coolheavy.h:13
ipxMg52855
long ipxMg52855
Definition: atmdat_readin.cpp:61
ipT2796
long ipT2796
Definition: atmdat_readin.cpp:59
t_CoolHeavy::Ar7263
double Ar7263
Definition: coolheavy.h:22
MAX2
#define MAX2
Definition: cddefines.h:782
ionbal
t_ionbal ionbal
Definition: ionbal.cpp:5
ipTMg625
long ipTMg625
Definition: atmdat_readin.cpp:62
ipfsCl233
long ipfsCl233
Definition: atmdat_readin.cpp:77
ipT1403
long ipT1403
Definition: atmdat_readin.cpp:70
t_CoolHeavy::c2972
double c2972
Definition: coolheavy.h:56
t_CoolHeavy::c6164
double c6164
Definition: coolheavy.h:27
ipTAl550
long ipTAl550
Definition: atmdat_readin.cpp:68
ipS4_1417
long ipS4_1417
Definition: atmdat_readin.cpp:90
t_CoolHeavy::c1806
double c1806
Definition: coolheavy.h:53
ipTSi35
long ipTSi35
Definition: atmdat_readin.cpp:73
t_dense::xIonDense
double xIonDense[LIMELM][LIMELM+1]
Definition: dense.h:125
t_sil::c1446
double c1446
Definition: sil.h:15
ipTAr8
long ipTAr8
Definition: atmdat_readin.cpp:80
TauLines
TransitionList TauLines("TauLines", &AnonStates)
ipTSi521
long ipTSi521
Definition: atmdat_readin.cpp:72
ipSi1_130m
long ipSi1_130m
Definition: atmdat_readin.cpp:97
t_sil::c1815
double c1815
Definition: sil.h:11
ipAl6366
long ipAl6366
Definition: atmdat_readin.cpp:66
sil.h
ipAr06453
long ipAr06453
Definition: atmdat_readin.cpp:80
PutLine
void PutLine(const TransitionProxy &t, const char *chComment, const char *chLabelTemp)
Definition: transition.cpp:449
ionbal.h
t_phycon::te10
double te10
Definition: phycon.h:55
t_CoolHeavy::Ar7237
double Ar7237
Definition: coolheavy.h:25
ipxMg72261
long ipxMg72261
Definition: atmdat_readin.cpp:61
t_CoolHeavy::c3109
double c3109
Definition: coolheavy.h:17
ipTSi3
long ipTSi3
Definition: atmdat_readin.cpp:74
t_embesq::em2669
realnum em2669
Definition: embesq.h:14
ipxMg71190
long ipxMg71190
Definition: atmdat_readin.cpp:61
ipAl8370
long ipAl8370
Definition: atmdat_readin.cpp:67
t_CoolHeavy::S10323
realnum S10323
Definition: coolheavy.h:107
taulines.h
ipT1527
long ipT1527
Definition: atmdat_readin.cpp:71
ipAl8575
long ipAl8575
Definition: atmdat_readin.cpp:66
CoolHeavy
t_CoolHeavy CoolHeavy
Definition: coolheavy.cpp:5
t_CoolHeavy::S10339
realnum S10339
Definition: coolheavy.h:110
ipT1863
long ipT1863
Definition: atmdat_readin.cpp:65
ipxNa6862
long ipxNa6862
Definition: atmdat_readin.cpp:58
phycon.h
ipTAr9
long ipTAr9
Definition: atmdat_readin.cpp:79
t_CoolHeavy::Cl5519
double Cl5519
Definition: coolheavy.h:33
t_phycon::te32
double te32
Definition: phycon.h:49
embesq
t_embesq embesq
Definition: embesq.cpp:5
ipS4_1424
long ipS4_1424
Definition: atmdat_readin.cpp:90
ipTMg610
long ipTMg610
Definition: atmdat_readin.cpp:62
t_CoolHeavy::Ar2868
double Ar2868
Definition: coolheavy.h:20
t_phycon::sqrte
double sqrte
Definition: phycon.h:48
ipT786
long ipT786
Definition: atmdat_readin.cpp:76
ipTSi4
long ipTSi4
Definition: atmdat_readin.cpp:74
ipSi1_68m
long ipSi1_68m
Definition: atmdat_readin.cpp:97
t_CoolHeavy::c8579
double c8579
Definition: coolheavy.h:26
ipTMg4
long ipTMg4
Definition: atmdat_readin.cpp:63
ipT1256
long ipT1256
Definition: atmdat_readin.cpp:75
ipP0318
long ipP0318
Definition: atmdat_readin.cpp:75
ipSi619
long ipSi619
Definition: atmdat_readin.cpp:71
ipCl04203
long ipCl04203
Definition: atmdat_readin.cpp:78
ipT639
long ipT639
Definition: atmdat_readin.cpp:67
ipxNa0746
long ipxNa0746
Definition: atmdat_readin.cpp:58
t_CoolHeavy::c2569
double c2569
Definition: coolheavy.h:54
lines_lv1_na_ar
void lines_lv1_na_ar(void)
Definition: prt_lines_lv1_na_ar.cpp:16
t_CoolHeavy::Cl8504
double Cl8504
Definition: coolheavy.h:36
t_phycon::te
double te
Definition: phycon.h:11
t_embesq::em1198
realnum em1198
Definition: embesq.h:16
t_CoolHeavy::c2428
double c2428
Definition: coolheavy.h:10
t_CoolHeavy::S6733
realnum S6733
Definition: coolheavy.h:103
ipSi2_2344
long ipSi2_2344
Definition: atmdat_readin.cpp:94
ipT1855
long ipT1855
Definition: atmdat_readin.cpp:65
t_sil::c949
double c949
Definition: sil.h:10
t_sil::c691
double c691
Definition: sil.h:13
t_CoolHeavy::c7136
double c7136
Definition: coolheavy.h:12
ipT5895
long ipT5895
Definition: atmdat_readin.cpp:57
t_CoolHeavy::S6718
realnum S6718
Definition: coolheavy.h:104
DEBUG_ENTRY
#define DEBUG_ENTRY(funcname)
Definition: cddefines.h:684
ipxMg51325
long ipxMg51325
Definition: atmdat_readin.cpp:60
ipTS1720
long ipTS1720
Definition: atmdat_readin.cpp:76
ipT1394
long ipT1394
Definition: atmdat_readin.cpp:70
t_CoolHeavy::S10289
realnum S10289
Definition: coolheavy.h:108
ipTS34
long ipTS34
Definition: atmdat_readin.cpp:77
t_trace::lgTrace
bool lgTrace
Definition: trace.h:12
ipSi2_2329
long ipSi2_2329
Definition: atmdat_readin.cpp:94
t_CoolHeavy::S4078
realnum S4078
Definition: coolheavy.h:106
t_CoolHeavy::Cl3344
double Cl3344
Definition: coolheavy.h:35
ipfsCl214
long ipfsCl214
Definition: atmdat_readin.cpp:77
TransitionList::Emis
EmissionList & Emis()
Definition: transition.h:329