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 /*CoolVana compute vanadium 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 #include "phycon.h" 00012 00013 void CoolVana(void) 00014 { 00015 realnum a21, 00016 a31, 00017 a32, 00018 g1, 00019 g2, 00020 g3, 00021 p2, 00022 p3, 00023 cs; 00024 00025 DEBUG_ENTRY( "CoolVana()" ); 00026 00027 /* V Vanadium cooling - element 23 00028 * 00029 * V III 8823 */ 00030 a21 = 0.05f; 00031 a31 = 0.10f; 00032 a32 = 0.00; 00033 g1 = 28.; 00034 g2 = 12.; 00035 g3 = 18.; 00036 00037 /* [V III] 8823, multiplet average */ 00038 p3 = (realnum)atom_pop3(g1,g2,g3,g1,g2,g3,a21,a31,a32,16303.,606.,&p2, 00039 dense.xIonDense[ipVANADIUM][2], 0.,0.,0.); 00040 00041 CoolHeavy.V38830 = p2*a21*2.25e-12; 00042 CoolHeavy.V38507 = p3*a31*2.34e-12; 00043 CoolAdd("V 3",8823,CoolHeavy.V38830); 00044 CoolAdd("V 3",8507,CoolHeavy.V38507); 00045 00046 /* V IV */ 00047 a21 = 0.054f; 00048 a31 = 0.039f; 00049 a32 = 0.007f; 00050 g1 = 21.; 00051 g2 = 5.; 00052 g3 = 9.; 00053 /* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2) 00054 * energies are in kelvin */ 00055 p3 = (realnum)atom_pop3(g1,g2,g3,g1,g2,g3,a21,a31,a32,15159.,3437.,&p2, 00056 dense.xIonDense[ipVANADIUM][3], 0.,0.,0.); 00057 /* 7735 ang - 3=>1 */ 00058 CoolHeavy.V47741 = p3*a31*2.57e-12; 00059 /* 9489 - 2=>1 */ 00060 CoolHeavy.V49496 = p2*a21*2.09e-12; 00061 /* 4.19 microns 3=>2 */ 00062 CoolHeavy.V44p2m = p3*a32*4.74e-13; 00063 CoolAdd("V 4",7735,CoolHeavy.V47741); 00064 CoolAdd("V 4",9489,CoolHeavy.V49496); 00065 CoolAdd("V 4",42,CoolHeavy.V44p2m); 00066 00067 /* [V VII] 1.3038 mic 00068 * Y(ik) from 00069 * >>refer v7 cs Pelan, J., & Berrington, K.A. 1995, A&A Suppl, 110, 209 */ 00070 PutCS(2.39,&TauLines[ipVa07130]); 00071 atom_level2(&TauLines[ipVa07130]); 00072 00073 /* [V 15] 1721.38, cs from 00074 * >>referold v15 cs Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */ 00075 /* >>refer v15 cs Berrington,K.A.,Saraph, H.E. & Tully, J.A. 1998, A&AS, 129, 161 */ 00076 /*>>chng 06 jul 19 Changes made-Humeshkar Nemala*/ 00077 if(phycon.te < 3.566E6) 00078 { 00079 cs = (realnum)(0.0149*phycon.te10*phycon.te05*phycon.te004*phycon.te0003); 00080 } 00081 else 00082 { 00083 cs = (realnum)(47.350653/((phycon.te40/phycon.te02)*phycon.te0002)); 00084 } 00085 00086 /*PutCS(0.10,&TauLines[ipVa15172]);*/ 00087 PutCS(cs,&TauLines[ipVa15172]); 00088 atom_level2(&TauLines[ipVa15172]); 00089 return; 00090 }