cloudy
trunk
|
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 }