cloudy  trunk
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
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,'c',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,'c',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,'c',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,'c',true,
60  " [Na VI] 2568.9 " );
61 
62  PntForLine(1357.,"Na 6",&ipnt);
63  lindst(CoolHeavy.c1357,1357,"Na 6",ipnt,'c',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,'c',true ,
68  " [Na VI] 2971.9 ");
69 
70  PntForLine(2872.,"Na 6",&ipnt);
71  lindst(CoolHeavy.c2972/(1.+3.02),2872,"Na 6",ipnt,'c',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,'c',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,'c',true,
201  " [Al VI] 2428.4 " );
202 
203  PntForLine(2601.,"Al 6",&ipnt);
204  lindst(CoolHeavy.c2428/(1.+3.73),2601,"Al 6",ipnt,'c',true,
205  " [Al VI] 2601.0");
206 
207  PntForLine(1170.,"Al 6",&ipnt);
208  lindst(CoolHeavy.c1170,1170,"Al 6",ipnt,'c',true,
209  " [Al VI] 1169.86 ");
210 
211  PntForLine(2125.,"Al 6",&ipnt);
212  lindst(CoolHeavy.c2125,2125,"Al 6",ipnt,'c',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 
308  PutLine(TauLines[ipT1207],
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,'c',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,'c',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,'c',true,
359  " SI VIII 1446, OIII like, collisionally excited" );
360 
361  PntForLine(1985.,"Si 9",&ipnt);
362  lindst(sil.c1985,1985,"Si 9",ipnt,'c',true,
363  " SI IX 1985, 2150, collisionally excited" );
364 
365 
366  PntForLine(949.,"Si 9",&ipnt);
367  lindst(sil.c949,949,"Si 9",ipnt,'c',true,
368  " collisionally excited" );
369 
370  PntForLine(1815.,"Si 9",&ipnt);
371  lindst(sil.c1815,1815,"Si 9",ipnt,'c',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,'c',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,'c',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,'c',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,'c',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,'c',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,'c',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,'c',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,'i',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,'c',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,'c',true,
461  " individual line from five level atom" );
462 
463  PntForLine(6716.,"S II",&ipnt);
464  lindst(CoolHeavy.S6718,6716,"S II",ipnt,'c',true,
465  " individual line from five level atom" );
466 
467  PntForLine(4070.,"S II",&ipnt);
468  lindst(CoolHeavy.S4070,4070,"S II",ipnt,'c',true,
469  " individual line from five level atom " );
470 
471  PntForLine(4078.,"S II",&ipnt);
472  lindst(CoolHeavy.S4078,4078,"S II",ipnt,'c',true,
473  " individual line from five level atom" );
474 
475  PntForLine(10330.,"S 2",&ipnt);
476  lindst(CoolHeavy.S10323,10323,"S II",ipnt,'c',false,
477  " individual line from five level atom " );
478 
479  lindst(CoolHeavy.S10289,10289,"S II",ipnt,'c',false,
480  " individual line from five level atom" );
481 
482  lindst(CoolHeavy.S10373,10373,"S II",ipnt,'c' ,false,
483  " individual line from five level atom ");
484 
485  lindst(CoolHeavy.S10339,10339,"S II",ipnt,'c',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,'c',true ,
502  " [S III] 9532 alone ");
503 
504  PntForLine(9069.,"S 3",&ipnt);
505  lindst(CoolHeavy.c9532/(1.+2.48),9069,"S 3",ipnt,'c',true,
506  " [S III] 9069 alone" );
507 
508  PntForLine(6312.,"S 3",&ipnt);
509  lindst(CoolHeavy.c6312,6312,"S 3",ipnt,'c',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,'c',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,'c',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,'c',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,'c',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,'c',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,'c',true,
589  " Cl III 5538 " );
590 
591  PntForLine(5518.,"Cl 3",&ipnt);
592  lindst(CoolHeavy.Cl5519,5518,"Cl 3",ipnt,'c',true,
593  " Cl III 5518" );
594 
595  PntForLine(3354.,"Cl 3",&ipnt);
596  lindst(CoolHeavy.Cl3354,3354,"Cl 3",ipnt,'c',true,
597  " Cl III 3354 " );
598 
599  PntForLine(3344.,"Cl 3",&ipnt);
600  lindst(CoolHeavy.Cl3344,3344,"Cl 3",ipnt,'c',true,
601  " Cl III 3344 " );
602 
603  PntForLine(8504.,"Cl 3",&ipnt);
604  lindst(CoolHeavy.Cl8504,8504,"Cl 3",ipnt,'c',true ,
605  " Cl III 8504 ");
606 
607  PntForLine(8436.,"Cl 3",&ipnt);
608  lindst(CoolHeavy.Cl8436,8436,"Cl 3",ipnt,'c',true,
609  " Cl III 8436" );
610 
611  PntForLine(8552.,"Cl 3",&ipnt);
612  lindst(CoolHeavy.Cl8552,8552,"Cl 3",ipnt,'c',true,
613  " Cl III 8552 " );
614 
615  PntForLine(8483.,"Cl 3",&ipnt);
616  lindst(CoolHeavy.Cl8483,8483,"Cl 3",ipnt,'c',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,'c',true,
627  " ClIV 8047" );
628 
629  PntForLine(7532.,"Cl 4",&ipnt);
630  lindst(CoolHeavy.c8047*0.333,7532,"Cl 4",ipnt,'c',true,
631  " ClIV 7532" );
632 
633  PntForLine(3119.,"Cl 4",&ipnt);
634  lindst(CoolHeavy.c3119,3119,"Cl 4",ipnt,'c',true,
635  " ClIV 3119" );
636 
637  PntForLine(5324.,"Cl 4",&ipnt);
638  lindst(CoolHeavy.c5324,5324,"Cl 4",ipnt,'c',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,'c',true ,
649  " Argon III 7135");
650 
651  PntForLine(7751.,"Ar 3",&ipnt);
652  lindst(CoolHeavy.c7136/(1.+4.144),7751,"Ar 3",ipnt,'c',true,
653  " Argon III 7751" );
654 
655  PntForLine(5192.,"Ar 3",&ipnt);
656  lindst(CoolHeavy.c5192,5192,"Ar 3",ipnt,'c',true,
657  " Argon III 5192" );
658 
659  PntForLine(3109.,"Ar 3",&ipnt);
660  lindst(CoolHeavy.c3109*0.9894,3109,"Ar 3",ipnt,'c',true,
661  " Argon III 3109" );
662 
663  PntForLine(3005.,"Ar 3",&ipnt);
664  lindst(CoolHeavy.c3109*(1.-0.9894),3005,"Ar 3",ipnt,'c',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,'c',true,
684  " [Ar IV] 4740" );
685 
686  PntForLine(4711.,"Ar 4",&ipnt);
687  lindst(CoolHeavy.Ar4711,4711,"Ar 4",ipnt,'c',true,
688  " [Ar IV] 4711" );
689 
690  PntForLine(2868.,"Ar 4",&ipnt);
691  lindst(CoolHeavy.Ar2868,2868,"Ar 4",ipnt,'c',true,
692  " [Ar IV] 2868" );
693 
694  PntForLine(2854.,"Ar 4",&ipnt);
695  lindst(CoolHeavy.Ar2854,2854,"Ar 4",ipnt,'c',true,
696  " [Ar IV] 2854" );
697 
698 
699  PntForLine(7263.,"Ar 4",&ipnt);
700  lindst(CoolHeavy.Ar7263,7263,"Ar 4",ipnt,'c',true,
701  " [Ar IV] 7263" );
702 
703  PntForLine(7171.,"Ar 4",&ipnt);
704  lindst(CoolHeavy.Ar7171,7171,"Ar 4",ipnt,'c',true,
705  " [Ar IV] 7171" );
706 
707  PntForLine(7331.,"Ar 4",&ipnt);
708  lindst(CoolHeavy.Ar7331,7331,"Ar 4",ipnt,'c',true,
709  " [Ar IV] 7331" );
710 
711  PntForLine(7237.,"Ar 4",&ipnt);
712  lindst(CoolHeavy.Ar7237,7237,"Ar 4",ipnt,'c',true,
713  " [Ar IV] 7237" );
714 
715  PntForLine(7005.,"Ar 5",&ipnt);
716  lindst(CoolHeavy.c7007/(1.+1./2.143),7005,"Ar 5",ipnt,'c',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,'c',true ,
721  " Argon V, 3P lines, 6435, collisionally excited");
722 
723  PntForLine(4626.,"Ar 5",&ipnt);
724  lindst( CoolHeavy.c4626 , 4626 ,"Ar 5",ipnt,'c',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,'c',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 }
double c583
Definition: sil.h:10
double c5192
Definition: coolheavy.h:9
double c1357
Definition: coolheavy.h:9
long ipAr06453
long ipT2796
double c1446
Definition: sil.h:10
long ipS4_1417
double Ar2854
Definition: coolheavy.h:9
long ipAl6366
long ipTAl550
double Ar7331
Definition: coolheavy.h:9
long ipSi10143
long ipTSi499
long ipTS1720
double p2_31
Definition: coolheavy.h:9
void lindst(double xInten, realnum wavelength, const char *chLab, long int ipnt, char chInfo, bool lgOutToo, const char *chComment)
double te03
Definition: phycon.h:48
double c691
Definition: sil.h:10
long ipTS11
long ipMgI2853
long ipxMg52417
double emit_frac(const TransitionProxy &t)
Definition: transition.cpp:84
long ipTSi35
double c2428
Definition: coolheavy.h:9
long ipxNa6143
long ipTSi65
t_sil sil
Definition: sil.cpp:5
double Ar4711
Definition: coolheavy.h:9
double c3679
Definition: coolheavy.h:9
long ipP713
t_embesq embesq
Definition: embesq.cpp:5
long ipxNa6862
double Ar7171
Definition: coolheavy.h:9
double c6312
Definition: coolheavy.h:9
long ipTMg14
#define MAX2
Definition: cddefines.h:786
realnum S6733
Definition: coolheavy.h:103
long ipMgI2026
long ipCl1_11m
long ipP0318
long ipxMg71190
t_phycon phycon
Definition: phycon.cpp:6
double Cl8504
Definition: coolheavy.h:9
long ipSi619
t_dense dense
Definition: dense.cpp:24
long ipfsNa490
realnum S4078
Definition: coolheavy.h:103
long ipT2670
sys_float sexp(sys_float x)
Definition: service.cpp:914
long ipSii2518
double c7007
Definition: coolheavy.h:9
t_CoolHeavy CoolHeavy
Definition: coolheavy.cpp:5
long ipT1808
double c4017
Definition: coolheavy.h:9
FILE * ioQQQ
Definition: cddefines.cpp:7
double Ar7263
Definition: coolheavy.h:9
double p2_21
Definition: coolheavy.h:9
long ipSi10_606
long ipT58
long ipS4_1424
long ipTMg610
long ipxMg51325
realnum S10323
Definition: coolheavy.h:103
long ipxMg72569
double c1815
Definition: sil.h:10
long ipS4_1398
long ipT705
const int ipSULPHUR
Definition: cddefines.h:324
void PntForLine(double wavelength, const char *chLabel, long int *ipnt)
double xIonDense[LIMELM][LIMELM+1]
Definition: dense.h:125
long ipTSi521
double c5324
Definition: coolheavy.h:9
long ipfsNa421
long ipT1895
long ipxNa0746
long ipT1527
long ipSi2_2334
long ipT1198
double c2972
Definition: coolheavy.h:9
long ipSi2_2344
long ipTSi3
double c1985
Definition: sil.h:10
t_trace trace
Definition: trace.cpp:5
long ipxMg72261
long ipTSi4
long ipT5895
double c10330
Definition: coolheavy.h:9
t_ionbal ionbal
Definition: ionbal.cpp:5
long ipCl04117
realnum em2669
Definition: embesq.h:9
realnum S10373
Definition: coolheavy.h:103
long ipfsMg790
long ipS4_1407
long ipTS34
realnum S6718
Definition: coolheavy.h:103
double c1365
Definition: coolheavy.h:9
long ipTAr8
long ipfsCl214
realnum S10289
Definition: coolheavy.h:103
long ipS4_1405
long ipT1305
long ipAlI3090
bool lgTrace
Definition: trace.h:12
double Cl8436
Definition: coolheavy.h:9
long ipTMg6
double c2148
Definition: sil.h:10
long ipAl09204
long ipT1207
long ipTSi41
long ipSi1_68m
long ipT4561
double Cl8552
Definition: coolheavy.h:9
long ipTAr13
double c5525
Definition: coolheavy.h:9
double c2569
Definition: coolheavy.h:9
long ipCl973
long ipfsNa373
long ipSi1_130m
double c3119
Definition: coolheavy.h:9
double Cl3354
Definition: coolheavy.h:9
long ipTMg625
double Ar4740
Definition: coolheavy.h:9
long ipT1855
long ipT1863
long ipAl529
long ipT2804
long ipTAr22
long ipfsMg755
void PutLine(const TransitionProxy &t, const char *chComment, const char *chLabelTemp)
Definition: transition.cpp:449
double **** PhotoRate_Shell
Definition: ionbal.h:111
double c9532
Definition: coolheavy.h:9
long ipS1_56m
long ipSi2_2336
const int ipSILICON
Definition: cddefines.h:322
long ipTSi25
realnum em1198
Definition: embesq.h:9
double c8579
Definition: coolheavy.h:9
realnum S10339
Definition: coolheavy.h:103
double c7136
Definition: coolheavy.h:9
long ipfsCl233
double c8494
Definition: coolheavy.h:9
double c3109
Definition: coolheavy.h:9
TransitionList TauLines("TauLines",&AnonStates)
double c6164
Definition: coolheavy.h:9
double c1806
Definition: coolheavy.h:9
double Ar2868
Definition: coolheavy.h:9
#define DEBUG_ENTRY(funcname)
Definition: cddefines.h:688
long ipP0260
long ipT1256
double c6731
Definition: coolheavy.h:9
double c1170
Definition: coolheavy.h:9
double te10
Definition: phycon.h:48
long ipT1394
long ipTAl48
long ipTMg4
double Cl5519
Definition: coolheavy.h:9
long ipTAl568
double eden
Definition: dense.h:190
long ipSii2215
long ipTS19
long ipAl8575
double c8047
Definition: coolheavy.h:9
void lines_lv1_na_ar(void)
long ipTAr9
double p2_32
Definition: coolheavy.h:9
long ipT786
long ipT1260
long ipT1194
double c2691
Definition: coolheavy.h:9
double sqrte
Definition: phycon.h:48
double Ar7237
Definition: coolheavy.h:9
void linadd(double xInten, realnum wavelength, const char *chLab, char chInfo, const char *chComment)
long ipP0233
long ipCl04203
double c2067
Definition: coolheavy.h:9
long ipSi2_2329
long ipT639
long ipxMg08303
realnum em1895
Definition: embesq.h:9
double te
Definition: phycon.h:11
double Cl3344
Definition: coolheavy.h:9
double c949
Definition: sil.h:10
long ipAl8370
long ipxMg52855
long ipTAr7
double c3350
Definition: coolheavy.h:9
long ipAlI3957
void PutExtra(double Extra)
Definition: transition.cpp:468
long ipT1403
long ipSi2_2350
double Cl8483
Definition: coolheavy.h:9
double Cl5539
Definition: coolheavy.h:9
EmissionList & Emis()
Definition: transition.h:329
double te32
Definition: phycon.h:48
double c4626
Definition: coolheavy.h:9
double c2125
Definition: coolheavy.h:9
realnum S4070
Definition: coolheavy.h:103
long ipAl6912
long ipS1_25m