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",'i',
00068 " coll excited calcium k+h " );
00069
00070 linadd(ca.Cair,8579,"Ca 2",'i' ,
00071 " infrared triplet ");
00072
00073 linadd(ca.c7306,7306,"Ca 2",'i',
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,'c',true,
00084 " individual lines from five level atom");
00085
00086
00087 PntForLine(3969.,"Ca2H",&ipnt);
00088 lindst(ca.Cah,3969,"Ca2H",ipnt,'c',true,
00089 " individual lines from five level atom" );
00090
00091
00092 PntForLine(8498.,"Ca2X",&ipnt);
00093 lindst(ca.Cax,8498,"Ca2X",ipnt,'c',true,
00094 " individual lines from five level atom " );
00095
00096
00097 PntForLine(8542.,"Ca2Y",&ipnt);
00098 lindst(ca.Cay,8542,"Ca2Y",ipnt,'c',true,
00099 " individual lines from five level atom" );
00100
00101
00102 PntForLine(8662.,"Ca2Z",&ipnt);
00103 lindst(ca.Caz,8662,"Ca2Z",ipnt,'c',true,
00104 " individual lines from five level atom" );
00105
00106
00107 PntForLine(7291.,"CaF1",&ipnt);
00108 lindst(ca.Caf1,7291,"CaF1",ipnt,'c',true,
00109 " individual lines from five level atom" );
00110
00111
00112 PntForLine(7324.,"CaF2",&ipnt);
00113 lindst(ca.Caf2,7324,"CaF2",ipnt,'c',true,
00114 " individual lines from five level atom" );
00115
00116 eff = dense.eden*dense.xIonDense[ipCALCIUM][2]*5.4e-21/(phycon.te/
00117 phycon.te10/phycon.te10);
00118 linadd(eff,3933,"Rec ",'i',
00119 " recombination contribution to CaII emission" );
00120
00121 PutLine(&TauLines[ipTCa3],
00122 " Ca IV 3.2 micron ");
00123
00124 PutLine(&TauLines[ipTCa4],
00125 " Ca V 4.16, 11.47 micron");
00126
00127 PutLine(&TauLines[ipTCa12],
00128 " Ca V 4.16, 11.47 micron ");
00129
00130 PntForLine(6087.,"Ca 5",&ipnt);
00131 lindst(ca.Ca6087,6087,"Ca 5",ipnt,'c',true ,
00132 " Ca V optical and UV lines, collisional excitation, 3-level atom");
00133
00134 PntForLine(5311.,"Ca 5",&ipnt);
00135 lindst(ca.c5311,5311,"Ca 5",ipnt,'c',true ,
00136 " Ca V optical and UV lines, collisional excitation, 3-level atom");
00137
00138 PntForLine(2414.,"Ca 5",&ipnt);
00139 lindst(ca.c2414,2414,"Ca 5",ipnt,'c',true ,
00140 " Ca V optical and UV lines, collisional excitation, 3-level atom");
00141
00142
00143 PntForLine(3997.,"Ca 5",&ipnt);
00144 lindst(ca.c3997,3997,"Ca 5",ipnt,'c',true,
00145 " Ca V optical and UV lines, collisional excitation, 3-level atom" );
00146
00147 PutLine(&TauLines[ipCa0741],
00148 " [Ca VII] 4.09 microns" );
00149
00150 PutLine(&TauLines[ipCa0761],
00151 " [Ca VII] 6.15 microns " );
00152
00153
00154 PntForLine(5620.,"Ca 7",&ipnt);
00155 lindst(ca.Ca5620,5620,"Ca 7",ipnt,'c',true,
00156 " Ca VII optical and UV lines, collisional excitation, 3-level atom" );
00157
00158
00159 PntForLine(4941.,"Ca 7",&ipnt);
00160 lindst(ca.Ca4941,4941,"Ca 7",ipnt,'c',true,
00161 " Ca VII optical and UV lines, collisional excitation, 3-level atom" );
00162
00163
00164 PntForLine(2112.,"Ca 7",&ipnt);
00165 lindst(ca.Ca2112,2112,"Ca 7",ipnt,'c',true ,
00166 " Ca VII optical and UV lines, collisional excitation, 3-level atom");
00167
00168
00169 PntForLine(3688.,"Ca 7",&ipnt);
00170 lindst(ca.Ca3688,3688,"Ca 7",ipnt,'c',true,
00171 " Ca VII optical and UV lines, collisional excitation, 3-level atom" );
00172
00173 PutLine(&TauLines[ipCa08232],
00174 " [Ca VIII] 2.32 microns, A Saraph and Strey ");
00175
00176 PutLine(&TauLines[ipCa12333],
00177 " [Ca 12] 3328.78A ");
00178
00179 PutLine(&TauLines[ipTCa302],
00180 " Ca 18 Li seq 2s2p, 302, 345 separate ");
00181
00182 PutLine(&TauLines[ipTCa345],
00183 " Ca 18 Li seq 2s2p, 302, 345 separate ");
00184
00185 PutLine(&TauLines[ipTCa19],
00186 " Ca 18 Li seq 2s3p, ");
00187
00188
00189 PntForLine(22.08e4,"Sc 2",&ipnt);
00190 lindst(CoolHeavy.Sc22p08m,22.08e4,"Sc 2",ipnt,'c',true,
00191 " Sc II 2.08 (1-3) " );
00192
00193
00194 PntForLine(24.1e4,"Sc 2",&ipnt);
00195 lindst(CoolHeavy.Sc24p1m,24.1e4,"Sc 2",ipnt,'c',true,
00196 " Sc II 4.1 micron (1-2)" );
00197
00198
00199 PntForLine(24.2e4,"Sc 2",&ipnt);
00200 lindst(CoolHeavy.Sc24p2m,24.2e4,"Sc 2",ipnt,'c',true,
00201 " Sc II 4.22 (2-3)" );
00202
00203
00204 PntForLine(3933.,"Sc 3",&ipnt);
00205 lindst(CoolHeavy.Sc33936,3933,"Sc 3",ipnt,'c',true,
00206 " Sc III 3936" );
00207
00208 PutLine(&TauLines[ipSc05231],
00209 " [Sc V] 1.46 microns ");
00210
00211
00212 PntForLine(5054.,"Sc 6",&ipnt);
00213 lindst(CoolHeavy.Sc45058,5054,"Sc 6",ipnt,'c',true ,
00214 " Sc VI 5054 (1-2)");
00215
00216
00217 PntForLine(3592.,"Sc 6",&ipnt);
00218 lindst(CoolHeavy.Sc43595,3592,"Sc 6",ipnt,'c',true,
00219 " Sc VI 3595 (2-3)" );
00220
00221
00222 PntForLine(2100.,"Sc 6",&ipnt);
00223 lindst(CoolHeavy.Sc42100,2100,"Sc 6",ipnt,'c',true,
00224 " Sc VI 2100 (1-3)" );
00225
00226 PutLine(&TauLines[ipSc13264],
00227 " [Sc 13] 2637.97A");
00228
00229
00230 PntForLine(1.21e4,"Ti 3",&ipnt);
00231 lindst(CoolHeavy.Ti3l21,1.21e4,"Ti 3",ipnt,'c',true,
00232 " Ti III 1.21 micron, (actually multiplet) 2-1 transition from model atom " );
00233
00234
00235 PntForLine(9594.,"Ti 3",&ipnt);
00236 lindst(CoolHeavy.Ti3l31,9594,"Ti 3",ipnt,'c',true,
00237 " Ti III 9594, 3-1 transition, (actually multiplet) from model atom" );
00238
00239 PntForLine(4.57e4,"Ti 3",&ipnt);
00240 lindst(CoolHeavy.Ti3l32,4.57e4,"Ti 3",ipnt,'c',true,
00241 " Ti III 4.57 micron, 3-2 transition, (actually multiplet) from model atom" );
00242
00243 PutLine(&TauLines[ipTi06172],
00244 " [Ti VI] 1.72 microns ");
00245
00246 PutLine(&TauLines[ipTi14212],
00247 " [Ti XIV] 2117.79 ");
00248
00249
00250 PntForLine(8823.,"V 3",&ipnt);
00251 lindst(CoolHeavy.V38830,8823,"V 3",ipnt,'c',true ,
00252 " V III 8823 ");
00253
00254
00255 PntForLine(8507.,"V 3",&ipnt);
00256 lindst(CoolHeavy.V38507,8507,"V 3",ipnt,'c',true,
00257 " V III 8507" );
00258
00259
00260 PntForLine(7735.,"V 4",&ipnt);
00261 lindst(CoolHeavy.V47741,7735,"V 4",ipnt,'c',true,
00262 " V IV 7741 1-3" );
00263
00264
00265 PntForLine(9489.,"V 4",&ipnt);
00266 lindst(CoolHeavy.V49496,9489,"V 4",ipnt,'c',true,
00267 " V IV 9496 2-1 " );
00268
00269
00270 PntForLine(4.19e4,"V 4",&ipnt);
00271 lindst(CoolHeavy.V44p2m,4.19e4,"V 4",ipnt,'c',true,
00272 " V IV 4.19 micron 3-2" );
00273
00274 PutLine(&TauLines[ipVa07130],
00275 " [V VII] 1.304 microns ");
00276
00277 PutLine(&TauLines[ipVa15172],
00278 " [V 15] 1721.38 ");
00279
00280 PntForLine(5828.,"Cr 3",&ipnt);
00281 lindst(CoolHeavy.Cr3l21,5828,"Cr 3",ipnt,'c',true,
00282 " [CrIII] multiplet blend at 5828A" );
00283
00284 PntForLine(7267.,"Cr 4",&ipnt);
00285 lindst(CoolHeavy.Cr4l21,7267,"Cr 4",ipnt,'c',true,
00286 " [CrIV] 2 - 1 multiplet blend at 7272" );
00287
00288
00289 PntForLine(6801.,"Cr 4",&ipnt);
00290 lindst(CoolHeavy.Cr4l31,6801,"Cr 4",ipnt,'c',true,
00291 " [CrIV] 3 - 1 multiplet blend at 6806" );
00292
00293
00294 PntForLine(7979.,"Cr 5",&ipnt);
00295 lindst(CoolHeavy.Cr5l21,7979,"Cr 5",ipnt,'c',true,
00296 " [CrV] 2 - 1 multiplet blend at 7985" );
00297
00298 PntForLine(6577.,"Cr 5",&ipnt);
00299 lindst(CoolHeavy.Cr5l31,6577,"Cr 5",ipnt,'c',true,
00300 " [CrV] 3 - 1 multiplet blend at 6582" );
00301
00302
00303 PntForLine(3.75e4,"Cr 5",&ipnt);
00304 lindst(CoolHeavy.Cr5l32,3.75e4,"Cr 5",ipnt,'c',true,
00305 " [CrV] 3 - 2 multiplet blend at 3.75 microns " );
00306
00307 PutLine(&TauLines[ipCr08101],
00308 " [Cr VIII] 1.01 microns ");
00309
00310 PutLine(&TauLines[ipCr16141],
00311 " [Cr 16] 1410.60 ");
00312
00313 PutLine(&TauLines[ipxMn0979],
00314 " [Mn IX] 7968.5 A ");
00315
00316 PutLine(&TauLines[ipxMn1712],
00317 " [Mn 17] 1169.59 ");
00318
00319
00320
00321
00322
00323
00324
00325
00326
00327 PutLine(&TauLines[ipFe1_24m],
00328 " Fe 1 24m ");
00329
00330 PutLine(&TauLines[ipFe1_35m],
00331 " Fe 1 35m ");
00332
00333 PutLine(&TauLines[ipFe1_54m],
00334 " Fe 1 54m ");
00335
00336 PutLine(&TauLines[ipFe1_111m],
00337 " Fe 1 111m ");
00338
00339 PutLine(&TauLines[ipFeI3884],
00340 " Fe 1 3884 ");
00341
00342 PutLine(&TauLines[ipFeI3729],
00343 " Fe 1 3729 ");
00344
00345 PutLine(&TauLines[ipFeI3457],
00346 " Fe 1 3457 ");
00347
00348 PutLine(&TauLines[ipFeI3021],
00349 " Fe 1 3021 ");
00350
00351 PutLine(&TauLines[ipFeI2966],
00352 " Fe 1 2966 ");
00353
00354 linadd(MAX2(0.,FeII.Fe2_large_cool+FeII.Fe2_UVsimp_cool),0,"Fe2c",'c' ,
00355 "total of all Fe 2 cooling, both simple UV and large atom together ");
00356
00357 linadd(MAX2(0.,-FeII.Fe2_large_cool-FeII.Fe2_UVsimp_cool),0,"Fe2h",'h' ,
00358 "total of all Fe 2 heating, both simple UV and large atom together ");
00359
00360 linadd(FeII.for7,4300,"Fe 2",'i' ,
00361 " Fe 2 forbidden 2-1 transition from Netzer's atom ");
00362
00363 PutLine(&TauLines[ipTuv3],
00364 " 2400 in simple Wills, Netzer, Wills FeII");
00365 PutLine(&TauLines[ipTr48],
00366 " 6200 in simple Wills, Netzer, Wills FeII");
00367 PutLine(&TauLines[ipTFe16],
00368 " 1080 in simple Wills, Netzer, Wills FeII");
00369 PutLine(&TauLines[ipTFe26],
00370 " 1500 in simple Wills, Netzer, Wills FeII");
00371 PutLine(&TauLines[ipTFe34],
00372 " 11500 in simple Wills, Netzer, Wills FeII");
00373 PutLine(&TauLines[ipTFe35],
00374 " 2500 in simple Wills, Netzer, Wills FeII");
00375 PutLine(&TauLines[ipTFe46],
00376 " 2300 in simple Wills, Netzer, Wills FeII");
00377 PutLine(&TauLines[ipTFe56],
00378 " 8900 in simple Wills, Netzer, Wills FeII");
00379
00380
00381
00382 FeIIAddLines();
00383
00384 for( long ipLo=0; ipLo < (FeII.nFeIILevel_malloc - 1); ipLo++ )
00385 {
00386 for( long ipHi=ipLo + 1; ipHi < FeII.nFeIILevel_malloc; ipHi++ )
00387 {
00388
00389 if( Fe2LevN[ipHi][ipLo].ipCont > 0 )
00390 PutLine( &Fe2LevN[ipHi][ipLo] ," Fe II emission" );
00391 }
00392 }
00393
00394
00395 for( i=0; i < nFeIIBands; i++ )
00396 {
00397 double SumBandInward;
00398
00399
00400
00401 eff = FeIISumBand(FeII_Bands[i][1],FeII_Bands[i][2],
00402 &SumBandInward);
00403
00404 linadd(eff,FeII_Bands[i][0],"Fe2b",'i' ,
00405 " total Fe II emission in Fe II bands, as defined in FeII_bands.ini ");
00406 linadd(SumBandInward,FeII_Bands[i][0],"Inwd",'i' ,
00407 " inward Fe II emission in Fe II bands, as defined in FeII_bands.ini ");
00408 }
00409
00410
00411 if( LineSave.ipass > 0 )
00412 {
00413
00414 if( nzone == 1 )
00415 {
00416 for( i=0; i < nFeIIConBins; i++ )
00417 {
00418
00419 FeII_Cont[i][1] = 0.;
00420 FeII_Cont[i][2] = 0.;
00421 }
00422 }
00423
00424
00425 for( i=0; i < nFeIIConBins; i++ )
00426 {
00427 double SumBandInward;
00428
00429
00430
00431 double TotalFeII = FeIISumBand(FeII_Cont[i][0],FeII_Cont[i+1][0],
00432 &SumBandInward);
00433 FeII_Cont[i][1] += (realnum)(SumBandInward*radius.dVeffAper);
00434 FeII_Cont[i][2] += (realnum)(MAX2(0.,TotalFeII-SumBandInward)*radius.dVeffAper);
00435
00436 }
00437 }
00438
00439 PutLine(&TauLines[ipT191],
00440 " anomalous Fe 2 transition at 1787, RMT 191");
00441
00442 linadd(fe.Fe3CoolTot,0,"Fe3c",'c' ,
00443 " chng 05 dec 16, FeIII code created by Kevin Blagrave Fe3c 0 - total cooling due to 14-level Fe 3 atom ");
00444
00445
00446
00447
00448
00449
00450
00451
00452
00453
00454
00455
00456
00457
00458
00459
00460
00461
00462
00463
00464
00465
00466
00467
00468
00469
00470
00471
00472
00473
00474
00475
00476
00477
00478
00479
00480
00481
00482
00483
00484
00485
00486
00487
00488
00489
00490
00491
00492
00493
00494
00495
00496
00497
00498
00499
00500
00501
00502
00503
00504
00505
00506
00507
00508
00509
00510
00511
00512
00513
00514
00515
00516
00517
00518
00519
00520
00521
00522
00523
00524
00525
00526
00527
00528
00529
00530
00531
00532
00533
00534
00535
00536
00537
00538 for( ihi=1; ihi<NLFE3; ++ihi )
00539 {
00540 for( ilo=0; ilo<ihi; ++ilo )
00541 {
00542
00543 PntForLine(fe.Fe3_wl[ihi][ilo],"Fe 3",&ipnt);
00544 # if 0
00545 fprintf( ioQQQ,"\t/* FeIII ");
00546 prt_wl( ioQQQ , (realnum)(fe.Fe3_wl[ihi][ilo]+0.5) );
00547 fprintf( ioQQQ," from Blagrave 14-level atom */\n" );
00548 # endif
00549 lindst( fe.Fe3_emiss[ihi][ilo] , (realnum)(fe.Fe3_wl[ihi][ilo]+0.5) , "Fe 3",ipnt,'c',true,
00550 " " );
00551 }
00552 }
00553
00554
00555
00556
00557
00558
00559
00560
00561
00562
00563
00564
00565
00566 PutLine(&TauLines[ipT1122]," Fe 3 1122 entire multiplet");
00567
00568 linadd(fe.Fe4CoolTot,0,"Fe4c",'i',
00569 " Fe4c 0 - total cooling due to 12-level Fe 4 atom " );
00570
00571
00572 PntForLine(3096.,"Fe 4",&ipnt);
00573 lindst(fe.fe40401,3096,"Fe 4",ipnt,'c',true,
00574 " Fe 4 3096.A, 4-1 and 5-1 transitions together" );
00575
00576
00577 PntForLine(2836.,"Fe 4",&ipnt);
00578 lindst(fe.fe42836,2836,"Fe 4",ipnt,'c',true,
00579 " Fe 4 2835.7A, 6-1 transition, 4P5/2 - 6S5/2 " );
00580
00581
00582 PntForLine(2829.,"Fe 4",&ipnt);
00583 lindst(fe.fe42829,2829,"Fe 4",ipnt,'c',true,
00584 " Fe 4 2829.4A, 7-1 transition, 4P3/2 - 6S5/2" );
00585
00586
00587 PntForLine(2567.,"Fe 4",&ipnt);
00588 lindst(fe.fe42567,2567,"Fe 4",ipnt,'c',true,
00589 " Fe 4 2567.6+ 2567.4. 11-1 and 12-1 transitions" );
00590
00591
00592 PntForLine(2.774e4,"Fe 4",&ipnt);
00593 lindst(fe.fe41207,2.774e4,"Fe 4",ipnt,'c',true,
00594 " Fe 4 2.774 microns 12-7 transition " );
00595
00596
00597 PntForLine(2.714e4,"Fe 4",&ipnt);
00598 lindst(fe.fe41206,2.714e4,"Fe 4",ipnt,'c',true,
00599 " Fe 4 2.714 microns 12-6 transition " );
00600
00601
00602 PntForLine(2.716e4,"Fe 4",&ipnt);
00603 lindst(fe.fe41106,2.716e4,"Fe 4",ipnt,'c',true,
00604 " Fe 4 2.716 microns 11-6 transition" );
00605
00606
00607 PntForLine(2.806e4,"Fe 4",&ipnt);
00608 lindst(fe.fe41007,2.806e4,"Fe 4",ipnt,'c',true,
00609 " Fe 4 2.806 microns 10-7 transition " );
00610
00611
00612 PntForLine(2.865e4,"Fe 4",&ipnt);
00613 lindst(fe.fe41008,2.865e4,"Fe 4",ipnt,'c',true ,
00614 " Fe 4 2.865 microns 10-8 transition");
00615
00616
00617 PntForLine(2.836e4,"Fe 4",&ipnt);
00618 lindst(fe.fe40906,2.836e4,"Fe 4",ipnt,'c',true,
00619 " Fe 4 2.836 microns 9-6 transition" );
00620
00621
00622 PntForLine(3892.,"Fe 5",&ipnt);
00623 lindst(CoolHeavy.c3892,3892,"Fe 5",ipnt,'c',true,
00624 " Fe 5 3892+3839" );
00625
00626 linadd(CoolHeavy.c5177,0,"Fe 6",'c' ,
00627 " all of 2G lines together first ");
00628
00629
00630 PntForLine(5177.,"Fe 6",&ipnt);
00631 lindst(CoolHeavy.c5177*0.354,5177,"Fe 6",ipnt,'c',true,
00632 " Fe 6 5177, approximate correct " );
00633
00634 linadd(fe.Fe7CoolTot,0,"Fe7c",'c' ,
00635 " Fe7c 0 - total cooling due to n-level Fe 7 atom ");
00636
00637
00638 for( ilo=0; ilo<NLFE7-1; ++ilo )
00639 {
00640
00641 for( ihi=MAX2(3,ilo+1); ihi<NLFE7; ++ihi )
00642 {
00643
00644 PntForLine(fe.Fe7_wl[ihi][ilo],"Fe 7",&ipnt);
00645 lindst( fe.Fe7_emiss[ihi][ilo] , (realnum)(fe.Fe7_wl[ihi][ilo]+0.5) , "Fe 7",ipnt,'c',true,
00646 " emission in these lines" );
00647 }
00648 }
00649 # if 0
00650 PntForLine(5721.,"Fe 7",&ipnt);
00651 lindst( fe.Fe7_5721 , 5721 , "Fe 7",ipnt,'c',true,
00652 " " );
00653
00654 PntForLine(6601.,"Fe 7",&ipnt);
00655 lindst( fe.Fe7_6601 , 6601 , "Fe 7",ipnt,'c',true,
00656 " " );
00657
00658 PntForLine(3760.,"Fe 7",&ipnt);
00659 lindst( fe.Fe7_3760 , 3760 , "Fe 7",ipnt,'c',true,
00660 " " );
00661
00662 PntForLine(3588.,"Fe 7",&ipnt);
00663 lindst( fe.Fe7_3588 , 3588 , "Fe 7",ipnt,'c',true,
00664 " " );
00665 # endif
00666
00667 PutLine(&TauLines[ipFe0795],
00668 " [Fe 7] 9.51 micron ");
00669
00670 PutLine(&TauLines[ipFe0778],
00671 " [Fe 7] 7.81 micron ");
00672
00673
00674
00675
00676
00677
00678
00679
00680
00681 PutLine(&TauLines[ipT245],
00682 " Be-seq lines ");
00683
00684 PntForLine(242.,"Fe 7",&ipnt);
00685 lindst(CoolHeavy.c242,242,"Fe 7",ipnt,'c',true,
00686 " Fe 9 242 j=1 slower decay");
00687
00688 PutLine(&TauLines[ipT352],
00689 " the E1 transition that can pump [Fe X] ");
00690
00691
00692 eff = 69.4/(69.4 + 0.27*dense.cdsqte);
00693
00694
00695
00696 c10 = TauLines[ipT352].Emis->phots*eff*.01676*TauLines[ipT352].EnergyErg*352/6374.;
00697
00698
00699
00700
00701 PutExtra( c10 );
00702
00703
00704 PutLine(&TauLines[ipFe106375],
00705 " [Fe 10] 6375, collisions with pumping too ");
00706
00707
00708
00709
00710
00711
00712
00713
00725
00726
00727
00728
00729
00730
00731
00732
00733 linadd(fe.Fe11CoolTot,0,"Fe11",'c' ,
00734 " >>chng 05 dec 18, add Fe 11 Fe11 0 - total cooling due to 5-level Fe 11 atom ");
00735
00736
00737 for( ihi=1; ihi<NLFE11; ++ihi )
00738 {
00739 for( ilo=0; ilo<ihi; ++ilo )
00740 {
00741 PntForLine(fe.Fe11_wl[ihi][ilo],"Fe11",&ipnt);
00742 lindst( fe.Fe11_emiss[ihi][ilo] , (realnum)(fe.Fe11_wl[ihi][ilo]+0.5) , "Fe11",ipnt,'c',true,
00743 " emission in these lines" );
00744 }
00745 }
00746
00747 PntForLine(1242.,"Fe12",&ipnt);
00748 lindst(CoolHeavy.c1242,1242,"Fe12",ipnt,'c',true ,
00749 " Fe 12, 1242, 1349 together, collisional excitation");
00750
00751 PntForLine(2170.,"Fe12",&ipnt);
00752 lindst(CoolHeavy.c2170,2170,"Fe12",ipnt,'c',true ,
00753 " Fe 12, 2170, 2406 together, collisional excitation");
00754
00755
00756 PntForLine(2568.,"Fe12",&ipnt);
00757 lindst(CoolHeavy.c2568,2568,"Fe12",ipnt,'c',true,
00758 " Fe12 2904, 2567, 3567, 3073 together, collisional excitation" );
00759
00760
00761
00762 linadd(fe.Fe13CoolTot,0,"Fe13",'c' ,
00763 "total cooling due to Fe 13 model atom ");
00764
00765
00766 for( ihi=1; ihi<NLFE13; ++ihi )
00767 {
00768 for( ilo=0; ilo<ihi; ++ilo )
00769 {
00770
00771 PntForLine(fe.Fe13_wl[ihi][ilo],"Fe13",&ipnt);
00772 lindst( fe.Fe13_emiss[ihi][ilo] , (realnum)(fe.Fe13_wl[ihi][ilo]+0.5) , "Fe13",ipnt,'c',true ,
00773 " Fe 13 emission");
00774 }
00775 }
00776
00777
00778 eff = 60.3/(60.3 + 0.23*dense.cdsqte/4.);
00779 r14 = dense.xIonDense[ipIRON][14-1]*fe.pfe14*rfield.flux[0][fe.ipfe10-1]*
00780 3.75e-12*eff/(1. + 24./.63*
00781 TauLines[ipT347].Emis->Pesc)*esc_PRD_1side(TauLines[ipT347].Emis->TauIn,1e-4);
00782
00783 linadd(CoolHeavy.c5303+r14,5303,"Fe14",'i',
00784 " total emission in Fe 14 5304");
00785
00786 PntForLine(5303.,"Fe14",&ipnt);
00787 lindst(CoolHeavy.c5303,5303,"Coll",ipnt,'c' ,true ,
00788 " Fe 14 5304 contribution from collisional excitation ");
00789
00790 lindst(r14,5303,"Pump",ipnt,'r' ,true ,
00791 " Fe 14 5304 continuum fluorescense ");
00792
00795 c14 = 0.;
00796 linadd(c14,5303," 347",'c' ,
00797 " collisional excitation of E1 line ");
00798
00799 PutLine(&TauLines[ipFe17_17],
00800 " Fe 17 17.1A M2 line");
00801
00802 PutLine(&TauLines[ipFe18975],
00803 " Fe 18 974.86A ");
00804
00805 PntForLine(7047.,"Fe19",&ipnt);
00806 lindst(CoolHeavy.c7082,7047,"Fe19",ipnt,'c',true,
00807 " O-like Fe19, 3P ground term, 7046.72A vacuum wl, 1328.90A >>chng 01 aug 10, updated wavelengths Fe 19 7047 '85 " );
00808
00809
00810 PntForLine(1328.,"Fe19",&ipnt);
00811 lindst(CoolHeavy.c1328,1329,"Fe19",ipnt,'c',true,
00812 " Fe 19 1329" );
00813
00814
00815 PntForLine(592.,"Fe19",&ipnt);
00816 lindst(CoolHeavy.c592,592,"Fe19",ipnt,'c',true,
00817 " Fe 19 from loulergue et al '85" );
00818
00819
00820 PntForLine(1118.,"Fe19",&ipnt);
00821 lindst(CoolHeavy.c1118,1118,"Fe19",ipnt,'c',true ,
00822 " Fe 19 from loulergue et al '85");
00823
00824 PutLine(&TauLines[ipTFe13],
00825 " next two 3p ground state lines, collisional excitation ");
00826
00827 PutLine(&TauLines[ipTFe23],
00828 " collisional excitation ");
00829
00830 PutLine(&TauLines[ipTFe20_578],
00831 " Fe20 721.40A, 578");
00832 PutLine(&TauLines[ipTFe20_721],
00833 " ");
00834
00835 linadd(
00836 TauLines[ipFe22_247].Emis->xIntensity+
00837 TauLines[ipFe22_217].Emis->xIntensity+
00838 TauLines[ipFe22_348].Emis->xIntensity+
00839 TauLines[ipFe22_292].Emis->xIntensity+
00840 TauLines[ipFe22_253].Emis->xIntensity,
00841 260,"TOTL",'i',
00842 " Fe 22 845.6A total intensity of Fe22, all lines in the multiplet " );
00843 PutLine(&TauLines[ipFe22_247],
00844 "Fe 22 247");
00845 PutLine(&TauLines[ipFe22_217],
00846 "Fe 22 217");
00847 PutLine(&TauLines[ipFe22_348],
00848 "Fe 22 348");
00849 PutLine(&TauLines[ipFe22_292],
00850 "Fe 22 292");
00851 PutLine(&TauLines[ipFe22_253],
00852 "Fe 22 253");
00853
00854
00855 PntForLine(263.,"Fe23",&ipnt);
00856 lindst(CoolHeavy.c263,263,"Fe23",ipnt,'c',true,
00857 "Fe 23 1909-like 262.6" );
00858
00859
00860 PutLine(&TauLines[ipT192],
00861 " Fe 24 only 192 of 255, 192 Li seq doublet, collisional excitation ");
00862
00863 PutLine(&TauLines[ipT255],
00864 " Fe 24 255 of 255, 192 Li seq doublet, collisional excitation ");
00865
00866 PutLine(&TauLines[ipT11],
00867 " Fe 24 Li seq 2s3p collisional excitation ");
00868
00869
00870 if( dense.lgElmtOn[ipIRON] )
00871 {
00872
00873
00874 fela = Transitions[ipH_LIKE][ipIRON][ipH2p][ipH1s].Emis->xIntensity;
00875 }
00876 else
00877 {
00878 fela = 0.;
00879 }
00880
00881
00882
00883
00884 if( dense.lgElmtOn[ipIRON] )
00885 {
00886 lindst((fe.fekcld+fe.fegrain)*1.03e-8+(fe.fekhot+fela)*1.11e-8,1.78f,"FeKa",
00887 Transitions[ipH_LIKE][ipIRON][ipH2p][ipH1s].ipCont,'i',false,
00888 "total intensity of K-alpha line" );
00889 }
00890
00891 linadd(fela*1.11e-8,2,"FeLr",'i' ,
00892 " recombination from fully stripped ion ");
00893
00894
00895 linadd((fe.fekhot+fela)*1.11e-8,2,"AugH",'i' ,
00896 " Auger hot iron, assumes case b for H and He-like ");
00897
00898 linadd(fe.fekcld*1.03e-8,2,"AugC",'i',
00899 " Auger production of cold iron, less than or 17 times ionized " );
00900
00901 linadd(fe.fegrain*1.03e-8,2,"AugG",'i' ,
00902 " grain production of cold iron ");
00903
00904 PutLine(&TauLines[ipCo11527],
00905 " [Co XI] 5168. A ");
00906
00907 PutLine(&TauLines[ipNi1_7m],
00908 " nickel [Ni I] 7m ");
00909
00910
00911
00912
00913 PutLine(&TauLines[ipNi1_11m],
00914 " [Ni I] 11m ");
00915
00916 PutLine(&TauLines[ipxNi1242],
00917 " [Ni XII] 4230.8 A ");
00918
00919
00920
00921
00922 PutLine(&TauLines[ipZn04363],
00923 "zinc iv 3.625 microns, cs and A just made up ");
00924
00925 if( trace.lgTrace )
00926 {
00927 fprintf( ioQQQ, " lines_lv1_k_zn returns\n" );
00928 }
00929 return;
00930 }