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 /*CoolPota compute potassium cooling */ 00004 #include "cddefines.h" 00005 #include "taulines.h" 00006 #include "phycon.h" 00007 #include "lines_service.h" 00008 #include "atoms.h" 00009 #include "cooling.h" 00010 00011 void CoolPota(void) 00012 { 00013 double cs; 00014 00015 DEBUG_ENTRY( "CoolPota()" ); 00016 00017 /* potasium lines 00018 * KI 7745 */ 00019 cs = 7.231e-4*phycon.te*phycon.te03*phycon.te02; 00020 PutCS(cs,&TauLines[ipKI7745]); 00021 atom_level2(&TauLines[ipKI7745]); 00022 00023 /* [K III] 4.62 microns 00024 * Y(ik) from 00025 * >>refer k3 cs Pelan, J., & Berrington, K.A. 1995, A&A Suppl, 110, 209 */ 00026 PutCS(2.2,&TauLines[ipxK03462]); 00027 atom_level2(&TauLines[ipxK03462]); 00028 00029 /* [KIV] 5.983, 15.39 mic, cs from 00030 * >>refer k4 cs Galavis, M.E., Mendoza, C., & Zeippen, C.J. 1995, A&AS, 111, 347 */ 00031 PutCS(4.3,&TauLines[ipxK04598]); 00032 PutCS(1.13,&TauLines[ipxK04154]); 00033 PutCS(1.3,&TauDummy); 00034 /* atom_level3( t10,t21,t20) */ 00035 atom_level3(&TauLines[ipxK04598],&TauLines[ipxK04154],&TauDummy); 00036 00037 /* [KVI] 8.823, 5.575 mic, cs from 00038 * >>refer k6 cs Galavis, M.E., Mendoza, C., & Zeippen, C.J. 1995, A&AS, 111, 347 */ 00039 cs = MIN2(1.505,0.274*phycon.te10*phycon.te05/phycon.te001/ 00040 phycon.te001); 00041 PutCS(cs,&TauLines[ipxK06882]); 00042 00043 cs = MIN2(4.632,1.909*phycon.te10/phycon.te003); 00044 cs = MAX2(4.0,cs); 00045 PutCS(cs,&TauLines[ipxK06557]); 00046 PutCS(1.2,&TauDummy); 00047 00048 atom_level3(&TauLines[ipxK06882],&TauLines[ipxK06557],&TauDummy); 00049 00050 /* [K VII] 3.189 microns cs from 00051 * >>refer k7 cs Saraph, H.E., & Storey, P.J. A&AS, 115, 151 */ 00052 PutCS(4.5,&TauLines[ipxK07319]); 00053 atom_level2(&TauLines[ipxK07319]); 00054 00055 /* K 11 4249.99A, cs from 00056 * >>refer k11 cs Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */ 00057 cs = MIN2(0.172,0.0109*phycon.te20*phycon.te02/ 00058 phycon.te001/phycon.te001); 00059 cs = MAX2(0.111,cs); 00060 PutCS(0.115,&TauLines[ipxK11425]); 00061 00062 atom_level2(&TauLines[ipxK11425]); 00063 00064 return; 00065 }