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 /*CoolSodi compute sodium cooling */ 00004 #include "cddefines.h" 00005 #include "coolheavy.h" 00006 #include "taulines.h" 00007 #include "dense.h" 00008 #include "phycon.h" 00009 #include "lines_service.h" 00010 #include "atoms.h" 00011 #include "cooling.h" 00012 00013 void CoolSodi(void) 00014 { 00015 double a12, 00016 a13, 00017 a23, 00018 cs, 00019 cs13, 00020 cs23, 00021 p3; 00022 realnum p2; 00023 00024 DEBUG_ENTRY( "CoolSodi()" ); 00025 00026 /* NaD Na D Na1 Na 1 lines 00027 * transition is 3s-3p, ^2P^o - ^2P^o 00028 * change 26 Feb 96 cs to Verner's value */ 00029 cs = 2.12e-2*(phycon.te70/phycon.te02); 00030 PutCS(cs,&TauLines[ipT5895]); 00031 atom_level2(&TauLines[ipT5895]); 00032 00033 /* [NaIII] 7.319 microns 00034 * cs 00035 * >>refer na3 cs Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */ 00036 cs = MIN2(0.40,0.198*phycon.te05*phycon.te01); 00037 cs = MAX2(0.35,cs); 00038 PutCS(cs,&TauLines[ipfsNa373]); 00039 atom_level2(&TauLines[ipfsNa373]); 00040 00041 /* collision data from 00042 * >>refer na4 cs Butler, K., & Zeippen, C.J. 1994, A&AS, 108, 1 00043 * [NaIV] 9.048 microns */ 00044 PutCS(0.802,&TauLines[ipfsNa490]); 00045 /* [NaIV] 21.29 mic */ 00046 PutCS(0.273,&TauLines[ipfsNa421]); 00047 /* following from pradhan review */ 00048 PutCS(.111,&TauDummy); 00049 atom_level3(&TauLines[ipfsNa490],&TauLines[ipfsNa421],&TauDummy); 00050 00051 /* >>chng 97 mar 19, added NaV lines 00052 * NaV lines, collision data from 00053 * >>refer na5 cs Mendoza, C. 1982, in Planetary Nebulae, IAU Symp No. 103, 00054 * >>refercon ed by D.R. Flower, (D. Reidel: Holland), 143 00055 * A's from 00056 * >>refer na5 as Kaufman, V., & Sugar, J. 1986, J Phys Chem Ref Dat, 15, 321 00057 * >>chng 97 jul 25, cs had been 1, changed to Me */ 00058 cs = 0.919; 00059 cs13 = 0.359; 00060 cs23 = 1.915; 00061 a12 = 8.16e-3; 00062 a13 = 3.81; 00063 a23 = 2.68; 00064 /* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2) */ 00065 p3 = atom_pop3(4.,10.,6.,cs,cs13,cs23,a12,a13,a23,6.96e4,3.58e4,&p2, 00066 dense.xIonDense[ipSODIUM][4],0.,0.,0.); 00067 CoolHeavy.c1365 = p3*a13*1.46e-11; 00068 CoolHeavy.c4017 = p3*a23*4.95e-12; 00069 CoolHeavy.c2067 = p2*a12*9.63e-12; 00070 CoolAdd("Na 5",1365,CoolHeavy.c1365); 00071 CoolAdd("Na 5",4017,CoolHeavy.c4017); 00072 CoolAdd("Na 5",2067,CoolHeavy.c2067); 00073 00074 /* [Na VI] 14.32 mic, 8.62 mic */ 00075 cs = MIN2(0.77,2.346/(phycon.te10*phycon.te02*phycon.te001)); 00076 PutCS(cs,&TauLines[ipxNa6143]); 00077 cs = MIN2(2.15,6.934/(phycon.te10*phycon.te03/phycon.te001/ 00078 phycon.te001)); 00079 PutCS(cs,&TauLines[ipxNa6862]); 00080 cs = MIN2(0.53,1.518/(phycon.te10*phycon.te01*phycon.te003* 00081 phycon.te003)); 00082 PutCS(cs,&TauDummy); 00083 00084 atom_level3(&TauLines[ipxNa6143],&TauLines[ipxNa6862],&TauDummy); 00085 00086 /* [Na VI] UV lines, 2971.9, 2872.7 doublet, 2578.9, 1356.6 00087 * POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2) */ 00088 cs = MIN2(0.2876,2.603e-3/(phycon.sqrte/phycon.te10)); 00089 cs = MAX2(0.1,cs); 00090 00091 CoolHeavy.c2569 = atom_pop3(9.,5.,1.,1.38,0.173,cs,1.68,16.9,5.27,5.01e4, 00092 5.60e4,&p2,dense.xIonDense[ipSODIUM][5],0.,0.,0.)*5.27*6.70e-12; 00093 00094 CoolHeavy.c1357 = CoolHeavy.c2569*(16.9/5.27)*(2972./1357.); 00095 CoolHeavy.c2972 = p2*1.68*6.70e-12; 00096 00097 CoolAdd("Na 6",2569,CoolHeavy.c2569); 00098 CoolAdd("Na 6",1357,CoolHeavy.c1357); 00099 CoolAdd("Na 6",2972,CoolHeavy.c2972); 00100 00101 /* [Na VII] 4.675 microns, no CS (interpolated), A NIST 00102 * collision strength interpolated from 00103 * >>refer na7 cs Lennon, D.J. Burke, V.M. 1994, A&AS, 103, 273 */ 00104 PutCS(1.5,&TauLines[ipxNa0746]); 00105 atom_level2(&TauLines[ipxNa0746]); 00106 return; 00107 }