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*ionbal.DielSupprs[1][2];
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*ionbal.DielSupprs[1][2];
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[ipP848],
00438 " [P VIII] 4.85 micron ");
00439
00440 PutLine(&TauLines[ipP817],
00441 " [P VIII] 1.735 microns ");
00442
00443 PutLine(&TauLines[ipP1027],
00444 " [P X] 2.708 micron ");
00445
00446 PutLine(&TauLines[ipP1018],
00447 " [P X] 1.868 microns ");
00448
00449 PutLine(&TauLines[ipS1_25m],
00450 " sulphur S I 25m ");
00451
00452 PutLine(&TauLines[ipS1_56m],
00453 " S I 56m ");
00454
00455 linadd(1.75e-22*dense.eden*dense.xIonDense[ipSULPHUR][1]/phycon.sqrte/
00456 phycon.te10/phycon.te03,1807,"S 1R",'i',
00457 " 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 " );
00458
00459 linadd(CoolHeavy.c6731,6720,"S 2",'i',
00460 " S II 6731 + 6716 together " );
00461
00462 linadd(CoolHeavy.S4070+CoolHeavy.S4078,4074,
00463 "S 2",'i',"S II 4070 +4078 together" );
00464
00465 PntForLine(10330.,"S 2",&ipnt);
00466 lindst(CoolHeavy.c10330,10330,"S 2",ipnt,'c',true,
00467 " S II N=3 lines, all four lines together " );
00468
00469 PntForLine(6731.,"S II",&ipnt);
00470 lindst(CoolHeavy.S6733,6731,"S II",ipnt,'c',true,
00471 " individual line from five level atom" );
00472
00473 PntForLine(6716.,"S II",&ipnt);
00474 lindst(CoolHeavy.S6718,6716,"S II",ipnt,'c',true,
00475 " individual line from five level atom" );
00476
00477 PntForLine(4070.,"S II",&ipnt);
00478 lindst(CoolHeavy.S4070,4070,"S II",ipnt,'c',true,
00479 " individual line from five level atom " );
00480
00481 PntForLine(4078.,"S II",&ipnt);
00482 lindst(CoolHeavy.S4078,4078,"S II",ipnt,'c',true,
00483 " individual line from five level atom" );
00484
00485 linadd(CoolHeavy.S10323,10323,"S II",'c',
00486 " individual line from five level atom " );
00487
00488 linadd(CoolHeavy.S10289,10289,"S II",'c',
00489 " individual line from five level atom" );
00490
00491 linadd(CoolHeavy.S10373,10373,"S II",'c' ,
00492 " individual line from five level atom ");
00493
00494 linadd(CoolHeavy.S10339,10339,"S II",'c',
00495 " individual line from five level atom " );
00496
00497 PutLine(&TauLines[ipT1256],
00498 " resonance line near NV, collisionally excited ");
00499
00500 PutLine(&TauLines[ipTS19],
00501 " S III fine structure 18.7 ");
00502
00503 PutLine(&TauLines[ipTS34],
00504 " S III fine structure 34 ");
00505
00506 PutLine(&TauLines[ipTS1720],
00507 " S III] 1713.12, 1728.94 ");
00508
00509 PntForLine(9532.,"S 3",&ipnt);
00510 lindst(CoolHeavy.c9532/(1.+1./2.48),9532,"S 3",ipnt,'c',true ,
00511 " [S III] 9532 alone ");
00512
00513 PntForLine(9069.,"S 3",&ipnt);
00514 lindst(CoolHeavy.c9532/(1.+2.48),9069,"S 3",ipnt,'c',true,
00515 " [S III] 9069 alone" );
00516
00517 PntForLine(6312.,"S 3",&ipnt);
00518 lindst(CoolHeavy.c6312,6312,"S 3",ipnt,'c',true ,
00519 " [S III] 6312, trans-auroral temperature sensitive ");
00520
00521 PntForLine(3722.,"S 3",&ipnt);
00522 lindst(CoolHeavy.c6312*0.59,3722,"S 3",ipnt,'c',true,
00523 " [S III] 3722, same upper level as 6312" );
00524
00525 PutLine(&TauLines[ipT1194],
00526 " WL, other data, from Ho + Henry Ap.J. 1984 ");
00527
00528 PutLine(&TauLines[ipTS11],
00529 " S IV 10.5 micron, collisionally excited (label is 105) ");
00530
00531 linadd(
00532 TauLines[ipS4_1405].Emis->xIntensity+
00533 TauLines[ipS4_1398].Emis->xIntensity+
00534 TauLines[ipS4_1424].Emis->xIntensity+
00535 TauLines[ipS4_1417].Emis->xIntensity+
00536 TauLines[ipS4_1407].Emis->xIntensity,
00537 1406,"TOTL",'i',
00538 " total intensity of S IV] 1406, all lines in the multiplet " );
00539 PutLine(&TauLines[ipS4_1405],
00540 " ");
00541 PutLine(&TauLines[ipS4_1398],
00542 " ");
00543 PutLine(&TauLines[ipS4_1424],
00544 " ");
00545 PutLine(&TauLines[ipS4_1417],
00546 " ");
00547 PutLine(&TauLines[ipS4_1407],
00548 " ");
00549
00550 linadd(embesq.em1198+TauLines[ipT1198].Emis->xIntensity,1198,"TOTL",'i',
00551 " S V 1198] both lines together " );
00552
00553 PutLine(&TauLines[ipT1198],
00554 " S V 1198] the stronger transition ");
00555
00556 linadd(embesq.em1198,1188,"S 5",'i',
00557 " Be seq, weaker of the two transitions" );
00558
00559 PutLine(&TauLines[ipT786],
00560 " S V 786.5, collisionally excited ");
00561
00562 linadd(TauLines[ipT944].Emis->xIntensity+TauLines[ipT933].Emis->xIntensity,933,"TOTL",'i',
00563 " S VI 933+944, collisionally excited total S VI 933+944 " );
00564 PutLine(&TauLines[ipT944],
00565 " ");
00566 PutLine(&TauLines[ipT933],
00567 " ");
00568
00569 PutLine(&TauLines[ipfsS810],
00570 " S VIII 9913 ");
00571
00572 PutLine(&TauLines[ipfsS912],
00573 " S IX 1.25 microns ");
00574
00575 PutLine(&TauLines[ipfsS938],
00576 " S IX 3.78 microns ");
00577
00578 PutLine(&TauLines[ipfsS1119],
00579 " S XI 1.91 microns ");
00580
00581 PutLine(&TauLines[ipfsS1114],
00582 " S XI 1.40 microns ");
00583
00584 PntForLine(1715.,"S 9",&ipnt);
00585 lindst(CoolHeavy.c1715,1715,"S 9",ipnt,'c',true ,
00586 " S IX 1715, 1987, collisionally excited");
00587
00588 PntForLine(1213.,"S 10",&ipnt);
00589 lindst(CoolHeavy.c1213,1213,"S 10",ipnt,'c',true,
00590 " S X 1213, 1197, collisionally excited " );
00591
00592 PntForLine(1826.,"S 11",&ipnt);
00593 lindst(CoolHeavy.c1826,1826,"S 11",ipnt,'c',true,
00594 " S XI 1615, 1826, collisionally excited" );
00595
00596 PutLine(&TauLines[ipfsS1207],
00597 " S XII 7611.2 ");
00598
00599
00600
00601
00602 PutLine(&TauLines[ipS12_520],
00603 " >>chng 03 sep 27, move to atom_level2 from expion");
00604
00605 PntForLine(488.,"S 13",&ipnt);
00606 lindst(CoolHeavy.c488,488,"S 13",ipnt,'c',true,
00607 " S XIII 488.4, 1909 like, collisionally excited" );
00608
00609 linadd(TauLines[ipTSu418].Emis->xIntensity+TauLines[ipTSu446].Emis->xIntensity,427,"TOTL",'i',
00610 " S 14 506 li seq 2s2p " );
00611 PutLine(&TauLines[ipTSu418],
00612 " ");
00613 PutLine(&TauLines[ipTSu446],
00614 " ");
00615
00616 PutLine(&TauLines[ipTSu30],
00617 " S 14 30 li seq 2s3p ");
00618
00619 PutLine(&TauLines[ipCl1_11m],
00620 " chlorine lines [Cl I] 11 micron ");
00621
00622 PutLine(&TauLines[ipfsCl233],
00623 " [Cl II] 33.281 micron ");
00624
00625 PutLine(&TauLines[ipfsCl214],
00626 " [Cl II] 14.3678 micron ");
00627
00628 PntForLine(8578.7,"Cl 2",&ipnt);
00629 lindst(CoolHeavy.c8579*0.791,8579,"Cl 2",ipnt,'c',true,
00630 " Chlorine II 8578.7, 9123.6 doublet");
00631
00632 PntForLine(9123.6,"Cl 2",&ipnt);
00633 lindst(CoolHeavy.c8579*0.209,9124,"Cl 2",ipnt,'c',true,
00634 " Chlorine II 8578.7, 9123.6 doublet" );
00635
00636 PntForLine(6161.8,"Cl 2",&ipnt);
00637 lindst(CoolHeavy.c6164,6162,"Cl 2",ipnt,'c',true,
00638 " Chlorine II 6161.8 auroral line >>chng 03 feb 24, change wavelength from 6164 to correct 6161.8 " );
00639
00640 PntForLine(3677.9,"Cl 2",&ipnt);
00641 lindst(CoolHeavy.c3679,3678,"Cl 2",ipnt,'c',true,
00642 " Chlorine II 3677.9 auroral line >>chng 03 feb 24, to correct wavelength " );
00643
00644 linadd(CoolHeavy.c5525,5525,"TOTL",'i',
00645 " Cl III 5519, 5539 doublet, both together " );
00646
00647 linadd(CoolHeavy.c3350,3350,"TOTL",'i',
00648 " Cl III 3354, 3344 doublet, both together " );
00649
00650 linadd(CoolHeavy.c8494,8494,"TOTL",'i',
00651 " Cl III 8504, 8436, 8552, 8483 multiplet, all together " );
00652
00653 PntForLine(5538.,"Cl 3",&ipnt);
00654 lindst(CoolHeavy.Cl5539,5538,"Cl 3",ipnt,'c',true,
00655 " Cl III 5538 " );
00656
00657 PntForLine(5518.,"Cl 3",&ipnt);
00658 lindst(CoolHeavy.Cl5519,5518,"Cl 3",ipnt,'c',true,
00659 " Cl III 5518" );
00660
00661 PntForLine(3354.,"Cl 3",&ipnt);
00662 lindst(CoolHeavy.Cl3354,3354,"Cl 3",ipnt,'c',true,
00663 " Cl III 3354 " );
00664
00665 PntForLine(3344.,"Cl 3",&ipnt);
00666 lindst(CoolHeavy.Cl3344,3344,"Cl 3",ipnt,'c',true,
00667 " Cl III 3344 " );
00668
00669 PntForLine(8504.,"Cl 3",&ipnt);
00670 lindst(CoolHeavy.Cl8504,8504,"Cl 3",ipnt,'c',true ,
00671 " Cl III 8504 ");
00672
00673 PntForLine(8436.,"Cl 3",&ipnt);
00674 lindst(CoolHeavy.Cl8436,8436,"Cl 3",ipnt,'c',true,
00675 " Cl III 8436" );
00676
00677 PntForLine(8552.,"Cl 3",&ipnt);
00678 lindst(CoolHeavy.Cl8552,8552,"Cl 3",ipnt,'c',true,
00679 " Cl III 8552 " );
00680
00681 PntForLine(8483.,"Cl 3",&ipnt);
00682 lindst(CoolHeavy.Cl8483,8483,"Cl 3",ipnt,'c',true,
00683 " Cl III 8483" );
00684
00685 PutLine(&TauLines[ipCl04203],
00686 " [Cl IV] fine structure line 20.354 microns");
00687
00688 PutLine(&TauLines[ipCl04117],
00689 " [Cl IV] fine structure line 11.741 microns ");
00690
00691 PntForLine(8047.,"Cl 4",&ipnt);
00692 lindst(CoolHeavy.c8047*0.667,8047,"Cl 4",ipnt,'c',true,
00693 " ClIV 8047" );
00694
00695 PntForLine(7532.,"Cl 4",&ipnt);
00696 lindst(CoolHeavy.c8047*0.333,7532,"Cl 4",ipnt,'c',true,
00697 " ClIV 7532" );
00698
00699 PntForLine(3119.,"Cl 4",&ipnt);
00700 lindst(CoolHeavy.c3119,3119,"Cl 4",ipnt,'c',true,
00701 " ClIV 3119" );
00702
00703 PntForLine(5324.,"Cl 4",&ipnt);
00704 lindst(CoolHeavy.c5324,5324,"Cl 4",ipnt,'c',true,
00705 " ClIV 5324" );
00706
00707 PutLine(&TauLines[ipCl973],
00708 " Cl IX 7334A ");
00709
00710 PutLine(&TauLines[ipCl1030],
00711 " Cl X 3.05 microns ");
00712
00713 PutLine(&TauLines[ipCl1092],
00714 " Cl X 9223A ");
00715
00716 PutLine(&TauLines[ipTAr7],
00717 " Argon II 7 micron ");
00718
00719 PntForLine(7135.,"Ar 3",&ipnt);
00720 lindst(CoolHeavy.c7136/(1.+1./4.144),7135,"Ar 3",ipnt,'c',true ,
00721 " Argon III 7135");
00722
00723 PntForLine(7751.,"Ar 3",&ipnt);
00724 lindst(CoolHeavy.c7136/(1.+4.144),7751,"Ar 3",ipnt,'c',true,
00725 " Argon III 7751" );
00726
00727 PntForLine(5192.,"Ar 3",&ipnt);
00728 lindst(CoolHeavy.c5192,5192,"Ar 3",ipnt,'c',true,
00729 " Argon III 5192" );
00730
00731 PntForLine(3109.,"Ar 3",&ipnt);
00732 lindst(CoolHeavy.c3109*0.9894,3109,"Ar 3",ipnt,'c',true,
00733 " Argon III 3109" );
00734
00735 PntForLine(3005.,"Ar 3",&ipnt);
00736 lindst(CoolHeavy.c3109*(1.-0.9894),3005,"Ar 3",ipnt,'c',true ,
00737 " Argon III 3005 ");
00738
00739 PutLine(&TauLines[ipTAr22],
00740 " Argon III 21.8, 9 micron lines");
00741
00742 PutLine(&TauLines[ipTAr9],
00743 " Argon III 21.8, 9 micron lines ");
00744
00745 linadd(CoolHeavy.Ar4740+CoolHeavy.Ar4711,4725,"TOTL",'i',
00746 " Argon IV 4711 + 4740 together, 4740=90%" );
00747
00748 linadd(CoolHeavy.Ar2868+CoolHeavy.Ar2854,2860,"TOTL",'i',
00749 " [AvIV] 2868, 2854 together " );
00750
00751 linadd(CoolHeavy.Ar7237+CoolHeavy.Ar7331+CoolHeavy.Ar7171+CoolHeavy.Ar7263,7250,"TOTL",'i',
00752 " [AvIV] auroral lines, 7237, 7331, 7171, 7263 " );
00753
00754 PntForLine(4740.,"Ar 4",&ipnt);
00755 lindst(CoolHeavy.Ar4740,4740,"Ar 4",ipnt,'c',true,
00756 " [Ar IV] 4740" );
00757
00758 PntForLine(4711.,"Ar 4",&ipnt);
00759 lindst(CoolHeavy.Ar4711,4711,"Ar 4",ipnt,'c',true,
00760 " [Ar IV] 4711" );
00761
00762 PntForLine(2868.,"Ar 4",&ipnt);
00763 lindst(CoolHeavy.Ar2868,2868,"Ar 4",ipnt,'c',true,
00764 " [Ar IV] 2868" );
00765
00766 PntForLine(2854.,"Ar 4",&ipnt);
00767 lindst(CoolHeavy.Ar2854,2854,"Ar 4",ipnt,'c',true,
00768 " [Ar IV] 2854" );
00769
00770
00771 PntForLine(7263.,"Ar 4",&ipnt);
00772 lindst(CoolHeavy.Ar7263,7263,"Ar 4",ipnt,'c',true,
00773 " [Ar IV] 7263" );
00774
00775 PntForLine(7171.,"Ar 4",&ipnt);
00776 lindst(CoolHeavy.Ar7171,7171,"Ar 4",ipnt,'c',true,
00777 " [Ar IV] 7171" );
00778
00779 PntForLine(7331.,"Ar 4",&ipnt);
00780 lindst(CoolHeavy.Ar7331,7331,"Ar 4",ipnt,'c',true,
00781 " [Ar IV] 7331" );
00782
00783 PntForLine(7237.,"Ar 4",&ipnt);
00784 lindst(CoolHeavy.Ar7237,7237,"Ar 4",ipnt,'c',true,
00785 " [Ar IV] 7237" );
00786
00787 PntForLine(7005.,"Ar 5",&ipnt);
00788 lindst(CoolHeavy.c7007/(1.+1./2.143),7005,"Ar 5",ipnt,'c',true,
00789 " Argon V, 3P lines, 7005, collisionally excited" );
00790
00791 PntForLine(6435.,"Ar 5",&ipnt);
00792 lindst(CoolHeavy.c7007/(1.+2.143),6435,"Ar 5",ipnt,'c',true ,
00793 " Argon V, 3P lines, 6435, collisionally excited");
00794
00795 PntForLine(4626.,"Ar 5",&ipnt);
00796 lindst( CoolHeavy.c4626 , 4626 ,"Ar 5",ipnt,'c',true,
00797 " >>chng 01 mar 10, add following two lines Argon V, 4626" );
00798
00799 PntForLine(2691.,"Ar 5",&ipnt);
00800 lindst( CoolHeavy.c2691 , 2691 ,"Ar 5",ipnt,'c',true ,
00801 " Argon V, 2691");
00802
00803 PutLine(&TauLines[ipTAr13],
00804 " Argon V fine structure lines, 13.09, 7.903 micron line ");
00805
00806 PutLine(&TauLines[ipTAr8],
00807 " Argon V fine structure lines, 13.09, 7.903 micron line ");
00808
00809 PutLine(&TauLines[ipAr06453],
00810 " [Ar VI] 4.53 micron ");
00811
00812 PutLine(&TauLines[ipAr1055],
00813 " Ar X 5533.4A ");
00814
00815 PutLine(&TauLines[ipAr1126],
00816 " Ar XI 2.60 microns ");
00817
00818 PutLine(&TauLines[ipAr1178],
00819 " Ar XI 6917A ");
00820
00821 PntForLine(4413.,"Ar14",&ipnt);
00822 lindst(CoolHeavy.fs4413,4413,"Ar14",ipnt,'c',true,
00823 " Ar XIV 4413, predicted lambda, not observed(?)" );
00824
00825 PntForLine(409.,"Ar15",&ipnt);
00826 lindst(CoolHeavy.c409,409,"Ar15",ipnt,'c',true,
00827 " collisionally excited " );
00828
00829 PutLine(&TauLines[ipT354],
00830 " collisionally excited Ar16 354, Li seq 2s-2p");
00831
00832 PutLine(&TauLines[ipT389],
00833 " collisionally excited Ar16 389, Li seq 2s-2p ");
00834
00835 PutLine(&TauLines[ipT25],
00836 " collisionally excited Ar16 365, Li seq 2s-3p ");
00837
00838 if( trace.lgTrace )
00839 {
00840 fprintf( ioQQQ, " lines_lv1_na_ar returns\n" );
00841 }
00842 return;
00843 }