cloudy trunk
|
00001 /* This file is part of Cloudy and is copyright (C)1978-2008 by Gary J. Ferland and 00002 * others. For conditions of distribution and use see copyright notice in license.txt */ 00003 /*lines_lv1_na_ar place lines of elements sodium through argon into lines storage stack */ 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 /* sum of inward fracs of lines */ 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 ""); 00104 00105 PutLine(&TauLines[ipT2804], 00106 ""); 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 /* SI II 1260, rough guess of dielec contribution */ 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,'i',true, 00291 " fac = (1.-TauLines[ipT1260].ColOvTot);" ); 00292 00293 00294 PntForLine(1909.,"Si 2",&ipnt); 00295 lindst(drec*1260./1909.,1909,"diel",ipnt,'i',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 /* >>chng 96 jul 8, added dielectronic recombination contribution to 1207 */ 00302 /*rec *= dense.xIonDense[13][3]*dense.eden*(1.-TauLines[ipT1207].ColOvTot)**/ 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,'i',true , 00321 " Si III] 1883 by itself"); 00322 00323 /*fac = (1.-TauLines[ipT1895].ColOvTot);*/ 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",'c', 00460 " S II 6731 + 6716 together " ); 00461 00462 linadd(CoolHeavy.S4070+CoolHeavy.S4078,4074, 00463 "S 2",'c',"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,'i',true, 00471 " now output the real thing individual line from five level atom" ); 00472 00473 PntForLine(6716.,"S II",&ipnt); 00474 lindst(CoolHeavy.S6718,6716,"S II",ipnt,'i',true, 00475 " individual line from five level atom" ); 00476 00477 PntForLine(4070.,"S II",&ipnt); 00478 lindst(CoolHeavy.S4070,4070,"S II",ipnt,'i',true, 00479 " individual line from five level atom " ); 00480 00481 PntForLine(4078.,"S II",&ipnt); 00482 lindst(CoolHeavy.S4078,4078,"S II",ipnt,'i',true, 00483 " individual line from five level atom" ); 00484 00485 linadd(CoolHeavy.S10323,10323,"S II",'i', 00486 " individual line from five level atom " ); 00487 00488 linadd(CoolHeavy.S10289,10289,"S II",'i', 00489 " individual line from five level atom" ); 00490 00491 linadd(CoolHeavy.S10373,10373,"S II",'i' , 00492 " individual line from five level atom "); 00493 00494 linadd(CoolHeavy.S10339,10339,"S II",'i', 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 /* group of four intercombination lines 00600 PntForLine(520.,"S 12",&ipnt); 00601 lindst(CoolHeavy.c520,520,"S 12",ipnt,'c',true ); */ 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",'c', 00645 " Cl III 5519, 5539 doublet, both together " ); 00646 00647 linadd(CoolHeavy.c3350,3350,"TOTL",'c', 00648 " Cl III 3354, 3344 doublet, both together " ); 00649 00650 linadd(CoolHeavy.c8494,8494,"TOTL",'c', 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,'i',true, 00655 " Cl III 5538 " ); 00656 00657 PntForLine(5518.,"Cl 3",&ipnt); 00658 lindst(CoolHeavy.Cl5519,5518,"Cl 3",ipnt,'i',true, 00659 " Cl III 5518" ); 00660 00661 PntForLine(3354.,"Cl 3",&ipnt); 00662 lindst(CoolHeavy.Cl3354,3354,"Cl 3",ipnt,'i',true, 00663 " Cl III 3354 " ); 00664 00665 PntForLine(3344.,"Cl 3",&ipnt); 00666 lindst(CoolHeavy.Cl3344,3344,"Cl 3",ipnt,'i',true, 00667 " Cl III 3344 " ); 00668 00669 PntForLine(8504.,"Cl 3",&ipnt); 00670 lindst(CoolHeavy.Cl8504,8504,"Cl 3",ipnt,'i',true , 00671 " Cl III 8504 "); 00672 00673 PntForLine(8436.,"Cl 3",&ipnt); 00674 lindst(CoolHeavy.Cl8436,8436,"Cl 3",ipnt,'i',true, 00675 " Cl III 8436" ); 00676 00677 PntForLine(8552.,"Cl 3",&ipnt); 00678 lindst(CoolHeavy.Cl8552,8552,"Cl 3",ipnt,'i',true, 00679 " Cl III 8552 " ); 00680 00681 PntForLine(8483.,"Cl 3",&ipnt); 00682 lindst(CoolHeavy.Cl8483,8483,"Cl 3",ipnt,'i',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 }