cloudy
trunk
Main Page
Related Pages
Data Structures
Files
File List
Globals
All
Data Structures
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Pages
source
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] );
39
LineConvRate2CS
( &
TauLines
[
ipNi1_11m
] , rate );
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] );
45
TauDummy
.
Hi
->
g
=
TauLines
[
ipNi1_11m
].
Hi
->
g
;
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
1.8.1.1