cloudy trunk

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,21,"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,41,"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,42,"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,12,"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,45,"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,42,"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 
00285         PntForLine(5828.,"Cr 3",&ipnt);
00286         lindst(CoolHeavy.Cr3l21,5828,"Cr 3",ipnt,'c',true,
00287                 " [CrIII] multiplet blend at 5828A" );
00288 
00289         PntForLine(7267.,"Cr 4",&ipnt);
00290         lindst(CoolHeavy.Cr4l21,7267,"Cr 4",ipnt,'c',true,
00291                 " [CrIV] 2 - 1 multiplet blend at 7272" );
00292 
00293 
00294         PntForLine(6801.,"Cr 4",&ipnt);
00295         lindst(CoolHeavy.Cr4l31,6801,"Cr 4",ipnt,'c',true,
00296                 " [CrIV] 3 - 1 multiplet blend at 6806" );
00297 
00298 
00299         PntForLine(7979.,"Cr 5",&ipnt);
00300         lindst(CoolHeavy.Cr5l21,7979,"Cr 5",ipnt,'c',true,
00301                 " [CrV] 2 - 1 multiplet blend at 7985" );
00302 
00303         PntForLine(6577.,"Cr 5",&ipnt);
00304         lindst(CoolHeavy.Cr5l31,6577,"Cr 5",ipnt,'c',true,
00305                 "  [CrV] 3 - 1 multiplet blend at 6582" );
00306 
00307 
00308         PntForLine(3.75e4,"Cr 5",&ipnt);
00309         lindst(CoolHeavy.Cr5l32,37,"Cr 5",ipnt,'c',true,
00310                 " [CrV] 3 - 2 multiplet blend at 3.75 microns " );
00311 
00312         PutLine(&TauLines[ipCr08101],
00313                 "  [Cr VIII] 1.01 microns ");
00314 
00315         PutLine(&TauLines[ipCr16141],
00316                 "  [Cr 16] 1410.60 ");
00317 
00318         PutLine(&TauLines[ipxMn0979],
00319                 " [Mn IX] 7968.5 A ");
00320 
00321         PutLine(&TauLines[ipxMn1712],
00322                 " [Mn 17] 1169.59 ");
00323 
00324         /* bob Rubin's UV line
00325          * f2 = dense.xIonDense(26,4)*sexp(50 764./te)*0.45*cdsqte/6.*7.01e-12
00326          * call linadd( f2 , 2837 , 'BobR' , 'i')
00327          * f2 = dense.xIonDense(26,4)*sexp(55 989./te)*0.384*cdsqte/6.*7.74e-12
00328          * call linadd( f2 , 2568 , 'BobR' , 'i') */
00329 
00330         /* iron */
00331 
00332         PutLine(&TauLines[ipFe1_24m],
00333                 "  Fe 1 24m ");
00334 
00335         PutLine(&TauLines[ipFe1_35m],
00336                 "  Fe 1 35m ");
00337 
00338         PutLine(&TauLines[ipFe1_54m],
00339                 "  Fe 1 54m ");
00340 
00341         PutLine(&TauLines[ipFe1_111m],
00342                 "  Fe 1 111m ");
00343 
00344         PutLine(&TauLines[ipFeI3884],
00345                 "  Fe 1 3884 ");
00346 
00347         PutLine(&TauLines[ipFeI3729],
00348                 "  Fe 1 3729 ");
00349 
00350         PutLine(&TauLines[ipFeI3457],
00351                 "  Fe 1 3457 ");
00352 
00353         PutLine(&TauLines[ipFeI3021],
00354                 "  Fe 1 3021 ");
00355 
00356         PutLine(&TauLines[ipFeI2966],
00357                 "  Fe 1 2966 ");
00358 
00359         linadd(MAX2(0.,FeII.Fe2_large_cool+FeII.Fe2_UVsimp_cool),0,"Fe2c",'c' ,
00360                 "total of all Fe 2, both simple UV and large atom together ");
00361 
00362         linadd(MAX2(0.,-FeII.Fe2_large_cool-FeII.Fe2_UVsimp_cool),0,"Fe2h",'h' ,
00363                 "total of all Fe 2 heating, both simple UV and large atom together ");
00364 
00365         linadd(FeII.for7,4300,"Fe 2",'i' ,
00366                 " Fe 2 forbidden 2-1 transition from Netzer's atom ");
00367 
00368         PutLine(&TauLines[ipTuv3],
00369                 " 2400 ");
00370         PutLine(&TauLines[ipTr48],
00371                 " 6200 ");
00372         PutLine(&TauLines[ipTFe16],
00373                 " 1080 ");
00374         PutLine(&TauLines[ipTFe26],
00375                 " 1500  ");
00376         PutLine(&TauLines[ipTFe34],
00377                 " 11500  ");
00378         PutLine(&TauLines[ipTFe35],
00379                 " 2500 ");
00380         PutLine(&TauLines[ipTFe46],
00381                 " 2300  ");
00382         PutLine(&TauLines[ipTFe56],
00383                 " 8900  ");
00384 
00385         /* >>chng 03 nov 02, changed wl in lindst calls to be wl in microns,
00386          * this will change wl in printout */
00387         /* results of FeIILevelPops Fe 2 atom, 
00388          * not all are printed, first is */
00389         /* total cooling due to 16 level atom 
00390          * >>chng 05 dec 04, rm this since no longer meaningful 
00391         linadd(FeII.Fe2_16levl_cool,0,"Fe 2",'c' );*/
00392 
00393 
00394         PntForLine(16640.,"Fe 2",&ipnt);
00395         lindst(FeII.fe21308,16640,"Fe 2",ipnt,'i',true,
00396                 " Fe 2 1.664 microns 8-13 " );
00397 
00398 
00399         PntForLine(15990.,"Fe 2",&ipnt);
00400         lindst(FeII.fe21207,15990,"Fe 2",ipnt,'i',true,
00401                 " Fe 2 1.599 microns 7-12" );
00402 
00403 
00404         PntForLine(15340.,"Fe 2",&ipnt);
00405         lindst(FeII.fe21106,15340,"Fe 2",ipnt,'i',true ,
00406                 " Fe 2 1.534 microns 6-11 ");
00407 
00408 
00409         PntForLine(16440.,"Fe 2",&ipnt);
00410         lindst(FeII.fe21006,16440,"Fe 2",ipnt,'i',true,
00411                 " Fe 2 1.644 microns 6-10 " );
00412 
00413 
00414         PntForLine(12790.,"Fe 2",&ipnt);
00415         lindst(FeII.fe21204,12790,"Fe 2",ipnt,'i',true,
00416                 " Fe 2 1.279 microns 12-4" );
00417 
00418 
00419         PntForLine(12950.,"Fe 2",&ipnt);
00420         lindst(FeII.fe21103,12950,"Fe 2",ipnt,'i',true,
00421                 " Fe 2 1.295 microns 11-3 " );
00422 
00423 
00424         PntForLine(13280.,"Fe 2",&ipnt);
00425         lindst(FeII.fe21104,13280,"Fe 2",ipnt,'i',true,
00426                 "  Fe 2 1.328 microns 11-4" );
00427 
00428 
00429         PntForLine(12570.,"Fe 2",&ipnt);
00430         lindst(FeII.fe21001,12570,"Fe 2",ipnt,'i',true,
00431                 " Fe 2 1.257 microns 10-1 " );
00432 
00433 
00434         PntForLine(13210.,"Fe 2",&ipnt);
00435         lindst(FeII.fe21002,13210,"Fe 2",ipnt,'i',true,
00436                 " Fe 2 1.321 microns 10-2 " );
00437 
00438 
00439         PntForLine(259880.,"Fe 2",&ipnt);
00440         lindst(FeII.fe20201,259880,"Fe 2",ipnt,'i',true,
00441                 "  Fe 2 25.988 microns 2-1" );
00442 
00443         PntForLine(353480.,"Fe 2",&ipnt);
00444         lindst(FeII.fe20302,353480,"Fe 2",ipnt,'i',true,
00445                 " Fe 2 35.348 microns 3-2" );
00446 
00447 
00448         PntForLine(179360.,"Fe 2",&ipnt);
00449         lindst(FeII.fe20706,179360,"Fe 2",ipnt,'i',true,
00450                 "  Fe 2 17.936 microns 7-6, label is 178 to be unique " );
00451 
00452 
00453         PntForLine(245180.,"Fe 2",&ipnt);
00454         lindst(FeII.fe20807,245180,"Fe 2",ipnt,'i',true,
00455                 " Fe 2 24.518 microns 8-7" );
00456 
00457         PntForLine(357760.,"Fe 2",&ipnt);
00458         lindst(FeII.fe20908,357760,"Fe 2",ipnt,'i',true,
00459                 " Fe 2 35.776 microns 9-8" );
00460 
00461         PntForLine(18100.,"Fe 2",&ipnt);
00462         lindst(FeII.fe21007,18100,"Fe 2",ipnt,'i',true,
00463                 " Fe 2 1.810 microns 10-7" );
00464 
00465 
00466         PntForLine(16770.,"Fe 2",&ipnt);
00467         lindst(FeII.fe21107,16770,"Fe 2",ipnt,'i',true,
00468                 "  Fe 2 1.677 microns 11-7 " );
00469 
00470 
00471         PntForLine(18000.,"Fe 2",&ipnt);
00472         lindst(FeII.fe21108,18000,"Fe 2",ipnt,'i',true,
00473                 "  Fe 2 1.800 microns 11-8" );
00474 
00475 
00476         PntForLine(17120.,"Fe 2",&ipnt);
00477         lindst(FeII.fe21208,17120,"Fe 2",ipnt,'i',true,
00478                 " Fe 2 1.712 microns 12-8 " );
00479 
00480 
00481         PntForLine(17980.,"Fe 2",&ipnt);
00482         lindst(FeII.fe21209,17980,"Fe 2",ipnt,'i',true,
00483                 " Fe 2 1.798 microns 12-9" );
00484 
00485 
00486         PntForLine(229020.,"Fe 2",&ipnt);
00487         lindst(FeII.fe21110,229020,"Fe 2",ipnt,'i',true,
00488                 " Fe 2 22.902 microns 11-10 " );
00489 
00490 
00491         PntForLine(346600.,"Fe 2",&ipnt);
00492         lindst(FeII.fe21211,346600,"Fe 2",ipnt,'i',true,
00493                 "  Fe 2 34.660 microns 12-11" );
00494 
00495         /* >>chng 97 mar 14, added following 4 lines after consultation with Verner
00496          * these are typically within 5-9% of their large model atom.  These are
00497          * predicted with both large and small atom */
00498 
00499 
00500         PntForLine(8617.,"Fe 2",&ipnt);
00501         lindst(FeII.fe21406,8617,"Fe 2",ipnt,'i',true,
00502                 "  Fe 2 8617A 14-06," );
00503 
00504 
00505         PntForLine(8892.,"Fe 2",&ipnt);
00506         lindst(FeII.fe21507,8892,"Fe 2",ipnt,'i',true,
00507                 " Fe 2 8892A 15-07" );
00508 
00509 
00510         PntForLine(9227.,"Fe 2",&ipnt);
00511         lindst(FeII.fe21508,9227,"Fe 2",ipnt,'i',true,
00512                 " Fe 2 9227A 15-08 " );
00513 
00514 
00515         PntForLine(9268.,"Fe 2",&ipnt);
00516         lindst(FeII.fe21609,9268,"Fe 2",ipnt,'i',true,
00517                 " Fe 2 9268A 16-09" );
00518 
00519         /* following lines are only non-zero when large atom is enabled.  When
00520          * running with small FeII atoms these are set zero in routine  */
00521 
00522 
00523         PntForLine(5159.,"Fe 2",&ipnt);
00524         lindst(FeII.fe25to6,5159,"Fe 2",ipnt,'i',true,
00525                 " Fe 2 5158.81 25-6" );
00526 
00527 
00528         PntForLine(5262.,"Fe 2",&ipnt);
00529         lindst(FeII.fe27to7,5262,"Fe 2",ipnt,'i',true,
00530                 " Fe 2 5261.61 27-7" );
00531 
00532         PntForLine(5334.,"Fe 2",&ipnt);
00533         lindst(FeII.fe28to8,5334,"Fe 2",ipnt,'i',true,
00534                 " Fe 2 5333.65 28-8 " );
00535 
00536 
00537         PntForLine(5376.,"Fe 2",&ipnt);
00538         lindst(FeII.fe29to9,5376,"Fe 2",ipnt,'i',true,
00539                 " Fe 2 5376.47 29-9 " );
00540 
00541         PntForLine(4815.,"Fe 2",&ipnt);
00542         lindst(FeII.fe32to6,4815,"Fe 2",ipnt,'i',true,
00543                 "  Fe 2 4814.55 32-6 " );
00544 
00545         PntForLine(4905.,"Fe 2",&ipnt);
00546         lindst(FeII.fe33to7,4905,"Fe 2",ipnt,'i',true,
00547                 " Fe 2 4905.35 33-7" );
00548 
00549 
00550         PntForLine(4347.,"Fe 2",&ipnt);
00551         lindst(FeII.fe37to7,4347,"Fe 2",ipnt,'i',true,
00552                 " Fe 2 4346.85 37-7 " );
00553 
00554 
00555         PntForLine(4353.,"Fe 2",&ipnt);
00556         lindst(FeII.fe39to8,4353,"Fe 2",ipnt,'i',true,
00557                 " Fe 2 4352.78 39-8 " );
00558 
00559 
00560         PntForLine(4372.,"Fe 2",&ipnt);
00561         lindst(FeII.fe40to9,4372,"Fe 2",ipnt,'i',true,
00562                 " Fe 2 4372.43 40-9 " );
00563 
00564 
00565         PntForLine(4244.,"Fe 2",&ipnt);
00566         lindst(FeII.fe37to6,4244,"Fe 2",ipnt,'i',true,
00567                 " Fe 2 4243.98 37-6" );
00568 
00569 
00570         PntForLine(4277.,"Fe 2",&ipnt);
00571         lindst(FeII.fe39to7,4277,"Fe 2",ipnt,'i',true,
00572                 " Fe 2 4276.83 39-7" );
00573 
00574 
00575         PntForLine(4320.,"Fe 2",&ipnt);
00576         lindst(FeII.fe40to8,4320,"Fe 2",ipnt,'i',true,
00577                 "  Fe 2 4319.62 40-8" );
00578 
00579 
00580         PntForLine(4358.,"Fe 2",&ipnt);
00581         lindst(FeII.fe41to9,4358,"Fe 2",ipnt,'i',true,
00582                 " Fe 2 4358.37 41-9 " );
00583 
00584 
00585         PntForLine(4177.,"Fe 2",&ipnt);
00586         lindst(FeII.fe39to6,4177,"Fe 2",ipnt,'i',true,
00587                 " Fe 2 4177.21 39-6 " );
00588 
00589 
00590         PntForLine(4245.,"Fe 2",&ipnt);
00591         lindst(FeII.fe40to7,4245,"Fe 2",ipnt,'i',true,
00592                 " Fe 2 4244.81 40-7 " );
00593 
00594 
00595         PntForLine(4306.,"Fe 2",&ipnt);
00596         lindst(FeII.fe41to8,4306,"Fe 2",ipnt,'i',true,
00597                 " Fe 2 4305.90 41-8" );
00598 
00599         /* >>chng 99 jun 18 new lines added as per Katya email, so all orion lines are in */
00600 
00601         PntForLine(4114.,"Fe 2",&ipnt);
00602         lindst(FeII.fe42to6,4114,"Fe 2",ipnt,'i',true,
00603                 "  Fe 2 4114.48 42-6 " );
00604 
00605 
00606         PntForLine(4179.,"Fe 2",&ipnt);
00607         lindst(FeII.fe43to7,4179,"Fe 2",ipnt,'i',true,
00608                 " Fe 2 4178.95  7-      43" );
00609 
00610 
00611         PntForLine(4211.,"Fe 2",&ipnt);
00612         lindst(FeII.fe42to7,4211,"Fe 2",ipnt,'i',true,
00613                 " Fe 2 4211.10  7-      42 " );
00614 
00615 
00616         PntForLine(4314.,"Fe 2",&ipnt);
00617         lindst(FeII.fe80to28,4314,"Fe 2",ipnt,'i',true,
00618                 "  Fe 2 4314.29 28-     80" );
00619 
00620         PntForLine(4359.,"Fe 2",&ipnt);
00621         lindst(FeII.fe36to2,4359,"Fe 2",ipnt,'i',true,
00622                 " Fe 2 4359.34  2-      36" );
00623 
00624 
00625         PntForLine(4414.,"Fe 2",&ipnt);
00626         lindst(FeII.fe36to3,4414,"Fe 2",ipnt,'i',true,
00627                 "  Fe 2 4413.78 3-      36 " );
00628 
00629 
00630         PntForLine(4416.,"Fe 2",&ipnt);
00631         lindst(FeII.fe32to1,4416,"Fe 2",ipnt,'i',true ,
00632                 " Fe 2 4416.27  1-      32");
00633 
00634 
00635         PntForLine(4458.,"Fe 2",&ipnt);
00636         lindst(FeII.fe33to2,4458,"Fe 2",ipnt,'i',true,
00637                 " Fe 2 4457.95  2-      33 " );
00638 
00639 
00640         PntForLine(4475.,"Fe 2",&ipnt);
00641         lindst(FeII.fe36to5,4475,"Fe 2",ipnt,'i',true,
00642                 " Fe 2 4474.91  5-      36 " );
00643 
00644 
00645         PntForLine(4493.,"Fe 2",&ipnt);
00646         lindst(FeII.fe32to2,4493,"Fe 2",ipnt,'i',true,
00647                 "  Fe 2 4492.64 2-      32" );
00648 
00649 
00650         PntForLine(4515.,"Fe 2",&ipnt);
00651         lindst(FeII.fe33to3,4515,"Fe 2",ipnt,'i',true,
00652                 "  Fe 2 4514.90 3-      33" );
00653 
00654 
00655         PntForLine(4728.,"Fe 2",&ipnt);
00656         lindst(FeII.fe30to3,4728,"Fe 2",ipnt,'i',true,
00657                 "  Fe 2 4728.07 3-      30" );
00658 
00659 
00660         PntForLine(4775.,"Fe 2",&ipnt);
00661         lindst(FeII.fe33to6,4775,"Fe 2",ipnt,'i',true,
00662                 " Fe 2 4774.74  6-      33" );
00663 
00664         PntForLine(4890.,"Fe 2",&ipnt);
00665         lindst(FeII.fe24to2,4890,"Fe 2",ipnt,'i',true,
00666                 "  Fe 2 4889.63 2-      24" );
00667 
00668 
00669         PntForLine(4947.,"Fe 2",&ipnt);
00670         lindst(FeII.fe32to7,4947,"Fe 2",ipnt,'i',true,
00671                 "  Fe 2 4947.37 7-      32" );
00672 
00673 
00674         PntForLine(4951.,"Fe 2",&ipnt);
00675         lindst(FeII.fe35to8,4951,"Fe 2",ipnt,'i',true,
00676                 " Fe 2 4950.74  8-      35" );
00677 
00678 
00679         PntForLine(4973.,"Fe 2",&ipnt);
00680         lindst(FeII.fe34to8,4973,"Fe 2",ipnt,'i',true,
00681                 " Fe 2 4973.39  8-      34" );
00682 
00683 
00684         PntForLine(5112.,"Fe 2",&ipnt);
00685         lindst(FeII.fe27to6,5112,"Fe 2",ipnt,'i',true,
00686                 " Fe 2 5111.63  6-      27" );
00687 
00688 
00689         PntForLine(5220.,"Fe 2",&ipnt);
00690         lindst(FeII.fe28to7,5220,"Fe 2",ipnt,'i',true,
00691                 " Fe 2 5220.06  7-      28" );
00692 
00693 
00694         PntForLine(5269.,"Fe 2",&ipnt);
00695         lindst(FeII.fe30to8,5269,"Fe 2",ipnt,'i',true,
00696                 " Fe 2 5268.880 8-      30 " );
00697 
00698 
00699         PntForLine(5273.,"Fe 2",&ipnt);
00700         lindst(FeII.fe24to6,5273,"Fe 2",ipnt,'i',true,
00701                 " Fe 2 5273.35  6-      24" );
00702 
00703         PntForLine(5297.,"Fe 2",&ipnt);
00704         lindst(FeII.fe29to8,5297,"Fe 2",ipnt,'i',true,
00705                 " Fe 2 5296.840 8-      29 " );
00706 
00707 
00708         PntForLine(5433.,"Fe 2",&ipnt);
00709         lindst(FeII.fe24to7,5433,"Fe 2",ipnt,'i',true,
00710                 " Fe 2 5433.13  7-      24" );
00711 
00712 
00713         PntForLine(5527.,"Fe 2",&ipnt);
00714         lindst(FeII.fe22to7,5527,"Fe 2",ipnt,'i',true,
00715                 " Fe 2 5527.340 7-      22" );
00716 
00717 
00718         PntForLine(5747.,"Fe 2",&ipnt);
00719         lindst(FeII.fe38to11,5747,"Fe 2",ipnt,'i',true,
00720                 "  Fe 2 5746.96 11-     38" );
00721 
00722 
00723         PntForLine(6440.,"Fe 2",&ipnt);
00724         lindst(FeII.fe19to8,6440,"Fe 2",ipnt,'i',true,
00725                 " Fe 2 6440.400 8-      19" );
00726 
00727         PntForLine(7155.,"Fe 2",&ipnt);
00728         lindst(FeII.fe17to6,7155,"Fe 2",ipnt,'i',true,
00729                 " Fe 2 7155.140 6-      17 " );
00730 
00731 
00732         PntForLine(7172.,"Fe 2",&ipnt);
00733         lindst(FeII.fe18to7,7172,"Fe 2",ipnt,'i',true,
00734                 " Fe 2 7171.980 7-      18" );
00735 
00736 
00737         PntForLine(7388.,"Fe 2",&ipnt);
00738         lindst(FeII.fe18to8,7388,"Fe 2",ipnt,'i',true,
00739                 " Fe 2 7388.160 8-      18" );
00740 
00741 
00742         PntForLine(7453.,"Fe 2",&ipnt);
00743         lindst(FeII.fe17to7,7453,"Fe 2",ipnt,'i',true,
00744                 " Fe 2 7452.540 7-      17" );
00745 
00746         /* option to save all intensities predicted by large FeII atom,
00747          * code is in FeIILevelPops */
00748         FeIIAddLines();
00749 
00750         /* some bands from the large model Fe II atom
00751          * only done when full large atom is turned on */
00752         if( FeII.lgFeIILargeOn )
00753         {
00754                 for( i=0; i < nFeIIBands; i++ )
00755                 {
00756                         /* [i][0] is center wavelength, [i][1] and [i][2] are upper and
00757                          * lower bounds in Angstroms.  These are set in FeIIZero 
00758                          * units are erg s-1 cm-3 */
00759                         eff = FeIISumBand(FeII_Bands[i][1],FeII_Bands[i][2]);
00760 
00761                         linadd(eff,FeII_Bands[i][0],"Fe2b",'i' ,
00762                                 " total Fe II emission in Fe II bands, as defined in bands_Fe2.dat ");
00763                 }
00765                 if( LineSave.ipass > 0 )
00766                 {
00767                         for( i=0; i < nFeIIConBins; i++ )
00768                         {
00769                                 /* [i][0] is total intensity in cell, [i][1] and [i][2] are lower and
00770                                  * upper bounds in Angstroms.  these are set in FeIIZero *
00771                                  * find total emission from large FeII atom, integrated over band */
00772                                 FeII_Cont[i][0] += 
00773                                         (realnum)(FeIISumBand(FeII_Cont[i][1],FeII_Cont[i][2])*radius.dVeff);
00774                                 /*fprintf(ioQQQ,"DEBUG feii\t%li\t%.2e\n", i, FeII_Cont[i][0]);*/
00775                         }
00776                 }
00777                 else
00778                 {
00779                         for( i=0; i < nFeIIConBins; i++ )
00780                         {
00781                                 /* initialize arrays */
00782                                 FeII_Cont[i][0] = 0.;
00783                         }
00784                 }
00785         }
00786         else
00787         {
00788                 /* must still add the bands in case of grid where FeII is turned off and on */
00789                 for( i=0; i < nFeIIBands; i++ )
00790                 {
00791 
00792                         linadd(0., FeII_Bands[i][0] ,"Fe2b",'i' ,
00793                                 " total Fe II emission in Fe II bands, as defined in bands_Fe2.dat ");
00794                 }
00795         }
00796 
00797         PutLine(&TauLines[ipT191],
00798                 "  anomalous Fe 2 transition at 1787, RMT 191");
00799 
00800         linadd(fe.Fe3CoolTot,0,"Fe3c",'c' ,
00801                 " chng 05 dec 16, FeIII code created by Kevin Blagrave  Fe3c 0 - total cooling due to 14-level Fe 3 atom ");
00802 
00803         /* Fe 3 14-level atom 
00804          * following from print statements within loop */
00805         /* Fe 3 22.92m from Blagrave 14-level atom */
00806         /* Fe 3 13.53m from Blagrave 14-level atom */
00807         /* Fe 3 33.03m from Blagrave 14-level atom */
00808         /* Fe 3 10.72m from Blagrave 14-level atom */
00809         /* Fe 3 20.15m from Blagrave 14-level atom */
00810         /* Fe 3 51.67m from Blagrave 14-level atom */
00811         /* Fe 3 9.732m from Blagrave 14-level atom */
00812         /* Fe 3 16.91m from Blagrave 14-level atom */
00813         /* Fe 3 34.66m from Blagrave 14-level atom */
00814         /* Fe 3 105.3m from Blagrave 14-level atom */
00815         /* Fe 3  5152A from Blagrave 14-level atom */
00816         /* Fe 3  5271A from Blagrave 14-level atom */
00817         /* Fe 3  5356A from Blagrave 14-level atom */
00818         /* Fe 3  5412A from Blagrave 14-level atom */
00819         /* Fe 3  5440A from Blagrave 14-level atom */
00820         /* Fe 3  4986A from Blagrave 14-level atom */
00821         /* Fe 3  5097A from Blagrave 14-level atom */
00822         /* Fe 3  5177A from Blagrave 14-level atom */
00823         /* Fe 3  5230A from Blagrave 14-level atom */
00824         /* Fe 3  5256A from Blagrave 14-level atom */
00825         /* Fe 3 15.47m from Blagrave 14-level atom */
00826         /* Fe 3  4925A from Blagrave 14-level atom */
00827         /* Fe 3  5033A from Blagrave 14-level atom */
00828         /* Fe 3  5111A from Blagrave 14-level atom */
00829         /* Fe 3  5162A from Blagrave 14-level atom */
00830         /* Fe 3  5188A from Blagrave 14-level atom */
00831         /* Fe 3 11.16m from Blagrave 14-level atom */
00832         /* Fe 3 40.04m from Blagrave 14-level atom */
00833         /* Fe 3  4881A from Blagrave 14-level atom */
00834         /* Fe 3  4988A from Blagrave 14-level atom */
00835         /* Fe 3  5064A from Blagrave 14-level atom */
00836         /* Fe 3  5114A from Blagrave 14-level atom */
00837         /* Fe 3  5139A from Blagrave 14-level atom */
00838         /* Fe 3 9.282m from Blagrave 14-level atom */
00839         /* Fe 3 23.21m from Blagrave 14-level atom */
00840         /* Fe 3 55.20m from Blagrave 14-level atom */
00841         /* Fe 3  4833A from Blagrave 14-level atom */
00842         /* Fe 3  4937A from Blagrave 14-level atom */
00843         /* Fe 3  5012A from Blagrave 14-level atom */
00844         /* Fe 3  5061A from Blagrave 14-level atom */
00845         /* Fe 3  5085A from Blagrave 14-level atom */
00846         /* Fe 3 7.789m from Blagrave 14-level atom */
00847         /* Fe 3 15.69m from Blagrave 14-level atom */
00848         /* Fe 3 25.79m from Blagrave 14-level atom */
00849         /* Fe 3 48.41m from Blagrave 14-level atom */
00850         /* Fe 3  4714A from Blagrave 14-level atom */
00851         /* Fe 3  4813A from Blagrave 14-level atom */
00852         /* Fe 3  4884A from Blagrave 14-level atom */
00853         /* Fe 3  4931A from Blagrave 14-level atom */
00854         /* Fe 3  4954A from Blagrave 14-level atom */
00855         /* Fe 3 5.543m from Blagrave 14-level atom */
00856         /* Fe 3 8.638m from Blagrave 14-level atom */
00857         /* Fe 3 11.01m from Blagrave 14-level atom */
00858         /* Fe 3 13.76m from Blagrave 14-level atom */
00859         /* Fe 3 19.22m from Blagrave 14-level atom */
00860         /* Fe 3  4659A from Blagrave 14-level atom */
00861         /* Fe 3  4755A from Blagrave 14-level atom */
00862         /* Fe 3  4825A from Blagrave 14-level atom */
00863         /* Fe 3  4870A from Blagrave 14-level atom */
00864         /* Fe 3  4893A from Blagrave 14-level atom */
00865         /* Fe 3 4.859m from Blagrave 14-level atom */
00866         /* Fe 3 7.085m from Blagrave 14-level atom */
00867         /* Fe 3 8.608m from Blagrave 14-level atom */
00868         /* Fe 3 10.20m from Blagrave 14-level atom */
00869         /* Fe 3 12.92m from Blagrave 14-level atom */
00870         /* Fe 3 39.41m from Blagrave 14-level atom */
00871         /* Fe 3  4608A from Blagrave 14-level atom */
00872         /* Fe 3  4702A from Blagrave 14-level atom */
00873         /* Fe 3  4770A from Blagrave 14-level atom */
00874         /* Fe 3  4814A from Blagrave 14-level atom */
00875         /* Fe 3  4836A from Blagrave 14-level atom */
00876         /* Fe 3 4.356m from Blagrave 14-level atom */
00877         /* Fe 3 6.063m from Blagrave 14-level atom */
00878         /* Fe 3 7.146m from Blagrave 14-level atom */
00879         /* Fe 3 8.208m from Blagrave 14-level atom */
00880         /* Fe 3 9.884m from Blagrave 14-level atom */
00881         /* Fe 3 20.34m from Blagrave 14-level atom */
00882         /* Fe 3 42.06m from Blagrave 14-level atom */
00883         /* Fe 3  4574A from Blagrave 14-level atom */
00884         /* Fe 3  4668A from Blagrave 14-level atom */
00885         /* Fe 3  4734A from Blagrave 14-level atom */
00886         /* Fe 3  4778A from Blagrave 14-level atom */
00887         /* Fe 3  4800A from Blagrave 14-level atom */
00888         /* Fe 3 4.077m from Blagrave 14-level atom */
00889         /* Fe 3 5.535m from Blagrave 14-level atom */
00890         /* Fe 3 6.423m from Blagrave 14-level atom */
00891         /* Fe 3 7.269m from Blagrave 14-level atom */
00892         /* Fe 3 8.554m from Blagrave 14-level atom */
00893         /* Fe 3 15.41m from Blagrave 14-level atom */
00894         /* Fe 3 25.31m from Blagrave 14-level atom */
00895         /* Fe 3 63.56m from Blagrave 14-level atom */
00896         for( ihi=1; ihi<NLFE3; ++ihi )
00897         {
00898                 for( ilo=0; ilo<ihi; ++ilo )
00899                 {
00900                         /* emission in these lines */
00901                         PntForLine(fe.Fe3_wl[ihi][ilo],"Fe 3",&ipnt);
00902 #                       if 0
00903                         fprintf( ioQQQ,"\t/* FeIII ");
00904                         prt_wl( ioQQQ , (realnum)(fe.Fe3_wl[ihi][ilo]+0.5) );
00905                         fprintf( ioQQQ," from Blagrave 14-level atom */\n ");
00906 #                       endif
00907                         lindst( fe.Fe3_emiss[ihi][ilo] , (realnum)(fe.Fe3_wl[ihi][ilo]+0.5) , "Fe 3",ipnt,'c',true,
00908                                 " " );
00909                 }
00910         }
00911 
00912         /*>>chng 05 dec 18, following are now in the above */
00913         /* sum of 3p and 3g states together */
00914         /*      linadd(CoolHeavy.c5270,0,"Fe 3",'c' ); */
00915 
00916         /* Fe 3 5270, predictions from Garstang et al 78
00917         PntForLine(5270.,"Fe 3",&ipnt);
00918         lindst(CoolHeavy.c5270*0.2090,5270,"Fe 3",ipnt,'c',true );*/
00919 
00920         /* Fe 3 5270, predictions from Garstang et al 78 
00921         PntForLine(4658.,"Fe 3",&ipnt);
00922         lindst(CoolHeavy.c5270*0.3667,4658,"Fe 3",ipnt,'c',true ); */
00923 
00924         PutLine(&TauLines[ipT1122]," Fe 3 1122 entire multiplet");
00925 
00926         linadd(fe.Fe4CoolTot,0,"Fe4c",'c',
00927                 " Fe4c 0 - total cooling due to 12-level Fe 4 atom " );
00928 
00929 
00930         PntForLine(3096.,"Fe 4",&ipnt);
00931         lindst(fe.fe40401,3096,"Fe 4",ipnt,'i',true,
00932                 " Fe 4 3096.A, 4-1 and 5-1 transitions together"  );
00933 
00934 
00935         PntForLine(2836.,"Fe 4",&ipnt);
00936         lindst(fe.fe42836,2836,"Fe 4",ipnt,'i',true,
00937                 " Fe 4 2835.7A, 6-1 transition, 4P5/2 - 6S5/2 "  );
00938 
00939 
00940         PntForLine(2829.,"Fe 4",&ipnt);
00941         lindst(fe.fe42829,2829,"Fe 4",ipnt,'i',true,
00942                 "   Fe 4 2829.4A, 7-1 transition, 4P3/2 - 6S5/2"  );
00943 
00944 
00945         PntForLine(2567.,"Fe 4",&ipnt);
00946         lindst(fe.fe42567,2567,"Fe 4",ipnt,'i',true,
00947                 "  Fe 4 2567.6+ 2567.4. 11-1 and 12-1 transitions"  );
00948 
00949 
00950         PntForLine(2.774e4,"Fe 4",&ipnt);
00951         lindst(fe.fe41207,277,"Fe 4",ipnt,'i',true,
00952                 " Fe 4 2.774 microns 12-7 transition "  );
00953 
00954 
00955         PntForLine(2.714e4,"Fe 4",&ipnt);
00956         lindst(fe.fe41206,271,"Fe 4",ipnt,'i',true,
00957                 " Fe 4 2.714 microns 12-6 transition "  );
00958 
00959 
00960         PntForLine(2.716e4,"Fe 4",&ipnt);
00961         lindst(fe.fe41106,272,"Fe 4",ipnt,'i',true,
00962                 " Fe 4 2.716 microns 11-6 transition"  );
00963 
00964 
00965         PntForLine(2.806e4,"Fe 4",&ipnt);
00966         lindst(fe.fe41007,281,"Fe 4",ipnt,'i',true,
00967                 " Fe 4 2.806 microns 10-7 transition "  );
00968 
00969 
00970         PntForLine(2.865e4,"Fe 4",&ipnt);
00971         lindst(fe.fe41008,287,"Fe 4",ipnt,'i',true ,
00972                 "  Fe 4 2.865 microns 10-8 transition");
00973 
00974 
00975         PntForLine(2.836e4,"Fe 4",&ipnt);
00976         lindst(fe.fe40906,284,"Fe 4",ipnt,'i',true,
00977                 " Fe 4 2.836 microns 9-6 transition" );
00978 
00979 
00980         PntForLine(3892.,"Fe 5",&ipnt);
00981         lindst(CoolHeavy.c3892,3892,"Fe 5",ipnt,'c',true,
00982                 " Fe 5  3892+3839" );
00983 
00984         linadd(CoolHeavy.c5177,0,"Fe 6",'c' ,
00985                 " all of 2G lines together first ");
00986 
00987 
00988         PntForLine(5177.,"Fe 6",&ipnt);
00989         lindst(CoolHeavy.c5177*0.354,5177,"Fe 6",ipnt,'c',true,
00990                 " Fe 6 5177, approximate correct " );
00991 
00992         linadd(fe.Fe7CoolTot,0,"Fe7c",'c' ,
00993                 " Fe7c 0 - total cooling due to n-level Fe 7 atom ");
00994 
00995         /* >>chng 04 nov 04, move to multi-level system */
00996         for( ilo=0; ilo<NLFE7-1; ++ilo )
00997         {
00998                 /* must not do 1-0 or 2-1, which are transferred lines */
00999                 for( ihi=MAX2(3,ilo+1); ihi<NLFE7; ++ihi )
01000                 {
01001 
01002                         PntForLine(fe.Fe7_wl[ihi][ilo],"Fe 7",&ipnt);
01003                         lindst( fe.Fe7_emiss[ihi][ilo] , (realnum)(fe.Fe7_wl[ihi][ilo]+0.5) , "Fe 7",ipnt,'c',true,
01004                                 " emission in these lines" );
01005                 }
01006         }
01007 #       if 0
01008         PntForLine(5721.,"Fe 7",&ipnt);
01009         lindst( fe.Fe7_5721 , 5721 , "Fe 7",ipnt,'c',true,
01010                 " " );
01011 
01012         PntForLine(6601.,"Fe 7",&ipnt);
01013         lindst( fe.Fe7_6601 , 6601 , "Fe 7",ipnt,'c',true,
01014                 " " );
01015 
01016         PntForLine(3760.,"Fe 7",&ipnt);
01017         lindst( fe.Fe7_3760 , 3760 , "Fe 7",ipnt,'c',true,
01018                 " " );
01019 
01020         PntForLine(3588.,"Fe 7",&ipnt);
01021         lindst( fe.Fe7_3588 , 3588 , "Fe 7",ipnt,'c',true,
01022                 " " );
01023 #       endif
01024 
01025         PutLine(&TauLines[ipFe0795],
01026                 "   [Fe 7] 9.51 micron ");
01027 
01028         PutLine(&TauLines[ipFe0778],
01029                 "  [Fe 7] 7.81 micron ");
01030 
01031         /* [Fe 7] 6087 
01032         PntForLine(6087.,"Fe 7",&ipnt);
01033         lindst(CoolHeavy.c6087,6087,"Fe 7",ipnt,'c',true );*/
01034 
01035         /* [Fe 7] 5722 
01036         PntForLine(5722.,"Fe 7",&ipnt);
01037         lindst(CoolHeavy.Fe5722,5722,"Fe 7",ipnt,'c',true );*/
01038 
01039         PutLine(&TauLines[ipT245],
01040                 "   Be-seq lines ");
01041 
01042         PntForLine(242.,"Fe 7",&ipnt);
01043         lindst(CoolHeavy.c242,242,"Fe 7",ipnt,'c',true,
01044                 "  Fe 9 242 j=1 slower decay");
01045 
01046         PutLine(&TauLines[ipT352],
01047                 " the E1 transition that can pump [Fe X] ");
01048 
01049         /* optically thin Fe X pumping */
01050         eff = 69.4/(69.4 + 0.27*dense.cdsqte);
01051 
01052         /* coll excitation of 352 which decays to excited state of 6374
01053          * assumes 17/56 (ratio of A's) go to excited state */
01054         c10 = TauLines[ipT352].Emis->phots*eff*.01676*TauLines[ipT352].EnergyErg*352/6374.;
01055 
01056         /* Fe 10 and Fe 14 from Mason 75
01057          * total (coll, pumped)
01058          * call linadd( C6374+C10 , 6374 , 'Fe10','i') */
01059         PutExtra( c10 );
01060 
01061 
01062         PutLine(&TauLines[ipFe106375],
01063                 " [Fe 10] 6375, collisions with pumping too ");
01064 
01065         /* collisional contribution
01066          * call linadd( C6374 , 6374 , 'Coll','c')
01067          * collisions of E1 line, plus pumped by continuum fluorescence
01068          * call linadd( C10 , 6374 , ' 352','c')
01069          * Fe XI 7892, 6.08 micron 
01070         PutLine(&TauLines[ipTFe07]);*/
01071 
01083         /* Fe 11 2649 collisional excitation
01084         PntForLine(2649.,"Fe11",&ipnt);
01085         lindst(CoolHeavy.c2649,2649,"Fe11",ipnt,'c',true ); */
01086 
01087         /*  Fe 11 1467 collisional excitation 
01088         PntForLine(1467.,"Fe11",&ipnt);
01089         lindst(CoolHeavy.c1467,1467,"Fe11",ipnt,'c',true );*/
01090 
01091         linadd(fe.Fe11CoolTot,0,"Fe11",'c' ,
01092                 " >>chng 05 dec 18, add Fe 11  Fe11 0 - total cooling due to 5-level Fe 11 atom ");
01093 
01094         /* Fe 11 5-level atom */
01095         for( ihi=1; ihi<NLFE11; ++ihi )
01096         {
01097                 for( ilo=0; ilo<ihi; ++ilo )
01098                 {
01099                         PntForLine(fe.Fe11_wl[ihi][ilo],"Fe11",&ipnt);
01100                         lindst( fe.Fe11_emiss[ihi][ilo] , (realnum)(fe.Fe11_wl[ihi][ilo]+0.5) , "Fe11",ipnt,'c',true,
01101                                 " emission in these lines" );
01102                 }
01103         }
01104 
01105         PntForLine(1242.,"Fe12",&ipnt);
01106         lindst(CoolHeavy.c1242,1242,"Fe12",ipnt,'c',true ,
01107                 " Fe 12, 1242, 1349 together, collisional excitation");
01108 
01109         PntForLine(2170.,"Fe12",&ipnt);
01110         lindst(CoolHeavy.c2170,2170,"Fe12",ipnt,'c',true ,
01111                 " Fe 12, 2170, 2406 together, collisional excitation");
01112 
01113 
01114         PntForLine(2568.,"Fe12",&ipnt);
01115         lindst(CoolHeavy.c2568,2568,"Fe12",ipnt,'c',true,
01116                 "  Fe12 2904, 2567, 3567, 3073 together, collisional excitation" );
01117 
01118         /* >>chng 05 dec 18, add Fe 13  */
01119         /* Fe13 0 - total cooling due to 5-level Fe 13 atom */
01120         linadd(fe.Fe13CoolTot,0,"Fe13",'c' ,
01121                    "total cooling due to Fe 13 model atom ");
01122 
01123         /* Fe 13 5-level atom */
01124         for( ihi=1; ihi<NLFE13; ++ihi )
01125         {
01126                 for( ilo=0; ilo<ihi; ++ilo )
01127                 {
01128 
01129                         PntForLine(fe.Fe13_wl[ihi][ilo],"Fe13",&ipnt);
01130                         lindst( fe.Fe13_emiss[ihi][ilo] , (realnum)(fe.Fe13_wl[ihi][ilo]+0.5) , "Fe13",ipnt,'c',true ,
01131                                 " Fe 13 emission");
01132                 }
01133         }
01134 
01135         /* Fe 14 optically thin in line 344 */
01136         eff = 60.3/(60.3 + 0.23*dense.cdsqte/4.);
01137         r14 = dense.xIonDense[ipIRON][14-1]*fe.pfe14*rfield.flux[fe.ipfe10-1]*
01138           3.75e-12*eff/(1. + 24./.63*
01139           TauLines[ipT347].Emis->Pesc)*esc_PRD_1side(TauLines[ipT347].Emis->TauIn,1e-4);
01140         PntForLine(5303.,"Fe14",&ipnt);
01141 
01142         lindst(CoolHeavy.c5303+r14,5303,"Fe14",ipnt,'i',true ," total emission in Fe 14 5304");
01143 
01144         linadd(CoolHeavy.c5303,5303,"Coll",'c' ,
01145                 " Fe 14 5304  contribution from collisional excitation ");
01146 
01147         linadd(r14,5303,"Pump",'r' ,
01148                 " Fe 14 5304  continuum fluorescense ");
01149 
01152         c14 = 0.;
01153         linadd(c14,5303," 347",'c' ,
01154                 " collisional excitation of E1 line ");
01155 
01156         PutLine(&TauLines[ipFe18975],
01157                 " Fe 18 974.86A ");
01158 
01159 
01160         PntForLine(7047.,"Fe19",&ipnt);
01161         lindst(CoolHeavy.c7082,7047,"Fe19",ipnt,'c',true,
01162                 "  O-like Fe19, 3P ground term, 7046.72A vacuum wl, 1328.90A  >>chng 01 aug 10, updated wavelengths   Fe 19 7047 '85 " );
01163 
01164 
01165         PntForLine(1328.,"Fe19",&ipnt);
01166         lindst(CoolHeavy.c1328,1329,"Fe19",ipnt,'c',true,
01167                 " Fe 19 1329" );
01168 
01169 
01170         PntForLine(592.,"Fe19",&ipnt);
01171         lindst(CoolHeavy.c592,592,"Fe19",ipnt,'c',true,
01172                 "  Fe 19 from loulergue et al '85" );
01173 
01174 
01175         PntForLine(1118.,"Fe19",&ipnt);
01176         lindst(CoolHeavy.c1118,1118,"Fe19",ipnt,'c',true ,
01177                 " Fe 19 from loulergue et al '85");
01178 
01179         PutLine(&TauLines[ipTFe13],
01180                 "   next two 3p ground state lines, collisional excitation ");
01181 
01182         PutLine(&TauLines[ipTFe23],
01183                 "  collisional excitation ");
01184 
01185         PutLine(&TauLines[ipTFe20_578],
01186                 "  Fe20 721.40A, 578");
01187         PutLine(&TauLines[ipTFe20_721],
01188                 " ");
01189 
01190         linadd(
01191                 TauLines[ipFe22_247].Emis->xIntensity+
01192                 TauLines[ipFe22_217].Emis->xIntensity+
01193                 TauLines[ipFe22_348].Emis->xIntensity+
01194                 TauLines[ipFe22_292].Emis->xIntensity+
01195                 TauLines[ipFe22_253].Emis->xIntensity,
01196                 260,"TOTL",'i',
01197                 " Fe 22 845.6A  total intensity of Fe22, all lines in the multiplet " );
01198         PutLine(&TauLines[ipFe22_247],
01199                 "Fe 22 247");
01200         PutLine(&TauLines[ipFe22_217],
01201                 "Fe 22 217");
01202         PutLine(&TauLines[ipFe22_348],
01203                 "Fe 22 348");
01204         PutLine(&TauLines[ipFe22_292],
01205                 "Fe 22 292");
01206         PutLine(&TauLines[ipFe22_253],
01207                 "Fe 22 253");
01208 
01209         /*  Fe 23 1909-like 262.6 */
01210         PntForLine(263.,"Fe23",&ipnt);
01211         lindst(CoolHeavy.c263,263,"Fe23",ipnt,'c',true,
01212                 "Fe 23 1909-like 262.6" );
01213 
01214 
01215         PutLine(&TauLines[ipT192],
01216                 " Fe 24 only 192 of 255, 192 Li seq doublet, collisional excitation ");
01217 
01218         PutLine(&TauLines[ipT255],
01219                 "  Fe 24 255 of 255, 192 Li seq doublet, collisional excitation ");
01220 
01221         PutLine(&TauLines[ipT11],
01222                 "  Fe 24 Li seq 2s3p collisional excitation ");
01223 
01224         /* recombination Ka */
01225         if( dense.lgElmtOn[ipIRON] )
01226         {
01227                 /* these lines added to outlin in metdif - following must be false
01228                  * fela = xLyaHeavy(nelem,nelem)*dense.xIonDense(nelem,nelem+1) */
01229                 fela = Transitions[ipH_LIKE][ipIRON][ipH2p][ipH1s].Emis->xIntensity;
01230         }
01231         else
01232         {
01233                 fela = 0.;
01234         }
01235 
01236         /* >>chng 02 jan 14, add grain fe to this sum */
01237         /* total intensity of K-alpha line */
01238         /*linadd((fe.fekcld+fe.fegrain)*1.03e-8+(fe.fekhot+fela)*1.11e-8,2,"FeKa",'i' );*/
01239         if( dense.lgElmtOn[ipIRON] )
01240         {
01241                 lindst((fe.fekcld+fe.fegrain)*1.03e-8+(fe.fekhot+fela)*1.11e-8,2,"FeKa",
01242                         Transitions[ipH_LIKE][ipIRON][ipH2p][ipH1s].ipCont,'i',false,
01243                            "total intensity of K-alpha line" );
01244         }
01245 
01246         linadd(fela*1.11e-8,2,"FeLr",'i' ,
01247                 " recombination from fully stripped ion ");
01248 
01249         /* >>chng 03 aug 14, label changed from TotH to AugH to be like rest total hot iron Ka; */
01250         linadd((fe.fekhot+fela)*1.11e-8,2,"AugH",'i' ,
01251                 "  Auger hot iron, assumes case b for H and He-like ");
01252 
01253         linadd(fe.fekcld*1.03e-8,2,"AugC",'i',
01254                 " Auger production of cold iron, less than or 17 times ionized " );
01255 
01256         linadd(fe.fegrain*1.03e-8,2,"AugG",'i' ,
01257                 " grain production of cold iron ");
01258 
01259         PutLine(&TauLines[ipCo11527],
01260                 "  [Co XI] 5168. A ");
01261 
01262         PutLine(&TauLines[ipNi1_7m],
01263                 "  nickel  [Ni I] 7m ");
01264 
01265         /* nickel*/
01266 
01267 
01268         PutLine(&TauLines[ipNi1_11m],
01269                 "  [Ni I] 11m ");
01270 
01271         PutLine(&TauLines[ipxNi1242],
01272                 " [Ni XII] 4230.8 A ");
01273 
01274         /* copper */
01275 
01276         /* zinc */
01277          PutLine(&TauLines[ipZn04363],
01278                  "zinc iv 3.625 microns, cs and A just made up ");
01279 
01280         if( trace.lgTrace )
01281         {
01282                 fprintf( ioQQQ, "   lines_lv1_k_zn returns\n" );
01283         }
01284         return;
01285 }
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Friends Defines

Generated for cloudy by doxygen 1.7.3