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_scan.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
/*CoolScan compute scandium cooling */
4
#include "
cddefines.h
"
5
#include "
taulines.h
"
6
#include "
coolheavy.h
"
7
#include "
dense.h
"
8
#include "
lines_service.h
"
9
#include "
atoms.h
"
10
#include "
cooling.h
"
11
12
void
CoolScan
(
void
)
13
{
14
double
a21,
15
a31,
16
a32,
17
g1,
18
g2,
19
g3,
20
p3;
21
realnum
p2;
22
23
DEBUG_ENTRY
(
"CoolScan()"
);
24
25
/* Sc scandium cooling
26
*
27
* Sc II
28
* these are 3 lines estimated by Jim Kingdon
29
* a's are bad, collision strengths just one */
30
a21 = 0.004;
31
a31 = 6.95e-3;
32
a32 = 5.0e-4;
33
g1 = 15.;
34
g2 = 5.;
35
g3 = 21.;
36
/* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2)
37
* energies are in kelvin */
38
p3 =
atom_pop3
(g1,g2,g3,g1,g2,g3,a21,a31,a32,3504.,3407.,&p2,
39
dense
.
xIonDense
[
ipSCANDIUM
][1], 0.,0.,0.);
40
41
/* 2.08 microns */
42
CoolHeavy
.
Sc22p08m
= p3*a31*9.56e-13;
43
44
/* 4.1 micron */
45
CoolHeavy
.
Sc24p1m
= p2*a21*4.85e-13;
46
47
/* 4.22 micron */
48
CoolHeavy
.
Sc24p2m
= p3*a32*4.71e-13;
49
50
CoolAdd
(
"Sc 2"
,41,
CoolHeavy
.
Sc24p1m
);
51
CoolAdd
(
"Sc 2"
,42,
CoolHeavy
.
Sc24p2m
);
52
CoolAdd
(
"Sc 2"
,21,
CoolHeavy
.
Sc22p08m
);
53
54
/* Sc III 3933
55
* POPEXC( O12,g1,g2,A21,excit,abund); result already*a21
56
* [Sc III] 3933, multiplet average */
57
g1 = 10.;
58
g2 = 2.;
59
CoolHeavy
.
Sc33936
=
atom_pop2
(g1,g1,g2,0.03,3.66e4,
60
dense
.
xIonDense
[
ipSCANDIUM
][2])* 5.06e-12;
61
CoolAdd
(
"Sc 3"
,3933,
CoolHeavy
.
Sc33936
);
62
63
/* [Sc V] 2.31 mic
64
* Y(ik) from
65
* >>refer sc5 cs Pelan, J., & Berrington, K.A. 1995, A&A Suppl, 110, 209 */
66
PutCS
(6.00,&
TauLines
[
ipSc05231
]);
67
atom_level2
(&
TauLines
[ipSc05231]);
68
69
/* Sc VI */
70
a21 = 4.94;
71
a31 = 49.24;
72
a32 = 4.3;
73
g1 = 9.;
74
g2 = 5.;
75
g3 = 1.;
76
/* POP3(G1,G2,G3,O12,O13,O23,A21,A31,A32,E12,E23,P2,ABUND,GAM2)
77
* energies are in kelvin */
78
p3 =
atom_pop3
(g1,g2,g3,g1,g2,g3,a21,a31,a32,28464.,40045.,&p2,
79
dense
.
xIonDense
[
ipSCANDIUM
][5], 0.,0.,0.);
80
/* 2100 ang - 3=>1 */
81
CoolHeavy
.
Sc42100
= p3*a31*9.47e-12;
82
/* 5054 - 2=>1 */
83
CoolHeavy
.
Sc45058
= p2*a21*3.93e-12;
84
/* 3592 3=>2 */
85
CoolHeavy
.
Sc43595
= p3*a32*5.53e-12;
86
CoolAdd
(
"Sc 6"
,2100,
CoolHeavy
.
Sc42100
);
87
CoolAdd
(
"Sc 6"
,5054,
CoolHeavy
.
Sc45058
);
88
CoolAdd
(
"Sc 6"
,3592,
CoolHeavy
.
Sc43595
);
89
90
/* [Sc 13] 2637.97A, cs from
91
* >>refer sc13 cs Saraph, H.E. & Tully, J.A. 1994, A&AS, 107, 29 */
92
PutCS
(0.182,&
TauLines
[
ipSc13264
]);
93
atom_level2
(&
TauLines
[ipSc13264]);
94
95
return
;
96
}
Generated for cloudy by
1.8.3.1