00001
00002
00003
00004 #include "cddefines.h"
00005 #include "coolheavy.h"
00006 #include "sil.h"
00007 #include "phycon.h"
00008 #include "embesq.h"
00009 #include "taulines.h"
00010 #include "dense.h"
00011 #include "ionbal.h"
00012 #include "trace.h"
00013 #include "lines_service.h"
00014 #include "lines.h"
00015
00016 void lines_lv1_na_ar(void)
00017 {
00018 long int ipnt;
00019 double drec,
00020 fac,
00021 rec,
00022 sum,
00023 t4;
00024
00025 DEBUG_ENTRY( "lines_lv1_na_ar()" );
00026
00027 if( trace.lgTrace )
00028 {
00029 fprintf( ioQQQ, " lines_lv1_na_ar called\n" );
00030 }
00031
00032 t4 = phycon.te/1e4;
00033
00034
00035
00036 PutLine(&TauLines[ipT5895]);
00037
00038
00039 PutLine(&TauLines[ipfsNa373]);
00040
00041
00042 PutLine(&TauLines[ipfsNa490]);
00043
00044
00045 PutLine(&TauLines[ipfsNa421]);
00046
00047
00048 PntForLine(1365.,"Na 5",&ipnt);
00049 lindst(CoolHeavy.c1365,1365,"Na 5",ipnt,'c',true );
00050
00051
00052 PntForLine(2067.,"Na 5",&ipnt);
00053 lindst(CoolHeavy.c2067,2067,"Na 5",ipnt,'c',true );
00054
00055
00056 PntForLine(4017.,"Na 5",&ipnt);
00057 lindst(CoolHeavy.c4017,4017,"Na 5",ipnt,'c',true );
00058
00059
00060 PntForLine(2569.,"Na 6",&ipnt);
00061 lindst(CoolHeavy.c2569,2569,"Na 6",ipnt,'c',true );
00062
00063
00064 PntForLine(1357.,"Na 6",&ipnt);
00065 lindst(CoolHeavy.c1357,1357,"Na 6",ipnt,'c',true );
00066
00067
00068 PntForLine(2972.,"Na 6",&ipnt);
00069 lindst(CoolHeavy.c2972/(1.+1./3.02),2972,"Na 6",ipnt,'c',true );
00070
00071
00072 PntForLine(2872.,"Na 6",&ipnt);
00073 lindst(CoolHeavy.c2972/(1.+3.02),2872,"Na 6",ipnt,'c',true );
00074
00075
00076 PutLine(&TauLines[ipxNa6143]);
00077
00078
00079 PutLine(&TauLines[ipxNa6862]);
00080
00081
00082 PutLine(&TauLines[ipxNa0746]);
00083
00084
00085 PutLine(&TauLines[ipT4561]);
00086
00087
00088 PutLine(&TauLines[ipMgI2853]);
00089
00090
00091 PutLine(&TauLines[ipMgI2026]);
00092
00093
00094 linadd(TauLines[ipT2796].xIntensity+TauLines[ipT2804].xIntensity,2798,
00095 "TOTL",'i' );
00096
00097 sum = TauLines[ipT2796].xIntensity*TauLines[ipT2796].FracInwd +
00098 TauLines[ipT2804].xIntensity*TauLines[ipT2804].FracInwd;
00099
00100
00101 linadd(sum,2798,"Inwd",'i' );
00102 PutLine(&TauLines[ipT2796]);
00103 PutLine(&TauLines[ipT2804]);
00104
00105
00106 PutLine(&TauLines[ipTMg4]);
00107
00108
00109 PutLine(&TauLines[ipTMg14]);
00110
00111
00112 PutLine(&TauLines[ipTMg6]);
00113
00114
00115 PutLine(&TauLines[ipxMg52855]);
00116
00117
00118 PutLine(&TauLines[ipxMg52417]);
00119
00120
00121 PutLine(&TauLines[ipxMg51325]);
00122
00123
00124 PntForLine(1806.,"Mg 6",&ipnt);
00125 lindst(CoolHeavy.c1806,1806,"Mg 6",ipnt,'c',true );
00126
00127
00128 PutLine(&TauLines[ipxMg72569]);
00129
00130
00131 PutLine(&TauLines[ipxMg72261]);
00132
00133
00134 PutLine(&TauLines[ipxMg71190]);
00135
00136
00137 PutLine(&TauLines[ipfsMg755]);
00138
00139
00140 PutLine(&TauLines[ipfsMg790]);
00141
00142
00143 PutLine(&TauLines[ipxMg08303]);
00144
00145
00146 PutLine(&TauLines[ipT705]);
00147
00148
00149 linadd(TauLines[ipTMg610].xIntensity+TauLines[ipTMg625].xIntensity,615,"TOTL",'i' );
00150 PutLine(&TauLines[ipTMg610]);
00151 PutLine(&TauLines[ipTMg625]);
00152
00153
00154
00155 PutLine(&TauLines[ipT58]);
00156 linadd(TauLines[ipTMg610].ots*TauLines[ipTMg610].EnergyErg+
00157 TauLines[ipTMg625].ots*TauLines[ipTMg625].EnergyErg,615,"dest",'i' );
00158
00159
00160
00161 PutLine(&TauLines[ipAlI3957]);
00162
00163
00164 PutLine(&TauLines[ipAlI3090]);
00165
00166
00167
00168 linadd(embesq.em2669+TauLines[ipT2670].xIntensity,2665,"totl",'i' );
00169 PutLine(&TauLines[ipT2670]);
00170
00171
00172 linadd(embesq.em2669,2660,"Al 2",'i' );
00173
00174
00175 linadd(TauLines[ipT1855].xIntensity+TauLines[ipT1863].xIntensity,1860,"TOTL",'i' );
00176 sum = TauLines[ipT1855].xIntensity*TauLines[ipT1855].FracInwd +
00177 TauLines[ipT1863].xIntensity* TauLines[ipT1863].FracInwd;
00178
00179
00180 linadd(sum,1860,"Inwd",'i' );
00181 PutLine(&TauLines[ipT1855]);
00182 PutLine(&TauLines[ipT1863]);
00183
00184
00185 PutLine(&TauLines[ipAl529]);
00186
00187
00188 PutLine(&TauLines[ipAl6366]);
00189
00190
00191 PutLine(&TauLines[ipAl6912]);
00192
00193
00194 PntForLine(2428.,"Al 6",&ipnt);
00195 lindst(CoolHeavy.c2428/(1.+1./3.73),2428,"Al 6",ipnt,'c',true );
00196
00197
00198 PntForLine(2601.,"Al 6",&ipnt);
00199 lindst(CoolHeavy.c2428/(1.+3.73),2601,"Al 6",ipnt,'c',true );
00200
00201
00202 PntForLine(1170.,"Al 6",&ipnt);
00203 lindst(CoolHeavy.c1170,1170,"Al 6",ipnt,'c',true );
00204
00205
00206 PntForLine(2125.,"Al 6",&ipnt);
00207 lindst(CoolHeavy.c2125,2125,"Al 6",ipnt,'c',true );
00208
00209
00210 PutLine(&TauLines[ipAl8575]);
00211
00212
00213 PutLine(&TauLines[ipAl8370]);
00214
00215
00216 PutLine(&TauLines[ipAl09204]);
00217
00218
00219 PutLine(&TauLines[ipT639]);
00220
00221
00222 linadd(TauLines[ipTAl550].xIntensity+TauLines[ipTAl568].xIntensity,556,"TOTL",'i' );
00223 PutLine(&TauLines[ipTAl550]);
00224 PutLine(&TauLines[ipTAl568]);
00225
00226
00227 PutLine(&TauLines[ipTAl48]);
00228
00229
00230
00231
00232 PutLine(&TauLines[ipSi1_130m]);
00233
00234
00235 PutLine(&TauLines[ipSi1_68m]);
00236
00237
00238 PutLine(&TauLines[ipSii2518]);
00239
00240
00241 PutLine(&TauLines[ipSii2215]);
00242
00243
00244 PutLine(&TauLines[ipTSi35]);
00245
00246
00247 linadd(
00248 TauLines[ipSi2_2334].xIntensity+
00249 TauLines[ipSi2_2329].xIntensity+
00250 TauLines[ipSi2_2350].xIntensity+
00251 TauLines[ipSi2_2344].xIntensity+
00252 TauLines[ipSi2_2336].xIntensity,
00253 2335,"TOTL",'i' );
00254 PutLine(&TauLines[ipSi2_2334]);
00255 PutLine(&TauLines[ipSi2_2329]);
00256 PutLine(&TauLines[ipSi2_2350]);
00257 PutLine(&TauLines[ipSi2_2344]);
00258 PutLine(&TauLines[ipSi2_2336]);
00259
00260
00261 PutLine(&TauLines[ipT1808]);
00262
00263
00264 PutLine(&TauLines[ipT1527]);
00265
00266
00267 PutLine(&TauLines[ipT1305]);
00268
00269
00270 PutLine(&TauLines[ipT1260]);
00271
00272
00273 drec = dense.xIonDense[13][2]*dense.eden*7.6e-7/phycon.te32*
00274 1.57e-11*ionbal.DielSupprs[1][2];
00275
00276 fac = emit_frac(&TauLines[ipT1260]);
00277 PntForLine(1260.,"Si 2",&ipnt);
00278 lindst(drec*fac,1260,"diel",ipnt,'i',true );
00279
00280
00281 PntForLine(1909.,"Si 2",&ipnt);
00282 lindst(drec*1260./1909.,1909,"diel",ipnt,'i',true );
00283
00284 rec = 1e-12*(0.1152/t4 - 0.3082 + 4.4734*t4 + 0.0207*t4*t4)/pow(t4,1.5)*
00285 sexp(0.2981/t4);
00286
00287
00288
00289 rec *= dense.xIonDense[13][3]*dense.eden*emit_frac(&TauLines[ipT1207])*
00290 1.65e-11*ionbal.DielSupprs[1][2];
00291
00292 PutExtra(MAX2(0.,rec));
00293
00294
00295 PutLine(&TauLines[ipT1207]);
00296
00297
00298 linadd(MAX2(0.,rec),1207,"rec ",'i' );
00299
00300
00301 linadd(embesq.em1895+TauLines[ipT1895].xIntensity,1888,"TOTL",'i' );
00302 PutLine(&TauLines[ipT1895]);
00303
00304
00305 PntForLine(1883.,"Si 3",&ipnt);
00306 lindst(embesq.em1895,1883,"Si 3",ipnt,'i',true );
00307
00308
00309 fac = emit_frac(&TauLines[ipT1895]);
00310
00311
00312 linadd(sil.p1895*1.05e-11*fac,1895,"PHOT",'i' );
00313
00314
00315 linadd(TauLines[ipT1403].xIntensity+TauLines[ipT1394].xIntensity,1397,"TOTL",'i' );
00316
00317 sum = TauLines[ipT1403].xIntensity*TauLines[ipT1403].FracInwd +
00318 TauLines[ipT1394].xIntensity* TauLines[ipT1394].FracInwd;
00319
00320
00321 linadd(sum,1397,"Inwd",'i' );
00322 PutLine(&TauLines[ipT1403]);
00323 PutLine(&TauLines[ipT1394]);
00324
00325
00326 PutLine(&TauLines[ipSi619]);
00327
00328
00329 PntForLine(2148.,"Si 7",&ipnt);
00330 lindst(sil.c2148,2148,"Si 7",ipnt,'c',true );
00331
00332
00333 PutLine(&TauLines[ipTSi25]);
00334
00335
00336 PutLine(&TauLines[ipTSi65]);
00337
00338
00339 PntForLine(1446.,"Si 8",&ipnt);
00340 lindst(sil.c1446,1446,"Si 8",ipnt,'c',true );
00341
00342
00343 PntForLine(1985.,"Si 9",&ipnt);
00344 lindst(sil.c1985,1985,"Si 9",ipnt,'c',true );
00345
00346
00347 PntForLine(949.,"Si 9",&ipnt);
00348 lindst(sil.c949,949,"Si 9",ipnt,'c',true );
00349
00350
00351 PntForLine(1815.,"Si 9",&ipnt);
00352 lindst(sil.c1815,1815,"Si 9",ipnt,'c',true );
00353
00354
00355 PutLine(&TauLines[ipTSi4]);
00356
00357
00358 PutLine(&TauLines[ipTSi3]);
00359
00360
00361 PntForLine(691.,"Si 9",&ipnt);
00362 lindst(sil.c691,691,"Si 9",ipnt,'c',true );
00363
00364
00365
00366
00367 PutLine(&TauLines[ipSi10_606]);
00368
00369
00370 PutLine(&TauLines[ipSi10143]);
00371
00372
00373
00374 PntForLine(583.,"Si11",&ipnt);
00375 lindst(sil.c583,581,"Si11",ipnt,'c',true );
00376
00377
00378 linadd(TauLines[ipTSi499].xIntensity+TauLines[ipTSi521].xIntensity,506,"TOTL",'i' );
00379
00380
00381 PutLine(&TauLines[ipTSi499]);
00382 PutLine(&TauLines[ipTSi521]);
00383
00384
00385 PutLine(&TauLines[ipTSi41]);
00386
00387
00388
00389 PutLine(&TauLines[ipP0260]);
00390
00391
00392 PutLine(&TauLines[ipP0233]);
00393
00394
00395 PntForLine(16400.,"P 2",&ipnt);
00396 lindst(CoolHeavy.p2_32,16400,"P 2",ipnt,'c',true );
00397
00398
00399
00400 PntForLine(4669.,"P 2",&ipnt);
00401 lindst(CoolHeavy.p2_31*0.75,4669,"P 2",ipnt,'c',true );
00402
00403
00404 PntForLine(4737.,"P 2",&ipnt);
00405 lindst(CoolHeavy.p2_31*0.25,4737,"P 2",ipnt,'c',true );
00406
00407
00408 PntForLine(11890.,"P 2",&ipnt);
00409 lindst(CoolHeavy.p2_21*0.75,11890,"P 2",ipnt,'c',true );
00410
00411
00412 PntForLine(11470.,"P 2",&ipnt);
00413 lindst(CoolHeavy.p2_21*0.25,11470,"P 2",ipnt,'c',true );
00414
00415
00416 PutLine(&TauLines[ipP0318]);
00417
00418
00419 PutLine(&TauLines[ipP713]);
00420
00421
00422 PutLine(&TauLines[ipP848]);
00423
00424
00425 PutLine(&TauLines[ipP817]);
00426
00427
00428 PutLine(&TauLines[ipP1027]);
00429
00430
00431 PutLine(&TauLines[ipP1018]);
00432
00433
00434
00435 PutLine(&TauLines[ipS1_25m]);
00436
00437
00438 PutLine(&TauLines[ipS1_56m]);
00439
00440
00441
00442
00443 linadd(1.75e-22*dense.eden*dense.xIonDense[ipSULPHUR][1]/phycon.sqrte/
00444 phycon.te10/phycon.te03,1807,"S 1R",'i' );
00445
00446
00447 linadd(CoolHeavy.c6731,6720,"S 2",'c' );
00448
00449
00450 linadd(CoolHeavy.S4070+CoolHeavy.S4078,4074,"S 2",'c' );
00451
00452
00453 PntForLine(10330.,"S 2",&ipnt);
00454 lindst(CoolHeavy.c10330,10330,"S 2",ipnt,'c',true );
00455
00456
00457
00458 PntForLine(6731.,"S II",&ipnt);
00459 lindst(CoolHeavy.S6733,6731,"S II",ipnt,'i',true );
00460
00461
00462 PntForLine(6716.,"S II",&ipnt);
00463 lindst(CoolHeavy.S6718,6716,"S II",ipnt,'i',true );
00464
00465
00466 PntForLine(4070.,"S II",&ipnt);
00467 lindst(CoolHeavy.S4070,4070,"S II",ipnt,'i',true );
00468
00469
00470 PntForLine(4078.,"S II",&ipnt);
00471 lindst(CoolHeavy.S4078,4078,"S II",ipnt,'i',true );
00472
00473
00474 linadd(CoolHeavy.S10323,10323,"S II",'i' );
00475
00476
00477 linadd(CoolHeavy.S10289,10289,"S II",'i' );
00478
00479
00480 linadd(CoolHeavy.S10373,10373,"S II",'i' );
00481
00482
00483 linadd(CoolHeavy.S10339,10339,"S II",'i' );
00484
00485
00486 PutLine(&TauLines[ipT1256]);
00487
00488
00489 PutLine(&TauLines[ipTS19]);
00490
00491
00492 PutLine(&TauLines[ipTS34]);
00493
00494
00495 PutLine(&TauLines[ipTS1720]);
00496
00497
00498 PntForLine(9532.,"S 3",&ipnt);
00499 lindst(CoolHeavy.c9532/(1.+1./2.48),9532,"S 3",ipnt,'c',true );
00500
00501
00502 PntForLine(9069.,"S 3",&ipnt);
00503 lindst(CoolHeavy.c9532/(1.+2.48),9069,"S 3",ipnt,'c',true );
00504
00505
00506 PntForLine(6312.,"S 3",&ipnt);
00507 lindst(CoolHeavy.c6312,6312,"S 3",ipnt,'c',true );
00508
00509
00510 PntForLine(3722.,"S 3",&ipnt);
00511 lindst(CoolHeavy.c6312*0.59,3722,"S 3",ipnt,'c',true );
00512
00513
00514 PutLine(&TauLines[ipT1194]);
00515
00516
00517 PutLine(&TauLines[ipTS11]);
00518
00519
00520 linadd(
00521 TauLines[ipS4_1405].xIntensity+
00522 TauLines[ipS4_1398].xIntensity+
00523 TauLines[ipS4_1424].xIntensity+
00524 TauLines[ipS4_1417].xIntensity+
00525 TauLines[ipS4_1407].xIntensity,
00526 1406,"TOTL",'i' );
00527 PutLine(&TauLines[ipS4_1405]);
00528 PutLine(&TauLines[ipS4_1398]);
00529 PutLine(&TauLines[ipS4_1424]);
00530 PutLine(&TauLines[ipS4_1417]);
00531 PutLine(&TauLines[ipS4_1407]);
00532
00533
00534 linadd(embesq.em1198+TauLines[ipT1198].xIntensity,1198,"TOTL",'i' );
00535
00536
00537 PutLine(&TauLines[ipT1198]);
00538
00539
00540 linadd(embesq.em1198,1188,"S 5",'i' );
00541
00542
00543 PutLine(&TauLines[ipT786]);
00544
00545
00546
00547 linadd(TauLines[ipT944].xIntensity+TauLines[ipT933].xIntensity,933,"TOTL",'i' );
00548 PutLine(&TauLines[ipT944]);
00549 PutLine(&TauLines[ipT933]);
00550
00551
00552 PutLine(&TauLines[ipfsS810]);
00553
00554
00555 PutLine(&TauLines[ipfsS912]);
00556
00557
00558 PutLine(&TauLines[ipfsS938]);
00559
00560
00561 PutLine(&TauLines[ipfsS1119]);
00562
00563
00564 PutLine(&TauLines[ipfsS1114]);
00565
00566
00567 PntForLine(1715.,"S 9",&ipnt);
00568 lindst(CoolHeavy.c1715,1715,"S 9",ipnt,'c',true );
00569
00570
00571 PntForLine(1213.,"S 10",&ipnt);
00572 lindst(CoolHeavy.c1213,1213,"S 10",ipnt,'c',true );
00573
00574
00575 PntForLine(1826.,"S 11",&ipnt);
00576 lindst(CoolHeavy.c1826,1826,"S 11",ipnt,'c',true );
00577
00578
00579 PutLine(&TauLines[ipfsS1207]);
00580
00581
00582
00583
00584
00585 PutLine(&TauLines[ipS12_520]);
00586
00587
00588 PntForLine(488.,"S 13",&ipnt);
00589 lindst(CoolHeavy.c488,488,"S 13",ipnt,'c',true );
00590
00591
00592 linadd(TauLines[ipTSu418].xIntensity+TauLines[ipTSu446].xIntensity,427,"TOTL",'i' );
00593 PutLine(&TauLines[ipTSu418]);
00594 PutLine(&TauLines[ipTSu446]);
00595
00596
00597 PutLine(&TauLines[ipTSu30]);
00598
00599
00600
00601
00602 PutLine(&TauLines[ipCl1_11m]);
00603
00604
00605 PutLine(&TauLines[ipfsCl233]);
00606
00607
00608 PutLine(&TauLines[ipfsCl214]);
00609
00610
00611 PntForLine(8578.7,"Cl 2",&ipnt);
00612 lindst(CoolHeavy.c8579*0.791,8579,"Cl 2",ipnt,'c',true );
00613
00614
00615 PntForLine(9123.6,"Cl 2",&ipnt);
00616 lindst(CoolHeavy.c8579*0.209,9124,"Cl 2",ipnt,'c',true );
00617
00618
00619
00620
00621 PntForLine(6161.8,"Cl 2",&ipnt);
00622 lindst(CoolHeavy.c6164,6162,"Cl 2",ipnt,'c',true );
00623
00624
00625
00626 PntForLine(3677.9,"Cl 2",&ipnt);
00627 lindst(CoolHeavy.c3679,3678,"Cl 2",ipnt,'c',true );
00628
00629
00630 linadd(CoolHeavy.c5525,5525,"TOTL",'c' );
00631
00632
00633 linadd(CoolHeavy.c3350,3350,"TOTL",'c' );
00634
00635
00636 linadd(CoolHeavy.c8494,8494,"TOTL",'c' );
00637
00638
00639 PntForLine(5538.,"Cl 3",&ipnt);
00640 lindst(CoolHeavy.Cl5539,5538,"Cl 3",ipnt,'i',true );
00641
00642
00643 PntForLine(5518.,"Cl 3",&ipnt);
00644 lindst(CoolHeavy.Cl5519,5518,"Cl 3",ipnt,'i',true );
00645
00646
00647 PntForLine(3354.,"Cl 3",&ipnt);
00648 lindst(CoolHeavy.Cl3354,3354,"Cl 3",ipnt,'i',true );
00649
00650
00651 PntForLine(3344.,"Cl 3",&ipnt);
00652 lindst(CoolHeavy.Cl3344,3344,"Cl 3",ipnt,'i',true );
00653
00654
00655 PntForLine(8504.,"Cl 3",&ipnt);
00656 lindst(CoolHeavy.Cl8504,8504,"Cl 3",ipnt,'i',true );
00657
00658
00659 PntForLine(8436.,"Cl 3",&ipnt);
00660 lindst(CoolHeavy.Cl8436,8436,"Cl 3",ipnt,'i',true );
00661
00662
00663 PntForLine(8552.,"Cl 3",&ipnt);
00664 lindst(CoolHeavy.Cl8552,8552,"Cl 3",ipnt,'i',true );
00665
00666
00667 PntForLine(8483.,"Cl 3",&ipnt);
00668 lindst(CoolHeavy.Cl8483,8483,"Cl 3",ipnt,'i',true );
00669
00670
00671 PutLine(&TauLines[ipCl04203]);
00672
00673
00674 PutLine(&TauLines[ipCl04117]);
00675
00676
00677 PntForLine(8047.,"Cl 4",&ipnt);
00678 lindst(CoolHeavy.c8047*0.667,8047,"Cl 4",ipnt,'c',true );
00679
00680
00681 PntForLine(7532.,"Cl 4",&ipnt);
00682 lindst(CoolHeavy.c8047*0.333,7532,"Cl 4",ipnt,'c',true );
00683
00684
00685 PntForLine(3119.,"Cl 4",&ipnt);
00686 lindst(CoolHeavy.c3119,3119,"Cl 4",ipnt,'c',true );
00687
00688
00689 PntForLine(5324.,"Cl 4",&ipnt);
00690 lindst(CoolHeavy.c5324,5324,"Cl 4",ipnt,'c',true );
00691
00692
00693 PutLine(&TauLines[ipCl973]);
00694
00695
00696 PutLine(&TauLines[ipCl1030]);
00697
00698
00699 PutLine(&TauLines[ipCl1092]);
00700
00701
00702 PutLine(&TauLines[ipTAr7]);
00703
00704
00705 PntForLine(7135.,"Ar 3",&ipnt);
00706 lindst(CoolHeavy.c7136/(1.+1./4.144),7135,"Ar 3",ipnt,'c',true );
00707
00708
00709 PntForLine(7751.,"Ar 3",&ipnt);
00710 lindst(CoolHeavy.c7136/(1.+4.144),7751,"Ar 3",ipnt,'c',true );
00711
00712
00713 PntForLine(5192.,"Ar 3",&ipnt);
00714 lindst(CoolHeavy.c5192,5192,"Ar 3",ipnt,'c',true );
00715
00716
00717 PntForLine(3109.,"Ar 3",&ipnt);
00718 lindst(CoolHeavy.c3109*0.9894,3109,"Ar 3",ipnt,'c',true );
00719
00720
00721 PntForLine(3005.,"Ar 3",&ipnt);
00722 lindst(CoolHeavy.c3109*(1.-0.9894),3005,"Ar 3",ipnt,'c',true );
00723
00724
00725 PutLine(&TauLines[ipTAr22]);
00726
00727
00728 PutLine(&TauLines[ipTAr9]);
00729
00730
00731 linadd(CoolHeavy.Ar4740+CoolHeavy.Ar4711,4725,"TOTL",'i' );
00732
00733
00734 linadd(CoolHeavy.Ar2868+CoolHeavy.Ar2854,2860,"TOTL",'i' );
00735
00736
00737 linadd(CoolHeavy.Ar7237+CoolHeavy.Ar7331+CoolHeavy.Ar7171+CoolHeavy.Ar7263,
00738 7250,"TOTL",'i' );
00739
00740
00741 PntForLine(4740.,"Ar 4",&ipnt);
00742 lindst(CoolHeavy.Ar4740,4740,"Ar 4",ipnt,'c',true );
00743
00744
00745 PntForLine(4711.,"Ar 4",&ipnt);
00746 lindst(CoolHeavy.Ar4711,4711,"Ar 4",ipnt,'c',true );
00747
00748
00749 PntForLine(2868.,"Ar 4",&ipnt);
00750 lindst(CoolHeavy.Ar2868,2868,"Ar 4",ipnt,'c',true );
00751
00752
00753 PntForLine(2854.,"Ar 4",&ipnt);
00754 lindst(CoolHeavy.Ar2854,2854,"Ar 4",ipnt,'c',true );
00755
00756
00757 PntForLine(7263.,"Ar 4",&ipnt);
00758 lindst(CoolHeavy.Ar7263,7263,"Ar 4",ipnt,'c',true );
00759
00760
00761 PntForLine(7171.,"Ar 4",&ipnt);
00762 lindst(CoolHeavy.Ar7171,7171,"Ar 4",ipnt,'c',true );
00763
00764
00765 PntForLine(7331.,"Ar 4",&ipnt);
00766 lindst(CoolHeavy.Ar7331,7331,"Ar 4",ipnt,'c',true );
00767
00768
00769 PntForLine(7237.,"Ar 4",&ipnt);
00770 lindst(CoolHeavy.Ar7237,7237,"Ar 4",ipnt,'c',true );
00771
00772
00773 PntForLine(7005.,"Ar 5",&ipnt);
00774 lindst(CoolHeavy.c7007/(1.+1./2.143),7005,"Ar 5",ipnt,'c',true );
00775
00776
00777 PntForLine(6435.,"Ar 5",&ipnt);
00778 lindst(CoolHeavy.c7007/(1.+2.143),6435,"Ar 5",ipnt,'c',true );
00779
00780
00781
00782 PntForLine(4626.,"Ar 5",&ipnt);
00783 lindst( CoolHeavy.c4626 , 4626 ,"Ar 5",ipnt,'c',true );
00784
00785
00786 PntForLine(2691.,"Ar 5",&ipnt);
00787 lindst( CoolHeavy.c2691 , 2691 ,"Ar 5",ipnt,'c',true );
00788
00789
00790 PutLine(&TauLines[ipTAr13]);
00791
00792
00793 PutLine(&TauLines[ipTAr8]);
00794
00795
00796 PutLine(&TauLines[ipAr06453]);
00797
00798
00799 PutLine(&TauLines[ipAr1055]);
00800
00801
00802 PutLine(&TauLines[ipAr1126]);
00803
00804
00805 PutLine(&TauLines[ipAr1178]);
00806
00807
00808 PntForLine(4413.,"Ar14",&ipnt);
00809 lindst(CoolHeavy.fs4413,4413,"Ar14",ipnt,'c',true );
00810
00811
00812 PntForLine(409.,"Ar15",&ipnt);
00813 lindst(CoolHeavy.c409,409,"Ar15",ipnt,'c',true );
00814
00815
00816 PutLine(&TauLines[ipT354]);
00817
00818
00819 PutLine(&TauLines[ipT389]);
00820
00821
00822 PutLine(&TauLines[ipT25]);
00823
00824 if( trace.lgTrace )
00825 {
00826 fprintf( ioQQQ, " lines_lv1_na_ar returns\n" );
00827 }
00828
00829 DEBUG_EXIT( "lines_lv1_na_ar()" );
00830 return;
00831 }
00832