00001
00002
00003
00004
00005
00006 #include "cddefines.h"
00007 #include "carb.h"
00008 #include "nitro.h"
00009 #include "oxy.h"
00010 #include "coolheavy.h"
00011 #include "atmdat.h"
00012 #include "doppvel.h"
00013 #include "ionbal.h"
00014 #include "dense.h"
00015 #include "phycon.h"
00016 #include "physconst.h"
00017 #include "atoms.h"
00018 #include "mole.h"
00019 #include "embesq.h"
00020 #include "taulines.h"
00021 #include "trace.h"
00022 #include "lines_service.h"
00023 #include "lines.h"
00024 STATIC double GetLineRec(
00025
00026 long int ip,
00027
00028 long int lWl);
00029
00030 void lines_lv1_li_ne(void)
00031 {
00032 long int ipnt;
00033 double
00034 chem ,
00035 corr,
00036 ct4363,
00037 ctRate,
00038 efac,
00039 effec,
00040 efficn2,
00041 fac,
00042 HBeta ,
00043 p386,
00044 pump,
00045 r4363,
00046 r6584,
00047 raten3,
00048 rb,
00049 rec,
00050 rn3mor,
00051 rn3tot,
00052 rnii,
00053 rp300,
00054 rp386,
00055 r12 ,
00056 r13 ,
00057 sum,
00058 rate_OH_dissoc;
00059 double rec7323 , rec7332, rec3730 , rec3726 , rec2471,
00060 reco23tot , reco22tot;
00061
00062 DEBUG_ENTRY( "lines_lv1_li_ne()" );
00063
00064 if( trace.lgTrace )
00065 {
00066 fprintf( ioQQQ, " lines_lv1_li_ne called\n" );
00067 }
00068
00069
00070 ipnt = StuffComment( "level 1 lines" );
00071 linadd( 0., (realnum)ipnt , "####", 'i',
00072 " start level 1 ines" );
00073
00074 linadd(CoolHeavy.colmet,0,"Mion",'c',
00075 " cooling due to collisional ionization of heavy elements" );
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087 rec = GetLineRec(3,1657)*(1.-TauLines[ipT1656].Emis().ColOvTot());
00088 rec = GetLineRec(3,1657)*emit_frac(TauLines[ipT1656]);
00089 PutExtra(rec);
00090
00091 PutLine(TauLines[ipT1656],
00092 " C 1 1656, collision strength from van Regemoter");
00093
00094 linadd(rec,1656,"REC ",'i',
00095 " C 1 1656 recomb; n.b. coll deexcitation not in" );
00096
00097 PntForLine(9850.,"C 1",&ipnt);
00098 lindst(carb.c9850,9850,"C Ic",ipnt,'c',true,
00099 " C 1 9850, coll excit" );
00100
00101
00102
00103
00104
00105
00106
00107 rec = (GetLineRec(1,9088) + GetLineRec(2,9658))*carb.r9850;
00108 lindst(rec,9850,"C Ir",ipnt,'r',true,
00109 " C I 9850 recombination contribution" );
00110
00111
00112
00113 linadd(rec+carb.c9850,9850,"TOTL",'i',
00114 " total intensity, all processes, C I 9850");
00115
00116
00117 PntForLine(8727.,"C 1",&ipnt);
00118 lindst(carb.c8727,8727,"C 1",ipnt,'c',true,
00119 "C 1 8727; equivalent to 4363" );
00120
00121 linadd(carb.c8727*1.22e-6,4621,"C 1",'c',
00122 " 1S - 3P" );
00123
00124 PutLine(TauLines[ipT610],
00125 " C 1 610 micron ");
00126
00127 PutLine(TauLines[ipT370],
00128 " C 1 370 micron ");
00129
00130 PutLine(TauLines[ipT157],
00131 " C 2 158 micron, both e- and H0, H2 in excitation ");
00132
00133 linadd(
00134 TauLines[ipC2_2325].Emis().xIntensity()+
00135 TauLines[ipC2_2324].Emis().xIntensity()+
00136 TauLines[ipC2_2329].Emis().xIntensity()+
00137 TauLines[ipC2_2328].Emis().xIntensity()+
00138 TauLines[ipC2_2327].Emis().xIntensity()+
00139 ionbal.PhotoRate_Shell[ipCARBON][0][1][0]*dense.xIonDense[ipCARBON][0]*0.1*8.6e-12,
00140 2326,"TOTL",'i',
00141 " total intensity of C II] 2326, all lines in the multiplet " );
00142 PutLine(TauLines[ipC2_2325],
00143 " ");
00144 PutLine(TauLines[ipC2_2324],
00145 " ");
00146 PutLine(TauLines[ipC2_2329],
00147 " ");
00148 PutLine(TauLines[ipC2_2328],
00149 " ");
00150 PutLine(TauLines[ipC2_2327],
00151 " ");
00152
00153 linadd(ionbal.PhotoRate_Shell[ipCARBON][0][1][0]*dense.xIonDense[ipCARBON][0]*0.1*8.6e-12,2326,"Phot",'i' ,
00154 " photoproduction, Helfand and Trefftz");
00155
00156
00157
00158
00159 rec = GetLineRec(11,1335)*emit_frac(TauLines[ipT1335]);
00160
00161
00162 PutExtra(MAX2(0.,rec));
00163 PutLine(TauLines[ipT1335],
00164 " total intensity of C 2 1335");
00165
00166 linadd(rec,1335,"REC ",'i',
00167 " C 2 1335 recombination," );
00168
00169
00170
00171
00172
00173 if( nWindLine > 0 )
00174 {
00175 pump = TauLine2[186].Emis().pump()*TauLine2[186].Emis().PopOpc();
00176 }
00177 else
00178 {
00179 pump = 0.;
00180 }
00181
00182
00183
00184 PntForLine(3920.,"C 2",&ipnt);
00185 lindst(pump*0.387 * 5.08e-12/(1.+dense.eden/1e12) ,3920,"pump",ipnt,'r',true ,
00186 " CII 3918.98/3920.68 is only pumped, no recombination part");
00187
00188
00189 rec = GetLineRec(8, 6580 );
00190
00191 pump *= 0.305 * 0.387 * 3.02e-12;
00192 linadd(rec/(1.+dense.eden/1e12),6580,"C 2r",'i',
00193 " recombination part of C II 6580 line " );
00194 linadd(pump/(1.+dense.eden/1e12),6580,"C 2p",'i',
00195 " pumped part of line C II 6580" );
00196
00197
00198 PntForLine(6580.,"C 2",&ipnt);
00199 lindst((rec+pump)/(1.+dense.eden/1e12),6580,"TOTL",ipnt,'r',true ,
00200 " total intensity, all processes, C II 6580");
00201
00202
00203
00204
00205
00206 rec = GetLineRec(179,977)*emit_frac(TauLines[ipT977]);
00207
00208
00209 rp386 = TauLines[ipT386].Emis().pump()*TauLines[ipT386].Emis().PopOpc();
00210
00211
00212 rp300 = TauLines[ipT310].Emis().pump()*TauLines[ipT310].Emis().PopOpc() +
00213 TauLines[ipT291].Emis().pump()*TauLines[ipT291].Emis().PopOpc() +
00214 TauLines[ipT280].Emis().pump()*TauLines[ipT280].Emis().PopOpc() +
00215 TauLines[ipT274].Emis().pump()*TauLines[ipT274].Emis().PopOpc() +
00216 TauLines[ipT270].Emis().pump()*TauLines[ipT270].Emis().PopOpc();
00217
00218
00219
00220
00221 p386 = (rp386 + rp300)*2.03e-11*emit_frac(TauLines[ipT977]);
00222
00223
00224 PutExtra(p386+MAX2(0.,rec));
00225
00226 PutLine(TauLines[ipT977],
00227 " total C 3] 977, recombination + collisional + pumped excitation ");
00228
00229 linadd(rec,977,"C3 R",'i',
00230 " dielectronic recombination contribution to C 3 977 " );
00231
00232 linadd(p386,977,"P386",'r',
00233 " C 3 977 pumped by continuum near 386A" );
00234
00235
00236 fac = embesq.em1908 + TauLines[ipT1909].Emis().xIntensity();
00237 lindst(fac,1909,"TOTL",TauLines[ipT1909].ipCont(),'i',false,
00238 "C 3 1909 collision, both lines together");
00239
00240
00241 PutLine(TauLines[ipT1909],
00242 "C III 19091");
00243
00244
00245 PntForLine(1907.,"C 3",&ipnt);
00246 lindst(embesq.em1908,1907,"C 3",ipnt,'c',true,
00247 " C 3 1908 j-2 to ground" );
00248
00249 lindst(embesq.em13C1910,1910,"13C3",ipnt,'c',true,
00250 " the 13C forbidden line of C III " );
00251
00252
00253
00254 corr = emit_frac(TauLines[ipT1909]);
00255 fac = dense.eden*dense.xIonDense[ipCARBON][3]/(phycon.te/phycon.te10);
00256
00257 linadd(3.1e-19*fac*corr,1909,"C3 R",'i',
00258 " C 3 1909 recombination from Storey" );
00259
00260 linadd(ionbal.PhotoRate_Shell[ipCARBON][1][1][0]*dense.xIonDense[ipCARBON][1]*0.62*corr*1.05e-11,1909,"Phot",'i',
00261 " C 3 1909 following relax following inner shell photoionization" );
00262
00263
00264
00265
00266 rec = GetLineRec(178,1176)*emit_frac(TauLines[ipc31175]);
00267 PutExtra(MAX2(0.,rec));
00268
00269 PutLine(TauLines[ipc31175],
00270 " C 3* 1175, excited state line, above 1909 ");
00271
00272 linadd(MAX2(0.,rec),1175,"Rec ",'i',
00273 " dielectronic recombination contribution to C 3 1175 " );
00274
00275
00276
00277
00278
00279 rec = GetLineRec(25,1549)*emit_frac(TauLines[ipT1550]);
00280
00281 linadd(
00282 TauLines[ipT1550].Emis().xIntensity()+
00283 TauLines[ipT1548].Emis().xIntensity()+
00284 rec,1549,"TOTL",'i',"total intensity of C 4 1549, all processes " );
00285
00286 sum =
00287 TauLines[ipT1550].Emis().xIntensity()*TauLines[ipT1550].Emis().FracInwd() +
00288 TauLines[ipT1548].Emis().xIntensity()*TauLines[ipT1548].Emis().FracInwd();
00289
00290 linadd(sum+rec*TauLines[ipT1550].Emis().FracInwd(),1549,"Inwd",'i',
00291 "inward part of C 4 " );
00292
00293 PutExtra(rec*.3333);
00294 PutLine(TauLines[ipT1550],
00295 " ");
00296
00297 PutExtra(rec*.6666);
00298 PutLine(TauLines[ipT1548],
00299 " ");
00300
00301 linadd((TauLines[ipT1550].Emis().ots()+
00302 TauLines[ipT1548].Emis().ots())*TauLines[ipT1548].EnergyErg(),
00303 1549,"DEST",'i',
00304 " part of line destroyed by photoionization of Balmer continuum " );
00305
00306 linadd(rec,1549,"C4 r",'i',
00307 " recombination C 4 1549 from CV" );
00308
00309 PutLine(TauLines[ipT312],
00310 " Li seq 2s 3p Li seq transition");
00311
00312
00313
00314
00315
00316
00317 PutLine(TauLines[ipNI_pumpIndirect],"FUV indirect excitation");
00318 for( int i=0; i < NI_NDP; ++i )
00319 PutLine(TauLines[ipNI_pumpDirect[i]],"FUV direct excitation");
00320
00321 PutLine(TauLines[ipT1200]," N I 1200, full multiplet, all processes");
00322
00323 PntForLine(3466.52,"N 1",&ipnt);
00324 lindst( nitro.xN3466, 3466, "N 1", ipnt, 'c', true, " [N I] 3466.50 only" );
00325
00326 lindst( nitro.xN3467, 3468, "N 1", ipnt, 'c', true, " [N I] 3466.54 only" );
00327
00328
00329 lindst( nitro.xN3467+nitro.xN3466, 3467, "TOTL", ipnt, 'i', false,
00330 " [N I] 3466.54, 3466.50 together" );
00331
00332 PntForLine(5199.,"N 1",&ipnt);
00333
00334 lindst( nitro.rec5199, 5199, "TOTr", ipnt, 'r', true,
00335 " estimate of contribution to [N I] 5199 by recombination" );
00336
00337
00338
00339
00340
00341 chem = mole.dissoc_rate("N") * 3.82e-12 * nitro.quench_5200;
00342 lindst( chem, 5199, "chem", ipnt, 'r', true,
00343 " upper limit to [N I] 5199 produced by chemistry" );
00344
00345
00346
00347 ctRate = atmdat.HCharExcRecTo_N0_2D*dense.xIonDense[ipHYDROGEN][0]*dense.xIonDense[ipNITROGEN][1] *
00348 3.82e-12 * nitro.quench_5200;
00349 lindst( ctRate, 5199, "H CT", ipnt, 'r', true,
00350 " upper limit to [N I] 5199 produced by charge transfer" );
00351
00352
00353 lindst( nitro.pump5199, 5199, "pump", ipnt, 'r', true,
00354 " estimate of contribution to [N I] 5199 by FUV pumping" );
00355
00356
00357 lindst( nitro.xN5200 + (ctRate+chem)*0.6, 5200, "N 1", ipnt, 'c', true,
00358 " [N I] 5200 - all processes - stat weight is 6 - total in term is 10" );
00359
00360 lindst( nitro.xN5198 + (ctRate+chem)*0.4, 5198, "N 1", ipnt, 'c', true,
00361 " [N I] 5198 - all processes - stat weight is 4 - total in term is 10" );
00362
00363
00364 lindst( nitro.xN5200+nitro.xN5198 + (ctRate+chem), 5199, "TOTL", ipnt, 'i', false,
00365 " [N I] 5200 + 5198 together" );
00366
00367 PntForLine(10403.,"N 1",&ipnt);
00368 lindst( nitro.xN10397, 10397, "N 1", ipnt, 'c', true, " [N I] 10397.7 only");
00369 lindst( nitro.xN10398, 10398, "N 1", ipnt, 'c', true, " [N I] 10398.2 only");
00370 lindst( nitro.xN10407, 10407, "N 1", ipnt, 'c', true, " [N I] 10407.2 only");
00371 lindst( nitro.xN10408, 10408, "N 1", ipnt, 'c', true, " [N I] 10407.6 only");
00372
00373
00374 lindst( nitro.xN10398+nitro.xN10397+nitro.xN10408+nitro.xN10407, 10403, "TOTL", ipnt, 'i', false,
00375 " [N I] 10398.2, 10397.7, 10407.6, 10407.2 together");
00376
00377
00378
00379 PutLine(TauLines[ipT671] ," 6 lines with fake collision strengths" );
00380 PutLine(TauLines[ipT315] ," 6 lines with fake collision strengths" );
00381 PutLine(TauLines[ipT333] ," 6 lines with fake collision strengths" );
00382 PutLine(TauLines[ipT324] ," 6 lines with fake collision strengths" );
00383 PutLine(TauLines[ipT374g]," 6 lines with fake collision strengths" );
00384 PutLine(TauLines[ipT374x]," 6 lines with fake collision strengths" );
00385
00386 PntForLine(6584.,"N 2",&ipnt);
00387 lindst(nitro.c6584/(1.+1./2.951),6584,"N 2",ipnt,'c',true ,
00388 " N 2 6584 alone " );
00389
00390 PntForLine(6548.,"N 2",&ipnt);
00391 lindst(nitro.c6584/(1.+2.951),6548,"N 2",ipnt,'c',true,
00392 " N 2 6548 alone " );
00393
00394 efficn2 = 4e-3/(4e-3 + 5.18e-6*dense.eden/phycon.sqrte);
00395 r6584 = 8e-22/(phycon.te70/phycon.te03/phycon.te03)*efficn2;
00396 linadd(r6584*dense.xIonDense[ipNITROGEN][2]*dense.eden,6584,"REC ", 'i',
00397 " N 2 6584 alone, recombination contribution" );
00398
00399
00400
00401 ctRate = 1.8e-11*dense.xIonDense[ipHELIUM][0]*dense.xIonDense[ipNITROGEN][2]*1.146/(1.146 +
00402 0.87*dense.cdsqte)*3.46e-12;
00403
00404
00405
00406
00407
00408 if( dense.xIonDense[ipHYDROGEN][1] > SMALLFLOAT )
00409 {
00410
00411
00412
00413
00414 HBeta = (pow(10.,-20.89 - 0.10612*POW2(phycon.alogte - 4.4)))/phycon.te;
00415
00416
00417
00418
00419 HBeta *= dense.eden * dense.xIonDense[ipHYDROGEN][1];
00420
00421
00422
00423
00424
00425 rec = nitro.xN2_A3_tot * HBeta *
00426 3.19 * phycon.te30 / 15.84893 * dense.xIonDense[ipNITROGEN][2]/dense.xIonDense[ipHYDROGEN][1];
00427 }
00428 else
00429 {
00430 HBeta = 0.;
00431 rec = 0.;
00432 }
00433
00434 linadd(nitro.c5755+ctRate+rec ,5755,"N 2",'i',
00435 " N 2 5755 total, collisions plus charge transfer plus recombination" );
00436
00437 PntForLine(5755.,"N 2",&ipnt);
00438 lindst(nitro.c5755,5755,"Coll",ipnt,'c',true,
00439 " N 2 5755 collisional contribution" );
00440
00441 lindst(ctRate,5755,"C T ",ipnt,'r',true,
00442 " N 2 5755 charge transfer contribution " );
00443
00444 lindst( rec ,5755,"N 2r",ipnt,'r',true,
00445 " N 2 5755 recombination contribution" );
00446
00447 PutLine(TauLines[ipT122],
00448 " N 2 fine structure line ");
00449
00450 PutLine(TauLines[ipT205],
00451 " N 2 fine structure line " );
00452
00453 PutLine(TauLines[ipT2140],
00454 " N 2 2140 intercombination line " );
00455
00456
00457
00458
00459 rec = GetLineRec(201,1085)*emit_frac(TauLines[ipT1085]);
00460 PutExtra(MAX2(0.,rec));
00461
00462 PutLine(TauLines[ipT1085],
00463 " N 2 1084, CS guess from g-bar " );
00464
00465 linadd(MAX2(0.,rec),1085,"Rec ",'i',
00466 " dielectronic recombination contribution to N 2 1085" );
00467
00468
00469 rnii = TauLines[ipT671].Emis().pump()*TauLines[ipT671].Emis().PopOpc();
00470
00471 linadd(rnii*0.377*0.75*3.02e-12*efficn2,6584,"N2cn",'i',
00472 " continuum pumped N 2 6584 " );
00473
00474 efficn2 = 1./(1. + COLL_CONST*dense.eden/phycon.sqrte);
00475 linadd(rnii*0.0117*3.46e-12*efficn2,5755,"N2cn",'i',
00476 " continuum pumped N 2 5755" );
00477
00478
00479
00480
00481
00482
00483
00484 if( nWindLine > 0 )
00485 {
00486 pump = TauLine2[265].Emis().pump()*TauLine2[265].Emis().PopOpc();
00487 }
00488 else
00489 {
00490 pump = 0.;
00491 }
00492
00493 PntForLine(3311.,"N 2",&ipnt);
00494 lindst(pump*0.236 * 6.01e-12/(1.+dense.eden/1e12) ,3311,"pump",ipnt,'r',true,
00495 " NII 3311.42 - 3331.31 (6 lines) are only pumped, no recombination part" );
00496
00497 PntForLine(3840.,"N 2",&ipnt);
00498 lindst(pump*0.186 * 5.18e-12/(1.+dense.eden/1e12) ,3840,"pump",ipnt,'r',true,
00499 " NII 3829.8-3856.06 (6 lines) are only pumped, no recombination part" );
00500
00501 PntForLine(3609.,"N 2",&ipnt);
00502 lindst(pump*0.025 * 5.52e-12/(1.+dense.eden/1e12) ,3609,"pump",ipnt,'r',true,
00503 " NII 3593.60/3609.1/3615.86 (3 lines) are only pumped, no recombination part" );
00504
00505
00506 PntForLine(4640.,"N 2",&ipnt);
00507 lindst(pump*0.186*0.595 * 4.31e-12/(1.+dense.eden/1e12) ,4640,"pump",ipnt,'r',true ,
00508 " NII 4601.5-4643.1 (6 lines) are only pumped, no recombination part");
00509
00510 PntForLine(5010.,"N 2",&ipnt);
00511 lindst(pump*0.025*0.442 * 3.97e-12/(1.+dense.eden/1e12) ,5010,"pump",ipnt,'r',true,
00512 " NII 5002.7/5010.6/5045.1 (3 lines) are only pumped, no recombination part " );
00513
00514
00515 rec = GetLineRec(44, 5679 );
00516
00517 pump *= 0.236 * 0.626 * 3.50e-12;
00518
00519 linadd(rec/(1.+dense.eden/1e12),5679,"N 2r",'i',
00520 " recombination part of N II 5679 line" );
00521
00522 linadd(pump/(1.+dense.eden/1e12),5679,"N 2p",'i',
00523 " pumped part of line N II 5679 " );
00524
00525 PntForLine(5679.,"N 2",&ipnt);
00526 lindst((rec+pump)/(1.+dense.eden/1e12),5679,"TOTL",ipnt,'r',true,
00527 " total intensity, all processes, N II 5679 " );
00528
00529 PutLine(TauLines[ipT57],
00530 " [N 3] 57 micron fine structure line");
00531
00532 linadd(
00533 TauLines[ipN3_1749].Emis().xIntensity()+
00534 TauLines[ipN3_1747].Emis().xIntensity()+
00535 TauLines[ipN3_1754].Emis().xIntensity()+
00536 TauLines[ipN3_1752].Emis().xIntensity()+
00537 TauLines[ipN3_1751].Emis().xIntensity(),
00538 1750,"TOTL",'i',
00539 " total intensity of N III] 1750, all lines in the multiplet " );
00540 PutLine(TauLines[ipN3_1749],
00541 " ");
00542 PutLine(TauLines[ipN3_1747],
00543 " ");
00544 PutLine(TauLines[ipN3_1754],
00545 " ");
00546 PutLine(TauLines[ipN3_1752],
00547 " ");
00548 PutLine(TauLines[ipN3_1751],
00549 " ");
00550
00551
00552
00553 raten3 = TauLines[ipT374x].Emis().PopOpc()*TauLines[ipT374x].Emis().pump();
00554
00555
00556 if( DoppVel.TurbVel < 200. )
00557 {
00558 rb = TauLines[ipT374x].Emis().PopOpc()*TauLines[ipT374x].Emis().pump() +
00559 TauLines[ipT374g].Emis().PopOpc()*TauLines[ipT374g].Emis().pump();
00560 }
00561 else
00562 {
00563
00564 rb = TauLines[ipT374g].Emis().PopOpc()*TauLines[ipT374g].Emis().pump();
00565 }
00566
00567 rn3mor =
00568 TauLines[ipT315].Emis().PopOpc()*TauLines[ipT315].Emis().pump()*0.448 +
00569 TauLines[ipT324].Emis().PopOpc()*TauLines[ipT324].Emis().pump()*0.78 +
00570 TauLines[ipT333].Emis().PopOpc()*TauLines[ipT333].Emis().pump()*0.434;
00571
00572
00573 rn3tot = (rb + raten3)*0.439 + rn3mor;
00574
00575 sum = raten3*4.29e-12;
00576 PntForLine(4640.,"N3cn",&ipnt);
00577 lindst(sum,4640,"N3cn",ipnt,'r',true ,
00578 " continuum pumped \"Bowen\" N 3, optically thin excited line ");
00579
00580
00581 sum = rb*4.29e-12*0.834;
00582 PntForLine(4634.,"N3cn",&ipnt);
00583 lindst(sum,4634,"N3cn",ipnt,'r',true ,
00584 " continuum pumped \"Bowen\" N 3, optically thin excited line");
00585
00586 sum = rb*4.29e-12*(1. - 0.834);
00587 PntForLine(4642.,"N3cn",&ipnt);
00588 lindst(sum,4642,"N3cn",ipnt,'r',true,
00589 " continuum pumped \"Bowen\" N 3, optically thin excited line" );
00590
00591
00592
00593
00594
00595 fac = 1.-emit_frac(TauLines[ipT990]);
00596
00597
00598
00599
00600 rec = GetLineRec(216,991)*emit_frac(TauLines[ipT990] );
00601 PutExtra(MAX2(0.,rec)+rn3tot*2.01e-11*fac);
00602
00603 PutLine(TauLines[ipT990],
00604 " N 3 990, all processes ");
00605
00606 linadd(rec+rn3tot*2.01e-11*fac,990,"extr",'i',
00607 " total N 3 990, both electron excitation and continuum pumping" );
00608
00609 linadd(rec,990,"rec ",'i',
00610 " part of N 3 990 due to recombination " );
00611
00612 linadd(rn3tot*2.01e-11,990,"N 3p",'r',
00613 " N 3 989.8, continuum pumped" );
00614
00615 linadd(embesq.em1486+TauLines[ipT1486].Emis().xIntensity(),1486,"TOTL",'i',
00616 " N 4] 1486, total intensity of both lines" );
00617
00618 PutLine(TauLines[ipT1486],
00619 " ");
00620
00621 PntForLine(1485.,"N 4",&ipnt);
00622 lindst(embesq.em1486,1485,"N 4",ipnt,'c',true ,
00623 " the N IV] slow transition by itself " );
00624
00625
00626
00627
00628 rec = GetLineRec(287,765)*emit_frac(TauLines[ipT765] );
00629
00630
00631 PutExtra(rec);
00632
00633 PutLine(TauLines[ipT765],
00634 " N 4 765, collisionally excited");
00635
00636 linadd(MAX2(0.,rec),765,"rec ",'i',
00637 " N 4 765 recombination," );
00638
00639 linadd(TauLines[ipT1243].Emis().xIntensity()+TauLines[ipT1239].Emis().xIntensity(),1240,"TOTL",'i',
00640 " continuum pumping of NV 1240, N 5 1240, total emission, collisions plus pumping " );
00641 sum = TauLines[ipT1243].Emis().xIntensity()*TauLines[ipT1243].Emis().FracInwd() + TauLines[ipT1239].Emis().xIntensity()*
00642 TauLines[ipT1239].Emis().FracInwd();
00643
00644 linadd(sum,1240,"Inwd",'i',
00645 " inward part of N 5 " );
00646 PutLine(TauLines[ipT1243],
00647 " ");
00648 PutLine(TauLines[ipT1239],
00649 " ");
00650
00651 PutLine(TauLines[ipT209],
00652 " N 5 209, 2s-3p Li seq ");
00653
00654 PntForLine(6300.,"O 1",&ipnt);
00655 lindst(CoolHeavy.c6300,6300,"O 1",ipnt,'c',true ,
00656 " oxygen total Oxygen I 6300, including continuum optical depth ");
00657
00658
00659
00660
00661
00662
00663
00664
00665
00666
00667
00668 rate_OH_dissoc = mole.findrate("PHOTON,OH=>O,H");
00669 r12 = rate_OH_dissoc * 0.55 * 3.16e-12 * CoolHeavy.c6300_frac_emit;
00670
00671 lindst( r12*TauLines[ipT6300].Emis().Aul()/(TauLines[ipT6300].Emis().Aul()+TauLines[ipT6363].Emis().Aul()) ,
00672 6300., "OH p",ipnt , 'i' , false,
00673 " the intensity of [OI] 6300 line due to OH photodistruction");
00674
00675 lindst( r12*TauLines[ipT6363].Emis().Aul()/(TauLines[ipT6300].Emis().Aul()+TauLines[ipT6363].Emis().Aul()) ,
00676 6363., "OH p",ipnt , 'i' , false,
00677 " the intensity of [OI] 6363 line due to OH photodistruction ");
00678
00679 PntForLine(6363.,"O 1",&ipnt);
00680 lindst(CoolHeavy.c6363,6363,"O 1",ipnt,'c',true,
00681 " total Oxygen I 6363, including continuum optical depth " );
00682
00683 PntForLine(5577.,"O 1",&ipnt);
00684 lindst(CoolHeavy.c5577,5577,"O 1",ipnt,'c',true,
00685 " auroral OI " );
00686
00687 r13 = rate_OH_dissoc * 0.05 * 3.57e-12 * 0.94*CoolHeavy.c5577_frac_emit;
00688 lindst( r13 , 5577., "OH p",ipnt , 'i', false,
00689 " 94% of excitations to highest level decay via 5577" );
00690
00691 PutLine(TauLines[ipT63],
00692 " O I fine structure line ");
00693
00694 PutLine(TauLines[ipT146],
00695 " O I fine structure line ");
00696
00697 linadd(MAX2(0.,CoolHeavy.coolOi),0,"TOIc",'c',
00698 " total collisional cooling due to 6-level OI atom" );
00699
00700 linadd(MAX2(0.,-CoolHeavy.coolOi),0,"TOIh",'h',
00701 " total collisional heating due to 6-level OI atom " );
00702
00703
00704
00705
00706 sum = atoms.popoi[2]*TauLines[ipT8446].Emis().Pesc_total()*TauLines[ipT8446].Emis().Aul()*2.36e-12;
00707 PntForLine(8446.,"O 1",&ipnt);
00708
00709 lindst(sum,8446,"6lev",ipnt,'i',false,
00710 " OI 8446 from six level atom" );
00711
00712 PntForLine(1304.,"O 1",&ipnt);
00713 sum = atoms.popoi[1]*TauLines[ipT1304].Emis().Pesc_total()*TauLines[ipT1304].Emis().Aul()*1.53e-11;
00714 lindst(sum,1304,"6lev",ipnt,'i',false,
00715 " OI 1304 from six level atom " );
00716
00717 PntForLine(1039.,"O 1",&ipnt);
00718 sum = atoms.popoi[3]*TauLines[ipT1039].Emis().Pesc_total()*TauLines[ipT1039].Emis().Aul()*1.92e-11;
00719 lindst(sum,1039,"6lev",ipnt,'i',false ,
00720 " OI 1039 from six level atom");
00721
00722 PntForLine(4368.,"O 1",&ipnt);
00723 sum = atoms.popoi[5]*TauLines[ipT4368].Emis().Pesc_total()*TauLines[ipT4368].Emis().Aul()*4.55e-12;
00724 lindst(sum,4368,"6lev",ipnt,'i',false,
00725 " OI 4368 from six level atom" );
00726
00727 PntForLine(13100.,"O 1",&ipnt);
00728 sum = atoms.popoi[3]*TauLines[ipTOI13].Emis().Pesc_total()*TauLines[ipTOI13].Emis().Aul()*1.52e-12;
00729 lindst(sum,13100,"6lev",ipnt,'i',false ,
00730 " OI 1.3 micron from six level atom");
00731
00732 PntForLine(11300.,"O 1",&ipnt);
00733 sum = atoms.popoi[4]*TauLines[ipTOI11].Emis().Pesc_total()*TauLines[ipTOI11].Emis().Aul()*1.76e-12;
00734 lindst(sum,11300,"6lev",ipnt,'i',false ,
00735 " OI 1.1 micron from six level atom");
00736
00737 PntForLine(29000.,"O 1",&ipnt);
00738 sum = atoms.popoi[5]*TauLines[ipTOI29].Emis().Pesc_total()*TauLines[ipTOI29].Emis().Aul()*6.86e-13;
00739 lindst(sum,29000,"6lev",ipnt,'i',false ,
00740 " OI 2.9 micron from six level atom");
00741
00742 PntForLine(46000.,"O 1",&ipnt);
00743 sum = atoms.popoi[5]*TauLines[ipTOI46].Emis().Pesc_total()*TauLines[ipTOI46].Emis().Aul()*4.32e-13;
00744 lindst(sum,46000,"6lev",ipnt,'i',false ,
00745 " OI 4.6 micron from six level atom");
00746
00747
00748
00749
00750
00751
00752
00753
00754
00755
00756
00757
00758
00759 if( dense.xIonDense[ipHYDROGEN][1] > SMALLFLOAT )
00760 {
00761
00762
00763 reco23tot = CoolHeavy.O2_A3_tot * HBeta *
00764 9.36 * phycon.te40*phycon.te04 / 57.544 * dense.xIonDense[ipOXYGEN][2]/dense.xIonDense[ipHYDROGEN][1];
00765 }
00766 else
00767 {
00768 reco23tot = 0.;
00769 }
00770
00771 sum = CoolHeavy.O2471*2471./7325. + CoolHeavy.O7323 + CoolHeavy.O7332;
00772 if( sum > SMALLFLOAT )
00773 {
00774
00775 reco23tot /= sum;
00776 }
00777 else
00778 {
00779 reco23tot = 0.;
00780 }
00781
00782 rec7323 = reco23tot * CoolHeavy.O7323;
00783 rec7332 = reco23tot * CoolHeavy.O7332;
00784 rec2471 = reco23tot * CoolHeavy.O2471*2471./7325. * 8.05e-12/2.72e-12;
00785
00786
00787
00788 reco22tot = 1.660e-10 / ( phycon.sqrte * phycon.te03 * phycon.te005 ) *
00789 dense.eden * dense.xIonDense[ipOXYGEN][2] * CoolHeavy.O2_A2_tot;
00790
00791 sum = CoolHeavy.O3726 + CoolHeavy.O3730;
00792 if( sum > SMALLFLOAT )
00793 {
00794 reco22tot /= sum;
00795 }
00796 else
00797 {
00798 reco22tot = 0.;
00799 }
00800
00801 rec3726 = reco22tot * CoolHeavy.O3726 * 5.34e-12;
00802 rec3730 = reco22tot * CoolHeavy.O3730 * 5.34e-12;
00803
00804
00805 oxy.s3727 = (realnum)((oxy.s3727 + oxy.s7325*0.5)*5.34e-12*
00806 9.7e-5/(9.7e-5 + dense.eden*1.15e-6/phycon.sqrte));
00807
00808 PntForLine(3727.,"O 2",&ipnt);
00809 fac = CoolHeavy.c3727+oxy.s3727+rec3726+rec3730;
00810 lindst(fac ,3727,"TOTL",ipnt,'c',true,
00811 " O II 3727, all lines of multiplet together " );
00812
00813 PntForLine(7325.,"O 2",&ipnt);
00814 fac = CoolHeavy.c7325+rec7323+rec7332;
00815 lindst( fac ,7325,"TOTL",ipnt,'c',true ,
00816 " O II 7325, all lines of multiplet together");
00817
00818 linadd(oxy.s3727,3727,"IONZ",'i',
00819 " line produced by photoionization of Oo; already in TOTL" );
00820 oxy.s7325 = (realnum)(oxy.s7325*2.72e-12*0.34/(0.34 + dense.eden*
00821 6.04e-6/phycon.sqrte));
00822
00823 linadd(oxy.s7325,7325,"IONZ",'i',
00824 " line produced by photoionization of Oo; already in TOTL" );
00825
00826 linadd(CoolHeavy.c7325,7325,"Coll",'i',
00827 " collisional contribution to line " );
00828
00829 linadd(CoolHeavy.c3727,3727,"Coll",'i',
00830 " collisional contribution to line " );
00831
00832 linadd(CoolHeavy.O3730,3729,"O II",'i',
00833 " five level atom calculations; D5/2 - S3/2" );
00834
00835 linadd(CoolHeavy.O3726,3726,"O II",'i',
00836 " D3/2 - S3/2 transition" );
00837
00838 linadd(CoolHeavy.O2471,2471,"O II",'c',
00839 " both 2P 1/2 and 3/2 to ground " );
00840
00841 linadd(CoolHeavy.O7323,7323,"O II",'i',
00842 " P1/2-D5/2 and P3/2-D5/2 together" );
00843
00844 linadd(CoolHeavy.O7332,7332,"O II",'i',
00845 " P1/2-D3/2 and P3/2-D3/2 together " );
00846
00847 linadd( rec3730 ,3729,"O 2r",'i',
00848 " recombination contribution refer o2 rec Liu, X-W., Storey, P.J., Barlow, M.J., Danziger, I.J.,refercon Cohen, M., & Bryce, M., 2000, MNRAS, 312, 585 recombination contributions five level atom calculations; D5/2 - S3/2 " );
00849
00850 linadd( rec3726 ,3726,"O 2r",'i',
00851 " D3/2 - S3/2 transition" );
00852
00853 linadd(rec2471,2471,"O 2r",'i',
00854 " both 2P 1/2 and 3/2 to ground " );
00855 linadd(rec7323,7323,"O 2r",'i',
00856 " P1/2-D5/2 and P3/2-D5/2 together " );
00857
00858 linadd(rec7332,7332,"O 2r",'i',
00859 " P1/2-D3/2 and P3/2-D3/2 together " );
00860
00861 PutLine(TauLines[ipT834],
00862 " O II 833.8 coll excit ");
00863
00864
00865
00866
00867
00868 if( nWindLine > 0 )
00869 {
00870 pump = TauLine2[387].Emis().pump()*TauLine2[387].Emis().PopOpc();
00871 }
00872 else
00873 {
00874 pump = 0.;
00875 }
00876
00877 PntForLine(3120.,"O 2",&ipnt);
00878 lindst(pump*0.336 * 6.37e-12/(1.+dense.eden/1e12) ,3120,"pump",ipnt,'r',true,
00879 " OII 3113.62 - 3139.68 (8 lines) are only pumped, no recombination part" );
00880
00881 PntForLine(3300.,"O 2",&ipnt);
00882 lindst(pump*0.147 * 6.03e-12/(1.+dense.eden/1e12) ,3300,"pump",ipnt,'r',true,
00883 " OII 3277.56 - 3306.45 (6 lines) are only pumped, no recombination part" );
00884
00885 PntForLine(3762.,"O 2",&ipnt);
00886 lindst(pump*0.087 * 5.29e-12/(1.+dense.eden/1e12) ,3762,"pump",ipnt,'r',true,
00887 " OII 3739.76/3762.47/3777.42 (3 lines) are only pumped, no recombination part" );
00888
00889
00890 rec = GetLineRec(82, 4651 );
00891 PntForLine(4651.,"O 2",&ipnt);
00892 lindst(rec,4651,"O 2r",ipnt,'r',true,
00893 " O II 4651 total recombination, 4638.86-4696.35 (8 lines) " );
00894
00895
00896
00897 linadd(pump* 0.336 * 0.933 * 4.27e-12/(1.+dense.eden/1e12),4651,"O 2p",'i',
00898 " pumped part of line O II 4651 " );
00899
00900
00901 rec = GetLineRec(83, 4341 );
00902
00903 linadd(rec/(1.+dense.eden/1e12),4341,"O 2r",'i',
00904 " recombination contribution to O II 4341 line " );
00905
00906 linadd(pump* 0.147 * 0.661 * 4.58e-12/(1.+dense.eden/1e12),4341,"O 2p",'i',
00907 " pumped part of line O II 4341 " );
00908
00909 PntForLine(4341.,"O 2",&ipnt);
00910 lindst(rec+pump* 0.147 * 0.661 * 4.58e-12/(1.+dense.eden/1e12),4341,"TOTL",ipnt,'r',true,
00911 " total intensity, all processes, O II 4341" );
00912
00913
00914 rec = GetLineRec(84, 3736 );
00915
00916
00917 linadd(rec/(1.+dense.eden/1e12),3736,"O 2r",'i',"\n recombination part of O II 3736 line " );
00918 linadd(pump* 0.087 * 0.763 * 5.33e-12/(1.+dense.eden/1e12),3736,"O 2p",'i',
00919 " pumped part of line O II 3736" );
00920
00921 PntForLine(3736.,"O 2",&ipnt);
00922 lindst((rec+pump* 0.087 * 0.763 * 5.33e-12)/(1.+dense.eden/1e12),3736,"TOTL",ipnt,'r',true,
00923 " total intensity, all processes, O II 3736" );
00924
00925
00926
00927 efac = (emit_frac(TauLines[ipT1666]) + emit_frac(TauLines[ipT1661]))*0.5;
00928
00929 linadd(TauLines[ipT1666].Emis().xIntensity()+TauLines[ipT1661].Emis().xIntensity(),1665,"TOTL",'i',
00930 "total intensity of OIII] 1665, all processes " );
00931 PutLine(TauLines[ipT1661]," ");
00932
00933 PutLine(TauLines[ipT1666]," ");
00934
00935 linadd(ionbal.PhotoRate_Shell[ipOXYGEN][3][1][0]*dense.xIonDense[ipOXYGEN][1]*0.3*1.20e-11*efac,1665,"Phot",'i',
00936 " contribution to OIII 1665 due to inner shell (2s^2) ionization " );
00937
00938 linadd(oxy.AugerO3*1.20e-11*efac*0.27,1665,"Augr",'i',
00939 " contribution to OIII 1665 due to K-shell ionization " );
00940
00941 PntForLine(5007.,"O 3",&ipnt);
00942 lindst(CoolHeavy.c5007/(1.+1./3.01),5007,"O 3",ipnt,'c',true ,
00943 " O III 5007 alone, collisions, tot OIII is this times 1.333 fac = c5007/(1.+1./2.887) >>chng 01 may 04, branching ratio had been 2.887, revised to 3 as per refer o3 as Storey, P.J., & Zeippen, C.J., 2000, 312, 813-816 ");
00944
00945 PntForLine(4959.,"O 3",&ipnt);
00946 lindst(CoolHeavy.c5007/(1.+3.01),4959,"O 3",ipnt,'c',true,
00947 " O III 4959 alone, collisions, tot OIII is this times 4" );
00948
00949 PntForLine(4931.,"O 3",&ipnt);
00950 lindst(CoolHeavy.c5007/(1.+3.01)*4.09e-4 ,4931,"O 3",ipnt,'c',true ,
00951 " O III 4931 alone, collisions >>chng 01 jul 11, added this line >>refer o3 as Nussbaumer, H., & Storey, P., 1981, A&A, 99, 177 >>refer o3 as Mathis, J.S., & Liu, X.-W., 1999, ApJ, 521, 212-216 ");
00952
00953 linadd(oxy.d5007t/1.25,5007,"LOST",'i',
00954 " O III 5007 lost through excited state photo" );
00955
00956
00957 effec = 1.6/(1.6 + 0.9*dense.cdsqte);
00958
00959
00960 r4363 = 6.3e-21/(phycon.te70*phycon.te10)*dense.eden*dense.xIonDense[ipOXYGEN][3]*
00961 effec;
00962
00963
00964
00965
00966
00967 ct4363 = phycon.sqrte*1.3e-12*4.561e-12*dense.xIonDense[ipHYDROGEN][0]*dense.xIonDense[ipOXYGEN][3]*
00968 effec;
00969
00970 fac = CoolHeavy.c4363 + r4363 + ct4363;
00971 linadd(fac,4363,"TOTL",'i',
00972 " O III 4363, sum of rec, coll, ct excitation" );
00973
00974 PntForLine(4363.,"O 3",&ipnt);
00975 lindst(CoolHeavy.c4363,4363,"Coll",ipnt,'c',true,
00976 " O III 4363,collisions from five level atom " );
00977
00978 lindst(r4363,4363,"Rec ",ipnt,'r',true,
00979 " O III 4363 recombination, coefficient from Burgess and Seaton " );
00980
00981 PntForLine(2321.,"O 3",&ipnt);
00982 lindst(CoolHeavy.c4363*0.236,2321,"O 3",ipnt,'c',true ,
00983 " collisional excitation of 2321, 5-level atom");
00984 linadd(ct4363,4363,"C EX",'i' ,
00985 " call linadd( c4363*0.236 , 2321 , 'O 3','c') charge exchange, Dalgarno+Sternberg ApJ Let 257, L87. ");
00986
00987 linadd(dense.xIonDense[ipHYDROGEN][0]*dense.xIonDense[ipOXYGEN][3]*0.225*3.56e-12*1.34e-11*phycon.sqrte,
00988 5592,"C EX",'i'," charge exchange rate, D+S " );
00989
00990 PutLine(TauLines[ipTO88],
00991 " O III 88 micron, collisionally excited");
00992
00993 PutLine(TauLines[ipT52],
00994 " O III 52 micron, collisionally excited ");
00995
00996
00997
00998 rec = GetLineRec(331,835)*emit_frac(TauLines[ipT835]);
00999 PutExtra(MAX2(0.,rec));
01000
01001 PutLine(TauLines[ipT835],
01002 " O III 834A, collisions and dielectronic recombination ");
01003
01004 linadd(MAX2(0.,rec),835,"rec ",'i',
01005 " O III 834A, dielectronic recombination only" );
01006
01007 PutLine(TauLines[ipT26],
01008 " O IV 26 micron ");
01009
01010 linadd(
01011 TauLines[ipO4_1400].Emis().xIntensity()+
01012 TauLines[ipO4_1397].Emis().xIntensity()+
01013 TauLines[ipO4_1407].Emis().xIntensity()+
01014 TauLines[ipO4_1405].Emis().xIntensity()+
01015 TauLines[ipO4_1401].Emis().xIntensity(),
01016 1402,"TOTL",'i',
01017 " total intensity of O IV] 1402, all lines in the multiplet " );
01018
01019 PutLine(TauLines[ipO4_1400],
01020 " ");
01021 PutLine(TauLines[ipO4_1397],
01022 " ");
01023 PutLine(TauLines[ipO4_1407],
01024 " ");
01025 PutLine(TauLines[ipO4_1405],
01026 " ");
01027 PutLine(TauLines[ipO4_1401],
01028 " ");
01029
01030 linadd(ionbal.PhotoRate_Shell[ipOXYGEN][2][1][0]*dense.xIonDense[ipOXYGEN][2]*0.43*1.42e-11,1401,"InSh",'i',
01031 " inner shell photoionization, relaxation " );
01032
01033
01034 rec = GetLineRec(378,789)*(1.-TauLines[ipT789].Emis().ColOvTot());
01035 rec = GetLineRec(378,789)*emit_frac(TauLines[ipT789]);
01036 PutExtra(MAX2(0.,rec));
01037
01038 PutLine(TauLines[ipT789]," O IV 789A");
01039
01040 linadd(MAX2(0.,rec),789,"rec ",'i',
01041 " O IV 789A, dielectronic recombination only" );
01042
01043
01044 rec = GetLineRec(466,630);
01045 PutExtra(MAX2(0.,rec));
01046
01047 PutLine(TauLines[ipT630],"O V 630, collisional excitation and dielectronic recombination");
01048
01049 linadd(MAX2(0.,rec),630,"rec ",'i',
01050 " O V 630A, dielectronic recombination only" );
01051
01052 linadd(embesq.em1218+TauLines[ipT1214].Emis().xIntensity(),1218,"TOTL",'i',
01053 " O V 1218], total intensity of both lines " );
01054 PutLine(TauLines[ipT1214], " ");
01055
01056 linadd(embesq.em1218,1211,"O 5",'i',
01057 " the slow transition by itself" );
01058
01059 linadd(1.4e-21/phycon.te70*dense.eden*dense.xIonDense[ipOXYGEN][5]*
01060
01061 emit_frac(TauLines[ipT1214]),5112,"O 5",'i',
01062 " BS O V 5112, recombination " );
01063
01064 linadd(TauLines[ipT1032].Emis().xIntensity()+TauLines[ipT1037].Emis().xIntensity(),1035,"TOTL",'i',
01065 " O VI 1035, total of pumping and collisional excitation " );
01066 sum = TauLines[ipT1032].Emis().xIntensity()*TauLines[ipT1032].Emis().FracInwd() +
01067 TauLines[ipT1037].Emis().xIntensity()* TauLines[ipT1037].Emis().FracInwd();
01068
01069 linadd(sum,1035,"Inwd",'i',
01070 " inward part of OVI line" );
01071 PutLine(TauLines[ipT1032],
01072 " ");
01073 PutLine(TauLines[ipT1037],
01074 " ");
01075
01076 PutLine(TauLines[ipT150],
01077 "O VI 150, Li seq 2s 3p ");
01078
01079 PutLine(TauLines[ipTNe13],
01080 " neon Neon II 12.8 micron ");
01081
01082 PutLine(TauLines[ipTNe16],
01083 " Ne III fine structure line ");
01084
01085 PutLine(TauLines[ipTNe36],
01086 " Ne III fine structure line ");
01087
01088 PntForLine(3869.,"Ne 3",&ipnt);
01089 lindst(CoolHeavy.c3869/(1.+1./3.318),3869,"Ne 3",ipnt,'c',true,
01090 " Ne III 3869, of 3968+3869 doublet" );
01091
01092 PntForLine(3968.,"Ne 3",&ipnt);
01093 lindst(CoolHeavy.c3869/(1.+3.318),3968,"Ne 3",ipnt,'c',true,
01094 " Ne III 3968, of 3968+3869 doublet" );
01095
01096 PntForLine(3343.,"Ne 3",&ipnt);
01097 lindst(CoolHeavy.c3343,3343,"Ne 3",ipnt,'c',true,
01098 " NeIII auroral line " );
01099
01100 PntForLine(1815.,"Ne 3",&ipnt);
01101 lindst(CoolHeavy.c3343*1.38,1815,"Ne 3",ipnt,'c',true ,
01102 " NeIII auroral line");
01103
01104 PntForLine(2424.,"Ne 4",&ipnt);
01105 lindst(CoolHeavy.c2424,2424,"Ne 4",ipnt,'c',true,
01106 " Ne IV 2424, collisional excitation" );
01107
01108 PntForLine(4720.,"Ne 4",&ipnt);
01109 lindst(CoolHeavy.c4720,4720,"Ne 4",ipnt,'c',true,
01110 " Ne IV N=3-2 lines, three level atom approx, this is the sum of the 4714.5, 4724.2, 4725.5 lines" );
01111
01112 PntForLine(1602.,"Ne 4",&ipnt);
01113 lindst(CoolHeavy.c4720*4.34,1602,"Ne 4",ipnt,'c',true ,
01114 " Ne IV N=3 lines, three level atom approx");
01115
01116 PntForLine(3426.,"Ne 5",&ipnt);
01117 lindst(CoolHeavy.c3426/(1.+1./2.738),3426,"Ne 5",ipnt,'c',true,
01118 " Ne V 3426 of 3426, 3346 doublet" );
01119
01120 PntForLine(3346.,"Ne 5",&ipnt);
01121 lindst(CoolHeavy.c3426/(1.+2.738),3346,"Ne 5",ipnt,'c',true,
01122 " Ne V 3346 of 3426, 3346 doublet " );
01123
01124 PntForLine(2976.,"Ne 5",&ipnt);
01125 lindst(CoolHeavy.c2975,2976,"Ne 5",ipnt,'c',true,
01126 " auroral line " );
01127
01128 PntForLine(1575.,"Ne 5",&ipnt);
01129 lindst(CoolHeavy.c1565,1575,"Ne 5",ipnt,'c',true,
01130 " collisionally excited" );
01131
01132 PutLine(TauLines[ipTNe24],"\n Ne V 24.2, 14.3 micron ");
01133
01134 PutLine(TauLines[ipTNe14],"\n Ne V 24.2, 14.3 micron ");
01135
01136 PntForLine(1141.,"Ne 5",&ipnt);
01137 lindst(CoolHeavy.c1134,1141,"Ne 5",ipnt,'c',true," both components of 5S-3P 1146.1, 1137.0 doublet " );
01138
01139 PutLine(TauLines[ipxNe0676],"\n [Ne VI] 7.6 microns ");
01140
01141 linadd(embesq.em895+TauLines[ipT895].Emis().xIntensity(),895,"TOTL",'i',
01142 " Ne VII 895, collisionally excited, both lines " );
01143
01144 PutLine(TauLines[ipT895],
01145 " Ne VII 895, only fast transition ");
01146
01147 linadd(embesq.em895,890,"Ne 7",'i',
01148 " Ne VII 890, single line " );
01149
01150 linadd(TauLines[ipT770].Emis().xIntensity()+TauLines[ipT780].Emis().xIntensity(),774,"TOTL",'i',
01151 " Ne VIII 774, collisionally excited " );
01152
01153 sum = TauLines[ipT770].Emis().xIntensity()*TauLines[ipT770].Emis().FracInwd() +
01154 TauLines[ipT780].Emis().xIntensity()*TauLines[ipT780].Emis().FracInwd();
01155 linadd(sum,774,"Inwd",'i',
01156 " inward part of NeVIII 774 line" );
01157
01158 PutLine(TauLines[ipT770],
01159 " the NeVIII 770 780 doublet ");
01160 PutLine(TauLines[ipT780],
01161 " ");
01162
01163 PutLine(TauLines[ipT88],
01164 " Ne VIII 88 2s 3p, collisionally excited ");
01165
01166 if( trace.lgTrace )
01167 {
01168 fprintf( ioQQQ, " lines_lv1_li_ne returns\n" );
01169 }
01170 return;
01171 }
01172
01173
01174
01175 STATIC double GetLineRec(
01176
01177 long int ip,
01178
01179 long int lWl)
01180 {
01181 double GetLineRec_v;
01182
01183 DEBUG_ENTRY( "GetLineRec()" );
01184
01185 if( (long)(LineSave.RecCoefCNO[2][ip]+0.5) != lWl )
01186 {
01187 fprintf( ioQQQ, " GetLineRec called with incorrect wavelength.\n" );
01188 fprintf( ioQQQ, " index, call and get wl are %5ld%5ld%5ld\n",
01189 ip, lWl, (long)(LineSave.RecCoefCNO[2][ip]+0.5) );
01190 cdEXIT(EXIT_FAILURE);
01191 }
01192
01193
01194 GetLineRec_v = LineSave.RecCoefCNO[3][ip]*dense.eden*
01195 dense.xIonDense[(long)(LineSave.RecCoefCNO[0][ip])-1][(long)(LineSave.RecCoefCNO[0][ip]-LineSave.RecCoefCNO[1][ip]+2)-1]*
01196 1.99e-8/LineSave.RecCoefCNO[2][ip];
01197
01198
01199
01200 LineSave.RecCoefCNO[3][ip] = 0.;
01201 return( GetLineRec_v );
01202 }