00001
00002
00003
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
00325
00326
00327
00328
00329
00330
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
00386
00387
00388
00389
00390
00391
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
00496
00497
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
00520
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
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
00747
00748 FeIIAddLines();
00749
00750
00751
00752 if( FeII.lgFeIILargeOn )
00753 {
00754 for( i=0; i < nFeIIBands; i++ )
00755 {
00756
00757
00758
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
00770
00771
00772 FeII_Cont[i][0] +=
00773 (realnum)(FeIISumBand(FeII_Cont[i][1],FeII_Cont[i][2])*radius.dVeff);
00774
00775 }
00776 }
00777 else
00778 {
00779 for( i=0; i < nFeIIConBins; i++ )
00780 {
00781
00782 FeII_Cont[i][0] = 0.;
00783 }
00784 }
00785 }
00786 else
00787 {
00788
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
00804
00805
00806
00807
00808
00809
00810
00811
00812
00813
00814
00815
00816
00817
00818
00819
00820
00821
00822
00823
00824
00825
00826
00827
00828
00829
00830
00831
00832
00833
00834
00835
00836
00837
00838
00839
00840
00841
00842
00843
00844
00845
00846
00847
00848
00849
00850
00851
00852
00853
00854
00855
00856
00857
00858
00859
00860
00861
00862
00863
00864
00865
00866
00867
00868
00869
00870
00871
00872
00873
00874
00875
00876
00877
00878
00879
00880
00881
00882
00883
00884
00885
00886
00887
00888
00889
00890
00891
00892
00893
00894
00895
00896 for( ihi=1; ihi<NLFE3; ++ihi )
00897 {
00898 for( ilo=0; ilo<ihi; ++ilo )
00899 {
00900
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
00913
00914
00915
00916
00917
00918
00919
00920
00921
00922
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
00996 for( ilo=0; ilo<NLFE7-1; ++ilo )
00997 {
00998
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
01032
01033
01034
01035
01036
01037
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
01050 eff = 69.4/(69.4 + 0.27*dense.cdsqte);
01051
01052
01053
01054 c10 = TauLines[ipT352].Emis->phots*eff*.01676*TauLines[ipT352].EnergyErg*352/6374.;
01055
01056
01057
01058
01059 PutExtra( c10 );
01060
01061
01062 PutLine(&TauLines[ipFe106375],
01063 " [Fe 10] 6375, collisions with pumping too ");
01064
01065
01066
01067
01068
01069
01070
01071
01083
01084
01085
01086
01087
01088
01089
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
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
01119
01120 linadd(fe.Fe13CoolTot,0,"Fe13",'c' ,
01121 "total cooling due to Fe 13 model atom ");
01122
01123
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
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
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
01225 if( dense.lgElmtOn[ipIRON] )
01226 {
01227
01228
01229 fela = Transitions[ipH_LIKE][ipIRON][ipH2p][ipH1s].Emis->xIntensity;
01230 }
01231 else
01232 {
01233 fela = 0.;
01234 }
01235
01236
01237
01238
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
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
01266
01267
01268 PutLine(&TauLines[ipNi1_11m],
01269 " [Ni I] 11m ");
01270
01271 PutLine(&TauLines[ipxNi1242],
01272 " [Ni XII] 4230.8 A ");
01273
01274
01275
01276
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 }
