00001
00002
00003
00004 #include "cddefines.h"
00005 #include "coolheavy.h"
00006 #include "sil.h"
00007 #include "phycon.h"
00008 #include "embesq.h"
00009 #include "taulines.h"
00010 #include "dense.h"
00011 #include "ionbal.h"
00012 #include "trace.h"
00013 #include "lines_service.h"
00014 #include "lines.h"
00015
00016 void lines_lv1_na_ar(void)
00017 {
00018 long int ipnt;
00019 double drec,
00020 fac,
00021 rec,
00022 sum,
00023 t4;
00024
00025 DEBUG_ENTRY( "lines_lv1_na_ar()" );
00026
00027 if( trace.lgTrace )
00028 {
00029 fprintf( ioQQQ, " lines_lv1_na_ar called\n" );
00030 }
00031
00032 t4 = phycon.te/1e4;
00033
00034 PutLine(TauLines[ipT5895],
00035 " sodium sum of Na D lines");
00036
00037 PutLine(TauLines[ipfsNa373],
00038 " [NaIII] 7.319 micron ");
00039
00040 PutLine(TauLines[ipfsNa490],
00041 " [NaIV] 9.048 micron ");
00042
00043 PutLine(TauLines[ipfsNa421],
00044 " [NaIV] 21.29 micron ");
00045
00046 PntForLine(1365.,"Na 5",&ipnt);
00047 lindst(CoolHeavy.c1365,1365,"Na 5",ipnt,'c',true ,
00048 " [NaV] 1365, sum of 1365.1+1365.8; cs only guess ");
00049
00050 PntForLine(2067.,"Na 5",&ipnt);
00051 lindst(CoolHeavy.c2067,2067,"Na 5",ipnt,'c',true ,
00052 " [NaV] 2067, sum of 2066.9+2068.4; cs only guess ");
00053
00054 PntForLine(4017.,"Na 5",&ipnt);
00055 lindst(CoolHeavy.c4017,4017,"Na 5",ipnt,'c',true,
00056 " [NaV] 4017, sum of 4010.9+4016.7+4022.7; cs only guess " );
00057
00058 PntForLine(2569.,"Na 6",&ipnt);
00059 lindst(CoolHeavy.c2569,2569,"Na 6",ipnt,'c',true,
00060 " [Na VI] 2568.9 " );
00061
00062 PntForLine(1357.,"Na 6",&ipnt);
00063 lindst(CoolHeavy.c1357,1357,"Na 6",ipnt,'c',true ,
00064 " [Na VI] 1356.6 ");
00065
00066 PntForLine(2972.,"Na 6",&ipnt);
00067 lindst(CoolHeavy.c2972/(1.+1./3.02),2972,"Na 6",ipnt,'c',true ,
00068 " [Na VI] 2971.9 ");
00069
00070 PntForLine(2872.,"Na 6",&ipnt);
00071 lindst(CoolHeavy.c2972/(1.+3.02),2872,"Na 6",ipnt,'c',true,
00072 " [Na VI] 2872.7 ");
00073
00074 PutLine(TauLines[ipxNa6143],
00075 " [NaVI] 14.32 micron ");
00076
00077 PutLine(TauLines[ipxNa6862],
00078 " [NaVI] 8.62 micron ");
00079
00080 PutLine(TauLines[ipxNa0746],
00081 " [NaVII] 4.675 micron ");
00082
00083 PutLine(TauLines[ipT4561],
00084 "Magnesium I 4571, O I data for coll strength and trans prob ");
00085
00086 PutLine(TauLines[ipMgI2853],
00087 " Mg I 2853 ");
00088
00089 PutLine(TauLines[ipMgI2026],
00090 " Mg I 2026 ");
00091
00092 linadd(TauLines[ipT2796].Emis().xIntensity()+TauLines[ipT2804].Emis().xIntensity(),2798,
00093 "TOTL",'i',"Mg II 2798 add both lines of multiplet together " );
00094
00095
00096 sum = TauLines[ipT2796].Emis().xIntensity()*TauLines[ipT2796].Emis().FracInwd() +
00097 TauLines[ipT2804].Emis().xIntensity()*TauLines[ipT2804].Emis().FracInwd();
00098
00099 linadd(sum,2798,"Inwd",'i',
00100 "inward part of Mg II 2798" );
00101
00102 PutLine(TauLines[ipT2796],
00103 "one member of Mg II multiplet");
00104
00105 PutLine(TauLines[ipT2804],
00106 "one member of Mg II multiplet");
00107
00108 PutLine(TauLines[ipTMg4],
00109 " Mg IV 4.5 micron ");
00110
00111 PutLine(TauLines[ipTMg14],
00112 " Mg V 13.5 micron emission");
00113
00114 PutLine(TauLines[ipTMg6],
00115 " Mg V 5.6 micron emission ");
00116
00117 PutLine(TauLines[ipxMg52855],
00118 " [Mg 5] 2571, 2893 ");
00119
00120 PutLine(TauLines[ipxMg52417],
00121 " [Mg 5] 2417.5A, 3-2 in model atom" );
00122
00123 PutLine(TauLines[ipxMg51325],
00124 " [Mg 5] 1324.58A, 3-1 in model atom ");
00125
00126 PntForLine(1806.,"Mg 6",&ipnt);
00127 lindst(CoolHeavy.c1806,1806,"Mg 6",ipnt,'c',true ,
00128 " MG VI");
00129
00130 PutLine(TauLines[ipxMg72569],
00131 " [Mg 7] 2510, 2629, 2-1 transitions, together");
00132
00133 PutLine(TauLines[ipxMg72261],
00134 " [Mg 7] 3-2 transition, 2261.5 ");
00135
00136 PutLine(TauLines[ipxMg71190],
00137 " [Mg 7] 3-1 transition, 1189.82A ");
00138
00139 PutLine(TauLines[ipfsMg755],
00140 " Mg 7 IR line 5.50 microns ");
00141
00142 PutLine(TauLines[ipfsMg790],
00143 " Mg 7 IR line 9.03 microns ");
00144
00145 PutLine(TauLines[ipxMg08303],
00146 " [Mg 8] 3.03 micron ");
00147
00148 PutLine(TauLines[ipT705]," Mg 9 704.5 ");
00149
00150 linadd(TauLines[ipTMg610].Emis().xIntensity()+TauLines[ipTMg625].Emis().xIntensity(),615,"TOTL",'i',
00151 " Mg 10 614.9 both of doublet, li seq 2s 2p" );
00152
00153 PutLine(TauLines[ipTMg610],
00154 "");
00155 PutLine(TauLines[ipTMg625],
00156 "");
00157
00158 PutLine(TauLines[ipT58],
00159 " part of Mg 10 destroyed by background opacity Mg 10 58 li seq 2s 3p ");
00160 linadd(TauLines[ipTMg610].Emis().ots()*TauLines[ipTMg610].EnergyErg()+
00161 TauLines[ipTMg625].Emis().ots()*TauLines[ipTMg625].EnergyErg(),615,"dest",'i',
00162 "" );
00163
00164 PutLine(TauLines[ipAlI3957],
00165 " Aluminum Al I 3957 ");
00166
00167 PutLine(TauLines[ipAlI3090],
00168 " Al I 3090 ");
00169
00170 linadd(embesq.em2669+TauLines[ipT2670].Emis().xIntensity(),2665,"totl",'i',
00171 "Al II 1671 total emission in Al II] 2669.7, 2660 doublet" );
00172 PutLine(TauLines[ipT2670],
00173 " ");
00174
00175 linadd(embesq.em2669,2660,"Al 2",'i',
00176 "emission in Al II] 2669 alone" );
00177
00178 linadd(TauLines[ipT1855].Emis().xIntensity()+TauLines[ipT1863].Emis().xIntensity(),1860,"TOTL",'i',
00179 " Al III" );
00180 sum = TauLines[ipT1855].Emis().xIntensity()*TauLines[ipT1855].Emis().FracInwd() +
00181 TauLines[ipT1863].Emis().xIntensity()* TauLines[ipT1863].Emis().FracInwd();
00182
00183 linadd(sum,1860,"Inwd",'i',
00184 " inward part of AlIII line" );
00185 PutLine(TauLines[ipT1855],
00186 "");
00187 PutLine(TauLines[ipT1863],
00188 "");
00189
00190 PutLine(TauLines[ipAl529],
00191 " [Al V] 2.905 micron ");
00192
00193 PutLine(TauLines[ipAl6366],
00194 " [Al VI] 3.66 micron ");
00195
00196 PutLine(TauLines[ipAl6912],
00197 " [Al VI] 9.12 micron");
00198
00199 PntForLine(2428.,"Al 6",&ipnt);
00200 lindst(CoolHeavy.c2428/(1.+1./3.73),2428,"Al 6",ipnt,'c',true,
00201 " [Al VI] 2428.4 " );
00202
00203 PntForLine(2601.,"Al 6",&ipnt);
00204 lindst(CoolHeavy.c2428/(1.+3.73),2601,"Al 6",ipnt,'c',true,
00205 " [Al VI] 2601.0");
00206
00207 PntForLine(1170.,"Al 6",&ipnt);
00208 lindst(CoolHeavy.c1170,1170,"Al 6",ipnt,'c',true,
00209 " [Al VI] 1169.86 ");
00210
00211 PntForLine(2125.,"Al 6",&ipnt);
00212 lindst(CoolHeavy.c2125,2125,"Al 6",ipnt,'c',true ,
00213 " [Al VI] 2124.95" );
00214
00215 PutLine(TauLines[ipAl8575],
00216 " [Al VIII] 5.75 micron " );
00217
00218 PutLine(TauLines[ipAl8370],
00219 " [Al VIII] 3.70 micron");
00220
00221 PutLine(TauLines[ipAl09204],
00222 " [Al IX] 2.04 micron, no collision strength, A NIST ");
00223
00224 PutLine(TauLines[ipT639]," Al X ");
00225
00226 linadd(TauLines[ipTAl550].Emis().xIntensity()+TauLines[ipTAl568].Emis().xIntensity(),556,"TOTL",'i',
00227 " Al 11, Li seq 2s2p" );
00228
00229 PutLine(TauLines[ipTAl550],
00230 "");
00231
00232 PutLine(TauLines[ipTAl568],
00233 "");
00234
00235 PutLine(TauLines[ipTAl48],
00236 " Al 11, Li seq 2s3p ");
00237
00238 PutLine(TauLines[ipSi1_130m],
00239 " silicon Silicon Si I 130m ");
00240
00241 PutLine(TauLines[ipSi1_68m],
00242 " Si I 68m");
00243
00244 PutLine(TauLines[ipSii2518],
00245 " Si I 2518A ");
00246
00247 PutLine(TauLines[ipSii2215],
00248 " Si I 2215A ");
00249
00250 PutLine(TauLines[ipTSi35],
00251 " Silicon II 35 micron ");
00252
00253 linadd(
00254 TauLines[ipSi2_2334].Emis().xIntensity()+
00255 TauLines[ipSi2_2329].Emis().xIntensity()+
00256 TauLines[ipSi2_2350].Emis().xIntensity()+
00257 TauLines[ipSi2_2344].Emis().xIntensity()+
00258 TauLines[ipSi2_2336].Emis().xIntensity(),
00259 2335,"TOTL",'i',
00260 "total intensity of S IV] 1406, all lines in the multiplet" );
00261 PutLine(TauLines[ipSi2_2334],
00262 " ");
00263 PutLine(TauLines[ipSi2_2329],
00264 " ");
00265 PutLine(TauLines[ipSi2_2350],
00266 " ");
00267 PutLine(TauLines[ipSi2_2344],
00268 " ");
00269 PutLine(TauLines[ipSi2_2336],
00270 " ");
00271
00272 PutLine(TauLines[ipT1808],
00273 " SI II 1808, permitted resonance line, collisionally excited ");
00274
00275 PutLine(TauLines[ipT1527],
00276 " SI II 1527, permitted resonance line, collisionally excited ");
00277
00278 PutLine(TauLines[ipT1305],
00279 " SI II 1305, permitted resonance line, collisionally excited ");
00280
00281 PutLine(TauLines[ipT1260],
00282 " SI II 1260, permitted resonance line, collisionally excited ");
00283
00284
00285 drec = dense.xIonDense[13][2]*dense.eden*7.6e-7/phycon.te32*
00286 1.57e-11;
00287
00288 fac = emit_frac(TauLines[ipT1260]);
00289 PntForLine(1260.,"Si 2",&ipnt);
00290 lindst(drec*fac,1260,"diel",ipnt,'r',true,
00291 " fac = (1.-TauLines[ipT1260].ColOvTot());" );
00292
00293
00294 PntForLine(1909.,"Si 2",&ipnt);
00295 lindst(drec*1260./1909.,1909,"diel",ipnt,'r',true,
00296 " dielectronic recombination SiII 1909" );
00297
00298 rec = 1e-12*(0.1152/t4 - 0.3082 + 4.4734*t4 + 0.0207*t4*t4)/pow(t4,1.5)*
00299 sexp(0.2981/t4);
00300
00301
00302
00303 rec *= dense.xIonDense[13][3]*dense.eden*emit_frac(TauLines[ipT1207])*
00304 1.65e-11;
00305
00306 PutExtra(MAX2(0.,rec));
00307
00308 PutLine(TauLines[ipT1207],
00309 " SI III 1207, collisional excitation and dielectronic recombination ");
00310
00311 linadd(MAX2(0.,rec),1207,"rec ",'i',
00312 " Si III 1207, dielectronic recombination only" );
00313
00314 linadd(embesq.em1895+TauLines[ipT1895].Emis().xIntensity(),1888,"TOTL",'i',
00315 " Si III] 1892+1883, total intensity of both lines" );
00316 PutLine(TauLines[ipT1895],
00317 " ");
00318
00319 PntForLine(1883.,"Si 3",&ipnt);
00320 lindst(embesq.em1895,1883,"Si 3",ipnt,'c',true ,
00321 " Si III] 1883 by itself");
00322
00323
00324 fac = emit_frac(TauLines[ipT1895]);
00325
00326 double p1895 = ionbal.PhotoRate_Shell[ipSILICON][1][2][0]*
00327 dense.xIonDense[ipSILICON][1]*0.85;
00328 linadd( p1895*1.05e-11*fac,1895,"PHOT",'i',
00329 " photoproduction by inner shell removal" );
00330
00331 linadd(TauLines[ipT1403].Emis().xIntensity()+TauLines[ipT1394].Emis().xIntensity(),1397,"TOTL",'i',
00332 " Si IV 1397, collisionally excited " );
00333
00334 sum = TauLines[ipT1403].Emis().xIntensity()*TauLines[ipT1403].Emis().FracInwd() +
00335 TauLines[ipT1394].Emis().xIntensity()* TauLines[ipT1394].Emis().FracInwd();
00336
00337 linadd(sum,1397,"Inwd",'i',
00338 " inward part of SiIV 1397" );
00339 PutLine(TauLines[ipT1403],
00340 " ");
00341 PutLine(TauLines[ipT1394],
00342 " ");
00343
00344 PutLine(TauLines[ipSi619],
00345 " SI VI 1.9641 micron ");
00346
00347 PntForLine(2148.,"Si 7",&ipnt);
00348 lindst(sil.c2148,2148,"Si 7",ipnt,'c',true,
00349 " SI VII, 2148, O III like, collisionally excited" );
00350
00351 PutLine(TauLines[ipTSi25],
00352 " Si VII 2.48, 6.49 micron, collisionally excited ");
00353
00354 PutLine(TauLines[ipTSi65],
00355 " Si VII 2.48, 6.49 micron, collisionally excited ");
00356
00357 PntForLine(1446.,"Si 8",&ipnt);
00358 lindst(sil.c1446,1446,"Si 8",ipnt,'c',true,
00359 " SI VIII 1446, OIII like, collisionally excited" );
00360
00361 PntForLine(1985.,"Si 9",&ipnt);
00362 lindst(sil.c1985,1985,"Si 9",ipnt,'c',true,
00363 " SI IX 1985, 2150, collisionally excited" );
00364
00365
00366 PntForLine(949.,"Si 9",&ipnt);
00367 lindst(sil.c949,949,"Si 9",ipnt,'c',true,
00368 " collisionally excited" );
00369
00370 PntForLine(1815.,"Si 9",&ipnt);
00371 lindst(sil.c1815,1815,"Si 9",ipnt,'c',true,
00372 " collisionally excited " );
00373
00374 PutLine(TauLines[ipTSi4],
00375 " SI 9, 3.86, 2.84 3P fine structure lines ");
00376
00377 PutLine(TauLines[ipTSi3],
00378 " SI 9, 3.86, 2.84 3P fine structure lines ");
00379
00380 PntForLine(691.,"Si 9",&ipnt);
00381 lindst(sil.c691,691,"Si 9",ipnt,'c',true ,
00382 " both components of 5S-3P doublet");
00383
00384 PutLine(TauLines[ipSi10_606],
00385 " SI 10 606A, actually group of 4 intercombination lines ");
00386
00387 PutLine(TauLines[ipSi10143],
00388 " [Si 10] 1.43 micron, collisionally excited ");
00389
00390 PntForLine(583.,"Si11",&ipnt);
00391 lindst(sil.c583,581,"Si11",ipnt,'c',true,
00392 " Si 11 582.9, collisionally excited >>chng 01 may 23, wavelength from 583 to 581" );
00393
00394 linadd(TauLines[ipTSi499].Emis().xIntensity()+TauLines[ipTSi521].Emis().xIntensity(),506,"TOTL",'i' ,
00395 "emission total Si 12 506 + 499 ");
00396
00397 PutLine(TauLines[ipTSi499],
00398 " Si 12 506 li seq 2s 2p ");
00399 PutLine(TauLines[ipTSi521],
00400 " ");
00401
00402 PutLine(TauLines[ipTSi41],
00403 " Si 12 40.9A, li seq 2s 3p ");
00404
00405 PutLine(TauLines[ipP0260],
00406 " phosphorus [P II] 60.64 micron ");
00407
00408 PutLine(TauLines[ipP0233],
00409 " [P II] 32.87 micron ");
00410
00411 PntForLine(16400.,"P 2",&ipnt);
00412 lindst(CoolHeavy.p2_32,16400,"P 2",ipnt,'c',true ,
00413 " 3-2 1.64 micron ");
00414
00415 PntForLine(4669.,"P 2",&ipnt);
00416 lindst(CoolHeavy.p2_31*0.75,4669,"P 2",ipnt,'c',true ,
00417 " >>chng 01 may 15, add these lines P 2 3-1 4670, 4738 vac wl, 4669, 4737 air");
00418
00419 PntForLine(4737.,"P 2",&ipnt);
00420 lindst(CoolHeavy.p2_31*0.25,4737,"P 2",ipnt,'c',true ,
00421 " P 2 3-1 4670, 4738 vac wl, 4737, 4737 air");
00422
00423 PntForLine(11890.,"P 2",&ipnt);
00424 lindst(CoolHeavy.p2_21*0.75,11890,"P 2",ipnt,'c',true,
00425 " 2-1 1.147, 1.189 micron" );
00426
00427 PntForLine(11470.,"P 2",&ipnt);
00428 lindst(CoolHeavy.p2_21*0.25,11470,"P 2",ipnt,'c',true ,
00429 " [P II] 1.14 micron");
00430
00431 PutLine(TauLines[ipP0318],
00432 " [P III] 17.885 micron ");
00433
00434 PutLine(TauLines[ipP713],
00435 " [P VII] 1.3745 micron ");
00436
00437 PutLine(TauLines[ipS1_25m],
00438 " sulphur S I 25m ");
00439
00440 PutLine(TauLines[ipS1_56m],
00441 " S I 56m ");
00442
00443 linadd(1.75e-22*dense.eden*dense.xIonDense[ipSULPHUR][1]/phycon.sqrte/
00444 phycon.te10/phycon.te03,1807,"S 1R",'i',
00445 " 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 " );
00446
00447 linadd(CoolHeavy.c6731,6720,"S 2",'i',
00448 " S II 6731 + 6716 together " );
00449
00450 linadd(CoolHeavy.S4070+CoolHeavy.S4078,4074,
00451 "S 2",'i',"S II 4070 +4078 together" );
00452
00453 PntForLine(10330.,"S 2",&ipnt);
00454 lindst(CoolHeavy.c10330,10330,"S 2",ipnt,'c',true,
00455 " S II N=3 lines, all four lines together " );
00456
00457 PntForLine(6731.,"S II",&ipnt);
00458 lindst(CoolHeavy.S6733,6731,"S II",ipnt,'c',true,
00459 " individual line from five level atom" );
00460
00461 PntForLine(6716.,"S II",&ipnt);
00462 lindst(CoolHeavy.S6718,6716,"S II",ipnt,'c',true,
00463 " individual line from five level atom" );
00464
00465 PntForLine(4070.,"S II",&ipnt);
00466 lindst(CoolHeavy.S4070,4070,"S II",ipnt,'c',true,
00467 " individual line from five level atom " );
00468
00469 PntForLine(4078.,"S II",&ipnt);
00470 lindst(CoolHeavy.S4078,4078,"S II",ipnt,'c',true,
00471 " individual line from five level atom" );
00472
00473 linadd(CoolHeavy.S10323,10323,"S II",'c',
00474 " individual line from five level atom " );
00475
00476 linadd(CoolHeavy.S10289,10289,"S II",'c',
00477 " individual line from five level atom" );
00478
00479 linadd(CoolHeavy.S10373,10373,"S II",'c' ,
00480 " individual line from five level atom ");
00481
00482 linadd(CoolHeavy.S10339,10339,"S II",'c',
00483 " individual line from five level atom " );
00484
00485 PutLine(TauLines[ipT1256],
00486 " resonance line near NV, collisionally excited ");
00487
00488 PutLine(TauLines[ipTS19],
00489 " S III fine structure 18.7 ");
00490
00491 PutLine(TauLines[ipTS34],
00492 " S III fine structure 34 ");
00493
00494 PutLine(TauLines[ipTS1720],
00495 " S III] 1713.12, 1728.94 ");
00496
00497 PntForLine(9532.,"S 3",&ipnt);
00498 lindst(CoolHeavy.c9532/(1.+1./2.48),9532,"S 3",ipnt,'c',true ,
00499 " [S III] 9532 alone ");
00500
00501 PntForLine(9069.,"S 3",&ipnt);
00502 lindst(CoolHeavy.c9532/(1.+2.48),9069,"S 3",ipnt,'c',true,
00503 " [S III] 9069 alone" );
00504
00505 PntForLine(6312.,"S 3",&ipnt);
00506 lindst(CoolHeavy.c6312,6312,"S 3",ipnt,'c',true ,
00507 " [S III] 6312, trans-auroral temperature sensitive ");
00508
00509 PntForLine(3722.,"S 3",&ipnt);
00510 lindst(CoolHeavy.c6312*0.59,3722,"S 3",ipnt,'c',true,
00511 " [S III] 3722, same upper level as 6312" );
00512
00513 PutLine(TauLines[ipT1194],
00514 " WL, other data, from Ho + Henry Ap.J. 1984 ");
00515
00516 PutLine(TauLines[ipTS11],
00517 " S IV 10.5 micron, collisionally excited (label is 105) ");
00518
00519 linadd(
00520 TauLines[ipS4_1405].Emis().xIntensity()+
00521 TauLines[ipS4_1398].Emis().xIntensity()+
00522 TauLines[ipS4_1424].Emis().xIntensity()+
00523 TauLines[ipS4_1417].Emis().xIntensity()+
00524 TauLines[ipS4_1407].Emis().xIntensity(),
00525 1406,"TOTL",'i',
00526 " total intensity of S IV] 1406, all lines in the multiplet " );
00527 PutLine(TauLines[ipS4_1405],
00528 " ");
00529 PutLine(TauLines[ipS4_1398],
00530 " ");
00531 PutLine(TauLines[ipS4_1424],
00532 " ");
00533 PutLine(TauLines[ipS4_1417],
00534 " ");
00535 PutLine(TauLines[ipS4_1407],
00536 " ");
00537
00538 linadd(embesq.em1198+TauLines[ipT1198].Emis().xIntensity(),1198,"TOTL",'i',
00539 " S V 1198] both lines together " );
00540
00541 PutLine(TauLines[ipT1198],
00542 " S V 1198] the stronger transition ");
00543
00544 linadd(embesq.em1198,1188,"S 5",'i',
00545 " Be seq, weaker of the two transitions" );
00546
00547 PutLine(TauLines[ipT786],
00548 " S V 786.5, collisionally excited ");
00549
00550 PutLine(TauLines[ipCl1_11m],
00551 " chlorine lines [Cl I] 11 micron ");
00552
00553 PutLine(TauLines[ipfsCl233],
00554 " [Cl II] 33.281 micron ");
00555
00556 PutLine(TauLines[ipfsCl214],
00557 " [Cl II] 14.3678 micron ");
00558
00559 PntForLine(8578.7,"Cl 2",&ipnt);
00560 lindst(CoolHeavy.c8579*0.791,8579,"Cl 2",ipnt,'c',true,
00561 " Chlorine II 8578.7, 9123.6 doublet");
00562
00563 PntForLine(9123.6,"Cl 2",&ipnt);
00564 lindst(CoolHeavy.c8579*0.209,9124,"Cl 2",ipnt,'c',true,
00565 " Chlorine II 8578.7, 9123.6 doublet" );
00566
00567 PntForLine(6161.8,"Cl 2",&ipnt);
00568 lindst(CoolHeavy.c6164,6162,"Cl 2",ipnt,'c',true,
00569 " Chlorine II 6161.8 auroral line >>chng 03 feb 24, change wavelength from 6164 to correct 6161.8 " );
00570
00571 PntForLine(3677.9,"Cl 2",&ipnt);
00572 lindst(CoolHeavy.c3679,3678,"Cl 2",ipnt,'c',true,
00573 " Chlorine II 3677.9 auroral line >>chng 03 feb 24, to correct wavelength " );
00574
00575 linadd(CoolHeavy.c5525,5525,"TOTL",'i',
00576 " Cl III 5519, 5539 doublet, both together " );
00577
00578 linadd(CoolHeavy.c3350,3350,"TOTL",'i',
00579 " Cl III 3354, 3344 doublet, both together " );
00580
00581 linadd(CoolHeavy.c8494,8494,"TOTL",'i',
00582 " Cl III 8504, 8436, 8552, 8483 multiplet, all together " );
00583
00584 PntForLine(5538.,"Cl 3",&ipnt);
00585 lindst(CoolHeavy.Cl5539,5538,"Cl 3",ipnt,'c',true,
00586 " Cl III 5538 " );
00587
00588 PntForLine(5518.,"Cl 3",&ipnt);
00589 lindst(CoolHeavy.Cl5519,5518,"Cl 3",ipnt,'c',true,
00590 " Cl III 5518" );
00591
00592 PntForLine(3354.,"Cl 3",&ipnt);
00593 lindst(CoolHeavy.Cl3354,3354,"Cl 3",ipnt,'c',true,
00594 " Cl III 3354 " );
00595
00596 PntForLine(3344.,"Cl 3",&ipnt);
00597 lindst(CoolHeavy.Cl3344,3344,"Cl 3",ipnt,'c',true,
00598 " Cl III 3344 " );
00599
00600 PntForLine(8504.,"Cl 3",&ipnt);
00601 lindst(CoolHeavy.Cl8504,8504,"Cl 3",ipnt,'c',true ,
00602 " Cl III 8504 ");
00603
00604 PntForLine(8436.,"Cl 3",&ipnt);
00605 lindst(CoolHeavy.Cl8436,8436,"Cl 3",ipnt,'c',true,
00606 " Cl III 8436" );
00607
00608 PntForLine(8552.,"Cl 3",&ipnt);
00609 lindst(CoolHeavy.Cl8552,8552,"Cl 3",ipnt,'c',true,
00610 " Cl III 8552 " );
00611
00612 PntForLine(8483.,"Cl 3",&ipnt);
00613 lindst(CoolHeavy.Cl8483,8483,"Cl 3",ipnt,'c',true,
00614 " Cl III 8483" );
00615
00616 PutLine(TauLines[ipCl04203],
00617 " [Cl IV] fine structure line 20.354 microns");
00618
00619 PutLine(TauLines[ipCl04117],
00620 " [Cl IV] fine structure line 11.741 microns ");
00621
00622 PntForLine(8047.,"Cl 4",&ipnt);
00623 lindst(CoolHeavy.c8047*0.667,8047,"Cl 4",ipnt,'c',true,
00624 " ClIV 8047" );
00625
00626 PntForLine(7532.,"Cl 4",&ipnt);
00627 lindst(CoolHeavy.c8047*0.333,7532,"Cl 4",ipnt,'c',true,
00628 " ClIV 7532" );
00629
00630 PntForLine(3119.,"Cl 4",&ipnt);
00631 lindst(CoolHeavy.c3119,3119,"Cl 4",ipnt,'c',true,
00632 " ClIV 3119" );
00633
00634 PntForLine(5324.,"Cl 4",&ipnt);
00635 lindst(CoolHeavy.c5324,5324,"Cl 4",ipnt,'c',true,
00636 " ClIV 5324" );
00637
00638 PutLine(TauLines[ipCl973],
00639 " Cl IX 7334A ");
00640
00641 PutLine(TauLines[ipTAr7],
00642 " Argon II 7 micron ");
00643
00644 PntForLine(7135.,"Ar 3",&ipnt);
00645 lindst(CoolHeavy.c7136/(1.+1./4.144),7135,"Ar 3",ipnt,'c',true ,
00646 " Argon III 7135");
00647
00648 PntForLine(7751.,"Ar 3",&ipnt);
00649 lindst(CoolHeavy.c7136/(1.+4.144),7751,"Ar 3",ipnt,'c',true,
00650 " Argon III 7751" );
00651
00652 PntForLine(5192.,"Ar 3",&ipnt);
00653 lindst(CoolHeavy.c5192,5192,"Ar 3",ipnt,'c',true,
00654 " Argon III 5192" );
00655
00656 PntForLine(3109.,"Ar 3",&ipnt);
00657 lindst(CoolHeavy.c3109*0.9894,3109,"Ar 3",ipnt,'c',true,
00658 " Argon III 3109" );
00659
00660 PntForLine(3005.,"Ar 3",&ipnt);
00661 lindst(CoolHeavy.c3109*(1.-0.9894),3005,"Ar 3",ipnt,'c',true ,
00662 " Argon III 3005 ");
00663
00664 PutLine(TauLines[ipTAr22],
00665 " Argon III 21.8, 9 micron lines");
00666
00667 PutLine(TauLines[ipTAr9],
00668 " Argon III 21.8, 9 micron lines ");
00669
00670 linadd(CoolHeavy.Ar4740+CoolHeavy.Ar4711,4725,"TOTL",'i',
00671 " Argon IV 4711 + 4740 together, 4740=90%" );
00672
00673 linadd(CoolHeavy.Ar2868+CoolHeavy.Ar2854,2860,"TOTL",'i',
00674 " [AvIV] 2868, 2854 together " );
00675
00676 linadd(CoolHeavy.Ar7237+CoolHeavy.Ar7331+CoolHeavy.Ar7171+CoolHeavy.Ar7263,7250,"TOTL",'i',
00677 " [AvIV] auroral lines, 7237, 7331, 7171, 7263 " );
00678
00679 PntForLine(4740.,"Ar 4",&ipnt);
00680 lindst(CoolHeavy.Ar4740,4740,"Ar 4",ipnt,'c',true,
00681 " [Ar IV] 4740" );
00682
00683 PntForLine(4711.,"Ar 4",&ipnt);
00684 lindst(CoolHeavy.Ar4711,4711,"Ar 4",ipnt,'c',true,
00685 " [Ar IV] 4711" );
00686
00687 PntForLine(2868.,"Ar 4",&ipnt);
00688 lindst(CoolHeavy.Ar2868,2868,"Ar 4",ipnt,'c',true,
00689 " [Ar IV] 2868" );
00690
00691 PntForLine(2854.,"Ar 4",&ipnt);
00692 lindst(CoolHeavy.Ar2854,2854,"Ar 4",ipnt,'c',true,
00693 " [Ar IV] 2854" );
00694
00695
00696 PntForLine(7263.,"Ar 4",&ipnt);
00697 lindst(CoolHeavy.Ar7263,7263,"Ar 4",ipnt,'c',true,
00698 " [Ar IV] 7263" );
00699
00700 PntForLine(7171.,"Ar 4",&ipnt);
00701 lindst(CoolHeavy.Ar7171,7171,"Ar 4",ipnt,'c',true,
00702 " [Ar IV] 7171" );
00703
00704 PntForLine(7331.,"Ar 4",&ipnt);
00705 lindst(CoolHeavy.Ar7331,7331,"Ar 4",ipnt,'c',true,
00706 " [Ar IV] 7331" );
00707
00708 PntForLine(7237.,"Ar 4",&ipnt);
00709 lindst(CoolHeavy.Ar7237,7237,"Ar 4",ipnt,'c',true,
00710 " [Ar IV] 7237" );
00711
00712 PntForLine(7005.,"Ar 5",&ipnt);
00713 lindst(CoolHeavy.c7007/(1.+1./2.143),7005,"Ar 5",ipnt,'c',true,
00714 " Argon V, 3P lines, 7005, collisionally excited" );
00715
00716 PntForLine(6435.,"Ar 5",&ipnt);
00717 lindst(CoolHeavy.c7007/(1.+2.143),6435,"Ar 5",ipnt,'c',true ,
00718 " Argon V, 3P lines, 6435, collisionally excited");
00719
00720 PntForLine(4626.,"Ar 5",&ipnt);
00721 lindst( CoolHeavy.c4626 , 4626 ,"Ar 5",ipnt,'c',true,
00722 " >>chng 01 mar 10, add following two lines Argon V, 4626" );
00723
00724 PntForLine(2691.,"Ar 5",&ipnt);
00725 lindst( CoolHeavy.c2691 , 2691 ,"Ar 5",ipnt,'c',true ,
00726 " Argon V, 2691");
00727
00728 PutLine(TauLines[ipTAr13],
00729 " Argon V fine structure lines, 13.09, 7.903 micron line ");
00730
00731 PutLine(TauLines[ipTAr8],
00732 " Argon V fine structure lines, 13.09, 7.903 micron line ");
00733
00734 PutLine(TauLines[ipAr06453],
00735 " [Ar VI] 4.53 micron ");
00736
00737 if( trace.lgTrace )
00738 {
00739 fprintf( ioQQQ, " lines_lv1_na_ar returns\n" );
00740 }
00741 return;
00742 }
