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 /*CoolTita compute titanium 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 CoolTita(void) 00013 { 00014 realnum a21, 00015 a31, 00016 a32, 00017 p2, 00018 p3; 00019 00020 DEBUG_ENTRY( "CoolTita()" ); 00021 00022 /* Ti Titanium cooling 00023 * 00024 * these are 3 lines estimated by Jim Kingdon 00025 * a's are bad, collision strengths just one */ 00026 a21 = 0.015f; 00027 a31 = 0.032f; 00028 a32 = 0.002f; 00029 00030 /* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2) 00031 * energies are in kelvin */ 00032 p3 = (realnum)atom_pop3(21.,5.,9.,21.,5.0,9.0,a21,a31,a32,11844.,3148.4,&p2, 00033 dense.xIonDense[ipTITANIUM][2],0.,0.,0.); 00034 00035 /* multiplet at roughly 9594 */ 00036 CoolHeavy.Ti3l31 = p3*a31*2.07e-12; 00037 /* multiplet at roughly 4.57 microns */ 00038 00039 CoolHeavy.Ti3l32 = p3*a32*4.35e-13; 00040 /* multiplet at roughly 1.21 microns */ 00041 00042 CoolHeavy.Ti3l21 = p2*a21*1.64e-12; 00043 CoolAdd("Ti 3",9594,CoolHeavy.Ti3l31); 00044 CoolAdd("Ti 3",4,CoolHeavy.Ti3l32); 00045 CoolAdd("Ti 3",1,CoolHeavy.Ti3l21); 00046 00047 /* [Ti VI] 1.7150 mic 00048 * Y(ik) from 00049 * >>refer ti6 cs Pelan, J., & Berrington, K.A. 1995, A&A Suppl, 110, 209 */ 00050 PutCS(3.48,&TauLines[ipTi06172]); 00051 atom_level2(&TauLines[ipTi06172]); 00052 00053 /* [Ti XIV] 2117.79, cs from 00054 * >>refer ti14 cs Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */ 00055 PutCS(0.23,&TauLines[ipTi14212]); 00056 atom_level2(&TauLines[ipTi14212]); 00057 return; 00058 }