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_k_zn place lines of elements potassium and heavier into lines storage stack */ 00004 #include "cddefines.h" 00005 #include "cddrive.h" 00006 #include "coolheavy.h" 00007 #include "ca.h" 00008 #include "fe.h" 00009 #include "rfield.h" 00010 #include "dense.h" 00011 #include "phycon.h" 00012 #include "radius.h" 00013 #include "taulines.h" 00014 #include "trace.h" 00015 #include "lines_service.h" 00016 #include "rt.h" 00017 #include "atomfeii.h" 00018 #include "lines.h" 00019 00020 void lines_lv1_k_zn(void) 00021 { 00022 long int i, 00023 ipnt, 00024 ilo, 00025 ihi; 00026 00027 double c10, 00028 c14, 00029 eff, 00030 fela, 00031 r14; 00032 00033 DEBUG_ENTRY( "lines_lv1_k_zn()" ); 00034 00035 if( trace.lgTrace ) 00036 { 00037 fprintf( ioQQQ, " lines_lv1_k_zn called\n" ); 00038 } 00039 00040 PutLine(&TauLines[ipKI7745], 00041 " potassium K I 7745 "); 00042 00043 PutLine(&TauLines[ipxK03462], 00044 " [K III] 4.62 microns "); 00045 00046 PutLine(&TauLines[ipxK04598], 00047 " [KIV] 5.983 min "); 00048 00049 PutLine(&TauLines[ipxK04154], 00050 " [KIV] 15.39 micron "); 00051 00052 PutLine(&TauLines[ipxK06882], 00053 " [KVI] 8.823 micron "); 00054 00055 PutLine(&TauLines[ipxK06557], 00056 " [KVI] 5.575 micron "); 00057 00058 PutLine(&TauLines[ipxK07319], 00059 " [K VII] 3.189 microns "); 00060 00061 PutLine(&TauLines[ipxK11425], 00062 " K 11 4249.99A "); 00063 00064 PutLine(&TauLines[ipCaI4228], 00065 " calcium Ca I 4228 "); 00066 00067 linadd(ca.Cakh,3933,"Ca 2",'c', 00068 " coll excited calcium k+h " ); 00069 00070 linadd(ca.Cair,8579,"Ca 2",'c' , 00071 " infrared triplet "); 00072 00073 linadd(ca.c7306,7306,"Ca 2",'c', 00074 " forbidden lines, 7291+7324 together " ); 00075 00076 linadd(ca.dCakh,3933,"Phot",'i' , 00077 " fraction H Ly-alpha destruction of excited levels "); 00078 00079 linadd(ca.dCaf12,7306,"Phot",'i' , 00080 " fraction H Ly-alpha destruction of excited levels "); 00081 00082 PntForLine(3934.,"Ca2K",&ipnt); 00083 lindst(ca.Cak,3934,"Ca2K",ipnt,'i',true, 00084 " individual lines from five level atom"); 00085 00086 00087 PntForLine(3969.,"Ca2H",&ipnt); 00088 lindst(ca.Cah,3969,"Ca2H",ipnt,'i',true, 00089 " individual lines from five level atom" ); 00090 00091 00092 PntForLine(8498.,"Ca2X",&ipnt); 00093 lindst(ca.Cax,8498,"Ca2X",ipnt,'i',true, 00094 " individual lines from five level atom " ); 00095 00096 00097 PntForLine(8542.,"Ca2Y",&ipnt); 00098 lindst(ca.Cay,8542,"Ca2Y",ipnt,'i',true, 00099 " individual lines from five level atom" ); 00100 00101 00102 PntForLine(8662.,"Ca2Z",&ipnt); 00103 lindst(ca.Caz,8662,"Ca2Z",ipnt,'i',true, 00104 " individual lines from five level atom" ); 00105 00106 00107 PntForLine(7291.,"CaF1",&ipnt); 00108 lindst(ca.Caf1,7291,"CaF1",ipnt,'i',true, 00109 " individual lines from five level atom" ); 00110 00111 00112 PntForLine(7324.,"CaF2",&ipnt); 00113 lindst(ca.Caf2,7324,"CaF2",ipnt,'i',true, 00114 " individual lines from five level atom" ); 00115 00116 00117 eff = dense.eden*dense.xIonDense[ipCALCIUM][2]*5.4e-21/(phycon.te/ 00118 phycon.te10/phycon.te10); 00119 linadd(eff,3933,"Rec ",'i', 00120 " recombination contribution to CaII emission" ); 00121 00122 PutLine(&TauLines[ipTCa3], 00123 " Ca IV 3.2 micron "); 00124 00125 PutLine(&TauLines[ipTCa4], 00126 " Ca V 4.16, 11.47 micron"); 00127 00128 PutLine(&TauLines[ipTCa12], 00129 " Ca V 4.16, 11.47 micron "); 00130 00131 00132 PntForLine(6087.,"Ca 5",&ipnt); 00133 lindst(ca.Ca6087,6087,"Ca 5",ipnt,'c',true , 00134 " Ca V optical and UV lines, collisional excitation, 3-level atom"); 00135 00136 00137 PntForLine(5311.,"Ca 5",&ipnt); 00138 lindst(ca.c5311,5311,"Ca 5",ipnt,'c',true , 00139 " Ca V optical and UV lines, collisional excitation, 3-level atom"); 00140 00141 00142 PntForLine(2414.,"Ca 5",&ipnt); 00143 lindst(ca.c2414,2414,"Ca 5",ipnt,'c',true , 00144 " Ca V optical and UV lines, collisional excitation, 3-level atom"); 00145 00146 00147 PntForLine(3997.,"Ca 5",&ipnt); 00148 lindst(ca.c3997,3997,"Ca 5",ipnt,'c',true, 00149 " Ca V optical and UV lines, collisional excitation, 3-level atom" ); 00150 00151 PutLine(&TauLines[ipCa0741], 00152 " [Ca VII] 4.09 microns" ); 00153 00154 PutLine(&TauLines[ipCa0761], 00155 " [Ca VII] 6.15 microns " ); 00156 00157 00158 PntForLine(5620.,"Ca 7",&ipnt); 00159 lindst(ca.Ca5620,5620,"Ca 7",ipnt,'c',true, 00160 " Ca VII optical and UV lines, collisional excitation, 3-level atom" ); 00161 00162 00163 PntForLine(4941.,"Ca 7",&ipnt); 00164 lindst(ca.Ca4941,4941,"Ca 7",ipnt,'c',true, 00165 " Ca VII optical and UV lines, collisional excitation, 3-level atom" ); 00166 00167 00168 PntForLine(2112.,"Ca 7",&ipnt); 00169 lindst(ca.Ca2112,2112,"Ca 7",ipnt,'c',true , 00170 " Ca VII optical and UV lines, collisional excitation, 3-level atom"); 00171 00172 00173 PntForLine(3688.,"Ca 7",&ipnt); 00174 lindst(ca.Ca3688,3688,"Ca 7",ipnt,'c',true, 00175 " Ca VII optical and UV lines, collisional excitation, 3-level atom" ); 00176 00177 PutLine(&TauLines[ipCa08232], 00178 " [Ca VIII] 2.32 microns, A Saraph and Strey "); 00179 00180 PutLine(&TauLines[ipCa12333], 00181 " [Ca 12] 3328.78A "); 00182 00183 PutLine(&TauLines[ipTCa302], 00184 " Ca 18 Li seq 2s2p, 302, 345 separate "); 00185 00186 PutLine(&TauLines[ipTCa345], 00187 " Ca 18 Li seq 2s2p, 302, 345 separate "); 00188 00189 PutLine(&TauLines[ipTCa19], 00190 " Ca 18 Li seq 2s3p, "); 00191 00192 00193 PntForLine(22.08e4,"Sc 2",&ipnt); 00194 lindst(CoolHeavy.Sc22p08m,21,"Sc 2",ipnt,'c',true, 00195 " Sc II 2.08 (1-3) " ); 00196 00197 00198 PntForLine(24.1e4,"Sc 2",&ipnt); 00199 lindst(CoolHeavy.Sc24p1m,41,"Sc 2",ipnt,'c',true, 00200 " Sc II 4.1 micron (1-2)" ); 00201 00202 00203 PntForLine(24.2e4,"Sc 2",&ipnt); 00204 lindst(CoolHeavy.Sc24p2m,42,"Sc 2",ipnt,'c',true, 00205 " Sc II 4.22 (2-3)" ); 00206 00207 00208 PntForLine(3933.,"Sc 3",&ipnt); 00209 lindst(CoolHeavy.Sc33936,3933,"Sc 3",ipnt,'c',true, 00210 " Sc III 3936" ); 00211 00212 PutLine(&TauLines[ipSc05231], 00213 " [Sc V] 1.46 microns "); 00214 00215 00216 PntForLine(5054.,"Sc 6",&ipnt); 00217 lindst(CoolHeavy.Sc45058,5054,"Sc 6",ipnt,'c',true , 00218 " Sc VI 5054 (1-2)"); 00219 00220 00221 PntForLine(3592.,"Sc 6",&ipnt); 00222 lindst(CoolHeavy.Sc43595,3592,"Sc 6",ipnt,'c',true, 00223 " Sc VI 3595 (2-3)" ); 00224 00225 00226 PntForLine(2100.,"Sc 6",&ipnt); 00227 lindst(CoolHeavy.Sc42100,2100,"Sc 6",ipnt,'c',true, 00228 " Sc VI 2100 (1-3)" ); 00229 00230 PutLine(&TauLines[ipSc13264], 00231 " [Sc 13] 2637.97A"); 00232 00233 00234 PntForLine(1.21e4,"Ti 3",&ipnt); 00235 lindst(CoolHeavy.Ti3l21,12,"Ti 3",ipnt,'c',true, 00236 " Ti III 1.21 micron, (actually multiplet) 2-1 transition from model atom " ); 00237 00238 00239 PntForLine(9594.,"Ti 3",&ipnt); 00240 lindst(CoolHeavy.Ti3l31,9594,"Ti 3",ipnt,'c',true, 00241 " Ti III 9594, 3-1 transition, (actually multiplet) from model atom" ); 00242 00243 PntForLine(4.57e4,"Ti 3",&ipnt); 00244 lindst(CoolHeavy.Ti3l32,45,"Ti 3",ipnt,'c',true, 00245 " Ti III 4.57 micron, 3-2 transition, (actually multiplet) from model atom" ); 00246 00247 PutLine(&TauLines[ipTi06172], 00248 " [Ti VI] 1.72 microns "); 00249 00250 PutLine(&TauLines[ipTi14212], 00251 " [Ti XIV] 2117.79 "); 00252 00253 00254 PntForLine(8823.,"V 3",&ipnt); 00255 lindst(CoolHeavy.V38830,8823,"V 3",ipnt,'c',true , 00256 " V III 8823 "); 00257 00258 00259 PntForLine(8507.,"V 3",&ipnt); 00260 lindst(CoolHeavy.V38507,8507,"V 3",ipnt,'c',true, 00261 " V III 8507" ); 00262 00263 00264 PntForLine(7735.,"V 4",&ipnt); 00265 lindst(CoolHeavy.V47741,7735,"V 4",ipnt,'c',true, 00266 " V IV 7741 1-3" ); 00267 00268 00269 PntForLine(9489.,"V 4",&ipnt); 00270 lindst(CoolHeavy.V49496,9489,"V 4",ipnt,'c',true, 00271 " V IV 9496 2-1 " ); 00272 00273 00274 PntForLine(4.19e4,"V 4",&ipnt); 00275 lindst(CoolHeavy.V44p2m,42,"V 4",ipnt,'c',true, 00276 " V IV 4.19 micron 3-2" ); 00277 00278 PutLine(&TauLines[ipVa07130], 00279 " [V VII] 1.304 microns "); 00280 00281 PutLine(&TauLines[ipVa15172], 00282 " [V 15] 1721.38 "); 00283 00284 00285 PntForLine(5828.,"Cr 3",&ipnt); 00286 lindst(CoolHeavy.Cr3l21,5828,"Cr 3",ipnt,'c',true, 00287 " [CrIII] multiplet blend at 5828A" ); 00288 00289 PntForLine(7267.,"Cr 4",&ipnt); 00290 lindst(CoolHeavy.Cr4l21,7267,"Cr 4",ipnt,'c',true, 00291 " [CrIV] 2 - 1 multiplet blend at 7272" ); 00292 00293 00294 PntForLine(6801.,"Cr 4",&ipnt); 00295 lindst(CoolHeavy.Cr4l31,6801,"Cr 4",ipnt,'c',true, 00296 " [CrIV] 3 - 1 multiplet blend at 6806" ); 00297 00298 00299 PntForLine(7979.,"Cr 5",&ipnt); 00300 lindst(CoolHeavy.Cr5l21,7979,"Cr 5",ipnt,'c',true, 00301 " [CrV] 2 - 1 multiplet blend at 7985" ); 00302 00303 PntForLine(6577.,"Cr 5",&ipnt); 00304 lindst(CoolHeavy.Cr5l31,6577,"Cr 5",ipnt,'c',true, 00305 " [CrV] 3 - 1 multiplet blend at 6582" ); 00306 00307 00308 PntForLine(3.75e4,"Cr 5",&ipnt); 00309 lindst(CoolHeavy.Cr5l32,37,"Cr 5",ipnt,'c',true, 00310 " [CrV] 3 - 2 multiplet blend at 3.75 microns " ); 00311 00312 PutLine(&TauLines[ipCr08101], 00313 " [Cr VIII] 1.01 microns "); 00314 00315 PutLine(&TauLines[ipCr16141], 00316 " [Cr 16] 1410.60 "); 00317 00318 PutLine(&TauLines[ipxMn0979], 00319 " [Mn IX] 7968.5 A "); 00320 00321 PutLine(&TauLines[ipxMn1712], 00322 " [Mn 17] 1169.59 "); 00323 00324 /* bob Rubin's UV line 00325 * f2 = dense.xIonDense(26,4)*sexp(50 764./te)*0.45*cdsqte/6.*7.01e-12 00326 * call linadd( f2 , 2837 , 'BobR' , 'i') 00327 * f2 = dense.xIonDense(26,4)*sexp(55 989./te)*0.384*cdsqte/6.*7.74e-12 00328 * call linadd( f2 , 2568 , 'BobR' , 'i') */ 00329 00330 /* iron */ 00331 00332 PutLine(&TauLines[ipFe1_24m], 00333 " Fe 1 24m "); 00334 00335 PutLine(&TauLines[ipFe1_35m], 00336 " Fe 1 35m "); 00337 00338 PutLine(&TauLines[ipFe1_54m], 00339 " Fe 1 54m "); 00340 00341 PutLine(&TauLines[ipFe1_111m], 00342 " Fe 1 111m "); 00343 00344 PutLine(&TauLines[ipFeI3884], 00345 " Fe 1 3884 "); 00346 00347 PutLine(&TauLines[ipFeI3729], 00348 " Fe 1 3729 "); 00349 00350 PutLine(&TauLines[ipFeI3457], 00351 " Fe 1 3457 "); 00352 00353 PutLine(&TauLines[ipFeI3021], 00354 " Fe 1 3021 "); 00355 00356 PutLine(&TauLines[ipFeI2966], 00357 " Fe 1 2966 "); 00358 00359 linadd(MAX2(0.,FeII.Fe2_large_cool+FeII.Fe2_UVsimp_cool),0,"Fe2c",'c' , 00360 "total of all Fe 2, both simple UV and large atom together "); 00361 00362 linadd(MAX2(0.,-FeII.Fe2_large_cool-FeII.Fe2_UVsimp_cool),0,"Fe2h",'h' , 00363 "total of all Fe 2 heating, both simple UV and large atom together "); 00364 00365 linadd(FeII.for7,4300,"Fe 2",'i' , 00366 " Fe 2 forbidden 2-1 transition from Netzer's atom "); 00367 00368 PutLine(&TauLines[ipTuv3], 00369 " 2400 "); 00370 PutLine(&TauLines[ipTr48], 00371 " 6200 "); 00372 PutLine(&TauLines[ipTFe16], 00373 " 1080 "); 00374 PutLine(&TauLines[ipTFe26], 00375 " 1500 "); 00376 PutLine(&TauLines[ipTFe34], 00377 " 11500 "); 00378 PutLine(&TauLines[ipTFe35], 00379 " 2500 "); 00380 PutLine(&TauLines[ipTFe46], 00381 " 2300 "); 00382 PutLine(&TauLines[ipTFe56], 00383 " 8900 "); 00384 00385 /* >>chng 03 nov 02, changed wl in lindst calls to be wl in microns, 00386 * this will change wl in printout */ 00387 /* results of FeIILevelPops Fe 2 atom, 00388 * not all are printed, first is */ 00389 /* total cooling due to 16 level atom 00390 * >>chng 05 dec 04, rm this since no longer meaningful 00391 linadd(FeII.Fe2_16levl_cool,0,"Fe 2",'c' );*/ 00392 00393 00394 PntForLine(16640.,"Fe 2",&ipnt); 00395 lindst(FeII.fe21308,16640,"Fe 2",ipnt,'i',true, 00396 " Fe 2 1.664 microns 8-13 " ); 00397 00398 00399 PntForLine(15990.,"Fe 2",&ipnt); 00400 lindst(FeII.fe21207,15990,"Fe 2",ipnt,'i',true, 00401 " Fe 2 1.599 microns 7-12" ); 00402 00403 00404 PntForLine(15340.,"Fe 2",&ipnt); 00405 lindst(FeII.fe21106,15340,"Fe 2",ipnt,'i',true , 00406 " Fe 2 1.534 microns 6-11 "); 00407 00408 00409 PntForLine(16440.,"Fe 2",&ipnt); 00410 lindst(FeII.fe21006,16440,"Fe 2",ipnt,'i',true, 00411 " Fe 2 1.644 microns 6-10 " ); 00412 00413 00414 PntForLine(12790.,"Fe 2",&ipnt); 00415 lindst(FeII.fe21204,12790,"Fe 2",ipnt,'i',true, 00416 " Fe 2 1.279 microns 12-4" ); 00417 00418 00419 PntForLine(12950.,"Fe 2",&ipnt); 00420 lindst(FeII.fe21103,12950,"Fe 2",ipnt,'i',true, 00421 " Fe 2 1.295 microns 11-3 " ); 00422 00423 00424 PntForLine(13280.,"Fe 2",&ipnt); 00425 lindst(FeII.fe21104,13280,"Fe 2",ipnt,'i',true, 00426 " Fe 2 1.328 microns 11-4" ); 00427 00428 00429 PntForLine(12570.,"Fe 2",&ipnt); 00430 lindst(FeII.fe21001,12570,"Fe 2",ipnt,'i',true, 00431 " Fe 2 1.257 microns 10-1 " ); 00432 00433 00434 PntForLine(13210.,"Fe 2",&ipnt); 00435 lindst(FeII.fe21002,13210,"Fe 2",ipnt,'i',true, 00436 " Fe 2 1.321 microns 10-2 " ); 00437 00438 00439 PntForLine(259880.,"Fe 2",&ipnt); 00440 lindst(FeII.fe20201,259880,"Fe 2",ipnt,'i',true, 00441 " Fe 2 25.988 microns 2-1" ); 00442 00443 PntForLine(353480.,"Fe 2",&ipnt); 00444 lindst(FeII.fe20302,353480,"Fe 2",ipnt,'i',true, 00445 " Fe 2 35.348 microns 3-2" ); 00446 00447 00448 PntForLine(179360.,"Fe 2",&ipnt); 00449 lindst(FeII.fe20706,179360,"Fe 2",ipnt,'i',true, 00450 " Fe 2 17.936 microns 7-6, label is 178 to be unique " ); 00451 00452 00453 PntForLine(245180.,"Fe 2",&ipnt); 00454 lindst(FeII.fe20807,245180,"Fe 2",ipnt,'i',true, 00455 " Fe 2 24.518 microns 8-7" ); 00456 00457 PntForLine(357760.,"Fe 2",&ipnt); 00458 lindst(FeII.fe20908,357760,"Fe 2",ipnt,'i',true, 00459 " Fe 2 35.776 microns 9-8" ); 00460 00461 PntForLine(18100.,"Fe 2",&ipnt); 00462 lindst(FeII.fe21007,18100,"Fe 2",ipnt,'i',true, 00463 " Fe 2 1.810 microns 10-7" ); 00464 00465 00466 PntForLine(16770.,"Fe 2",&ipnt); 00467 lindst(FeII.fe21107,16770,"Fe 2",ipnt,'i',true, 00468 " Fe 2 1.677 microns 11-7 " ); 00469 00470 00471 PntForLine(18000.,"Fe 2",&ipnt); 00472 lindst(FeII.fe21108,18000,"Fe 2",ipnt,'i',true, 00473 " Fe 2 1.800 microns 11-8" ); 00474 00475 00476 PntForLine(17120.,"Fe 2",&ipnt); 00477 lindst(FeII.fe21208,17120,"Fe 2",ipnt,'i',true, 00478 " Fe 2 1.712 microns 12-8 " ); 00479 00480 00481 PntForLine(17980.,"Fe 2",&ipnt); 00482 lindst(FeII.fe21209,17980,"Fe 2",ipnt,'i',true, 00483 " Fe 2 1.798 microns 12-9" ); 00484 00485 00486 PntForLine(229020.,"Fe 2",&ipnt); 00487 lindst(FeII.fe21110,229020,"Fe 2",ipnt,'i',true, 00488 " Fe 2 22.902 microns 11-10 " ); 00489 00490 00491 PntForLine(346600.,"Fe 2",&ipnt); 00492 lindst(FeII.fe21211,346600,"Fe 2",ipnt,'i',true, 00493 " Fe 2 34.660 microns 12-11" ); 00494 00495 /* >>chng 97 mar 14, added following 4 lines after consultation with Verner 00496 * these are typically within 5-9% of their large model atom. These are 00497 * predicted with both large and small atom */ 00498 00499 00500 PntForLine(8617.,"Fe 2",&ipnt); 00501 lindst(FeII.fe21406,8617,"Fe 2",ipnt,'i',true, 00502 " Fe 2 8617A 14-06," ); 00503 00504 00505 PntForLine(8892.,"Fe 2",&ipnt); 00506 lindst(FeII.fe21507,8892,"Fe 2",ipnt,'i',true, 00507 " Fe 2 8892A 15-07" ); 00508 00509 00510 PntForLine(9227.,"Fe 2",&ipnt); 00511 lindst(FeII.fe21508,9227,"Fe 2",ipnt,'i',true, 00512 " Fe 2 9227A 15-08 " ); 00513 00514 00515 PntForLine(9268.,"Fe 2",&ipnt); 00516 lindst(FeII.fe21609,9268,"Fe 2",ipnt,'i',true, 00517 " Fe 2 9268A 16-09" ); 00518 00519 /* following lines are only non-zero when large atom is enabled. When 00520 * running with small FeII atoms these are set zero in routine */ 00521 00522 00523 PntForLine(5159.,"Fe 2",&ipnt); 00524 lindst(FeII.fe25to6,5159,"Fe 2",ipnt,'i',true, 00525 " Fe 2 5158.81 25-6" ); 00526 00527 00528 PntForLine(5262.,"Fe 2",&ipnt); 00529 lindst(FeII.fe27to7,5262,"Fe 2",ipnt,'i',true, 00530 " Fe 2 5261.61 27-7" ); 00531 00532 PntForLine(5334.,"Fe 2",&ipnt); 00533 lindst(FeII.fe28to8,5334,"Fe 2",ipnt,'i',true, 00534 " Fe 2 5333.65 28-8 " ); 00535 00536 00537 PntForLine(5376.,"Fe 2",&ipnt); 00538 lindst(FeII.fe29to9,5376,"Fe 2",ipnt,'i',true, 00539 " Fe 2 5376.47 29-9 " ); 00540 00541 PntForLine(4815.,"Fe 2",&ipnt); 00542 lindst(FeII.fe32to6,4815,"Fe 2",ipnt,'i',true, 00543 " Fe 2 4814.55 32-6 " ); 00544 00545 PntForLine(4905.,"Fe 2",&ipnt); 00546 lindst(FeII.fe33to7,4905,"Fe 2",ipnt,'i',true, 00547 " Fe 2 4905.35 33-7" ); 00548 00549 00550 PntForLine(4347.,"Fe 2",&ipnt); 00551 lindst(FeII.fe37to7,4347,"Fe 2",ipnt,'i',true, 00552 " Fe 2 4346.85 37-7 " ); 00553 00554 00555 PntForLine(4353.,"Fe 2",&ipnt); 00556 lindst(FeII.fe39to8,4353,"Fe 2",ipnt,'i',true, 00557 " Fe 2 4352.78 39-8 " ); 00558 00559 00560 PntForLine(4372.,"Fe 2",&ipnt); 00561 lindst(FeII.fe40to9,4372,"Fe 2",ipnt,'i',true, 00562 " Fe 2 4372.43 40-9 " ); 00563 00564 00565 PntForLine(4244.,"Fe 2",&ipnt); 00566 lindst(FeII.fe37to6,4244,"Fe 2",ipnt,'i',true, 00567 " Fe 2 4243.98 37-6" ); 00568 00569 00570 PntForLine(4277.,"Fe 2",&ipnt); 00571 lindst(FeII.fe39to7,4277,"Fe 2",ipnt,'i',true, 00572 " Fe 2 4276.83 39-7" ); 00573 00574 00575 PntForLine(4320.,"Fe 2",&ipnt); 00576 lindst(FeII.fe40to8,4320,"Fe 2",ipnt,'i',true, 00577 " Fe 2 4319.62 40-8" ); 00578 00579 00580 PntForLine(4358.,"Fe 2",&ipnt); 00581 lindst(FeII.fe41to9,4358,"Fe 2",ipnt,'i',true, 00582 " Fe 2 4358.37 41-9 " ); 00583 00584 00585 PntForLine(4177.,"Fe 2",&ipnt); 00586 lindst(FeII.fe39to6,4177,"Fe 2",ipnt,'i',true, 00587 " Fe 2 4177.21 39-6 " ); 00588 00589 00590 PntForLine(4245.,"Fe 2",&ipnt); 00591 lindst(FeII.fe40to7,4245,"Fe 2",ipnt,'i',true, 00592 " Fe 2 4244.81 40-7 " ); 00593 00594 00595 PntForLine(4306.,"Fe 2",&ipnt); 00596 lindst(FeII.fe41to8,4306,"Fe 2",ipnt,'i',true, 00597 " Fe 2 4305.90 41-8" ); 00598 00599 /* >>chng 99 jun 18 new lines added as per Katya email, so all orion lines are in */ 00600 00601 PntForLine(4114.,"Fe 2",&ipnt); 00602 lindst(FeII.fe42to6,4114,"Fe 2",ipnt,'i',true, 00603 " Fe 2 4114.48 42-6 " ); 00604 00605 00606 PntForLine(4179.,"Fe 2",&ipnt); 00607 lindst(FeII.fe43to7,4179,"Fe 2",ipnt,'i',true, 00608 " Fe 2 4178.95 7- 43" ); 00609 00610 00611 PntForLine(4211.,"Fe 2",&ipnt); 00612 lindst(FeII.fe42to7,4211,"Fe 2",ipnt,'i',true, 00613 " Fe 2 4211.10 7- 42 " ); 00614 00615 00616 PntForLine(4314.,"Fe 2",&ipnt); 00617 lindst(FeII.fe80to28,4314,"Fe 2",ipnt,'i',true, 00618 " Fe 2 4314.29 28- 80" ); 00619 00620 PntForLine(4359.,"Fe 2",&ipnt); 00621 lindst(FeII.fe36to2,4359,"Fe 2",ipnt,'i',true, 00622 " Fe 2 4359.34 2- 36" ); 00623 00624 00625 PntForLine(4414.,"Fe 2",&ipnt); 00626 lindst(FeII.fe36to3,4414,"Fe 2",ipnt,'i',true, 00627 " Fe 2 4413.78 3- 36 " ); 00628 00629 00630 PntForLine(4416.,"Fe 2",&ipnt); 00631 lindst(FeII.fe32to1,4416,"Fe 2",ipnt,'i',true , 00632 " Fe 2 4416.27 1- 32"); 00633 00634 00635 PntForLine(4458.,"Fe 2",&ipnt); 00636 lindst(FeII.fe33to2,4458,"Fe 2",ipnt,'i',true, 00637 " Fe 2 4457.95 2- 33 " ); 00638 00639 00640 PntForLine(4475.,"Fe 2",&ipnt); 00641 lindst(FeII.fe36to5,4475,"Fe 2",ipnt,'i',true, 00642 " Fe 2 4474.91 5- 36 " ); 00643 00644 00645 PntForLine(4493.,"Fe 2",&ipnt); 00646 lindst(FeII.fe32to2,4493,"Fe 2",ipnt,'i',true, 00647 " Fe 2 4492.64 2- 32" ); 00648 00649 00650 PntForLine(4515.,"Fe 2",&ipnt); 00651 lindst(FeII.fe33to3,4515,"Fe 2",ipnt,'i',true, 00652 " Fe 2 4514.90 3- 33" ); 00653 00654 00655 PntForLine(4728.,"Fe 2",&ipnt); 00656 lindst(FeII.fe30to3,4728,"Fe 2",ipnt,'i',true, 00657 " Fe 2 4728.07 3- 30" ); 00658 00659 00660 PntForLine(4775.,"Fe 2",&ipnt); 00661 lindst(FeII.fe33to6,4775,"Fe 2",ipnt,'i',true, 00662 " Fe 2 4774.74 6- 33" ); 00663 00664 PntForLine(4890.,"Fe 2",&ipnt); 00665 lindst(FeII.fe24to2,4890,"Fe 2",ipnt,'i',true, 00666 " Fe 2 4889.63 2- 24" ); 00667 00668 00669 PntForLine(4947.,"Fe 2",&ipnt); 00670 lindst(FeII.fe32to7,4947,"Fe 2",ipnt,'i',true, 00671 " Fe 2 4947.37 7- 32" ); 00672 00673 00674 PntForLine(4951.,"Fe 2",&ipnt); 00675 lindst(FeII.fe35to8,4951,"Fe 2",ipnt,'i',true, 00676 " Fe 2 4950.74 8- 35" ); 00677 00678 00679 PntForLine(4973.,"Fe 2",&ipnt); 00680 lindst(FeII.fe34to8,4973,"Fe 2",ipnt,'i',true, 00681 " Fe 2 4973.39 8- 34" ); 00682 00683 00684 PntForLine(5112.,"Fe 2",&ipnt); 00685 lindst(FeII.fe27to6,5112,"Fe 2",ipnt,'i',true, 00686 " Fe 2 5111.63 6- 27" ); 00687 00688 00689 PntForLine(5220.,"Fe 2",&ipnt); 00690 lindst(FeII.fe28to7,5220,"Fe 2",ipnt,'i',true, 00691 " Fe 2 5220.06 7- 28" ); 00692 00693 00694 PntForLine(5269.,"Fe 2",&ipnt); 00695 lindst(FeII.fe30to8,5269,"Fe 2",ipnt,'i',true, 00696 " Fe 2 5268.880 8- 30 " ); 00697 00698 00699 PntForLine(5273.,"Fe 2",&ipnt); 00700 lindst(FeII.fe24to6,5273,"Fe 2",ipnt,'i',true, 00701 " Fe 2 5273.35 6- 24" ); 00702 00703 PntForLine(5297.,"Fe 2",&ipnt); 00704 lindst(FeII.fe29to8,5297,"Fe 2",ipnt,'i',true, 00705 " Fe 2 5296.840 8- 29 " ); 00706 00707 00708 PntForLine(5433.,"Fe 2",&ipnt); 00709 lindst(FeII.fe24to7,5433,"Fe 2",ipnt,'i',true, 00710 " Fe 2 5433.13 7- 24" ); 00711 00712 00713 PntForLine(5527.,"Fe 2",&ipnt); 00714 lindst(FeII.fe22to7,5527,"Fe 2",ipnt,'i',true, 00715 " Fe 2 5527.340 7- 22" ); 00716 00717 00718 PntForLine(5747.,"Fe 2",&ipnt); 00719 lindst(FeII.fe38to11,5747,"Fe 2",ipnt,'i',true, 00720 " Fe 2 5746.96 11- 38" ); 00721 00722 00723 PntForLine(6440.,"Fe 2",&ipnt); 00724 lindst(FeII.fe19to8,6440,"Fe 2",ipnt,'i',true, 00725 " Fe 2 6440.400 8- 19" ); 00726 00727 PntForLine(7155.,"Fe 2",&ipnt); 00728 lindst(FeII.fe17to6,7155,"Fe 2",ipnt,'i',true, 00729 " Fe 2 7155.140 6- 17 " ); 00730 00731 00732 PntForLine(7172.,"Fe 2",&ipnt); 00733 lindst(FeII.fe18to7,7172,"Fe 2",ipnt,'i',true, 00734 " Fe 2 7171.980 7- 18" ); 00735 00736 00737 PntForLine(7388.,"Fe 2",&ipnt); 00738 lindst(FeII.fe18to8,7388,"Fe 2",ipnt,'i',true, 00739 " Fe 2 7388.160 8- 18" ); 00740 00741 00742 PntForLine(7453.,"Fe 2",&ipnt); 00743 lindst(FeII.fe17to7,7453,"Fe 2",ipnt,'i',true, 00744 " Fe 2 7452.540 7- 17" ); 00745 00746 /* option to save all intensities predicted by large FeII atom, 00747 * code is in FeIILevelPops */ 00748 FeIIAddLines(); 00749 00750 /* some bands from the large model Fe II atom 00751 * only done when full large atom is turned on */ 00752 if( FeII.lgFeIILargeOn ) 00753 { 00754 for( i=0; i < nFeIIBands; i++ ) 00755 { 00756 /* [i][0] is center wavelength, [i][1] and [i][2] are upper and 00757 * lower bounds in Angstroms. These are set in FeIIZero 00758 * units are erg s-1 cm-3 */ 00759 eff = FeIISumBand(FeII_Bands[i][1],FeII_Bands[i][2]); 00760 00761 linadd(eff,FeII_Bands[i][0],"Fe2b",'i' , 00762 " total Fe II emission in Fe II bands, as defined in bands_Fe2.dat "); 00763 } 00765 if( LineSave.ipass > 0 ) 00766 { 00767 for( i=0; i < nFeIIConBins; i++ ) 00768 { 00769 /* [i][0] is total intensity in cell, [i][1] and [i][2] are lower and 00770 * upper bounds in Angstroms. these are set in FeIIZero * 00771 * find total emission from large FeII atom, integrated over band */ 00772 FeII_Cont[i][0] += 00773 (realnum)(FeIISumBand(FeII_Cont[i][1],FeII_Cont[i][2])*radius.dVeff); 00774 /*fprintf(ioQQQ,"DEBUG feii\t%li\t%.2e\n", i, FeII_Cont[i][0]);*/ 00775 } 00776 } 00777 else 00778 { 00779 for( i=0; i < nFeIIConBins; i++ ) 00780 { 00781 /* initialize arrays */ 00782 FeII_Cont[i][0] = 0.; 00783 } 00784 } 00785 } 00786 else 00787 { 00788 /* must still add the bands in case of grid where FeII is turned off and on */ 00789 for( i=0; i < nFeIIBands; i++ ) 00790 { 00791 00792 linadd(0., FeII_Bands[i][0] ,"Fe2b",'i' , 00793 " total Fe II emission in Fe II bands, as defined in bands_Fe2.dat "); 00794 } 00795 } 00796 00797 PutLine(&TauLines[ipT191], 00798 " anomalous Fe 2 transition at 1787, RMT 191"); 00799 00800 linadd(fe.Fe3CoolTot,0,"Fe3c",'c' , 00801 " chng 05 dec 16, FeIII code created by Kevin Blagrave Fe3c 0 - total cooling due to 14-level Fe 3 atom "); 00802 00803 /* Fe 3 14-level atom 00804 * following from print statements within loop */ 00805 /* Fe 3 22.92m from Blagrave 14-level atom */ 00806 /* Fe 3 13.53m from Blagrave 14-level atom */ 00807 /* Fe 3 33.03m from Blagrave 14-level atom */ 00808 /* Fe 3 10.72m from Blagrave 14-level atom */ 00809 /* Fe 3 20.15m from Blagrave 14-level atom */ 00810 /* Fe 3 51.67m from Blagrave 14-level atom */ 00811 /* Fe 3 9.732m from Blagrave 14-level atom */ 00812 /* Fe 3 16.91m from Blagrave 14-level atom */ 00813 /* Fe 3 34.66m from Blagrave 14-level atom */ 00814 /* Fe 3 105.3m from Blagrave 14-level atom */ 00815 /* Fe 3 5152A from Blagrave 14-level atom */ 00816 /* Fe 3 5271A from Blagrave 14-level atom */ 00817 /* Fe 3 5356A from Blagrave 14-level atom */ 00818 /* Fe 3 5412A from Blagrave 14-level atom */ 00819 /* Fe 3 5440A from Blagrave 14-level atom */ 00820 /* Fe 3 4986A from Blagrave 14-level atom */ 00821 /* Fe 3 5097A from Blagrave 14-level atom */ 00822 /* Fe 3 5177A from Blagrave 14-level atom */ 00823 /* Fe 3 5230A from Blagrave 14-level atom */ 00824 /* Fe 3 5256A from Blagrave 14-level atom */ 00825 /* Fe 3 15.47m from Blagrave 14-level atom */ 00826 /* Fe 3 4925A from Blagrave 14-level atom */ 00827 /* Fe 3 5033A from Blagrave 14-level atom */ 00828 /* Fe 3 5111A from Blagrave 14-level atom */ 00829 /* Fe 3 5162A from Blagrave 14-level atom */ 00830 /* Fe 3 5188A from Blagrave 14-level atom */ 00831 /* Fe 3 11.16m from Blagrave 14-level atom */ 00832 /* Fe 3 40.04m from Blagrave 14-level atom */ 00833 /* Fe 3 4881A from Blagrave 14-level atom */ 00834 /* Fe 3 4988A from Blagrave 14-level atom */ 00835 /* Fe 3 5064A from Blagrave 14-level atom */ 00836 /* Fe 3 5114A from Blagrave 14-level atom */ 00837 /* Fe 3 5139A from Blagrave 14-level atom */ 00838 /* Fe 3 9.282m from Blagrave 14-level atom */ 00839 /* Fe 3 23.21m from Blagrave 14-level atom */ 00840 /* Fe 3 55.20m from Blagrave 14-level atom */ 00841 /* Fe 3 4833A from Blagrave 14-level atom */ 00842 /* Fe 3 4937A from Blagrave 14-level atom */ 00843 /* Fe 3 5012A from Blagrave 14-level atom */ 00844 /* Fe 3 5061A from Blagrave 14-level atom */ 00845 /* Fe 3 5085A from Blagrave 14-level atom */ 00846 /* Fe 3 7.789m from Blagrave 14-level atom */ 00847 /* Fe 3 15.69m from Blagrave 14-level atom */ 00848 /* Fe 3 25.79m from Blagrave 14-level atom */ 00849 /* Fe 3 48.41m from Blagrave 14-level atom */ 00850 /* Fe 3 4714A from Blagrave 14-level atom */ 00851 /* Fe 3 4813A from Blagrave 14-level atom */ 00852 /* Fe 3 4884A from Blagrave 14-level atom */ 00853 /* Fe 3 4931A from Blagrave 14-level atom */ 00854 /* Fe 3 4954A from Blagrave 14-level atom */ 00855 /* Fe 3 5.543m from Blagrave 14-level atom */ 00856 /* Fe 3 8.638m from Blagrave 14-level atom */ 00857 /* Fe 3 11.01m from Blagrave 14-level atom */ 00858 /* Fe 3 13.76m from Blagrave 14-level atom */ 00859 /* Fe 3 19.22m from Blagrave 14-level atom */ 00860 /* Fe 3 4659A from Blagrave 14-level atom */ 00861 /* Fe 3 4755A from Blagrave 14-level atom */ 00862 /* Fe 3 4825A from Blagrave 14-level atom */ 00863 /* Fe 3 4870A from Blagrave 14-level atom */ 00864 /* Fe 3 4893A from Blagrave 14-level atom */ 00865 /* Fe 3 4.859m from Blagrave 14-level atom */ 00866 /* Fe 3 7.085m from Blagrave 14-level atom */ 00867 /* Fe 3 8.608m from Blagrave 14-level atom */ 00868 /* Fe 3 10.20m from Blagrave 14-level atom */ 00869 /* Fe 3 12.92m from Blagrave 14-level atom */ 00870 /* Fe 3 39.41m from Blagrave 14-level atom */ 00871 /* Fe 3 4608A from Blagrave 14-level atom */ 00872 /* Fe 3 4702A from Blagrave 14-level atom */ 00873 /* Fe 3 4770A from Blagrave 14-level atom */ 00874 /* Fe 3 4814A from Blagrave 14-level atom */ 00875 /* Fe 3 4836A from Blagrave 14-level atom */ 00876 /* Fe 3 4.356m from Blagrave 14-level atom */ 00877 /* Fe 3 6.063m from Blagrave 14-level atom */ 00878 /* Fe 3 7.146m from Blagrave 14-level atom */ 00879 /* Fe 3 8.208m from Blagrave 14-level atom */ 00880 /* Fe 3 9.884m from Blagrave 14-level atom */ 00881 /* Fe 3 20.34m from Blagrave 14-level atom */ 00882 /* Fe 3 42.06m from Blagrave 14-level atom */ 00883 /* Fe 3 4574A from Blagrave 14-level atom */ 00884 /* Fe 3 4668A from Blagrave 14-level atom */ 00885 /* Fe 3 4734A from Blagrave 14-level atom */ 00886 /* Fe 3 4778A from Blagrave 14-level atom */ 00887 /* Fe 3 4800A from Blagrave 14-level atom */ 00888 /* Fe 3 4.077m from Blagrave 14-level atom */ 00889 /* Fe 3 5.535m from Blagrave 14-level atom */ 00890 /* Fe 3 6.423m from Blagrave 14-level atom */ 00891 /* Fe 3 7.269m from Blagrave 14-level atom */ 00892 /* Fe 3 8.554m from Blagrave 14-level atom */ 00893 /* Fe 3 15.41m from Blagrave 14-level atom */ 00894 /* Fe 3 25.31m from Blagrave 14-level atom */ 00895 /* Fe 3 63.56m from Blagrave 14-level atom */ 00896 for( ihi=1; ihi<NLFE3; ++ihi ) 00897 { 00898 for( ilo=0; ilo<ihi; ++ilo ) 00899 { 00900 /* emission in these lines */ 00901 PntForLine(fe.Fe3_wl[ihi][ilo],"Fe 3",&ipnt); 00902 # if 0 00903 fprintf( ioQQQ,"\t/* FeIII "); 00904 prt_wl( ioQQQ , (realnum)(fe.Fe3_wl[ihi][ilo]+0.5) ); 00905 fprintf( ioQQQ," from Blagrave 14-level atom */\n "); 00906 # endif 00907 lindst( fe.Fe3_emiss[ihi][ilo] , (realnum)(fe.Fe3_wl[ihi][ilo]+0.5) , "Fe 3",ipnt,'c',true, 00908 " " ); 00909 } 00910 } 00911 00912 /*>>chng 05 dec 18, following are now in the above */ 00913 /* sum of 3p and 3g states together */ 00914 /* linadd(CoolHeavy.c5270,0,"Fe 3",'c' ); */ 00915 00916 /* Fe 3 5270, predictions from Garstang et al 78 00917 PntForLine(5270.,"Fe 3",&ipnt); 00918 lindst(CoolHeavy.c5270*0.2090,5270,"Fe 3",ipnt,'c',true );*/ 00919 00920 /* Fe 3 5270, predictions from Garstang et al 78 00921 PntForLine(4658.,"Fe 3",&ipnt); 00922 lindst(CoolHeavy.c5270*0.3667,4658,"Fe 3",ipnt,'c',true ); */ 00923 00924 PutLine(&TauLines[ipT1122]," Fe 3 1122 entire multiplet"); 00925 00926 linadd(fe.Fe4CoolTot,0,"Fe4c",'c', 00927 " Fe4c 0 - total cooling due to 12-level Fe 4 atom " ); 00928 00929 00930 PntForLine(3096.,"Fe 4",&ipnt); 00931 lindst(fe.fe40401,3096,"Fe 4",ipnt,'i',true, 00932 " Fe 4 3096.A, 4-1 and 5-1 transitions together" ); 00933 00934 00935 PntForLine(2836.,"Fe 4",&ipnt); 00936 lindst(fe.fe42836,2836,"Fe 4",ipnt,'i',true, 00937 " Fe 4 2835.7A, 6-1 transition, 4P5/2 - 6S5/2 " ); 00938 00939 00940 PntForLine(2829.,"Fe 4",&ipnt); 00941 lindst(fe.fe42829,2829,"Fe 4",ipnt,'i',true, 00942 " Fe 4 2829.4A, 7-1 transition, 4P3/2 - 6S5/2" ); 00943 00944 00945 PntForLine(2567.,"Fe 4",&ipnt); 00946 lindst(fe.fe42567,2567,"Fe 4",ipnt,'i',true, 00947 " Fe 4 2567.6+ 2567.4. 11-1 and 12-1 transitions" ); 00948 00949 00950 PntForLine(2.774e4,"Fe 4",&ipnt); 00951 lindst(fe.fe41207,277,"Fe 4",ipnt,'i',true, 00952 " Fe 4 2.774 microns 12-7 transition " ); 00953 00954 00955 PntForLine(2.714e4,"Fe 4",&ipnt); 00956 lindst(fe.fe41206,271,"Fe 4",ipnt,'i',true, 00957 " Fe 4 2.714 microns 12-6 transition " ); 00958 00959 00960 PntForLine(2.716e4,"Fe 4",&ipnt); 00961 lindst(fe.fe41106,272,"Fe 4",ipnt,'i',true, 00962 " Fe 4 2.716 microns 11-6 transition" ); 00963 00964 00965 PntForLine(2.806e4,"Fe 4",&ipnt); 00966 lindst(fe.fe41007,281,"Fe 4",ipnt,'i',true, 00967 " Fe 4 2.806 microns 10-7 transition " ); 00968 00969 00970 PntForLine(2.865e4,"Fe 4",&ipnt); 00971 lindst(fe.fe41008,287,"Fe 4",ipnt,'i',true , 00972 " Fe 4 2.865 microns 10-8 transition"); 00973 00974 00975 PntForLine(2.836e4,"Fe 4",&ipnt); 00976 lindst(fe.fe40906,284,"Fe 4",ipnt,'i',true, 00977 " Fe 4 2.836 microns 9-6 transition" ); 00978 00979 00980 PntForLine(3892.,"Fe 5",&ipnt); 00981 lindst(CoolHeavy.c3892,3892,"Fe 5",ipnt,'c',true, 00982 " Fe 5 3892+3839" ); 00983 00984 linadd(CoolHeavy.c5177,0,"Fe 6",'c' , 00985 " all of 2G lines together first "); 00986 00987 00988 PntForLine(5177.,"Fe 6",&ipnt); 00989 lindst(CoolHeavy.c5177*0.354,5177,"Fe 6",ipnt,'c',true, 00990 " Fe 6 5177, approximate correct " ); 00991 00992 linadd(fe.Fe7CoolTot,0,"Fe7c",'c' , 00993 " Fe7c 0 - total cooling due to n-level Fe 7 atom "); 00994 00995 /* >>chng 04 nov 04, move to multi-level system */ 00996 for( ilo=0; ilo<NLFE7-1; ++ilo ) 00997 { 00998 /* must not do 1-0 or 2-1, which are transferred lines */ 00999 for( ihi=MAX2(3,ilo+1); ihi<NLFE7; ++ihi ) 01000 { 01001 01002 PntForLine(fe.Fe7_wl[ihi][ilo],"Fe 7",&ipnt); 01003 lindst( fe.Fe7_emiss[ihi][ilo] , (realnum)(fe.Fe7_wl[ihi][ilo]+0.5) , "Fe 7",ipnt,'c',true, 01004 " emission in these lines" ); 01005 } 01006 } 01007 # if 0 01008 PntForLine(5721.,"Fe 7",&ipnt); 01009 lindst( fe.Fe7_5721 , 5721 , "Fe 7",ipnt,'c',true, 01010 " " ); 01011 01012 PntForLine(6601.,"Fe 7",&ipnt); 01013 lindst( fe.Fe7_6601 , 6601 , "Fe 7",ipnt,'c',true, 01014 " " ); 01015 01016 PntForLine(3760.,"Fe 7",&ipnt); 01017 lindst( fe.Fe7_3760 , 3760 , "Fe 7",ipnt,'c',true, 01018 " " ); 01019 01020 PntForLine(3588.,"Fe 7",&ipnt); 01021 lindst( fe.Fe7_3588 , 3588 , "Fe 7",ipnt,'c',true, 01022 " " ); 01023 # endif 01024 01025 PutLine(&TauLines[ipFe0795], 01026 " [Fe 7] 9.51 micron "); 01027 01028 PutLine(&TauLines[ipFe0778], 01029 " [Fe 7] 7.81 micron "); 01030 01031 /* [Fe 7] 6087 01032 PntForLine(6087.,"Fe 7",&ipnt); 01033 lindst(CoolHeavy.c6087,6087,"Fe 7",ipnt,'c',true );*/ 01034 01035 /* [Fe 7] 5722 01036 PntForLine(5722.,"Fe 7",&ipnt); 01037 lindst(CoolHeavy.Fe5722,5722,"Fe 7",ipnt,'c',true );*/ 01038 01039 PutLine(&TauLines[ipT245], 01040 " Be-seq lines "); 01041 01042 PntForLine(242.,"Fe 7",&ipnt); 01043 lindst(CoolHeavy.c242,242,"Fe 7",ipnt,'c',true, 01044 " Fe 9 242 j=1 slower decay"); 01045 01046 PutLine(&TauLines[ipT352], 01047 " the E1 transition that can pump [Fe X] "); 01048 01049 /* optically thin Fe X pumping */ 01050 eff = 69.4/(69.4 + 0.27*dense.cdsqte); 01051 01052 /* coll excitation of 352 which decays to excited state of 6374 01053 * assumes 17/56 (ratio of A's) go to excited state */ 01054 c10 = TauLines[ipT352].Emis->phots*eff*.01676*TauLines[ipT352].EnergyErg*352/6374.; 01055 01056 /* Fe 10 and Fe 14 from Mason 75 01057 * total (coll, pumped) 01058 * call linadd( C6374+C10 , 6374 , 'Fe10','i') */ 01059 PutExtra( c10 ); 01060 01061 01062 PutLine(&TauLines[ipFe106375], 01063 " [Fe 10] 6375, collisions with pumping too "); 01064 01065 /* collisional contribution 01066 * call linadd( C6374 , 6374 , 'Coll','c') 01067 * collisions of E1 line, plus pumped by continuum fluorescence 01068 * call linadd( C10 , 6374 , ' 352','c') 01069 * Fe XI 7892, 6.08 micron 01070 PutLine(&TauLines[ipTFe07]);*/ 01071 01083 /* Fe 11 2649 collisional excitation 01084 PntForLine(2649.,"Fe11",&ipnt); 01085 lindst(CoolHeavy.c2649,2649,"Fe11",ipnt,'c',true ); */ 01086 01087 /* Fe 11 1467 collisional excitation 01088 PntForLine(1467.,"Fe11",&ipnt); 01089 lindst(CoolHeavy.c1467,1467,"Fe11",ipnt,'c',true );*/ 01090 01091 linadd(fe.Fe11CoolTot,0,"Fe11",'c' , 01092 " >>chng 05 dec 18, add Fe 11 Fe11 0 - total cooling due to 5-level Fe 11 atom "); 01093 01094 /* Fe 11 5-level atom */ 01095 for( ihi=1; ihi<NLFE11; ++ihi ) 01096 { 01097 for( ilo=0; ilo<ihi; ++ilo ) 01098 { 01099 PntForLine(fe.Fe11_wl[ihi][ilo],"Fe11",&ipnt); 01100 lindst( fe.Fe11_emiss[ihi][ilo] , (realnum)(fe.Fe11_wl[ihi][ilo]+0.5) , "Fe11",ipnt,'c',true, 01101 " emission in these lines" ); 01102 } 01103 } 01104 01105 PntForLine(1242.,"Fe12",&ipnt); 01106 lindst(CoolHeavy.c1242,1242,"Fe12",ipnt,'c',true , 01107 " Fe 12, 1242, 1349 together, collisional excitation"); 01108 01109 PntForLine(2170.,"Fe12",&ipnt); 01110 lindst(CoolHeavy.c2170,2170,"Fe12",ipnt,'c',true , 01111 " Fe 12, 2170, 2406 together, collisional excitation"); 01112 01113 01114 PntForLine(2568.,"Fe12",&ipnt); 01115 lindst(CoolHeavy.c2568,2568,"Fe12",ipnt,'c',true, 01116 " Fe12 2904, 2567, 3567, 3073 together, collisional excitation" ); 01117 01118 /* >>chng 05 dec 18, add Fe 13 */ 01119 /* Fe13 0 - total cooling due to 5-level Fe 13 atom */ 01120 linadd(fe.Fe13CoolTot,0,"Fe13",'c' , 01121 "total cooling due to Fe 13 model atom "); 01122 01123 /* Fe 13 5-level atom */ 01124 for( ihi=1; ihi<NLFE13; ++ihi ) 01125 { 01126 for( ilo=0; ilo<ihi; ++ilo ) 01127 { 01128 01129 PntForLine(fe.Fe13_wl[ihi][ilo],"Fe13",&ipnt); 01130 lindst( fe.Fe13_emiss[ihi][ilo] , (realnum)(fe.Fe13_wl[ihi][ilo]+0.5) , "Fe13",ipnt,'c',true , 01131 " Fe 13 emission"); 01132 } 01133 } 01134 01135 /* Fe 14 optically thin in line 344 */ 01136 eff = 60.3/(60.3 + 0.23*dense.cdsqte/4.); 01137 r14 = dense.xIonDense[ipIRON][14-1]*fe.pfe14*rfield.flux[fe.ipfe10-1]* 01138 3.75e-12*eff/(1. + 24./.63* 01139 TauLines[ipT347].Emis->Pesc)*esc_PRD_1side(TauLines[ipT347].Emis->TauIn,1e-4); 01140 PntForLine(5303.,"Fe14",&ipnt); 01141 01142 lindst(CoolHeavy.c5303+r14,5303,"Fe14",ipnt,'i',true ," total emission in Fe 14 5304"); 01143 01144 linadd(CoolHeavy.c5303,5303,"Coll",'c' , 01145 " Fe 14 5304 contribution from collisional excitation "); 01146 01147 linadd(r14,5303,"Pump",'r' , 01148 " Fe 14 5304 continuum fluorescense "); 01149 01152 c14 = 0.; 01153 linadd(c14,5303," 347",'c' , 01154 " collisional excitation of E1 line "); 01155 01156 PutLine(&TauLines[ipFe18975], 01157 " Fe 18 974.86A "); 01158 01159 01160 PntForLine(7047.,"Fe19",&ipnt); 01161 lindst(CoolHeavy.c7082,7047,"Fe19",ipnt,'c',true, 01162 " O-like Fe19, 3P ground term, 7046.72A vacuum wl, 1328.90A >>chng 01 aug 10, updated wavelengths Fe 19 7047 '85 " ); 01163 01164 01165 PntForLine(1328.,"Fe19",&ipnt); 01166 lindst(CoolHeavy.c1328,1329,"Fe19",ipnt,'c',true, 01167 " Fe 19 1329" ); 01168 01169 01170 PntForLine(592.,"Fe19",&ipnt); 01171 lindst(CoolHeavy.c592,592,"Fe19",ipnt,'c',true, 01172 " Fe 19 from loulergue et al '85" ); 01173 01174 01175 PntForLine(1118.,"Fe19",&ipnt); 01176 lindst(CoolHeavy.c1118,1118,"Fe19",ipnt,'c',true , 01177 " Fe 19 from loulergue et al '85"); 01178 01179 PutLine(&TauLines[ipTFe13], 01180 " next two 3p ground state lines, collisional excitation "); 01181 01182 PutLine(&TauLines[ipTFe23], 01183 " collisional excitation "); 01184 01185 PutLine(&TauLines[ipTFe20_578], 01186 " Fe20 721.40A, 578"); 01187 PutLine(&TauLines[ipTFe20_721], 01188 " "); 01189 01190 linadd( 01191 TauLines[ipFe22_247].Emis->xIntensity+ 01192 TauLines[ipFe22_217].Emis->xIntensity+ 01193 TauLines[ipFe22_348].Emis->xIntensity+ 01194 TauLines[ipFe22_292].Emis->xIntensity+ 01195 TauLines[ipFe22_253].Emis->xIntensity, 01196 260,"TOTL",'i', 01197 " Fe 22 845.6A total intensity of Fe22, all lines in the multiplet " ); 01198 PutLine(&TauLines[ipFe22_247], 01199 "Fe 22 247"); 01200 PutLine(&TauLines[ipFe22_217], 01201 "Fe 22 217"); 01202 PutLine(&TauLines[ipFe22_348], 01203 "Fe 22 348"); 01204 PutLine(&TauLines[ipFe22_292], 01205 "Fe 22 292"); 01206 PutLine(&TauLines[ipFe22_253], 01207 "Fe 22 253"); 01208 01209 /* Fe 23 1909-like 262.6 */ 01210 PntForLine(263.,"Fe23",&ipnt); 01211 lindst(CoolHeavy.c263,263,"Fe23",ipnt,'c',true, 01212 "Fe 23 1909-like 262.6" ); 01213 01214 01215 PutLine(&TauLines[ipT192], 01216 " Fe 24 only 192 of 255, 192 Li seq doublet, collisional excitation "); 01217 01218 PutLine(&TauLines[ipT255], 01219 " Fe 24 255 of 255, 192 Li seq doublet, collisional excitation "); 01220 01221 PutLine(&TauLines[ipT11], 01222 " Fe 24 Li seq 2s3p collisional excitation "); 01223 01224 /* recombination Ka */ 01225 if( dense.lgElmtOn[ipIRON] ) 01226 { 01227 /* these lines added to outlin in metdif - following must be false 01228 * fela = xLyaHeavy(nelem,nelem)*dense.xIonDense(nelem,nelem+1) */ 01229 fela = Transitions[ipH_LIKE][ipIRON][ipH2p][ipH1s].Emis->xIntensity; 01230 } 01231 else 01232 { 01233 fela = 0.; 01234 } 01235 01236 /* >>chng 02 jan 14, add grain fe to this sum */ 01237 /* total intensity of K-alpha line */ 01238 /*linadd((fe.fekcld+fe.fegrain)*1.03e-8+(fe.fekhot+fela)*1.11e-8,2,"FeKa",'i' );*/ 01239 if( dense.lgElmtOn[ipIRON] ) 01240 { 01241 lindst((fe.fekcld+fe.fegrain)*1.03e-8+(fe.fekhot+fela)*1.11e-8,2,"FeKa", 01242 Transitions[ipH_LIKE][ipIRON][ipH2p][ipH1s].ipCont,'i',false, 01243 "total intensity of K-alpha line" ); 01244 } 01245 01246 linadd(fela*1.11e-8,2,"FeLr",'i' , 01247 " recombination from fully stripped ion "); 01248 01249 /* >>chng 03 aug 14, label changed from TotH to AugH to be like rest total hot iron Ka; */ 01250 linadd((fe.fekhot+fela)*1.11e-8,2,"AugH",'i' , 01251 " Auger hot iron, assumes case b for H and He-like "); 01252 01253 linadd(fe.fekcld*1.03e-8,2,"AugC",'i', 01254 " Auger production of cold iron, less than or 17 times ionized " ); 01255 01256 linadd(fe.fegrain*1.03e-8,2,"AugG",'i' , 01257 " grain production of cold iron "); 01258 01259 PutLine(&TauLines[ipCo11527], 01260 " [Co XI] 5168. A "); 01261 01262 PutLine(&TauLines[ipNi1_7m], 01263 " nickel [Ni I] 7m "); 01264 01265 /* nickel*/ 01266 01267 01268 PutLine(&TauLines[ipNi1_11m], 01269 " [Ni I] 11m "); 01270 01271 PutLine(&TauLines[ipxNi1242], 01272 " [Ni XII] 4230.8 A "); 01273 01274 /* copper */ 01275 01276 /* zinc */ 01277 PutLine(&TauLines[ipZn04363], 01278 "zinc iv 3.625 microns, cs and A just made up "); 01279 01280 if( trace.lgTrace ) 01281 { 01282 fprintf( ioQQQ, " lines_lv1_k_zn returns\n" ); 01283 } 01284 return; 01285 }