cloudy: /home66/gary/public_html/cloudy/c08_branch/source/cool

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 /*CoolScan compute scandium cooling */
00004 #include "cddefines.h"
00005 #include "taulines.h"
00006 #include "coolheavy.h"
00007 #include "dense.h"
00008 #include "lines_service.h"
00009 #include "atoms.h"
00010 #include "cooling.h"
00011 
00012 void CoolScan(void)
00013 {
00014         double a21, 
00015           a31, 
00016           a32, 
00017           g1, 
00018           g2, 
00019           g3, 
00020           p3;
00021         realnum p2;
00022 
00023         DEBUG_ENTRY( "CoolScan()" );
00024 
00025         /* Sc scandium cooling
00026          *
00027          * Sc II
00028          * these are 3 lines estimated by Jim Kingdon
00029          * a's are bad, collision strengths just one */
00030         a21 = 0.004;
00031         a31 = 6.95e-3;
00032         a32 = 5.0e-4;
00033         g1 = 15.;
00034         g2 = 5.;
00035         g3 = 21.;
00036         /* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2)
00037          * energies are in kelvin */
00038         p3 = atom_pop3(g1,g2,g3,g1,g2,g3,a21,a31,a32,3504.,3407.,&p2,
00039                 dense.xIonDense[ipSCANDIUM][1], 0.,0.,0.);
00040 
00041         /* 2.08 microns */
00042         CoolHeavy.Sc22p08m = p3*a31*9.56e-13;
00043 
00044         /* 4.1 micron */
00045         CoolHeavy.Sc24p1m = p2*a21*4.85e-13;
00046 
00047         /* 4.22 micron */
00048         CoolHeavy.Sc24p2m = p3*a32*4.71e-13;
00049 
00050         CoolAdd("Sc 2",41,CoolHeavy.Sc24p1m);
00051         CoolAdd("Sc 2",42,CoolHeavy.Sc24p2m);
00052         CoolAdd("Sc 2",21,CoolHeavy.Sc22p08m);
00053 
00054         /* Sc III 3933
00055          * POPEXC( O12,g1,g2,A21,excit,abund); result already*a21
00056          * [Sc III] 3933, multiplet average */
00057         g1 = 10.;
00058         g2 = 2.;
00059         CoolHeavy.Sc33936 = atom_pop2(g1,g1,g2,0.03,3.66e4,
00060                 dense.xIonDense[ipSCANDIUM][2])* 5.06e-12;
00061         CoolAdd("Sc 3",3933,CoolHeavy.Sc33936);
00062 
00063         /* [Sc V] 2.31 mic
00064          * Y(ik) from 
00065          * >>refer      sc5     cs      Pelan, J., & Berrington, K.A. 1995, A&A Suppl, 110, 209 */
00066         PutCS(6.00,&TauLines[ipSc05231]);
00067         atom_level2(&TauLines[ipSc05231]);
00068 
00069         /* Sc VI */
00070         a21 = 4.94;
00071         a31 = 49.24;
00072         a32 = 4.3;
00073         g1 = 9.;
00074         g2 = 5.;
00075         g3 = 1.;
00076         /* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2)
00077          * energies are in kelvin */
00078         p3 = atom_pop3(g1,g2,g3,g1,g2,g3,a21,a31,a32,28464.,40045.,&p2,
00079                 dense.xIonDense[ipSCANDIUM][5], 0.,0.,0.);
00080         /* 2100 ang - 3=>1 */
00081         CoolHeavy.Sc42100 = p3*a31*9.47e-12;
00082         /* 5054 - 2=>1 */
00083         CoolHeavy.Sc45058 = p2*a21*3.93e-12;
00084         /* 3592 3=>2 */
00085         CoolHeavy.Sc43595 = p3*a32*5.53e-12;
00086         CoolAdd("Sc 6",2100,CoolHeavy.Sc42100);
00087         CoolAdd("Sc 6",5054,CoolHeavy.Sc45058);
00088         CoolAdd("Sc 6",3592,CoolHeavy.Sc43595);
00089 
00090         /* [Sc 13] 2637.97A, cs from 
00091          * >>refer      sc13    cs      Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */
00092         PutCS(0.182,&TauLines[ipSc13264]);
00093         atom_level2(&TauLines[ipSc13264]);
00094 
00095         return;
00096 }
/home66/gary/public_html/cloudy/c08_branch/source/cool_scan.cpp
