59 double rec7323 , rec7332, rec3730 , rec3726 , rec2471,
60 reco23tot , reco22tot;
66 fprintf(
ioQQQ,
" lines_lv1_li_ne called\n" );
72 " start level 1 ines" );
75 " cooling due to collisional ionization of heavy elements" );
92 " C 1 1656, collision strength from van Regemoter");
94 linadd(rec,1656,
"REC ",
'i',
95 " C 1 1656 recomb; n.b. coll deexcitation not in" );
98 " C 1 9850, coll excit" );
107 linadd(rec,9850,
"C Ir",
'i',
108 " C I 9850 recombination contribution" );
112 " total intensity, all processes, C I 9850");
117 "C 1 8727; equivalent to 4363" );
129 " C 2 158 micron, both e- and H0, H2 in excitation ");
138 " total intensity of C II] 2326, all lines in the multiplet " );
151 " photoproduction, Helfand and Trefftz");
161 " total intensity of C 2 1335");
163 linadd(rec,1335,
"REC ",
'i',
164 " C 2 1335 recombination," );
181 lindst(pump*0.387 * 5.08e-12/(1.+
dense.
eden/1e12) ,3920,
"pump",ipnt,
'i',
true ,
182 " CII 3918.98/3920.68 is only pumped, no recombination part");
187 pump *= 0.305 * 0.387 * 3.02e-12;
189 " recombination part of C II 6580 line " );
191 " pumped part of line C II 6580" );
195 " total intensity, all processes, C II 6580");
222 " total C 3] 977, recombination + collisional + pumped excitation ");
224 linadd(rec,977,
"C3 R",
'i',
225 " dielectronic recombination contribution to C 3 977 " );
227 linadd(p386,977,
"P386",
'r',
228 " C 3 977 pumped by continuum near 386A" );
233 "C 3 1909 collision, both lines together");
240 " C 3 1908 j-2 to ground" );
243 " the 13C forbidden line of C III " );
251 linadd(3.1e-19*fac*corr,1909,
"C3 R",
'i',
252 " C 3 1909 recombination from Storey" );
255 " C 3 1909 following relax following inner shell photoionization" );
264 " C 3* 1175, excited state line, above 1909 ");
267 " dielectronic recombination contribution to C 3 1175 " );
278 rec,1549,
"TOTL",
'i',
"total intensity of C 4 1549, all processes " );
285 "inward part of C 4 " );
298 " part of line destroyed by photoionization of Balmer continuum " );
300 linadd(rec,1549,
"C4 r",
'i',
301 " recombination C 4 1549 from CV" );
304 " Li seq 2s 3p Li seq transition");
319 double recrate_new = 1.108e-13 * pow( (
phycon.
te*1e-4) , -0.6085) /
320 (1. - 0.0041 * pow( (
phycon.
te*1e-4) , -0.3975 ) ) *
333 linadd(rec,5199,
"TOTr",
'i',
334 " recombination Nitrogen I 5200 + 5198 together " );
340 linadd(chem,5200,
"chem",
'i',
341 " upper limit to [N I] 5200 produced by chemistry" );
345 linadd(pump,5199,
"pump",
'i',
346 " guess of contribution from pumping by FUV lines, evaluated in CoolNitr assumes 100% of triplet pumps goes to singlets - this is surely a large overestimate so will only print intensity, do not add to total line " );
353 linadd(ctRate,5200,
"H CT",
'i',
354 " upper limit to [N I] 5200 produced by charge transfer" );
357 " [N I] 5200 - just collisions" );
360 " [N I] 5198 - just collisions" );
364 "[N I] 5200 - all processes - stat weight is 6 - total in term is 10" );
367 "[N I] 5198 - all processes - stat weight is 4 - total in term is 10 " );
371 " [N I] 5200 + 5198 together " );
374 "[N I] 3466 together " );
377 " [N I] 3467 together " );
380 " [N I] 3467, 3466 together " );
386 " Nitrogen I 10397, 10408, 10407 together");
389 " Nitrogen I 10398, 10397, 10408, 10407 together");
395 " N 1 5200, both 5198, 5200, collisions and recombination ");
398 " N 1 5200, both 5198, 5200, collisions and recombination " );
400 linadd(rec,5200,
"REC ",
'i',
401 " recombination contribution to [NI] 5200" );
405 " [N 1] 3466, 3 - 1 transition, whole multiplet" );
409 " [N 1] 10400 3 - 2 transition, whole multiplet");
412 " collisionally excited part ");
415 " 6 lines with fake collision strengths" );
429 " N 2 6584 alone " );
433 " N 2 6548 alone " );
438 " N 2 6584 alone, recombination contribution" );
477 " N 2 5755 total, collisions plus charge transfer plus recombination" );
480 " N 2 5755 collisional contribution" );
482 linadd(ctRate,5755,
"C T ",
'i',
483 " N 2 5755 charge transfer contribution " );
485 linadd( rec ,5755,
"N 2r",
'i',
486 " N 2 5755 recombination contribution" );
489 " N 2 fine structure line ");
492 " N 2 fine structure line " );
495 " N 2 2140 intercombination line " );
504 " N 2 1084, CS guess from g-bar " );
507 " dielectronic recombination contribution to N 2 1085" );
512 linadd(rnii*0.377*0.75*3.02e-12*efficn2,6584,
"N2cn",
'i',
513 " continuum pumped N 2 6584 " );
516 linadd(rnii*0.0117*3.46e-12*efficn2,5755,
"N2cn",
'i',
517 " continuum pumped N 2 5755" );
535 lindst(pump*0.236 * 6.01e-12/(1.+
dense.
eden/1e12) ,3311,
"pump",ipnt,
'i',
true,
536 " NII 3311.42 - 3331.31 (6 lines) are only pumped, no recombination part" );
539 lindst(pump*0.186 * 5.18e-12/(1.+
dense.
eden/1e12) ,3840,
"pump",ipnt,
'i',
true,
540 " NII 3829.8-3856.06 (6 lines) are only pumped, no recombination part" );
543 lindst(pump*0.025 * 5.52e-12/(1.+
dense.
eden/1e12) ,3609,
"pump",ipnt,
'i',
true,
544 " NII 3593.60/3609.1/3615.86 (3 lines) are only pumped, no recombination part" );
548 lindst(pump*0.186*0.595 * 4.31e-12/(1.+
dense.
eden/1e12) ,4640,
"pump",ipnt,
'i',
true ,
549 " NII 4601.5-4643.1 (6 lines) are only pumped, no recombination part");
552 lindst(pump*0.025*0.442 * 3.97e-12/(1.+
dense.
eden/1e12) ,5010,
"pump",ipnt,
'i',
true,
553 " NII 5002.7/5010.6/5045.1 (3 lines) are only pumped, no recombination part " );
558 pump *= 0.236 * 0.626 * 3.50e-12;
561 " recombination part of N II 5679 line" );
564 " pumped part of line N II 5679 " );
568 " total intensity, all processes, N II 5679 " );
571 " [N 3] 57 micron fine structure line");
580 " total intensity of N III] 1750, all lines in the multiplet " );
614 rn3tot = (rb + raten3)*0.439 + rn3mor;
616 sum = raten3*4.29e-12;
618 lindst(sum,4640,
"N3cn",ipnt,
'i',
true ,
619 " continuum pumped \"Bowen\" N 3, optically thin excited line ");
622 sum = rb*4.29e-12*0.834;
624 lindst(sum,4634,
"N3cn",ipnt,
'i',
true ,
625 " continuum pumped \"Bowen\" N 3, optically thin excited line");
627 sum = rb*4.29e-12*(1. - 0.834);
629 lindst(sum,4642,
"N3cn",ipnt,
'i',
true,
630 " continuum pumped \"Bowen\" N 3, optically thin excited line" );
645 " N 3 990, all processes ");
647 linadd(rec+rn3tot*2.01e-11*fac,990,
"extr",
'i',
648 " total N 3 990, both electron excitation and continuum pumping" );
650 linadd(rec,990,
"rec ",
'i',
651 " part of N 3 990 due to recombination " );
653 linadd(rn3tot*2.01e-11,990,
"N 3p",
'r',
654 " N 3 989.8, continuum pumped" );
657 " N 4] 1486, total intensity of both lines" );
663 " the N 4] slow transition by itself " );
674 " N 4 765, collisionally excited");
677 " N 4 765 recombination," );
680 " continuum pumping of NV 1240, N 5 1240, total emission, collisions plus pumping " );
684 linadd(sum,1240,
"Inwd",
'i',
685 " inward part of N 5 " );
692 " N 5 209, 2s-3p Li seq ");
696 " oxygen total Oxygen I 6300, including continuum optical depth ");
712 6300.,
"OH p",ipnt ,
'i' ,
false,
713 " the intensity of [OI] 6300 line due to OH photodistruction");
716 6363.,
"OH p",ipnt ,
'i' ,
false,
717 " the intensity of [OI] 6363 line due to OH photodistruction ");
721 " total Oxygen I 6363, including continuum optical depth " );
728 lindst( r13 , 5577.,
"OH p",ipnt ,
'i',
false,
729 " 94% of excitations to highest level decay via 5577" );
732 " O I fine structure line ");
735 " O I fine structure line ");
738 " total collisional cooling due to 6-level OI atom" );
741 " total collisional heating due to 6-level OI atom " );
750 lindst(sum,8446,
"6lev",ipnt,
'i',
false,
751 " * \todo 2 these are actually in the optical depth arrays, should be moved to call PutLine - but some needed vars are not defined in structure " );
755 lindst(sum,1304,
"6lev",ipnt,
'i',
false,
756 " OI 1304 from six level atom " );
760 lindst(sum,1039,
"6lev",ipnt,
'i',
false ,
761 " OI 1039 from six level atom");
765 lindst(sum,4368,
"6lev",ipnt,
'i',
false,
766 " OI 4368 from six level atom" );
770 lindst(sum,13100,
"6lev",ipnt,
'i',
false ,
771 "OI 1.3 micron from six level atom");
775 lindst(sum,11300,
"6lev",ipnt,
'i',
false ,
776 " OI 1.1 micron from six level atom");
780 lindst(sum,29000,
"6lev",ipnt,
'i',
false ,
781 " OI 2.9 micron from six level atom");
785 lindst(sum,46000,
"6lev",ipnt,
'i',
false ,
786 " OI 4.6 micron from six level atom");
825 rec2471 = reco23tot *
CoolHeavy.
O2471*2471./7325. * 8.05e-12/2.72e-12;
851 lindst(fac ,3727,
"TOTL",ipnt,
'c',
true,
852 " O II 3727, all lines of multiplet together " );
856 lindst( fac ,7325,
"TOTL",ipnt,
'c',
true ,
857 " O II 7325, all lines of multiplet together");
860 " line produced by photoionization of Oo; already in TOTL" );
865 " line produced by photoionization of Oo; already in TOTL" );
868 " collisional contribution to line " );
871 " collisional contribution to line " );
874 " five level atom calculations; D5/2 - S3/2" );
877 " D3/2 - S3/2 transition" );
880 " both 2P 1/2 and 3/2 to ground " );
883 " P1/2-D5/2 and P3/2-D5/2 together" );
886 " P1/2-D3/2 and P3/2-D3/2 together " );
888 linadd( rec3730 ,3729,
"O 2r",
'i',
889 " chng 01 jul 08, add 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 " );
891 linadd( rec3726 ,3726,
"O 2r",
'i',
892 " D3/2 - S3/2 transition" );
894 linadd(rec2471,2471,
"O 2r",
'i',
895 " both 2P 1/2 and 3/2 to ground " );
896 linadd(rec7323,7323,
"O 2r",
'i',
897 " P1/2-D5/2 and P3/2-D5/2 together " );
899 linadd(rec7332,7332,
"O 2r",
'i',
900 " P1/2-D3/2 and P3/2-D3/2 together " );
903 " O II 833.8 coll excit ");
919 lindst(pump*0.336 * 6.37e-12/(1.+
dense.
eden/1e12) ,3120,
"pump",ipnt,
'i',
true,
920 " OII 3113.62 - 3139.68 (8 lines) are only pumped, no recombination part" );
923 lindst(pump*0.147 * 6.03e-12/(1.+
dense.
eden/1e12) ,3300,
"pump",ipnt,
'i',
true,
924 " OII 3277.56 - 3306.45 (6 lines) are only pumped, no recombination part" );
927 lindst(pump*0.087 * 5.29e-12/(1.+
dense.
eden/1e12) ,3762,
"pump",ipnt,
'i',
true,
928 " OII 3739.76/3762.47/3777.42 (3 lines) are only pumped, no recombination part" );
933 lindst(rec,4651,
"O 2r",ipnt,
'i',
true,
934 " O II 4651 total recombination, 4638.86-4696.35 (8 lines) " );
938 linadd(pump* 0.336 * 0.933 * 4.27e-12/(1.+
dense.
eden/1e12),4651,
"O 2p",
'i',
939 " pumped part of line O II 4651 " );
945 " recombination contribution to O II 4341 line " );
947 linadd(pump* 0.147 * 0.661 * 4.58e-12/(1.+
dense.
eden/1e12),4341,
"O 2p",
'i',
948 " pumped part of line O II 4341 " );
951 lindst(rec+pump* 0.147 * 0.661 * 4.58e-12/(1.+
dense.
eden/1e12),4341,
"TOTL",ipnt,
'i',
true,
952 " total intensity, all processes, O II 4341" );
958 linadd(rec/(1.+
dense.
eden/1e12),3736,
"O 2r",
'i',
"\n recombination part of O II 3736 line " );
959 linadd(pump* 0.087 * 0.763 * 5.33e-12/(1.+
dense.
eden/1e12),3736,
"O 2p",
'i',
960 " pumped part of line O II 3736" );
963 lindst((rec+pump* 0.087 * 0.763 * 5.33e-12)/(1.+
dense.
eden/1e12),3736,
"TOTL",ipnt,
'i',
true,
964 " total intensity, all processes, O II 3736" );
971 "total intensity of OIII] 1665, all processes " );
977 " contribution to OIII 1665 due to inner shell (2s^2) ionization " );
980 " contribution to OIII 1665 due to K-shell ionization " );
984 " 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 ");
988 " O III 4959 alone, collisions, tot OIII is this times 4" );
992 " 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 ");
995 " O III 5007 lost through excited state photo" );
1013 lindst(fac,4363,
"TOTL",ipnt,
'i',
true,
1014 " O III 4363, sum of rec, coll, ct excitation" );
1017 " O III 4363,collisions from five level atom " );
1019 linadd(r4363,4363,
"Rec ",
'i',
1020 " O III 4363 recombination, coefficient from Burgess and Seaton " );
1024 " collisional excitation of 2321, 5-level atom");
1025 linadd(ct4363,4363,
"C EX",
'i' ,
1026 " call linadd( c4363*0.236 , 2321 , 'O 3','c') charge exchange, Dalgarno+Sternberg ApJ Let 257, L87. ");
1029 5592,
"C EX",
'i',
" charge exchange rate, D+S " );
1032 " O III 88 micron, collisionally excited");
1035 " O III 52 micron, collisionally excited ");
1043 " O III 834A, collisions and dielectronic recombination ");
1046 " O III 834A, dielectronic recombination only" );
1049 " O IV 26 micron ");
1058 " total intensity of O IV] 1402, all lines in the multiplet " );
1072 " inner shell photoionization, relaxation " );
1082 " O IV 789A, dielectronic recombination only" );
1091 " O V 630A, dielectronic recombination only" );
1094 " O V 1218], total intensity of both lines " );
1098 " the slow transition by itself" );
1103 " BS O V 5112, recombination " );
1106 " O VI 1035, total of pumping and collisional excitation " );
1110 linadd(sum,1035,
"Inwd",
'i',
1111 " inward part of OVI line" );
1118 "O VI 150, Li seq 2s 3p ");
1121 " fluorine * [F II] 29.33 microns");
1124 " [F II] 67.2 microns ");
1127 " F IV 44.07 micron ");
1130 " F IV 25.83 micron");
1133 " neon Neon II 12.8 micron ");
1136 " Ne III fine structure line ");
1139 " Ne III fine structure line ");
1143 " Ne III 3869, of 3968+3869 doublet" );
1147 " Ne III 3968, of 3968+3869 doublet" );
1151 " NeIII auroral line " );
1155 " NeIII auroral line");
1159 " Ne IV 2424, collisional excitation" );
1163 " Ne IV N=3-2 lines, three level atom approx, this is the sum of the 4714.5, 4724.2, 4725.5 lines" );
1167 " Ne IV N=3 lines, three level atom approx");
1171 " Ne V 3426 of 3426, 3346 doublet" );
1175 " Ne V 3346 of 3426, 3346 doublet " );
1183 " collisionally excited" );
1190 lindst(
CoolHeavy.
c1134,1141,
"Ne 5",ipnt,
'c',
true,
" both components of 5S-3P 1146.1, 1137.0 doublet " );
1195 " Ne VII 895, collisionally excited, both lines " );
1198 " Ne VII 895, only fast transition ");
1201 " Ne VII 890, single line " );
1204 " Ne VIII 774, collisionally excited " );
1208 linadd(sum,774,
"Inwd",
'i',
1209 " inward part of NeVIII 774 line" );
1212 " the NeVIII 770 780 doublet ");
1217 " Ne VIII 88 2s 3p, collisionally excited ");
1221 fprintf(
ioQQQ,
" lines_lv1_li_ne returns\n" );
1234 double GetLineRec_v;
1240 fprintf(
ioQQQ,
" GetLineRec called with incorrect wavelength.\n" );
1241 fprintf(
ioQQQ,
" index, call and get wl are %5ld%5ld%5ld\n",
1254 return( GetLineRec_v );