cloudy  trunk
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
cool_nick.cpp
Go to the documentation of this file.
1 /* This file is part of Cloudy and is copyright (C)1978-2008 by Gary J. Ferland and
2  * others. For conditions of distribution and use see copyright notice in license.txt */
3 /*CoolNick compute nickel cooling */
4 #include "cddefines.h"
5 #include "taulines.h"
6 #include "phycon.h"
7 #include "lines_service.h"
8 #include "dense.h"
9 #include "atoms.h"
10 #include "cooling.h"
11 
12 void CoolNick(void)
13 {
14  double cs,
15  tused;
16  realnum rate;
17 
18  DEBUG_ENTRY( "CoolNick()" );
19 
20  /*>>refer Ni1 cs Hollenbach, D. & McKee, C.F. 1989, ApJ, 342, 306 */
21  /* rates are said to be ok over range 30 - 3000K */
22  tused = MAX2( 30. , phycon.te );
23  tused = MIN2( 3000. , phycon.te );
24  tused /= 100.;
25 
26  /* the 7.5 micron line */
27  /* >>chng 03 nov 15, add these lines */
28  rate = (realnum)(1.2e-7 * dense.eden +
29  /*8.0e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0]) / dense.eden);*/
30  /* >>chng 05 jul 05, eden to cdsqte */
31  8.0e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0] );
32  LineConvRate2CS( &TauLines[ipNi1_7m] , rate );
33 
34  /* the 11.3 micron line */
35  rate = (realnum)(9.3e-8 * dense.eden +
36  /* >>chng 05 jul 05, eden to cdsqte */
37  /*5.3e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0]) / dense.eden);*/
38  5.3e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0] );
40 
41  rate = (realnum)(1.2e-7 * dense.eden +
42  /* >>chng 05 jul 05, eden to cdsqte */
43  /*6.9e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0]) / dense.eden);*/
44  6.9e-10*pow(tused, 0.17 )*dense.xIonDense[ipHYDROGEN][0] );
46  LineConvRate2CS( &TauDummy , rate );
47  /* this says that line is a dummy, not real one */
48  TauDummy.Hi->g = 0.;
49 
50  atom_level3(&TauLines[ipNi1_7m],&TauLines[ipNi1_11m],&TauDummy);
51 
52  /* [Ni XII] 4230.8A
53  * Y(ik) from
54  * >>refer ni12 cs Pelan, J., & Berrington, K.A. 1995, A&A Suppl, 110, 209
55  * call PutCS( 1.68 , xNi1242 )
56  * >>chng 98 jan 09, update collision strength to
57  * >>refer ni12 cs Mathews, A., Ramsbottom, C.A., Bell, K.L., & Keenan, F.P., 1998,
58  * >>refercon ApJ 492, 415 */
59  if( phycon.alogte < 5.4 )
60  {
61  cs = MAX2(0.6,0.6+0.8182*(phycon.alogte-3.2));
62  }
63  else
64  {
65  cs = MAX2(0.6,2.4-1.5*(phycon.alogte-5.4));
66  }
67  PutCS(cs,&TauLines[ipxNi1242]);
68  atom_level2(&TauLines[ipxNi1242]);
69  return;
70 }

Generated for cloudy by doxygen 1.8.3.1