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_sodi.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
/*CoolSodi compute sodium cooling */
4
#include "
cddefines.h
"
5
#include "
coolheavy.h
"
6
#include "
taulines.h
"
7
#include "
dense.h
"
8
#include "
phycon.h
"
9
#include "
lines_service.h
"
10
#include "
atoms.h
"
11
#include "
cooling.h
"
12
13
void
CoolSodi
(
void
)
14
{
15
double
a12,
16
a13,
17
a23,
18
cs,
19
cs13,
20
cs23,
21
p3;
22
realnum
p2;
23
24
DEBUG_ENTRY
(
"CoolSodi()"
);
25
26
/* NaD Na D Na1 Na 1 lines
27
* transition is 3s-3p, ^2P^o - ^2P^o
28
* change 26 Feb 96 cs to Verner's value */
29
cs = 2.12e-2*(
phycon
.
te70
/
phycon
.
te02
);
30
PutCS
(cs,&
TauLines
[
ipT5895
]);
31
atom_level2
(&
TauLines
[ipT5895]);
32
33
/* [NaIII] 7.319 microns
34
* cs
35
* >>refer na3 cs Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */
36
cs =
MIN2
(0.40,0.198*
phycon
.
te05
*
phycon
.
te01
);
37
cs =
MAX2
(0.35,cs);
38
PutCS
(cs,&
TauLines
[
ipfsNa373
]);
39
atom_level2
(&
TauLines
[ipfsNa373]);
40
41
/* collision data from
42
* >>refer na4 cs Butler, K., & Zeippen, C.J. 1994, A&AS, 108, 1
43
* [NaIV] 9.048 microns */
44
PutCS
(0.802,&
TauLines
[
ipfsNa490
]);
45
/* [NaIV] 21.29 mic */
46
PutCS
(0.273,&
TauLines
[
ipfsNa421
]);
47
/* following from pradhan review */
48
PutCS
(.111,&
TauDummy
);
49
atom_level3
(&
TauLines
[ipfsNa490],&
TauLines
[ipfsNa421],&
TauDummy
);
50
51
/* >>chng 97 mar 19, added NaV lines
52
* NaV lines, collision data from
53
* >>refer na5 cs Mendoza, C. 1982, in Planetary Nebulae, IAU Symp No. 103,
54
* >>refercon ed by D.R. Flower, (D. Reidel: Holland), 143
55
* A's from
56
* >>refer na5 as Kaufman, V., & Sugar, J. 1986, J Phys Chem Ref Dat, 15, 321
57
* >>chng 97 jul 25, cs had been 1, changed to Me */
58
cs = 0.919;
59
cs13 = 0.359;
60
cs23 = 1.915;
61
a12 = 8.16e-3;
62
a13 = 3.81;
63
a23 = 2.68;
64
/* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2) */
65
p3 =
atom_pop3
(4.,10.,6.,cs,cs13,cs23,a12,a13,a23,6.96e4,3.58e4,&p2,
66
dense
.
xIonDense
[
ipSODIUM
][4],0.,0.,0.);
67
CoolHeavy
.
c1365
= p3*a13*1.46e-11;
68
CoolHeavy
.
c4017
= p3*a23*4.95e-12;
69
CoolHeavy
.
c2067
= p2*a12*9.63e-12;
70
CoolAdd
(
"Na 5"
,1365,
CoolHeavy
.
c1365
);
71
CoolAdd
(
"Na 5"
,4017,
CoolHeavy
.
c4017
);
72
CoolAdd
(
"Na 5"
,2067,
CoolHeavy
.
c2067
);
73
74
/* [Na VI] 14.32 mic, 8.62 mic */
75
cs =
MIN2
(0.77,2.346/(
phycon
.
te10
*
phycon
.
te02
*
phycon
.
te001
));
76
PutCS
(cs,&
TauLines
[
ipxNa6143
]);
77
cs =
MIN2
(2.15,6.934/(
phycon
.
te10
*
phycon
.
te03
/
phycon
.
te001
/
78
phycon
.
te001
));
79
PutCS
(cs,&
TauLines
[
ipxNa6862
]);
80
cs =
MIN2
(0.53,1.518/(
phycon
.
te10
*
phycon
.
te01
*
phycon
.
te003
*
81
phycon
.
te003
));
82
PutCS
(cs,&
TauDummy
);
83
84
atom_level3
(&
TauLines
[ipxNa6143],&
TauLines
[ipxNa6862],&
TauDummy
);
85
86
/* [Na VI] UV lines, 2971.9, 2872.7 doublet, 2578.9, 1356.6
87
* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2) */
88
cs =
MIN2
(0.2876,2.603e-3/(
phycon
.
sqrte
/
phycon
.
te10
));
89
cs =
MAX2
(0.1,cs);
90
91
CoolHeavy
.
c2569
=
atom_pop3
(9.,5.,1.,1.38,0.173,cs,1.68,16.9,5.27,5.01e4,
92
5.60e4,&p2,
dense
.
xIonDense
[
ipSODIUM
][5],0.,0.,0.)*5.27*6.70e-12;
93
94
CoolHeavy
.
c1357
=
CoolHeavy
.
c2569
*(16.9/5.27)*(2972./1357.);
95
CoolHeavy
.
c2972
= p2*1.68*6.70e-12;
96
97
CoolAdd
(
"Na 6"
,2569,
CoolHeavy
.
c2569
);
98
CoolAdd
(
"Na 6"
,1357,
CoolHeavy
.
c1357
);
99
CoolAdd
(
"Na 6"
,2972,
CoolHeavy
.
c2972
);
100
101
/* [Na VII] 4.675 microns, no CS (interpolated), A NIST
102
* collision strength interpolated from
103
* >>refer na7 cs Lennon, D.J. Burke, V.M. 1994, A&AS, 103, 273 */
104
PutCS
(1.5,&
TauLines
[
ipxNa0746
]);
105
atom_level2
(&
TauLines
[ipxNa0746]);
106
return
;
107
}
Generated for cloudy by
1.8.4