/home66/gary/public_html/cloudy/c08_branch/source/prt_lines_lv1_k_zn.cpp

Go to the documentation of this file.
00001 /* This file is part of Cloudy and is copyright (C)1978-2008 by Gary J. Ferland and
00002  * others.  For conditions of distribution and use see copyright notice in license.txt */
00003 /*lines_lv1_k_zn place lines of elements potassium and heavier into lines storage stack */
00004 #include "cddefines.h"
00005 #include "cddrive.h"
00006 #include "coolheavy.h"
00007 #include "ca.h"
00008 #include "fe.h"
00009 #include "rfield.h"
00010 #include "dense.h"
00011 #include "phycon.h"
00012 #include "radius.h"
00013 #include "taulines.h"
00014 #include "trace.h"
00015 #include "lines_service.h"
00016 #include "rt.h"
00017 #include "atomfeii.h"
00018 #include "lines.h"
00019 
00020 void lines_lv1_k_zn(void)
00021 {
00022         long int i, 
00023           ipnt,
00024           ilo,
00025           ihi;
00026 
00027         double c10, 
00028           c14, 
00029           eff, 
00030           fela, 
00031           r14;
00032 
00033         DEBUG_ENTRY( "lines_lv1_k_zn()" );
00034 
00035         if( trace.lgTrace )
00036         {
00037                 fprintf( ioQQQ, "   lines_lv1_k_zn called\n" );
00038         }
00039 
00040         PutLine(&TauLines[ipKI7745],
00041                 "  potassium K I 7745 ");
00042 
00043         PutLine(&TauLines[ipxK03462],
00044                 "  [K III] 4.62 microns ");
00045 
00046         PutLine(&TauLines[ipxK04598],
00047                 "  [KIV] 5.983 min  ");
00048 
00049         PutLine(&TauLines[ipxK04154],
00050                 "  [KIV] 15.39 micron ");
00051 
00052         PutLine(&TauLines[ipxK06882],
00053                 "  [KVI] 8.823 micron ");
00054 
00055         PutLine(&TauLines[ipxK06557],
00056                 "  [KVI]  5.575 micron ");
00057 
00058         PutLine(&TauLines[ipxK07319],
00059                 "  [K VII] 3.189 microns ");
00060 
00061         PutLine(&TauLines[ipxK11425],
00062                 " K 11 4249.99A ");
00063 
00064         PutLine(&TauLines[ipCaI4228],
00065                 "  calcium Ca I 4228 ");
00066 
00067         linadd(ca.Cakh,3933,"Ca 2",'c',
00068                 " coll excited calcium k+h " );
00069 
00070         linadd(ca.Cair,8579,"Ca 2",'c' ,
00071                 " infrared triplet ");
00072 
00073         linadd(ca.c7306,7306,"Ca 2",'c',
00074                 " forbidden lines, 7291+7324 together " );
00075 
00076         linadd(ca.dCakh,3933,"Phot",'i' ,
00077                 " fraction H Ly-alpha destruction of excited levels ");
00078 
00079         linadd(ca.dCaf12,7306,"Phot",'i' ,
00080                 " fraction H Ly-alpha destruction of excited levels ");
00081 
00082         PntForLine(3934.,"Ca2K",&ipnt);
00083         lindst(ca.Cak,3934,"Ca2K",ipnt,'i',true,
00084                 " individual lines from five level atom");
00085 
00086 
00087         PntForLine(3969.,"Ca2H",&ipnt);
00088         lindst(ca.Cah,3969,"Ca2H",ipnt,'i',true,
00089                 " individual lines from five level atom" );
00090 
00091 
00092         PntForLine(8498.,"Ca2X",&ipnt);
00093         lindst(ca.Cax,8498,"Ca2X",ipnt,'i',true,
00094                 " individual lines from five level atom " );
00095 
00096 
00097         PntForLine(8542.,"Ca2Y",&ipnt);
00098         lindst(ca.Cay,8542,"Ca2Y",ipnt,'i',true,
00099                 "  individual lines from five level atom" );
00100 
00101 
00102         PntForLine(8662.,"Ca2Z",&ipnt);
00103         lindst(ca.Caz,8662,"Ca2Z",ipnt,'i',true,
00104                 " individual lines from five level atom" );
00105 
00106 
00107         PntForLine(7291.,"CaF1",&ipnt);
00108         lindst(ca.Caf1,7291,"CaF1",ipnt,'i',true,
00109                 " individual lines from five level atom" );
00110 
00111 
00112         PntForLine(7324.,"CaF2",&ipnt);
00113         lindst(ca.Caf2,7324,"CaF2",ipnt,'i',true,
00114                 " individual lines from five level atom" );
00115 
00116 
00117         eff = dense.eden*dense.xIonDense[ipCALCIUM][2]*5.4e-21/(phycon.te/
00118           phycon.te10/phycon.te10);
00119         linadd(eff,3933,"Rec ",'i',
00120                 " recombination contribution to CaII emission" );
00121 
00122         PutLine(&TauLines[ipTCa3],
00123                 "  Ca IV 3.2 micron ");
00124 
00125         PutLine(&TauLines[ipTCa4],
00126                 "  Ca V 4.16, 11.47 micron");
00127 
00128         PutLine(&TauLines[ipTCa12],
00129                 "  Ca V 4.16, 11.47 micron ");
00130 
00131 
00132         PntForLine(6087.,"Ca 5",&ipnt);
00133         lindst(ca.Ca6087,6087,"Ca 5",ipnt,'c',true ,
00134                 "  Ca V optical and UV lines, collisional excitation, 3-level atom");
00135 
00136 
00137         PntForLine(5311.,"Ca 5",&ipnt);
00138         lindst(ca.c5311,5311,"Ca 5",ipnt,'c',true ,
00139                 " Ca V optical and UV lines, collisional excitation, 3-level atom");
00140 
00141 
00142         PntForLine(2414.,"Ca 5",&ipnt);
00143         lindst(ca.c2414,2414,"Ca 5",ipnt,'c',true ,
00144                 "  Ca V optical and UV lines, collisional excitation, 3-level atom");
00145 
00146 
00147         PntForLine(3997.,"Ca 5",&ipnt);
00148         lindst(ca.c3997,3997,"Ca 5",ipnt,'c',true,
00149                 " Ca V optical and UV lines, collisional excitation, 3-level atom" );
00150 
00151         PutLine(&TauLines[ipCa0741],
00152                 "  [Ca VII] 4.09 microns" );
00153 
00154         PutLine(&TauLines[ipCa0761],
00155                 "  [Ca VII] 6.15 microns " );
00156 
00157 
00158         PntForLine(5620.,"Ca 7",&ipnt);
00159         lindst(ca.Ca5620,5620,"Ca 7",ipnt,'c',true,
00160                 " Ca VII optical and UV lines, collisional excitation, 3-level atom" );
00161 
00162 
00163         PntForLine(4941.,"Ca 7",&ipnt);
00164         lindst(ca.Ca4941,4941,"Ca 7",ipnt,'c',true,
00165                 " Ca VII optical and UV lines, collisional excitation, 3-level atom" );
00166 
00167 
00168         PntForLine(2112.,"Ca 7",&ipnt);
00169         lindst(ca.Ca2112,2112,"Ca 7",ipnt,'c',true ,
00170                 " Ca VII optical and UV lines, collisional excitation, 3-level atom");
00171 
00172 
00173         PntForLine(3688.,"Ca 7",&ipnt);
00174         lindst(ca.Ca3688,3688,"Ca 7",ipnt,'c',true,
00175                 " Ca VII optical and UV lines, collisional excitation, 3-level atom" );
00176 
00177         PutLine(&TauLines[ipCa08232],
00178                 "  [Ca VIII]  2.32 microns, A Saraph and Strey ");
00179 
00180         PutLine(&TauLines[ipCa12333],
00181                 "  [Ca 12] 3328.78A ");
00182 
00183         PutLine(&TauLines[ipTCa302],
00184                 "  Ca 18 Li seq 2s2p, 302, 345 separate ");
00185 
00186         PutLine(&TauLines[ipTCa345],
00187                 "  Ca 18 Li seq 2s2p, 302, 345 separate ");
00188 
00189         PutLine(&TauLines[ipTCa19],
00190                 "  Ca 18 Li seq 2s3p,  ");
00191 
00192 
00193         PntForLine(22.08e4,"Sc 2",&ipnt);
00194         lindst(CoolHeavy.Sc22p08m,22.08e4,"Sc 2",ipnt,'c',true,
00195                 " Sc II 2.08 (1-3) " );
00196 
00197 
00198         PntForLine(24.1e4,"Sc 2",&ipnt);
00199         lindst(CoolHeavy.Sc24p1m,24.1e4,"Sc 2",ipnt,'c',true,
00200                 " Sc II 4.1 micron (1-2)" );
00201 
00202 
00203         PntForLine(24.2e4,"Sc 2",&ipnt);
00204         lindst(CoolHeavy.Sc24p2m,24.2e4,"Sc 2",ipnt,'c',true,
00205                 "  Sc II 4.22 (2-3)" );
00206 
00207 
00208         PntForLine(3933.,"Sc 3",&ipnt);
00209         lindst(CoolHeavy.Sc33936,3933,"Sc 3",ipnt,'c',true,
00210                 " Sc III 3936" );
00211 
00212         PutLine(&TauLines[ipSc05231],
00213                 "  [Sc V] 1.46 microns ");
00214 
00215 
00216         PntForLine(5054.,"Sc 6",&ipnt);
00217         lindst(CoolHeavy.Sc45058,5054,"Sc 6",ipnt,'c',true ,
00218                 " Sc VI 5054 (1-2)");
00219 
00220 
00221         PntForLine(3592.,"Sc 6",&ipnt);
00222         lindst(CoolHeavy.Sc43595,3592,"Sc 6",ipnt,'c',true,
00223                 " Sc VI 3595 (2-3)" );
00224 
00225 
00226         PntForLine(2100.,"Sc 6",&ipnt);
00227         lindst(CoolHeavy.Sc42100,2100,"Sc 6",ipnt,'c',true,
00228                 "  Sc VI 2100 (1-3)" );
00229 
00230         PutLine(&TauLines[ipSc13264],
00231                 "  [Sc 13] 2637.97A");
00232 
00233 
00234         PntForLine(1.21e4,"Ti 3",&ipnt);
00235         lindst(CoolHeavy.Ti3l21,1.21e4,"Ti 3",ipnt,'c',true,
00236                 " Ti III 1.21 micron, (actually multiplet) 2-1 transition from model atom " );
00237 
00238 
00239         PntForLine(9594.,"Ti 3",&ipnt);
00240         lindst(CoolHeavy.Ti3l31,9594,"Ti 3",ipnt,'c',true,
00241                 " Ti III 9594, 3-1 transition, (actually multiplet) from model atom" );
00242 
00243         PntForLine(4.57e4,"Ti 3",&ipnt);
00244         lindst(CoolHeavy.Ti3l32,4.57e4,"Ti 3",ipnt,'c',true,
00245                 " Ti III 4.57 micron, 3-2 transition, (actually multiplet) from model atom" );
00246 
00247         PutLine(&TauLines[ipTi06172],
00248                 "  [Ti VI] 1.72 microns ");
00249 
00250         PutLine(&TauLines[ipTi14212],
00251                 "  [Ti XIV] 2117.79 ");
00252 
00253 
00254         PntForLine(8823.,"V  3",&ipnt);
00255         lindst(CoolHeavy.V38830,8823,"V  3",ipnt,'c',true ,
00256                 "  V III 8823 ");
00257 
00258 
00259         PntForLine(8507.,"V  3",&ipnt);
00260         lindst(CoolHeavy.V38507,8507,"V  3",ipnt,'c',true,
00261                 "  V III 8507" );
00262 
00263 
00264         PntForLine(7735.,"V  4",&ipnt);
00265         lindst(CoolHeavy.V47741,7735,"V  4",ipnt,'c',true,
00266                 "  V IV 7741 1-3" );
00267 
00268 
00269         PntForLine(9489.,"V  4",&ipnt);
00270         lindst(CoolHeavy.V49496,9489,"V  4",ipnt,'c',true,
00271                 " V IV 9496 2-1 " );
00272 
00273 
00274         PntForLine(4.19e4,"V  4",&ipnt);
00275         lindst(CoolHeavy.V44p2m,4.19e4,"V  4",ipnt,'c',true,
00276                 "  V IV 4.19 micron 3-2" );
00277 
00278         PutLine(&TauLines[ipVa07130],
00279                 "  [V VII] 1.304 microns ");
00280 
00281         PutLine(&TauLines[ipVa15172],
00282                 " [V 15] 1721.38 ");
00283 
00284         PntForLine(5828.,"Cr 3",&ipnt);
00285         lindst(CoolHeavy.Cr3l21,5828,"Cr 3",ipnt,'c',true,
00286                 " [CrIII] multiplet blend at 5828A" );
00287 
00288         PntForLine(7267.,"Cr 4",&ipnt);
00289         lindst(CoolHeavy.Cr4l21,7267,"Cr 4",ipnt,'c',true,
00290                 " [CrIV] 2 - 1 multiplet blend at 7272" );
00291 
00292 
00293         PntForLine(6801.,"Cr 4",&ipnt);
00294         lindst(CoolHeavy.Cr4l31,6801,"Cr 4",ipnt,'c',true,
00295                 " [CrIV] 3 - 1 multiplet blend at 6806" );
00296 
00297 
00298         PntForLine(7979.,"Cr 5",&ipnt);
00299         lindst(CoolHeavy.Cr5l21,7979,"Cr 5",ipnt,'c',true,
00300                 " [CrV] 2 - 1 multiplet blend at 7985" );
00301 
00302         PntForLine(6577.,"Cr 5",&ipnt);
00303         lindst(CoolHeavy.Cr5l31,6577,"Cr 5",ipnt,'c',true,
00304                 "  [CrV] 3 - 1 multiplet blend at 6582" );
00305 
00306 
00307         PntForLine(3.75e4,"Cr 5",&ipnt);
00308         lindst(CoolHeavy.Cr5l32,3.75e4,"Cr 5",ipnt,'c',true,
00309                 " [CrV] 3 - 2 multiplet blend at 3.75 microns " );
00310 
00311         PutLine(&TauLines[ipCr08101],
00312                 "  [Cr VIII] 1.01 microns ");
00313 
00314         PutLine(&TauLines[ipCr16141],
00315                 "  [Cr 16] 1410.60 ");
00316 
00317         PutLine(&TauLines[ipxMn0979],
00318                 " [Mn IX] 7968.5 A ");
00319 
00320         PutLine(&TauLines[ipxMn1712],
00321                 " [Mn 17] 1169.59 ");
00322 
00323         /* bob Rubin's UV line
00324          * f2 = dense.xIonDense(26,4)*sexp(50 764./te)*0.45*cdsqte/6.*7.01e-12
00325          * call linadd( f2 , 2837 , 'BobR' , 'i')
00326          * f2 = dense.xIonDense(26,4)*sexp(55 989./te)*0.384*cdsqte/6.*7.74e-12
00327          * call linadd( f2 , 2568 , 'BobR' , 'i') */
00328 
00329         /* iron */
00330 
00331         PutLine(&TauLines[ipFe1_24m],
00332                 "  Fe 1 24m ");
00333 
00334         PutLine(&TauLines[ipFe1_35m],
00335                 "  Fe 1 35m ");
00336 
00337         PutLine(&TauLines[ipFe1_54m],
00338                 "  Fe 1 54m ");
00339 
00340         PutLine(&TauLines[ipFe1_111m],
00341                 "  Fe 1 111m ");
00342 
00343         PutLine(&TauLines[ipFeI3884],
00344                 "  Fe 1 3884 ");
00345 
00346         PutLine(&TauLines[ipFeI3729],
00347                 "  Fe 1 3729 ");
00348 
00349         PutLine(&TauLines[ipFeI3457],
00350                 "  Fe 1 3457 ");
00351 
00352         PutLine(&TauLines[ipFeI3021],
00353                 "  Fe 1 3021 ");
00354 
00355         PutLine(&TauLines[ipFeI2966],
00356                 "  Fe 1 2966 ");
00357 
00358         linadd(MAX2(0.,FeII.Fe2_large_cool+FeII.Fe2_UVsimp_cool),0,"Fe2c",'c' ,
00359                 "total of all Fe 2, both simple UV and large atom together ");
00360 
00361         linadd(MAX2(0.,-FeII.Fe2_large_cool-FeII.Fe2_UVsimp_cool),0,"Fe2h",'h' ,
00362                 "total of all Fe 2 heating, both simple UV and large atom together ");
00363 
00364         linadd(FeII.for7,4300,"Fe 2",'i' ,
00365                 " Fe 2 forbidden 2-1 transition from Netzer's atom ");
00366 
00367         PutLine(&TauLines[ipTuv3],
00368                 " 2400 ");
00369         PutLine(&TauLines[ipTr48],
00370                 " 6200 ");
00371         PutLine(&TauLines[ipTFe16],
00372                 " 1080 ");
00373         PutLine(&TauLines[ipTFe26],
00374                 " 1500  ");
00375         PutLine(&TauLines[ipTFe34],
00376                 " 11500  ");
00377         PutLine(&TauLines[ipTFe35],
00378                 " 2500 ");
00379         PutLine(&TauLines[ipTFe46],
00380                 " 2300  ");
00381         PutLine(&TauLines[ipTFe56],
00382                 " 8900  ");
00383 
00384         /* >>chng 03 nov 02, changed wl in lindst calls to be wl in microns,
00385          * this will change wl in printout */
00386         /* results of FeIILevelPops Fe 2 atom, 
00387          * not all are printed, first is */
00388         /* total cooling due to 16 level atom 
00389          * >>chng 05 dec 04, rm this since no longer meaningful 
00390         linadd(FeII.Fe2_16levl_cool,0,"Fe 2",'c' );*/
00391 
00392 
00393         PntForLine(16640.,"Fe 2",&ipnt);
00394         lindst(FeII.fe21308,16640,"Fe 2",ipnt,'i',true,
00395                 " Fe 2 1.664 microns 8-13 " );
00396 
00397 
00398         PntForLine(15990.,"Fe 2",&ipnt);
00399         lindst(FeII.fe21207,15990,"Fe 2",ipnt,'i',true,
00400                 " Fe 2 1.599 microns 7-12" );
00401 
00402 
00403         PntForLine(15340.,"Fe 2",&ipnt);
00404         lindst(FeII.fe21106,15340,"Fe 2",ipnt,'i',true ,
00405                 " Fe 2 1.534 microns 6-11 ");
00406 
00407 
00408         PntForLine(16440.,"Fe 2",&ipnt);
00409         lindst(FeII.fe21006,16440,"Fe 2",ipnt,'i',true,
00410                 " Fe 2 1.644 microns 6-10 " );
00411 
00412 
00413         PntForLine(12790.,"Fe 2",&ipnt);
00414         lindst(FeII.fe21204,12790,"Fe 2",ipnt,'i',true,
00415                 " Fe 2 1.279 microns 12-4" );
00416 
00417 
00418         PntForLine(12950.,"Fe 2",&ipnt);
00419         lindst(FeII.fe21103,12950,"Fe 2",ipnt,'i',true,
00420                 " Fe 2 1.295 microns 11-3 " );
00421 
00422 
00423         PntForLine(13280.,"Fe 2",&ipnt);
00424         lindst(FeII.fe21104,13280,"Fe 2",ipnt,'i',true,
00425                 "  Fe 2 1.328 microns 11-4" );
00426 
00427 
00428         PntForLine(12570.,"Fe 2",&ipnt);
00429         lindst(FeII.fe21001,12570,"Fe 2",ipnt,'i',true,
00430                 " Fe 2 1.257 microns 10-1 " );
00431 
00432 
00433         PntForLine(13210.,"Fe 2",&ipnt);
00434         lindst(FeII.fe21002,13210,"Fe 2",ipnt,'i',true,
00435                 " Fe 2 1.321 microns 10-2 " );
00436 
00437 
00438         PntForLine(259880.,"Fe 2",&ipnt);
00439         lindst(FeII.fe20201,259880,"Fe 2",ipnt,'i',true,
00440                 "  Fe 2 25.988 microns 2-1" );
00441 
00442         PntForLine(353480.,"Fe 2",&ipnt);
00443         lindst(FeII.fe20302,353480,"Fe 2",ipnt,'i',true,
00444                 " Fe 2 35.348 microns 3-2" );
00445 
00446 
00447         PntForLine(179360.,"Fe 2",&ipnt);
00448         lindst(FeII.fe20706,179360,"Fe 2",ipnt,'i',true,
00449                 "  Fe 2 17.936 microns 7-6, label is 178 to be unique " );
00450 
00451 
00452         PntForLine(245180.,"Fe 2",&ipnt);
00453         lindst(FeII.fe20807,245180,"Fe 2",ipnt,'i',true,
00454                 " Fe 2 24.518 microns 8-7" );
00455 
00456         PntForLine(357760.,"Fe 2",&ipnt);
00457         lindst(FeII.fe20908,357760,"Fe 2",ipnt,'i',true,
00458                 " Fe 2 35.776 microns 9-8" );
00459 
00460         PntForLine(18100.,"Fe 2",&ipnt);
00461         lindst(FeII.fe21007,18100,"Fe 2",ipnt,'i',true,
00462                 " Fe 2 1.810 microns 10-7" );
00463 
00464 
00465         PntForLine(16770.,"Fe 2",&ipnt);
00466         lindst(FeII.fe21107,16770,"Fe 2",ipnt,'i',true,
00467                 "  Fe 2 1.677 microns 11-7 " );
00468 
00469 
00470         PntForLine(18000.,"Fe 2",&ipnt);
00471         lindst(FeII.fe21108,18000,"Fe 2",ipnt,'i',true,
00472                 "  Fe 2 1.800 microns 11-8" );
00473 
00474 
00475         PntForLine(17120.,"Fe 2",&ipnt);
00476         lindst(FeII.fe21208,17120,"Fe 2",ipnt,'i',true,
00477                 " Fe 2 1.712 microns 12-8 " );
00478 
00479 
00480         PntForLine(17980.,"Fe 2",&ipnt);
00481         lindst(FeII.fe21209,17980,"Fe 2",ipnt,'i',true,
00482                 " Fe 2 1.798 microns 12-9" );
00483 
00484 
00485         PntForLine(229020.,"Fe 2",&ipnt);
00486         lindst(FeII.fe21110,229020,"Fe 2",ipnt,'i',true,
00487                 " Fe 2 22.902 microns 11-10 " );
00488 
00489 
00490         PntForLine(346600.,"Fe 2",&ipnt);
00491         lindst(FeII.fe21211,346600,"Fe 2",ipnt,'i',true,
00492                 "  Fe 2 34.660 microns 12-11" );
00493 
00494         /* >>chng 97 mar 14, added following 4 lines after consultation with Verner
00495          * these are typically within 5-9% of their large model atom.  These are
00496          * predicted with both large and small atom */
00497 
00498 
00499         PntForLine(8617.,"Fe 2",&ipnt);
00500         lindst(FeII.fe21406,8617,"Fe 2",ipnt,'i',true,
00501                 "  Fe 2 8617A 14-06," );
00502 
00503 
00504         PntForLine(8892.,"Fe 2",&ipnt);
00505         lindst(FeII.fe21507,8892,"Fe 2",ipnt,'i',true,
00506                 " Fe 2 8892A 15-07" );
00507 
00508 
00509         PntForLine(9227.,"Fe 2",&ipnt);
00510         lindst(FeII.fe21508,9227,"Fe 2",ipnt,'i',true,
00511                 " Fe 2 9227A 15-08 " );
00512 
00513 
00514         PntForLine(9268.,"Fe 2",&ipnt);
00515         lindst(FeII.fe21609,9268,"Fe 2",ipnt,'i',true,
00516                 " Fe 2 9268A 16-09" );
00517 
00518         /* following lines are only non-zero when large atom is enabled.  When
00519          * running with small FeII atoms these are set zero in routine  */
00520 
00521 
00522         PntForLine(5159.,"Fe 2",&ipnt);
00523         lindst(FeII.fe25to6,5159,"Fe 2",ipnt,'i',true,
00524                 " Fe 2 5158.81 25-6" );
00525 
00526 
00527         PntForLine(5262.,"Fe 2",&ipnt);
00528         lindst(FeII.fe27to7,5262,"Fe 2",ipnt,'i',true,
00529                 " Fe 2 5261.61 27-7" );
00530 
00531         PntForLine(5334.,"Fe 2",&ipnt);
00532         lindst(FeII.fe28to8,5334,"Fe 2",ipnt,'i',true,
00533                 " Fe 2 5333.65 28-8 " );
00534 
00535 
00536         PntForLine(5376.,"Fe 2",&ipnt);
00537         lindst(FeII.fe29to9,5376,"Fe 2",ipnt,'i',true,
00538                 " Fe 2 5376.47 29-9 " );
00539 
00540         PntForLine(4815.,"Fe 2",&ipnt);
00541         lindst(FeII.fe32to6,4815,"Fe 2",ipnt,'i',true,
00542                 "  Fe 2 4814.55 32-6 " );
00543 
00544         PntForLine(4905.,"Fe 2",&ipnt);
00545         lindst(FeII.fe33to7,4905,"Fe 2",ipnt,'i',true,
00546                 " Fe 2 4905.35 33-7" );
00547 
00548 
00549         PntForLine(4347.,"Fe 2",&ipnt);
00550         lindst(FeII.fe37to7,4347,"Fe 2",ipnt,'i',true,
00551                 " Fe 2 4346.85 37-7 " );
00552 
00553 
00554         PntForLine(4353.,"Fe 2",&ipnt);
00555         lindst(FeII.fe39to8,4353,"Fe 2",ipnt,'i',true,
00556                 " Fe 2 4352.78 39-8 " );
00557 
00558 
00559         PntForLine(4372.,"Fe 2",&ipnt);
00560         lindst(FeII.fe40to9,4372,"Fe 2",ipnt,'i',true,
00561                 " Fe 2 4372.43 40-9 " );
00562 
00563 
00564         PntForLine(4244.,"Fe 2",&ipnt);
00565         lindst(FeII.fe37to6,4244,"Fe 2",ipnt,'i',true,
00566                 " Fe 2 4243.98 37-6" );
00567 
00568 
00569         PntForLine(4277.,"Fe 2",&ipnt);
00570         lindst(FeII.fe39to7,4277,"Fe 2",ipnt,'i',true,
00571                 " Fe 2 4276.83 39-7" );
00572 
00573 
00574         PntForLine(4320.,"Fe 2",&ipnt);
00575         lindst(FeII.fe40to8,4320,"Fe 2",ipnt,'i',true,
00576                 "  Fe 2 4319.62 40-8" );
00577 
00578 
00579         PntForLine(4358.,"Fe 2",&ipnt);
00580         lindst(FeII.fe41to9,4358,"Fe 2",ipnt,'i',true,
00581                 " Fe 2 4358.37 41-9 " );
00582 
00583 
00584         PntForLine(4177.,"Fe 2",&ipnt);
00585         lindst(FeII.fe39to6,4177,"Fe 2",ipnt,'i',true,
00586                 " Fe 2 4177.21 39-6 " );
00587 
00588 
00589         PntForLine(4245.,"Fe 2",&ipnt);
00590         lindst(FeII.fe40to7,4245,"Fe 2",ipnt,'i',true,
00591                 " Fe 2 4244.81 40-7 " );
00592 
00593 
00594         PntForLine(4306.,"Fe 2",&ipnt);
00595         lindst(FeII.fe41to8,4306,"Fe 2",ipnt,'i',true,
00596                 " Fe 2 4305.90 41-8" );
00597 
00598         /* >>chng 99 jun 18 new lines added as per Katya email, so all orion lines are in */
00599 
00600         PntForLine(4114.,"Fe 2",&ipnt);
00601         lindst(FeII.fe42to6,4114,"Fe 2",ipnt,'i',true,
00602                 "  Fe 2 4114.48 42-6 " );
00603 
00604 
00605         PntForLine(4179.,"Fe 2",&ipnt);
00606         lindst(FeII.fe43to7,4179,"Fe 2",ipnt,'i',true,
00607                 " Fe 2 4178.95  7-      43" );
00608 
00609 
00610         PntForLine(4211.,"Fe 2",&ipnt);
00611         lindst(FeII.fe42to7,4211,"Fe 2",ipnt,'i',true,
00612                 " Fe 2 4211.10  7-      42 " );
00613 
00614 
00615         PntForLine(4314.,"Fe 2",&ipnt);
00616         lindst(FeII.fe80to28,4314,"Fe 2",ipnt,'i',true,
00617                 "  Fe 2 4314.29 28-     80" );
00618 
00619         PntForLine(4359.,"Fe 2",&ipnt);
00620         lindst(FeII.fe36to2,4359,"Fe 2",ipnt,'i',true,
00621                 " Fe 2 4359.34  2-      36" );
00622 
00623 
00624         PntForLine(4414.,"Fe 2",&ipnt);
00625         lindst(FeII.fe36to3,4414,"Fe 2",ipnt,'i',true,
00626                 "  Fe 2 4413.78 3-      36 " );
00627 
00628 
00629         PntForLine(4416.,"Fe 2",&ipnt);
00630         lindst(FeII.fe32to1,4416,"Fe 2",ipnt,'i',true ,
00631                 " Fe 2 4416.27  1-      32");
00632 
00633 
00634         PntForLine(4458.,"Fe 2",&ipnt);
00635         lindst(FeII.fe33to2,4458,"Fe 2",ipnt,'i',true,
00636                 " Fe 2 4457.95  2-      33 " );
00637 
00638 
00639         PntForLine(4475.,"Fe 2",&ipnt);
00640         lindst(FeII.fe36to5,4475,"Fe 2",ipnt,'i',true,
00641                 " Fe 2 4474.91  5-      36 " );
00642 
00643 
00644         PntForLine(4493.,"Fe 2",&ipnt);
00645         lindst(FeII.fe32to2,4493,"Fe 2",ipnt,'i',true,
00646                 "  Fe 2 4492.64 2-      32" );
00647 
00648 
00649         PntForLine(4515.,"Fe 2",&ipnt);
00650         lindst(FeII.fe33to3,4515,"Fe 2",ipnt,'i',true,
00651                 "  Fe 2 4514.90 3-      33" );
00652 
00653 
00654         PntForLine(4728.,"Fe 2",&ipnt);
00655         lindst(FeII.fe30to3,4728,"Fe 2",ipnt,'i',true,
00656                 "  Fe 2 4728.07 3-      30" );
00657 
00658 
00659         PntForLine(4775.,"Fe 2",&ipnt);
00660         lindst(FeII.fe33to6,4775,"Fe 2",ipnt,'i',true,
00661                 " Fe 2 4774.74  6-      33" );
00662 
00663         PntForLine(4890.,"Fe 2",&ipnt);
00664         lindst(FeII.fe24to2,4890,"Fe 2",ipnt,'i',true,
00665                 "  Fe 2 4889.63 2-      24" );
00666 
00667 
00668         PntForLine(4947.,"Fe 2",&ipnt);
00669         lindst(FeII.fe32to7,4947,"Fe 2",ipnt,'i',true,
00670                 "  Fe 2 4947.37 7-      32" );
00671 
00672 
00673         PntForLine(4951.,"Fe 2",&ipnt);
00674         lindst(FeII.fe35to8,4951,"Fe 2",ipnt,'i',true,
00675                 " Fe 2 4950.74  8-      35" );
00676 
00677 
00678         PntForLine(4973.,"Fe 2",&ipnt);
00679         lindst(FeII.fe34to8,4973,"Fe 2",ipnt,'i',true,
00680                 " Fe 2 4973.39  8-      34" );
00681 
00682 
00683         PntForLine(5112.,"Fe 2",&ipnt);
00684         lindst(FeII.fe27to6,5112,"Fe 2",ipnt,'i',true,
00685                 " Fe 2 5111.63  6-      27" );
00686 
00687 
00688         PntForLine(5220.,"Fe 2",&ipnt);
00689         lindst(FeII.fe28to7,5220,"Fe 2",ipnt,'i',true,
00690                 " Fe 2 5220.06  7-      28" );
00691 
00692 
00693         PntForLine(5269.,"Fe 2",&ipnt);
00694         lindst(FeII.fe30to8,5269,"Fe 2",ipnt,'i',true,
00695                 " Fe 2 5268.880 8-      30 " );
00696 
00697 
00698         PntForLine(5273.,"Fe 2",&ipnt);
00699         lindst(FeII.fe24to6,5273,"Fe 2",ipnt,'i',true,
00700                 " Fe 2 5273.35  6-      24" );
00701 
00702         PntForLine(5297.,"Fe 2",&ipnt);
00703         lindst(FeII.fe29to8,5297,"Fe 2",ipnt,'i',true,
00704                 " Fe 2 5296.840 8-      29 " );
00705 
00706 
00707         PntForLine(5433.,"Fe 2",&ipnt);
00708         lindst(FeII.fe24to7,5433,"Fe 2",ipnt,'i',true,
00709                 " Fe 2 5433.13  7-      24" );
00710 
00711 
00712         PntForLine(5527.,"Fe 2",&ipnt);
00713         lindst(FeII.fe22to7,5527,"Fe 2",ipnt,'i',true,
00714                 " Fe 2 5527.340 7-      22" );
00715 
00716 
00717         PntForLine(5747.,"Fe 2",&ipnt);
00718         lindst(FeII.fe38to11,5747,"Fe 2",ipnt,'i',true,
00719                 "  Fe 2 5746.96 11-     38" );
00720 
00721 
00722         PntForLine(6440.,"Fe 2",&ipnt);
00723         lindst(FeII.fe19to8,6440,"Fe 2",ipnt,'i',true,
00724                 " Fe 2 6440.400 8-      19" );
00725 
00726         PntForLine(7155.,"Fe 2",&ipnt);
00727         lindst(FeII.fe17to6,7155,"Fe 2",ipnt,'i',true,
00728                 " Fe 2 7155.140 6-      17 " );
00729 
00730 
00731         PntForLine(7172.,"Fe 2",&ipnt);
00732         lindst(FeII.fe18to7,7172,"Fe 2",ipnt,'i',true,
00733                 " Fe 2 7171.980 7-      18" );
00734 
00735 
00736         PntForLine(7388.,"Fe 2",&ipnt);
00737         lindst(FeII.fe18to8,7388,"Fe 2",ipnt,'i',true,
00738                 " Fe 2 7388.160 8-      18" );
00739 
00740 
00741         PntForLine(7453.,"Fe 2",&ipnt);
00742         lindst(FeII.fe17to7,7453,"Fe 2",ipnt,'i',true,
00743                 " Fe 2 7452.540 7-      17" );
00744 
00745         /* option to save all intensities predicted by large FeII atom,
00746          * code is in FeIILevelPops */
00747         FeIIAddLines();
00748 
00749         /* some bands from the large model Fe II atom
00750          * only done when full large atom is turned on */
00751         if( FeII.lgFeIILargeOn )
00752         {
00753                 for( i=0; i < nFeIIBands; i++ )
00754                 {
00755                         /* [i][0] is center wavelength, [i][1] and [i][2] are upper and
00756                          * lower bounds in Angstroms.  These are set in FeIIZero 
00757                          * units are erg s-1 cm-3 */
00758                         eff = FeIISumBand(FeII_Bands[i][1],FeII_Bands[i][2]);
00759 
00760                         linadd(eff,FeII_Bands[i][0],"Fe2b",'i' ,
00761                                 " total Fe II emission in Fe II bands, as defined in bands_Fe2.dat ");
00762                 }
00764                 if( LineSave.ipass > 0 )
00765                 {
00766                         for( i=0; i < nFeIIConBins; i++ )
00767                         {
00768                                 /* [i][0] is total intensity in cell, [i][1] and [i][2] are lower and
00769                                  * upper bounds in Angstroms.  these are set in FeIIZero *
00770                                  * find total emission from large FeII atom, integrated over band */
00771                                 FeII_Cont[i][0] += 
00772                                         (realnum)(FeIISumBand(FeII_Cont[i][1],FeII_Cont[i][2])*radius.dVeff);
00773                                 /*fprintf(ioQQQ,"DEBUG feii\t%li\t%.2e\n", i, FeII_Cont[i][0]);*/
00774                         }
00775                 }
00776                 else
00777                 {
00778                         for( i=0; i < nFeIIConBins; i++ )
00779                         {
00780                                 /* initialize arrays */
00781                                 FeII_Cont[i][0] = 0.;
00782                         }
00783                 }
00784         }
00785         else
00786         {
00787                 /* must still add the bands in case of grid where FeII is turned off and on */
00788                 for( i=0; i < nFeIIBands; i++ )
00789                 {
00790 
00791                         linadd(0., FeII_Bands[i][0] ,"Fe2b",'i' ,
00792                                 " total Fe II emission in Fe II bands, as defined in bands_Fe2.dat ");
00793                 }
00794         }
00795 
00796         PutLine(&TauLines[ipT191],
00797                 "  anomalous Fe 2 transition at 1787, RMT 191");
00798 
00799         linadd(fe.Fe3CoolTot,0,"Fe3c",'c' ,
00800                 " chng 05 dec 16, FeIII code created by Kevin Blagrave  Fe3c 0 - total cooling due to 14-level Fe 3 atom ");
00801 
00802         /* Fe 3 14-level atom 
00803          * following from print statements within loop */
00804         /* Fe 3 22.92m from Blagrave 14-level atom */
00805         /* Fe 3 13.53m from Blagrave 14-level atom */
00806         /* Fe 3 33.03m from Blagrave 14-level atom */
00807         /* Fe 3 10.72m from Blagrave 14-level atom */
00808         /* Fe 3 20.15m from Blagrave 14-level atom */
00809         /* Fe 3 51.67m from Blagrave 14-level atom */
00810         /* Fe 3 9.732m from Blagrave 14-level atom */
00811         /* Fe 3 16.91m from Blagrave 14-level atom */
00812         /* Fe 3 34.66m from Blagrave 14-level atom */
00813         /* Fe 3 105.3m from Blagrave 14-level atom */
00814         /* Fe 3  5152A from Blagrave 14-level atom */
00815         /* Fe 3  5271A from Blagrave 14-level atom */
00816         /* Fe 3  5356A from Blagrave 14-level atom */
00817         /* Fe 3  5412A from Blagrave 14-level atom */
00818         /* Fe 3  5440A from Blagrave 14-level atom */
00819         /* Fe 3  4986A from Blagrave 14-level atom */
00820         /* Fe 3  5097A from Blagrave 14-level atom */
00821         /* Fe 3  5177A from Blagrave 14-level atom */
00822         /* Fe 3  5230A from Blagrave 14-level atom */
00823         /* Fe 3  5256A from Blagrave 14-level atom */
00824         /* Fe 3 15.47m from Blagrave 14-level atom */
00825         /* Fe 3  4925A from Blagrave 14-level atom */
00826         /* Fe 3  5033A from Blagrave 14-level atom */
00827         /* Fe 3  5111A from Blagrave 14-level atom */
00828         /* Fe 3  5162A from Blagrave 14-level atom */
00829         /* Fe 3  5188A from Blagrave 14-level atom */
00830         /* Fe 3 11.16m from Blagrave 14-level atom */
00831         /* Fe 3 40.04m from Blagrave 14-level atom */
00832         /* Fe 3  4881A from Blagrave 14-level atom */
00833         /* Fe 3  4988A from Blagrave 14-level atom */
00834         /* Fe 3  5064A from Blagrave 14-level atom */
00835         /* Fe 3  5114A from Blagrave 14-level atom */
00836         /* Fe 3  5139A from Blagrave 14-level atom */
00837         /* Fe 3 9.282m from Blagrave 14-level atom */
00838         /* Fe 3 23.21m from Blagrave 14-level atom */
00839         /* Fe 3 55.20m from Blagrave 14-level atom */
00840         /* Fe 3  4833A from Blagrave 14-level atom */
00841         /* Fe 3  4937A from Blagrave 14-level atom */
00842         /* Fe 3  5012A from Blagrave 14-level atom */
00843         /* Fe 3  5061A from Blagrave 14-level atom */
00844         /* Fe 3  5085A from Blagrave 14-level atom */
00845         /* Fe 3 7.789m from Blagrave 14-level atom */
00846         /* Fe 3 15.69m from Blagrave 14-level atom */
00847         /* Fe 3 25.79m from Blagrave 14-level atom */
00848         /* Fe 3 48.41m from Blagrave 14-level atom */
00849         /* Fe 3  4714A from Blagrave 14-level atom */
00850         /* Fe 3  4813A from Blagrave 14-level atom */
00851         /* Fe 3  4884A from Blagrave 14-level atom */
00852         /* Fe 3  4931A from Blagrave 14-level atom */
00853         /* Fe 3  4954A from Blagrave 14-level atom */
00854         /* Fe 3 5.543m from Blagrave 14-level atom */
00855         /* Fe 3 8.638m from Blagrave 14-level atom */
00856         /* Fe 3 11.01m from Blagrave 14-level atom */
00857         /* Fe 3 13.76m from Blagrave 14-level atom */
00858         /* Fe 3 19.22m from Blagrave 14-level atom */
00859         /* Fe 3  4659A from Blagrave 14-level atom */
00860         /* Fe 3  4755A from Blagrave 14-level atom */
00861         /* Fe 3  4825A from Blagrave 14-level atom */
00862         /* Fe 3  4870A from Blagrave 14-level atom */
00863         /* Fe 3  4893A from Blagrave 14-level atom */
00864         /* Fe 3 4.859m from Blagrave 14-level atom */
00865         /* Fe 3 7.085m from Blagrave 14-level atom */
00866         /* Fe 3 8.608m from Blagrave 14-level atom */
00867         /* Fe 3 10.20m from Blagrave 14-level atom */
00868         /* Fe 3 12.92m from Blagrave 14-level atom */
00869         /* Fe 3 39.41m from Blagrave 14-level atom */
00870         /* Fe 3  4608A from Blagrave 14-level atom */
00871         /* Fe 3  4702A from Blagrave 14-level atom */
00872         /* Fe 3  4770A from Blagrave 14-level atom */
00873         /* Fe 3  4814A from Blagrave 14-level atom */
00874         /* Fe 3  4836A from Blagrave 14-level atom */
00875         /* Fe 3 4.356m from Blagrave 14-level atom */
00876         /* Fe 3 6.063m from Blagrave 14-level atom */
00877         /* Fe 3 7.146m from Blagrave 14-level atom */
00878         /* Fe 3 8.208m from Blagrave 14-level atom */
00879         /* Fe 3 9.884m from Blagrave 14-level atom */
00880         /* Fe 3 20.34m from Blagrave 14-level atom */
00881         /* Fe 3 42.06m from Blagrave 14-level atom */
00882         /* Fe 3  4574A from Blagrave 14-level atom */
00883         /* Fe 3  4668A from Blagrave 14-level atom */
00884         /* Fe 3  4734A from Blagrave 14-level atom */
00885         /* Fe 3  4778A from Blagrave 14-level atom */
00886         /* Fe 3  4800A from Blagrave 14-level atom */
00887         /* Fe 3 4.077m from Blagrave 14-level atom */
00888         /* Fe 3 5.535m from Blagrave 14-level atom */
00889         /* Fe 3 6.423m from Blagrave 14-level atom */
00890         /* Fe 3 7.269m from Blagrave 14-level atom */
00891         /* Fe 3 8.554m from Blagrave 14-level atom */
00892         /* Fe 3 15.41m from Blagrave 14-level atom */
00893         /* Fe 3 25.31m from Blagrave 14-level atom */
00894         /* Fe 3 63.56m from Blagrave 14-level atom */
00895         for( ihi=1; ihi<NLFE3; ++ihi )
00896         {
00897                 for( ilo=0; ilo<ihi; ++ilo )
00898                 {
00899                         /* emission in these lines */
00900                         PntForLine(fe.Fe3_wl[ihi][ilo],"Fe 3",&ipnt);
00901 #                       if 0
00902                         fprintf( ioQQQ,"\t/* FeIII ");
00903                         prt_wl( ioQQQ , (realnum)(fe.Fe3_wl[ihi][ilo]+0.5) );
00904                         fprintf( ioQQQ," from Blagrave 14-level atom */\n ");
00905 #                       endif
00906                         lindst( fe.Fe3_emiss[ihi][ilo] , (realnum)(fe.Fe3_wl[ihi][ilo]+0.5) , "Fe 3",ipnt,'c',true,
00907                                 " " );
00908                 }
00909         }
00910 
00911         /*>>chng 05 dec 18, following are now in the above */
00912         /* sum of 3p and 3g states together */
00913         /*      linadd(CoolHeavy.c5270,0,"Fe 3",'c' ); */
00914 
00915         /* Fe 3 5270, predictions from Garstang et al 78
00916         PntForLine(5270.,"Fe 3",&ipnt);
00917         lindst(CoolHeavy.c5270*0.2090,5270,"Fe 3",ipnt,'c',true );*/
00918 
00919         /* Fe 3 5270, predictions from Garstang et al 78 
00920         PntForLine(4658.,"Fe 3",&ipnt);
00921         lindst(CoolHeavy.c5270*0.3667,4658,"Fe 3",ipnt,'c',true ); */
00922 
00923         PutLine(&TauLines[ipT1122]," Fe 3 1122 entire multiplet");
00924 
00925         linadd(fe.Fe4CoolTot,0,"Fe4c",'c',
00926                 " Fe4c 0 - total cooling due to 12-level Fe 4 atom " );
00927 
00928 
00929         PntForLine(3096.,"Fe 4",&ipnt);
00930         lindst(fe.fe40401,3096,"Fe 4",ipnt,'i',true,
00931                 " Fe 4 3096.A, 4-1 and 5-1 transitions together"  );
00932 
00933 
00934         PntForLine(2836.,"Fe 4",&ipnt);
00935         lindst(fe.fe42836,2836,"Fe 4",ipnt,'i',true,
00936                 " Fe 4 2835.7A, 6-1 transition, 4P5/2 - 6S5/2 "  );
00937 
00938 
00939         PntForLine(2829.,"Fe 4",&ipnt);
00940         lindst(fe.fe42829,2829,"Fe 4",ipnt,'i',true,
00941                 "   Fe 4 2829.4A, 7-1 transition, 4P3/2 - 6S5/2"  );
00942 
00943 
00944         PntForLine(2567.,"Fe 4",&ipnt);
00945         lindst(fe.fe42567,2567,"Fe 4",ipnt,'i',true,
00946                 "  Fe 4 2567.6+ 2567.4. 11-1 and 12-1 transitions"  );
00947 
00948 
00949         PntForLine(2.774e4,"Fe 4",&ipnt);
00950         lindst(fe.fe41207,2.774e4,"Fe 4",ipnt,'i',true,
00951                 " Fe 4 2.774 microns 12-7 transition "  );
00952 
00953 
00954         PntForLine(2.714e4,"Fe 4",&ipnt);
00955         lindst(fe.fe41206,2.714e4,"Fe 4",ipnt,'i',true,
00956                 " Fe 4 2.714 microns 12-6 transition "  );
00957 
00958 
00959         PntForLine(2.716e4,"Fe 4",&ipnt);
00960         lindst(fe.fe41106,2.716e4,"Fe 4",ipnt,'i',true,
00961                 " Fe 4 2.716 microns 11-6 transition"  );
00962 
00963 
00964         PntForLine(2.806e4,"Fe 4",&ipnt);
00965         lindst(fe.fe41007,2.806e4,"Fe 4",ipnt,'i',true,
00966                 " Fe 4 2.806 microns 10-7 transition "  );
00967 
00968 
00969         PntForLine(2.865e4,"Fe 4",&ipnt);
00970         lindst(fe.fe41008,2.865e4,"Fe 4",ipnt,'i',true ,
00971                 "  Fe 4 2.865 microns 10-8 transition");
00972 
00973 
00974         PntForLine(2.836e4,"Fe 4",&ipnt);
00975         lindst(fe.fe40906,2.836e4,"Fe 4",ipnt,'i',true,
00976                 " Fe 4 2.836 microns 9-6 transition" );
00977 
00978 
00979         PntForLine(3892.,"Fe 5",&ipnt);
00980         lindst(CoolHeavy.c3892,3892,"Fe 5",ipnt,'c',true,
00981                 " Fe 5  3892+3839" );
00982 
00983         linadd(CoolHeavy.c5177,0,"Fe 6",'c' ,
00984                 " all of 2G lines together first ");
00985 
00986 
00987         PntForLine(5177.,"Fe 6",&ipnt);
00988         lindst(CoolHeavy.c5177*0.354,5177,"Fe 6",ipnt,'c',true,
00989                 " Fe 6 5177, approximate correct " );
00990 
00991         linadd(fe.Fe7CoolTot,0,"Fe7c",'c' ,
00992                 " Fe7c 0 - total cooling due to n-level Fe 7 atom ");
00993 
00994         /* >>chng 04 nov 04, move to multi-level system */
00995         for( ilo=0; ilo<NLFE7-1; ++ilo )
00996         {
00997                 /* must not do 1-0 or 2-1, which are transferred lines */
00998                 for( ihi=MAX2(3,ilo+1); ihi<NLFE7; ++ihi )
00999                 {
01000 
01001                         PntForLine(fe.Fe7_wl[ihi][ilo],"Fe 7",&ipnt);
01002                         lindst( fe.Fe7_emiss[ihi][ilo] , (realnum)(fe.Fe7_wl[ihi][ilo]+0.5) , "Fe 7",ipnt,'c',true,
01003                                 " emission in these lines" );
01004                 }
01005         }
01006 #       if 0
01007         PntForLine(5721.,"Fe 7",&ipnt);
01008         lindst( fe.Fe7_5721 , 5721 , "Fe 7",ipnt,'c',true,
01009                 " " );
01010 
01011         PntForLine(6601.,"Fe 7",&ipnt);
01012         lindst( fe.Fe7_6601 , 6601 , "Fe 7",ipnt,'c',true,
01013                 " " );
01014 
01015         PntForLine(3760.,"Fe 7",&ipnt);
01016         lindst( fe.Fe7_3760 , 3760 , "Fe 7",ipnt,'c',true,
01017                 " " );
01018 
01019         PntForLine(3588.,"Fe 7",&ipnt);
01020         lindst( fe.Fe7_3588 , 3588 , "Fe 7",ipnt,'c',true,
01021                 " " );
01022 #       endif
01023 
01024         PutLine(&TauLines[ipFe0795],
01025                 "   [Fe 7] 9.51 micron ");
01026 
01027         PutLine(&TauLines[ipFe0778],
01028                 "  [Fe 7] 7.81 micron ");
01029 
01030         /* [Fe 7] 6087 
01031         PntForLine(6087.,"Fe 7",&ipnt);
01032         lindst(CoolHeavy.c6087,6087,"Fe 7",ipnt,'c',true );*/
01033 
01034         /* [Fe 7] 5722 
01035         PntForLine(5722.,"Fe 7",&ipnt);
01036         lindst(CoolHeavy.Fe5722,5722,"Fe 7",ipnt,'c',true );*/
01037 
01038         PutLine(&TauLines[ipT245],
01039                 "   Be-seq lines ");
01040 
01041         PntForLine(242.,"Fe 7",&ipnt);
01042         lindst(CoolHeavy.c242,242,"Fe 7",ipnt,'c',true,
01043                 "  Fe 9 242 j=1 slower decay");
01044 
01045         PutLine(&TauLines[ipT352],
01046                 " the E1 transition that can pump [Fe X] ");
01047 
01048         /* optically thin Fe X pumping */
01049         eff = 69.4/(69.4 + 0.27*dense.cdsqte);
01050 
01051         /* coll excitation of 352 which decays to excited state of 6374
01052          * assumes 17/56 (ratio of A's) go to excited state */
01053         c10 = TauLines[ipT352].Emis->phots*eff*.01676*TauLines[ipT352].EnergyErg*352/6374.;
01054 
01055         /* Fe 10 and Fe 14 from Mason 75
01056          * total (coll, pumped)
01057          * call linadd( C6374+C10 , 6374 , 'Fe10','i') */
01058         PutExtra( c10 );
01059 
01060 
01061         PutLine(&TauLines[ipFe106375],
01062                 " [Fe 10] 6375, collisions with pumping too ");
01063 
01064         /* collisional contribution
01065          * call linadd( C6374 , 6374 , 'Coll','c')
01066          * collisions of E1 line, plus pumped by continuum fluorescence
01067          * call linadd( C10 , 6374 , ' 352','c')
01068          * Fe XI 7892, 6.08 micron 
01069         PutLine(&TauLines[ipTFe07]);*/
01070 
01082         /* Fe 11 2649 collisional excitation
01083         PntForLine(2649.,"Fe11",&ipnt);
01084         lindst(CoolHeavy.c2649,2649,"Fe11",ipnt,'c',true ); */
01085 
01086         /*  Fe 11 1467 collisional excitation 
01087         PntForLine(1467.,"Fe11",&ipnt);
01088         lindst(CoolHeavy.c1467,1467,"Fe11",ipnt,'c',true );*/
01089 
01090         linadd(fe.Fe11CoolTot,0,"Fe11",'c' ,
01091                 " >>chng 05 dec 18, add Fe 11  Fe11 0 - total cooling due to 5-level Fe 11 atom ");
01092 
01093         /* Fe 11 5-level atom */
01094         for( ihi=1; ihi<NLFE11; ++ihi )
01095         {
01096                 for( ilo=0; ilo<ihi; ++ilo )
01097                 {
01098                         PntForLine(fe.Fe11_wl[ihi][ilo],"Fe11",&ipnt);
01099                         lindst( fe.Fe11_emiss[ihi][ilo] , (realnum)(fe.Fe11_wl[ihi][ilo]+0.5) , "Fe11",ipnt,'c',true,
01100                                 " emission in these lines" );
01101                 }
01102         }
01103 
01104         PntForLine(1242.,"Fe12",&ipnt);
01105         lindst(CoolHeavy.c1242,1242,"Fe12",ipnt,'c',true ,
01106                 " Fe 12, 1242, 1349 together, collisional excitation");
01107 
01108         PntForLine(2170.,"Fe12",&ipnt);
01109         lindst(CoolHeavy.c2170,2170,"Fe12",ipnt,'c',true ,
01110                 " Fe 12, 2170, 2406 together, collisional excitation");
01111 
01112 
01113         PntForLine(2568.,"Fe12",&ipnt);
01114         lindst(CoolHeavy.c2568,2568,"Fe12",ipnt,'c',true,
01115                 "  Fe12 2904, 2567, 3567, 3073 together, collisional excitation" );
01116 
01117         /* >>chng 05 dec 18, add Fe 13  */
01118         /* Fe13 0 - total cooling due to 5-level Fe 13 atom */
01119         linadd(fe.Fe13CoolTot,0,"Fe13",'c' ,
01120                    "total cooling due to Fe 13 model atom ");
01121 
01122         /* Fe 13 5-level atom */
01123         for( ihi=1; ihi<NLFE13; ++ihi )
01124         {
01125                 for( ilo=0; ilo<ihi; ++ilo )
01126                 {
01127 
01128                         PntForLine(fe.Fe13_wl[ihi][ilo],"Fe13",&ipnt);
01129                         lindst( fe.Fe13_emiss[ihi][ilo] , (realnum)(fe.Fe13_wl[ihi][ilo]+0.5) , "Fe13",ipnt,'c',true ,
01130                                 " Fe 13 emission");
01131                 }
01132         }
01133 
01134         /* Fe 14 optically thin in line 344 */
01135         eff = 60.3/(60.3 + 0.23*dense.cdsqte/4.);
01136         r14 = dense.xIonDense[ipIRON][14-1]*fe.pfe14*rfield.flux[fe.ipfe10-1]*
01137           3.75e-12*eff/(1. + 24./.63*
01138           TauLines[ipT347].Emis->Pesc)*esc_PRD_1side(TauLines[ipT347].Emis->TauIn,1e-4);
01139         PntForLine(5303.,"Fe14",&ipnt);
01140 
01141         lindst(CoolHeavy.c5303+r14,5303,"Fe14",ipnt,'i',true ," total emission in Fe 14 5304");
01142 
01143         linadd(CoolHeavy.c5303,5303,"Coll",'c' ,
01144                 " Fe 14 5304  contribution from collisional excitation ");
01145 
01146         linadd(r14,5303,"Pump",'r' ,
01147                 " Fe 14 5304  continuum fluorescense ");
01148 
01151         c14 = 0.;
01152         linadd(c14,5303," 347",'c' ,
01153                 " collisional excitation of E1 line ");
01154 
01155         PutLine(&TauLines[ipFe18975],
01156                 " Fe 18 974.86A ");
01157 
01158 
01159         PntForLine(7047.,"Fe19",&ipnt);
01160         lindst(CoolHeavy.c7082,7047,"Fe19",ipnt,'c',true,
01161                 "  O-like Fe19, 3P ground term, 7046.72A vacuum wl, 1328.90A  >>chng 01 aug 10, updated wavelengths   Fe 19 7047 '85 " );
01162 
01163 
01164         PntForLine(1328.,"Fe19",&ipnt);
01165         lindst(CoolHeavy.c1328,1329,"Fe19",ipnt,'c',true,
01166                 " Fe 19 1329" );
01167 
01168 
01169         PntForLine(592.,"Fe19",&ipnt);
01170         lindst(CoolHeavy.c592,592,"Fe19",ipnt,'c',true,
01171                 "  Fe 19 from loulergue et al '85" );
01172 
01173 
01174         PntForLine(1118.,"Fe19",&ipnt);
01175         lindst(CoolHeavy.c1118,1118,"Fe19",ipnt,'c',true ,
01176                 " Fe 19 from loulergue et al '85");
01177 
01178         PutLine(&TauLines[ipTFe13],
01179                 "   next two 3p ground state lines, collisional excitation ");
01180 
01181         PutLine(&TauLines[ipTFe23],
01182                 "  collisional excitation ");
01183 
01184         PutLine(&TauLines[ipTFe20_578],
01185                 "  Fe20 721.40A, 578");
01186         PutLine(&TauLines[ipTFe20_721],
01187                 " ");
01188 
01189         linadd(
01190                 TauLines[ipFe22_247].Emis->xIntensity+
01191                 TauLines[ipFe22_217].Emis->xIntensity+
01192                 TauLines[ipFe22_348].Emis->xIntensity+
01193                 TauLines[ipFe22_292].Emis->xIntensity+
01194                 TauLines[ipFe22_253].Emis->xIntensity,
01195                 260,"TOTL",'i',
01196                 " Fe 22 845.6A  total intensity of Fe22, all lines in the multiplet " );
01197         PutLine(&TauLines[ipFe22_247],
01198                 "Fe 22 247");
01199         PutLine(&TauLines[ipFe22_217],
01200                 "Fe 22 217");
01201         PutLine(&TauLines[ipFe22_348],
01202                 "Fe 22 348");
01203         PutLine(&TauLines[ipFe22_292],
01204                 "Fe 22 292");
01205         PutLine(&TauLines[ipFe22_253],
01206                 "Fe 22 253");
01207 
01208         /*  Fe 23 1909-like 262.6 */
01209         PntForLine(263.,"Fe23",&ipnt);
01210         lindst(CoolHeavy.c263,263,"Fe23",ipnt,'c',true,
01211                 "Fe 23 1909-like 262.6" );
01212 
01213 
01214         PutLine(&TauLines[ipT192],
01215                 " Fe 24 only 192 of 255, 192 Li seq doublet, collisional excitation ");
01216 
01217         PutLine(&TauLines[ipT255],
01218                 "  Fe 24 255 of 255, 192 Li seq doublet, collisional excitation ");
01219 
01220         PutLine(&TauLines[ipT11],
01221                 "  Fe 24 Li seq 2s3p collisional excitation ");
01222 
01223         /* recombination Ka */
01224         if( dense.lgElmtOn[ipIRON] )
01225         {
01226                 /* these lines added to outlin in metdif - following must be false
01227                  * fela = xLyaHeavy(nelem,nelem)*dense.xIonDense(nelem,nelem+1) */
01228                 fela = Transitions[ipH_LIKE][ipIRON][ipH2p][ipH1s].Emis->xIntensity;
01229         }
01230         else
01231         {
01232                 fela = 0.;
01233         }
01234 
01235         /* >>chng 02 jan 14, add grain fe to this sum */
01236         /* total intensity of K-alpha line */
01237         /*linadd((fe.fekcld+fe.fegrain)*1.03e-8+(fe.fekhot+fela)*1.11e-8,2,"FeKa",'i' );*/
01238         if( dense.lgElmtOn[ipIRON] )
01239         {
01240                 lindst((fe.fekcld+fe.fegrain)*1.03e-8+(fe.fekhot+fela)*1.11e-8,1.78f,"FeKa",
01241                         Transitions[ipH_LIKE][ipIRON][ipH2p][ipH1s].ipCont,'i',false,
01242                            "total intensity of K-alpha line" );
01243         }
01244 
01245         linadd(fela*1.11e-8,2,"FeLr",'i' ,
01246                 " recombination from fully stripped ion ");
01247 
01248         /* >>chng 03 aug 14, label changed from TotH to AugH to be like rest total hot iron Ka; */
01249         linadd((fe.fekhot+fela)*1.11e-8,2,"AugH",'i' ,
01250                 "  Auger hot iron, assumes case b for H and He-like ");
01251 
01252         linadd(fe.fekcld*1.03e-8,2,"AugC",'i',
01253                 " Auger production of cold iron, less than or 17 times ionized " );
01254 
01255         linadd(fe.fegrain*1.03e-8,2,"AugG",'i' ,
01256                 " grain production of cold iron ");
01257 
01258         PutLine(&TauLines[ipCo11527],
01259                 "  [Co XI] 5168. A ");
01260 
01261         PutLine(&TauLines[ipNi1_7m],
01262                 "  nickel  [Ni I] 7m ");
01263 
01264         /* nickel*/
01265 
01266 
01267         PutLine(&TauLines[ipNi1_11m],
01268                 "  [Ni I] 11m ");
01269 
01270         PutLine(&TauLines[ipxNi1242],
01271                 " [Ni XII] 4230.8 A ");
01272 
01273         /* copper */
01274 
01275         /* zinc */
01276          PutLine(&TauLines[ipZn04363],
01277                  "zinc iv 3.625 microns, cs and A just made up ");
01278 
01279         if( trace.lgTrace )
01280         {
01281                 fprintf( ioQQQ, "   lines_lv1_k_zn returns\n" );
01282         }
01283         return;
01284 }
Generated on Mon Feb 16 12:01:26 2009 for cloudy by  doxygen 1.4.7