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
00041
00042 PutLine(&TauLines[ipKI7745]);
00043
00044
00045 PutLine(&TauLines[ipxK03462]);
00046
00047
00048 PutLine(&TauLines[ipxK04598]);
00049
00050
00051 PutLine(&TauLines[ipxK04154]);
00052
00053
00054 PutLine(&TauLines[ipxK06882]);
00055
00056
00057 PutLine(&TauLines[ipxK06557]);
00058
00059
00060 PutLine(&TauLines[ipxK07319]);
00061
00062
00063 PutLine(&TauLines[ipxK11425]);
00064
00065
00066
00067 PutLine(&TauLines[ipCaI4228]);
00068
00069
00070 linadd(ca.Cakh,3933,"Ca 2",'c' );
00071
00072
00073 linadd(ca.Cair,8579,"Ca 2",'c' );
00074
00075
00076 linadd(ca.c7306,7306,"Ca 2",'c' );
00077
00078
00079 linadd(ca.dCakh,3933,"Phot",'i' );
00080
00081
00082 linadd(ca.dCaf12,7306,"Phot",'i' );
00083
00084
00085 PntForLine(3934.,"Ca2K",&ipnt);
00086 lindst(ca.Cak,3934,"Ca2K",ipnt,'i',true );
00087
00088
00089 PntForLine(3969.,"Ca2H",&ipnt);
00090 lindst(ca.Cah,3969,"Ca2H",ipnt,'i',true );
00091
00092
00093 PntForLine(8498.,"Ca2X",&ipnt);
00094 lindst(ca.Cax,8498,"Ca2X",ipnt,'i',true );
00095
00096
00097 PntForLine(8542.,"Ca2Y",&ipnt);
00098 lindst(ca.Cay,8542,"Ca2Y",ipnt,'i',true );
00099
00100
00101 PntForLine(8662.,"Ca2Z",&ipnt);
00102 lindst(ca.Caz,8662,"Ca2Z",ipnt,'i',true );
00103
00104
00105 PntForLine(7291.,"CaF1",&ipnt);
00106 lindst(ca.Caf1,7291,"CaF1",ipnt,'i',true );
00107
00108
00109 PntForLine(7324.,"CaF2",&ipnt);
00110 lindst(ca.Caf2,7324,"CaF2",ipnt,'i',true );
00111
00112
00113 eff = dense.eden*dense.xIonDense[ipCALCIUM][2]*5.4e-21/(phycon.te/
00114 phycon.te10/phycon.te10);
00115 linadd(eff,3933,"Rec ",'i' );
00116
00117
00118 PutLine(&TauLines[ipTCa3]);
00119
00120
00121 PutLine(&TauLines[ipTCa4]);
00122
00123
00124 PutLine(&TauLines[ipTCa12]);
00125
00126
00127 PntForLine(6087.,"Ca 5",&ipnt);
00128 lindst(ca.Ca6087,6087,"Ca 5",ipnt,'c',true );
00129
00130
00131 PntForLine(5311.,"Ca 5",&ipnt);
00132 lindst(ca.c5311,5311,"Ca 5",ipnt,'c',true );
00133
00134
00135 PntForLine(2414.,"Ca 5",&ipnt);
00136 lindst(ca.c2414,2414,"Ca 5",ipnt,'c',true );
00137
00138
00139 PntForLine(3997.,"Ca 5",&ipnt);
00140 lindst(ca.c3997,3997,"Ca 5",ipnt,'c',true );
00141
00142
00143 PutLine(&TauLines[ipCa0741]);
00144
00145
00146 PutLine(&TauLines[ipCa0761]);
00147
00148
00149 PntForLine(5620.,"Ca 7",&ipnt);
00150 lindst(ca.Ca5620,5620,"Ca 7",ipnt,'c',true );
00151
00152
00153 PntForLine(4941.,"Ca 7",&ipnt);
00154 lindst(ca.Ca4941,4941,"Ca 7",ipnt,'c',true );
00155
00156
00157 PntForLine(2112.,"Ca 7",&ipnt);
00158 lindst(ca.Ca2112,2112,"Ca 7",ipnt,'c',true );
00159
00160
00161 PntForLine(3688.,"Ca 7",&ipnt);
00162 lindst(ca.Ca3688,3688,"Ca 7",ipnt,'c',true );
00163
00164
00165 PutLine(&TauLines[ipCa08232]);
00166
00167
00168 PutLine(&TauLines[ipCa12333]);
00169
00170
00171 PutLine(&TauLines[ipTCa302]);
00172
00173
00174 PutLine(&TauLines[ipTCa345]);
00175
00176
00177 PutLine(&TauLines[ipTCa19]);
00178
00179
00180 PntForLine(22.08e4,"Sc 2",&ipnt);
00181 lindst(CoolHeavy.Sc22p08m,21,"Sc 2",ipnt,'c',true );
00182
00183
00184 PntForLine(24.1e4,"Sc 2",&ipnt);
00185 lindst(CoolHeavy.Sc24p1m,41,"Sc 2",ipnt,'c',true );
00186
00187
00188 PntForLine(24.2e4,"Sc 2",&ipnt);
00189 lindst(CoolHeavy.Sc24p2m,42,"Sc 2",ipnt,'c',true );
00190
00191
00192 PntForLine(3933.,"Sc 3",&ipnt);
00193 lindst(CoolHeavy.Sc33936,3933,"Sc 3",ipnt,'c',true );
00194
00195
00196 PutLine(&TauLines[ipSc05231]);
00197
00198
00199 PntForLine(5054.,"Sc 6",&ipnt);
00200 lindst(CoolHeavy.Sc45058,5054,"Sc 6",ipnt,'c',true );
00201
00202
00203 PntForLine(3592.,"Sc 6",&ipnt);
00204 lindst(CoolHeavy.Sc43595,3592,"Sc 6",ipnt,'c',true );
00205
00206
00207 PntForLine(2100.,"Sc 6",&ipnt);
00208 lindst(CoolHeavy.Sc42100,2100,"Sc 6",ipnt,'c',true );
00209
00210
00211 PutLine(&TauLines[ipSc13264]);
00212
00213
00214 PntForLine(1.21e4,"Ti 3",&ipnt);
00215 lindst(CoolHeavy.Ti3l21,12,"Ti 3",ipnt,'c',true );
00216
00217
00218 PntForLine(9594.,"Ti 3",&ipnt);
00219 lindst(CoolHeavy.Ti3l31,9594,"Ti 3",ipnt,'c',true );
00220
00221
00222 PntForLine(4.57e4,"Ti 3",&ipnt);
00223 lindst(CoolHeavy.Ti3l32,45,"Ti 3",ipnt,'c',true );
00224
00225
00226 PutLine(&TauLines[ipTi06172]);
00227
00228
00229 PutLine(&TauLines[ipTi14212]);
00230
00231
00232 PntForLine(8823.,"V 3",&ipnt);
00233 lindst(CoolHeavy.V38830,8823,"V 3",ipnt,'c',true );
00234
00235
00236 PntForLine(8507.,"V 3",&ipnt);
00237 lindst(CoolHeavy.V38507,8507,"V 3",ipnt,'c',true );
00238
00239
00240 PntForLine(7735.,"V 4",&ipnt);
00241 lindst(CoolHeavy.V47741,7735,"V 4",ipnt,'c',true );
00242
00243
00244 PntForLine(9489.,"V 4",&ipnt);
00245 lindst(CoolHeavy.V49496,9489,"V 4",ipnt,'c',true );
00246
00247
00248 PntForLine(4.19e4,"V 4",&ipnt);
00249 lindst(CoolHeavy.V44p2m,42,"V 4",ipnt,'c',true );
00250
00251
00252 PutLine(&TauLines[ipVa07130]);
00253
00254
00255 PutLine(&TauLines[ipVa15172]);
00256
00257
00258 PntForLine(5828.,"Cr 3",&ipnt);
00259 lindst(CoolHeavy.Cr3l21,5828,"Cr 3",ipnt,'c',true );
00260
00261
00262 PntForLine(7267.,"Cr 4",&ipnt);
00263 lindst(CoolHeavy.Cr4l21,7267,"Cr 4",ipnt,'c',true );
00264
00265
00266 PntForLine(6801.,"Cr 4",&ipnt);
00267 lindst(CoolHeavy.Cr4l31,6801,"Cr 4",ipnt,'c',true );
00268
00269
00270 PntForLine(7979.,"Cr 5",&ipnt);
00271 lindst(CoolHeavy.Cr5l21,7979,"Cr 5",ipnt,'c',true );
00272
00273
00274 PntForLine(6577.,"Cr 5",&ipnt);
00275 lindst(CoolHeavy.Cr5l31,6577,"Cr 5",ipnt,'c',true );
00276
00277
00278 PntForLine(3.75e4,"Cr 5",&ipnt);
00279 lindst(CoolHeavy.Cr5l32,37,"Cr 5",ipnt,'c',true );
00280
00281
00282 PutLine(&TauLines[ipCr08101]);
00283
00284
00285 PutLine(&TauLines[ipCr16141]);
00286
00287
00288 PutLine(&TauLines[ipxMn0979]);
00289
00290
00291 PutLine(&TauLines[ipxMn1712]);
00292
00293
00294
00295
00296
00297
00298
00299
00300
00301
00302 PutLine(&TauLines[ipFe1_24m]);
00303
00304
00305 PutLine(&TauLines[ipFe1_35m]);
00306
00307
00308 PutLine(&TauLines[ipFe1_54m]);
00309
00310
00311 PutLine(&TauLines[ipFe1_111m]);
00312
00313
00314 PutLine(&TauLines[ipFeI3884]);
00315
00316
00317 PutLine(&TauLines[ipFeI3729]);
00318
00319
00320 PutLine(&TauLines[ipFeI3457]);
00321
00322
00323 PutLine(&TauLines[ipFeI3021]);
00324
00325
00326 PutLine(&TauLines[ipFeI2966]);
00327
00328
00329
00330 linadd(MAX2(0.,FeII.Fe2_large_cool+FeII.Fe2_UVsimp_cool),0,"Fe2c",'c' );
00331
00332
00333
00334 linadd(MAX2(0.,-FeII.Fe2_large_cool-FeII.Fe2_UVsimp_cool),0,"Fe2h",'h' );
00335
00336
00337 linadd(FeII.for7,4300,"Fe 2",'i' );
00338
00339 PutLine(&TauLines[ipTuv3]);
00340 PutLine(&TauLines[ipTr48]);
00341 PutLine(&TauLines[ipTFe16]);
00342 PutLine(&TauLines[ipTFe26]);
00343 PutLine(&TauLines[ipTFe34]);
00344 PutLine(&TauLines[ipTFe35]);
00345 PutLine(&TauLines[ipTFe46]);
00346 PutLine(&TauLines[ipTFe56]);
00347
00348
00349
00350
00351
00352
00353
00354
00355
00356
00357 PntForLine(16640.,"Fe 2",&ipnt);
00358 lindst(FeII.fe21308,16640,"Fe 2",ipnt,'i',true );
00359
00360
00361 PntForLine(15990.,"Fe 2",&ipnt);
00362 lindst(FeII.fe21207,15990,"Fe 2",ipnt,'i',true );
00363
00364
00365 PntForLine(15340.,"Fe 2",&ipnt);
00366 lindst(FeII.fe21106,15340,"Fe 2",ipnt,'i',true );
00367
00368
00369 PntForLine(16440.,"Fe 2",&ipnt);
00370 lindst(FeII.fe21006,16440,"Fe 2",ipnt,'i',true );
00371
00372
00373 PntForLine(12790.,"Fe 2",&ipnt);
00374 lindst(FeII.fe21204,12790,"Fe 2",ipnt,'i',true );
00375
00376
00377 PntForLine(12950.,"Fe 2",&ipnt);
00378 lindst(FeII.fe21103,12950,"Fe 2",ipnt,'i',true );
00379
00380
00381 PntForLine(13280.,"Fe 2",&ipnt);
00382 lindst(FeII.fe21104,13280,"Fe 2",ipnt,'i',true );
00383
00384
00385 PntForLine(12570.,"Fe 2",&ipnt);
00386 lindst(FeII.fe21001,12570,"Fe 2",ipnt,'i',true );
00387
00388
00389 PntForLine(13210.,"Fe 2",&ipnt);
00390 lindst(FeII.fe21002,13210,"Fe 2",ipnt,'i',true );
00391
00392
00393 PntForLine(259880.,"Fe 2",&ipnt);
00394 lindst(FeII.fe20201,259880,"Fe 2",ipnt,'i',true );
00395
00396
00397 PntForLine(353480.,"Fe 2",&ipnt);
00398 lindst(FeII.fe20302,353480,"Fe 2",ipnt,'i',true );
00399
00400
00401 PntForLine(179360.,"Fe 2",&ipnt);
00402 lindst(FeII.fe20706,179360,"Fe 2",ipnt,'i',true );
00403
00404
00405 PntForLine(245180.,"Fe 2",&ipnt);
00406 lindst(FeII.fe20807,245180,"Fe 2",ipnt,'i',true );
00407
00408
00409 PntForLine(357760.,"Fe 2",&ipnt);
00410 lindst(FeII.fe20908,357760,"Fe 2",ipnt,'i',true );
00411
00412
00413 PntForLine(18100.,"Fe 2",&ipnt);
00414 lindst(FeII.fe21007,18100,"Fe 2",ipnt,'i',true );
00415
00416
00417 PntForLine(16770.,"Fe 2",&ipnt);
00418 lindst(FeII.fe21107,16770,"Fe 2",ipnt,'i',true );
00419
00420
00421 PntForLine(18000.,"Fe 2",&ipnt);
00422 lindst(FeII.fe21108,18000,"Fe 2",ipnt,'i',true );
00423
00424
00425 PntForLine(17120.,"Fe 2",&ipnt);
00426 lindst(FeII.fe21208,17120,"Fe 2",ipnt,'i',true );
00427
00428
00429 PntForLine(17980.,"Fe 2",&ipnt);
00430 lindst(FeII.fe21209,17980,"Fe 2",ipnt,'i',true );
00431
00432
00433 PntForLine(229020.,"Fe 2",&ipnt);
00434 lindst(FeII.fe21110,229020,"Fe 2",ipnt,'i',true );
00435
00436
00437 PntForLine(346600.,"Fe 2",&ipnt);
00438 lindst(FeII.fe21211,346600,"Fe 2",ipnt,'i',true );
00439
00440
00441
00442
00443
00444
00445 PntForLine(8617.,"Fe 2",&ipnt);
00446 lindst(FeII.fe21406,8617,"Fe 2",ipnt,'i',true );
00447
00448
00449 PntForLine(8892.,"Fe 2",&ipnt);
00450 lindst(FeII.fe21507,8892,"Fe 2",ipnt,'i',true );
00451
00452
00453 PntForLine(9227.,"Fe 2",&ipnt);
00454 lindst(FeII.fe21508,9227,"Fe 2",ipnt,'i',true );
00455
00456
00457 PntForLine(9268.,"Fe 2",&ipnt);
00458 lindst(FeII.fe21609,9268,"Fe 2",ipnt,'i',true );
00459
00460
00461
00462
00463
00464 PntForLine(5159.,"Fe 2",&ipnt);
00465 lindst(FeII.fe25to6,5159,"Fe 2",ipnt,'i',true );
00466
00467
00468 PntForLine(5262.,"Fe 2",&ipnt);
00469 lindst(FeII.fe27to7,5262,"Fe 2",ipnt,'i',true );
00470
00471
00472 PntForLine(5334.,"Fe 2",&ipnt);
00473 lindst(FeII.fe28to8,5334,"Fe 2",ipnt,'i',true );
00474
00475
00476 PntForLine(5376.,"Fe 2",&ipnt);
00477 lindst(FeII.fe29to9,5376,"Fe 2",ipnt,'i',true );
00478
00479
00480 PntForLine(4815.,"Fe 2",&ipnt);
00481 lindst(FeII.fe32to6,4815,"Fe 2",ipnt,'i',true );
00482
00483
00484 PntForLine(4905.,"Fe 2",&ipnt);
00485 lindst(FeII.fe33to7,4905,"Fe 2",ipnt,'i',true );
00486
00487
00488 PntForLine(4347.,"Fe 2",&ipnt);
00489 lindst(FeII.fe37to7,4347,"Fe 2",ipnt,'i',true );
00490
00491
00492 PntForLine(4353.,"Fe 2",&ipnt);
00493 lindst(FeII.fe39to8,4353,"Fe 2",ipnt,'i',true );
00494
00495
00496 PntForLine(4372.,"Fe 2",&ipnt);
00497 lindst(FeII.fe40to9,4372,"Fe 2",ipnt,'i',true );
00498
00499
00500 PntForLine(4244.,"Fe 2",&ipnt);
00501 lindst(FeII.fe37to6,4244,"Fe 2",ipnt,'i',true );
00502
00503
00504 PntForLine(4277.,"Fe 2",&ipnt);
00505 lindst(FeII.fe39to7,4277,"Fe 2",ipnt,'i',true );
00506
00507
00508 PntForLine(4320.,"Fe 2",&ipnt);
00509 lindst(FeII.fe40to8,4320,"Fe 2",ipnt,'i',true );
00510
00511
00512 PntForLine(4358.,"Fe 2",&ipnt);
00513 lindst(FeII.fe41to9,4358,"Fe 2",ipnt,'i',true );
00514
00515
00516 PntForLine(4177.,"Fe 2",&ipnt);
00517 lindst(FeII.fe39to6,4177,"Fe 2",ipnt,'i',true );
00518
00519
00520 PntForLine(4245.,"Fe 2",&ipnt);
00521 lindst(FeII.fe40to7,4245,"Fe 2",ipnt,'i',true );
00522
00523
00524 PntForLine(4306.,"Fe 2",&ipnt);
00525 lindst(FeII.fe41to8,4306,"Fe 2",ipnt,'i',true );
00526
00527
00528
00529
00530 PntForLine(4114.,"Fe 2",&ipnt);
00531 lindst(FeII.fe42to6,4114,"Fe 2",ipnt,'i',true );
00532
00533
00534 PntForLine(4179.,"Fe 2",&ipnt);
00535 lindst(FeII.fe43to7,4179,"Fe 2",ipnt,'i',true );
00536
00537
00538 PntForLine(4211.,"Fe 2",&ipnt);
00539 lindst(FeII.fe42to7,4211,"Fe 2",ipnt,'i',true );
00540
00541
00542 PntForLine(4314.,"Fe 2",&ipnt);
00543 lindst(FeII.fe80to28,4314,"Fe 2",ipnt,'i',true );
00544
00545
00546 PntForLine(4359.,"Fe 2",&ipnt);
00547 lindst(FeII.fe36to2,4359,"Fe 2",ipnt,'i',true );
00548
00549
00550 PntForLine(4414.,"Fe 2",&ipnt);
00551 lindst(FeII.fe36to3,4414,"Fe 2",ipnt,'i',true );
00552
00553
00554 PntForLine(4416.,"Fe 2",&ipnt);
00555 lindst(FeII.fe32to1,4416,"Fe 2",ipnt,'i',true );
00556
00557
00558 PntForLine(4458.,"Fe 2",&ipnt);
00559 lindst(FeII.fe33to2,4458,"Fe 2",ipnt,'i',true );
00560
00561
00562 PntForLine(4475.,"Fe 2",&ipnt);
00563 lindst(FeII.fe36to5,4475,"Fe 2",ipnt,'i',true );
00564
00565
00566 PntForLine(4493.,"Fe 2",&ipnt);
00567 lindst(FeII.fe32to2,4493,"Fe 2",ipnt,'i',true );
00568
00569
00570 PntForLine(4515.,"Fe 2",&ipnt);
00571 lindst(FeII.fe33to3,4515,"Fe 2",ipnt,'i',true );
00572
00573
00574 PntForLine(4728.,"Fe 2",&ipnt);
00575 lindst(FeII.fe30to3,4728,"Fe 2",ipnt,'i',true );
00576
00577
00578 PntForLine(4775.,"Fe 2",&ipnt);
00579 lindst(FeII.fe33to6,4775,"Fe 2",ipnt,'i',true );
00580
00581
00582 PntForLine(4890.,"Fe 2",&ipnt);
00583 lindst(FeII.fe24to2,4890,"Fe 2",ipnt,'i',true );
00584
00585
00586 PntForLine(4947.,"Fe 2",&ipnt);
00587 lindst(FeII.fe32to7,4947,"Fe 2",ipnt,'i',true );
00588
00589
00590 PntForLine(4951.,"Fe 2",&ipnt);
00591 lindst(FeII.fe35to8,4951,"Fe 2",ipnt,'i',true );
00592
00593
00594 PntForLine(4973.,"Fe 2",&ipnt);
00595 lindst(FeII.fe34to8,4973,"Fe 2",ipnt,'i',true );
00596
00597
00598 PntForLine(5112.,"Fe 2",&ipnt);
00599 lindst(FeII.fe27to6,5112,"Fe 2",ipnt,'i',true );
00600
00601
00602 PntForLine(5220.,"Fe 2",&ipnt);
00603 lindst(FeII.fe28to7,5220,"Fe 2",ipnt,'i',true );
00604
00605
00606 PntForLine(5269.,"Fe 2",&ipnt);
00607 lindst(FeII.fe30to8,5269,"Fe 2",ipnt,'i',true );
00608
00609
00610 PntForLine(5273.,"Fe 2",&ipnt);
00611 lindst(FeII.fe24to6,5273,"Fe 2",ipnt,'i',true );
00612
00613
00614 PntForLine(5297.,"Fe 2",&ipnt);
00615 lindst(FeII.fe29to8,5297,"Fe 2",ipnt,'i',true );
00616
00617
00618 PntForLine(5433.,"Fe 2",&ipnt);
00619 lindst(FeII.fe24to7,5433,"Fe 2",ipnt,'i',true );
00620
00621
00622 PntForLine(5527.,"Fe 2",&ipnt);
00623 lindst(FeII.fe22to7,5527,"Fe 2",ipnt,'i',true );
00624
00625
00626 PntForLine(5747.,"Fe 2",&ipnt);
00627 lindst(FeII.fe38to11,5747,"Fe 2",ipnt,'i',true );
00628
00629
00630 PntForLine(6440.,"Fe 2",&ipnt);
00631 lindst(FeII.fe19to8,6440,"Fe 2",ipnt,'i',true );
00632
00633
00634 PntForLine(7155.,"Fe 2",&ipnt);
00635 lindst(FeII.fe17to6,7155,"Fe 2",ipnt,'i',true );
00636
00637
00638 PntForLine(7172.,"Fe 2",&ipnt);
00639 lindst(FeII.fe18to7,7172,"Fe 2",ipnt,'i',true );
00640
00641
00642 PntForLine(7388.,"Fe 2",&ipnt);
00643 lindst(FeII.fe18to8,7388,"Fe 2",ipnt,'i',true );
00644
00645
00646 PntForLine(7453.,"Fe 2",&ipnt);
00647 lindst(FeII.fe17to7,7453,"Fe 2",ipnt,'i',true );
00648
00649
00650
00651 FeIIAddLines();
00652
00653
00654
00655 if( FeII.lgFeIILargeOn )
00656 {
00657 for( i=0; i < nFeIIBands; i++ )
00658 {
00659
00660
00661
00662 eff = FeIISumBand(FeII_Bands[i][1],FeII_Bands[i][2]);
00663
00664 linadd(eff,FeII_Bands[i][0],"Fe2b",'i' );
00665 }
00667 if( LineSave.ipass > 0 )
00668 {
00669 for( i=0; i < nFeIIConBins; i++ )
00670 {
00671
00672
00673
00674 FeII_Cont[i][0] +=
00675 (float)(FeIISumBand(FeII_Cont[i][1],FeII_Cont[i][2])*radius.dVeff);
00676
00677 }
00678 }
00679 else
00680 {
00681 for( i=0; i < nFeIIConBins; i++ )
00682 {
00683
00684 FeII_Cont[i][0] = 0.;
00685 }
00686 }
00687
00688 # if 0
00689
00690 for( i=0; i < nFeIIConBins; i++ )
00691 {
00692
00693 linadd(0., FeII_Cont[i][0] ,"Fe2c",'i' );
00694 }
00695 # endif
00696 }
00697 else
00698 {
00699
00700 for( i=0; i < nFeIIBands; i++ )
00701 {
00702
00703 linadd(0., FeII_Bands[i][0] ,"Fe2b",'i' );
00704 }
00705
00706 # if 0
00707 for( i=0; i < nFeIIConBins; i++ )
00708 {
00709
00710 linadd(0., FeII_Cont[i][0] ,"Fe2c",'i' );
00711 }
00712 # endif
00713 }
00714
00715
00716 PutLine(&TauLines[ipT191]);
00717
00718
00719
00720 linadd(fe.Fe3CoolTot,0,"Fe3c",'c' );
00721
00722
00723
00724
00725
00726
00727
00728
00729
00730
00731
00732
00733
00734
00735
00736
00737
00738
00739
00740
00741
00742
00743
00744
00745
00746
00747
00748
00749
00750
00751
00752
00753
00754
00755
00756
00757
00758
00759
00760
00761
00762
00763
00764
00765
00766
00767
00768
00769
00770
00771
00772
00773
00774
00775
00776
00777
00778
00779
00780
00781
00782
00783
00784
00785
00786
00787
00788
00789
00790
00791
00792
00793
00794
00795
00796
00797
00798
00799
00800
00801
00802
00803
00804
00805
00806
00807
00808
00809
00810
00811
00812
00813
00814
00815 for( ihi=1; ihi<NLFE3; ++ihi )
00816 {
00817 for( ilo=0; ilo<ihi; ++ilo )
00818 {
00819
00820 PntForLine(fe.Fe3_wl[ihi][ilo],"Fe 3",&ipnt);
00821 # if 0
00822 fprintf( ioQQQ,"\t/* FeIII ");
00823 prt_wl( ioQQQ , (float)(fe.Fe3_wl[ihi][ilo]+0.5) );
00824 fprintf( ioQQQ," from Blagrave 14-level atom */\n ");
00825 # endif
00826 lindst( fe.Fe3_emiss[ihi][ilo] , (float)(fe.Fe3_wl[ihi][ilo]+0.5) , "Fe 3",ipnt,'c',true );
00827 }
00828 }
00829
00830
00831
00832
00833
00834
00835
00836
00837
00838
00839
00840
00841
00842
00843 PutLine(&TauLines[ipT1122]);
00844
00845
00846 linadd(fe.Fe4CoolTot,0,"Fe4c",'c' );
00847
00848
00849 PntForLine(3096.,"Fe 4",&ipnt);
00850 lindst(fe.fe40401,3096,"Fe 4",ipnt,'i',true );
00851
00852
00853 PntForLine(2836.,"Fe 4",&ipnt);
00854 lindst(fe.fe42836,2836,"Fe 4",ipnt,'i',true );
00855
00856
00857 PntForLine(2829.,"Fe 4",&ipnt);
00858 lindst(fe.fe42829,2829,"Fe 4",ipnt,'i',true );
00859
00860
00861 PntForLine(2567.,"Fe 4",&ipnt);
00862 lindst(fe.fe42567,2567,"Fe 4",ipnt,'i',true );
00863
00864
00865 PntForLine(2.774e4,"Fe 4",&ipnt);
00866 lindst(fe.fe41207,277,"Fe 4",ipnt,'i',true );
00867
00868
00869 PntForLine(2.714e4,"Fe 4",&ipnt);
00870 lindst(fe.fe41206,271,"Fe 4",ipnt,'i',true );
00871
00872
00873 PntForLine(2.716e4,"Fe 4",&ipnt);
00874 lindst(fe.fe41106,272,"Fe 4",ipnt,'i',true );
00875
00876
00877 PntForLine(2.806e4,"Fe 4",&ipnt);
00878 lindst(fe.fe41007,281,"Fe 4",ipnt,'i',true );
00879
00880
00881 PntForLine(2.865e4,"Fe 4",&ipnt);
00882 lindst(fe.fe41008,287,"Fe 4",ipnt,'i',true );
00883
00884
00885 PntForLine(2.836e4,"Fe 4",&ipnt);
00886 lindst(fe.fe40906,284,"Fe 4",ipnt,'i',true );
00887
00888
00889 PntForLine(3892.,"Fe 5",&ipnt);
00890 lindst(CoolHeavy.c3892,3892,"Fe 5",ipnt,'c',true );
00891
00892
00893 linadd(CoolHeavy.c5177,0,"Fe 6",'c' );
00894
00895
00896 PntForLine(5177.,"Fe 6",&ipnt);
00897 lindst(CoolHeavy.c5177*0.354,5177,"Fe 6",ipnt,'c',true );
00898
00899
00900 linadd(fe.Fe7CoolTot,0,"Fe7c",'c' );
00901
00902
00903 for( ilo=0; ilo<NLFE7-1; ++ilo )
00904 {
00905
00906 for( ihi=MAX2(3,ilo+1); ihi<NLFE7; ++ihi )
00907 {
00908
00909 PntForLine(fe.Fe7_wl[ihi][ilo],"Fe 7",&ipnt);
00910 lindst( fe.Fe7_emiss[ihi][ilo] , (float)(fe.Fe7_wl[ihi][ilo]+0.5) , "Fe 7",ipnt,'c',true );
00911 }
00912 }
00913 # if 0
00914 PntForLine(5721.,"Fe 7",&ipnt);
00915 lindst( fe.Fe7_5721 , 5721 , "Fe 7",ipnt,'c',true );
00916
00917 PntForLine(6601.,"Fe 7",&ipnt);
00918 lindst( fe.Fe7_6601 , 6601 , "Fe 7",ipnt,'c',true );
00919
00920 PntForLine(3760.,"Fe 7",&ipnt);
00921 lindst( fe.Fe7_3760 , 3760 , "Fe 7",ipnt,'c',true );
00922
00923 PntForLine(3588.,"Fe 7",&ipnt);
00924 lindst( fe.Fe7_3588 , 3588 , "Fe 7",ipnt,'c',true );
00925 # endif
00926
00927
00928 PutLine(&TauLines[ipFe0795]);
00929
00930
00931 PutLine(&TauLines[ipFe0778]);
00932
00933
00934
00935
00936
00937
00938
00939
00940
00941
00942 PutLine(&TauLines[ipT245]);
00943
00944
00945 PntForLine(242.,"Fe 7",&ipnt);
00946 lindst(CoolHeavy.c242,242,"Fe 7",ipnt,'c',true );
00947
00948
00949 PutLine(&TauLines[ipT352]);
00950
00951
00952 eff = 69.4/(69.4 + 0.27*dense.cdsqte);
00953
00954
00955
00956 c10 = TauLines[ipT352].phots*eff*.01676*TauLines[ipT352].EnergyErg*352/6374.;
00957
00958
00959
00960
00961 PutExtra( c10 );
00962
00963
00964 PutLine(&TauLines[ipFe106375]);
00965
00966
00967
00968
00969
00970
00971
00972
00984
00985
00986
00987
00988
00989
00990
00991
00992
00993
00994 linadd(fe.Fe11CoolTot,0,"Fe11",'c' );
00995
00996
00997 for( ihi=1; ihi<NLFE11; ++ihi )
00998 {
00999 for( ilo=0; ilo<ihi; ++ilo )
01000 {
01001
01002 PntForLine(fe.Fe11_wl[ihi][ilo],"Fe11",&ipnt);
01003 lindst( fe.Fe11_emiss[ihi][ilo] , (float)(fe.Fe11_wl[ihi][ilo]+0.5) , "Fe11",ipnt,'c',true );
01004 }
01005 }
01006
01007
01008 PntForLine(1242.,"Fe12",&ipnt);
01009 lindst(CoolHeavy.c1242,1242,"Fe12",ipnt,'c',true );
01010
01011
01012 PntForLine(2170.,"Fe12",&ipnt);
01013 lindst(CoolHeavy.c2170,2170,"Fe12",ipnt,'c',true );
01014
01015
01016 PntForLine(2568.,"Fe12",&ipnt);
01017 lindst(CoolHeavy.c2568,2568,"Fe12",ipnt,'c',true );
01018
01019
01020
01021
01022
01023
01024
01025
01026
01027
01028
01029
01030
01031
01032
01033
01034
01035
01036
01037
01038
01039
01040
01041
01042
01043 linadd(fe.Fe13CoolTot,0,"Fe13",'c' );
01044
01045
01046 for( ihi=1; ihi<NLFE13; ++ihi )
01047 {
01048 for( ilo=0; ilo<ihi; ++ilo )
01049 {
01050
01051 PntForLine(fe.Fe13_wl[ihi][ilo],"Fe13",&ipnt);
01052 lindst( fe.Fe13_emiss[ihi][ilo] , (float)(fe.Fe13_wl[ihi][ilo]+0.5) , "Fe13",ipnt,'c',true );
01053 }
01054 }
01055
01056
01057 eff = 60.3/(60.3 + 0.23*dense.cdsqte/4.);
01058 r14 = dense.xIonDense[ipIRON][14-1]*fe.pfe14*rfield.flux[fe.ipfe10-1]*
01059 3.75e-12*eff/(1. + 24./.63*
01060 TauLines[ipT347].Pesc)*esc_PRD_1side(TauLines[ipT347].TauIn,1e-4);
01061 PntForLine(5303.,"Fe14",&ipnt);
01062
01063
01064 lindst(CoolHeavy.c5303+r14,5303,"Fe14",ipnt,'i',true );
01065
01066
01067 linadd(CoolHeavy.c5303,5303,"Coll",'c' );
01068
01069
01070 linadd(r14,5303,"Pump",'r' );
01071
01073
01074 c14 = 0.;
01075 linadd(c14,5303," 347",'c' );
01076
01077
01078 PutLine(&TauLines[ipFe18975]);
01079
01080
01081
01082
01083 PntForLine(7047.,"Fe19",&ipnt);
01084 lindst(CoolHeavy.c7082,7047,"Fe19",ipnt,'c',true );
01085
01086
01087 PntForLine(1328.,"Fe19",&ipnt);
01088 lindst(CoolHeavy.c1328,1329,"Fe19",ipnt,'c',true );
01089
01090
01091 PntForLine(592.,"Fe19",&ipnt);
01092 lindst(CoolHeavy.c592,592,"Fe19",ipnt,'c',true );
01093
01094
01095 PntForLine(1118.,"Fe19",&ipnt);
01096 lindst(CoolHeavy.c1118,1118,"Fe19",ipnt,'c',true );
01097
01098
01099 PutLine(&TauLines[ipTFe13]);
01100
01101
01102 PutLine(&TauLines[ipTFe23]);
01103
01104
01105 PutLine(&TauLines[ipTFe20_578]);
01106 PutLine(&TauLines[ipTFe20_721]);
01107
01108
01109
01110 linadd(
01111 TauLines[ipFe22_247].xIntensity+
01112 TauLines[ipFe22_217].xIntensity+
01113 TauLines[ipFe22_348].xIntensity+
01114 TauLines[ipFe22_292].xIntensity+
01115 TauLines[ipFe22_253].xIntensity,
01116 260,"TOTL",'i' );
01117 PutLine(&TauLines[ipFe22_247]);
01118 PutLine(&TauLines[ipFe22_217]);
01119 PutLine(&TauLines[ipFe22_348]);
01120 PutLine(&TauLines[ipFe22_292]);
01121 PutLine(&TauLines[ipFe22_253]);
01122
01123
01124 PntForLine(263.,"Fe23",&ipnt);
01125 lindst(CoolHeavy.c263,263,"Fe23",ipnt,'c',true );
01126
01127
01128 PutLine(&TauLines[ipT192]);
01129
01130
01131 PutLine(&TauLines[ipT255]);
01132
01133
01134 PutLine(&TauLines[ipT11]);
01135
01136
01137 if( dense.lgElmtOn[ipIRON] )
01138 {
01139
01140
01141 fela = EmisLines[ipH_LIKE][ipIRON][ipH2p][ipH1s].xIntensity;
01142 }
01143 else
01144 {
01145 fela = 0.;
01146 }
01147
01148
01149
01150
01151 if( dense.lgElmtOn[ipIRON] )
01152 {
01153 lindst((fe.fekcld+fe.fegrain)*1.03e-8+(fe.fekhot+fela)*1.11e-8,2,"FeKa",
01154 EmisLines[ipH_LIKE][ipIRON][ipH2p][ipH1s].ipCont,'i',false );
01155 }
01156
01157
01158 linadd(fela*1.11e-8,2,"FeLr",'i' );
01159
01160
01161
01162 linadd((fe.fekhot+fela)*1.11e-8,2,"AugH",'i' );
01163
01164
01165 linadd(fe.fekcld*1.03e-8,2,"AugC",'i' );
01166
01167
01168 linadd(fe.fegrain*1.03e-8,2,"AugG",'i' );
01169
01170
01171 PutLine(&TauLines[ipCo11527]);
01172
01173
01174
01175
01176 PutLine(&TauLines[ipNi1_7m]);
01177
01178
01179 PutLine(&TauLines[ipNi1_11m]);
01180
01181
01182 PutLine(&TauLines[ipxNi1242]);
01183
01184
01185
01186
01187
01188 PutLine(&TauLines[ipZn04363]);
01189
01190 if( trace.lgTrace )
01191 {
01192 fprintf( ioQQQ, " lines_lv1_k_zn returns\n" );
01193 }
01194
01195 DEBUG_EXIT( "lines_lv1_k_zn()" );
01196 return;
01197 }