MPI-AMRVAC 3.2
The MPI - Adaptive Mesh Refinement - Versatile Advection Code (development version)
Loading...
Searching...
No Matches
mod_radloss_tables.t
Go to the documentation of this file.
1!> module containing all optically thin radiative cooling tables
2!>
3!> Used by mod_radiative_cooling.t for HD and ffhd and MHD and twofl
4!> Assumptions: full ionized plasma dominated by H and He, ionization equilibrium
5!> Formula: Q=-n_H*n_e*f(T), positive f(T) function is pre-computed and tabulated or a piecewise power law
6!>
7module mod_radloss_tables
8!
9! For interpolatable tables: contain log_10 temperature values and corresponding
10! log_10 luminosity values.
11!
12 implicit none
13
14 double precision :: t_hildner(1:6), t_fm(1:5), t_rosner(1:10), t_klimchuk(1:8), &
15 t_spex_dm_rough(1:8), t_spex_dm_fine(1:15)
16
17 double precision :: x_hildner(1:5), x_fm(1:4), x_rosner(1:9), x_klimchuk(1:7), &
18 x_spex_dm_rough(1:7), x_spex_dm_fine(1:14)
19
20 double precision :: a_hildner(1:5), a_fm(1:4), a_rosner(1:9), a_klimchuk(1:7), &
21 a_spex_dm_rough(1:7), a_spex_dm_fine(1:14)
22
23 integer :: n_hildner, n_fm, n_rosner, n_klimchuk, n_spex_dm_rough, n_spex_dm_fine
24
25 data n_hildner / 5 /
26
27 data t_hildner / 3.00000d0, 4.17609d0, 4.90309d0, 5.47712d0, 5.90309d0, 10.00000d0 /
28
29 data x_hildner / -53.30803d0, -29.92082d0, -21.09691d0, -7.40450d0, -16.25885d0 /
30
31 data a_hildner / 7.4d0, 1.8d0, 0.0d0, -2.5d0, -1.0d0 /
32
33 data n_fm / 4 /
34
35 data t_fm / 3.00000d0, 4.30103d0, 5.60206d0, 7.00000d0, 10.00000d0 /
36
37 data x_fm / -31.15813d0, -27.50227d0, -18.25885d0, -35.75893d0 /
38
39 data a_fm / 2.0d0, 1.15d0, -0.5d0, 2.0d0 /
40
41 data n_rosner / 9 /
42
43 data t_rosner / 3.00000d0, 3.89063d0, 4.30195d0, 4.57500d0, 4.90000d0, &
44 5.40000d0, 5.77000d0, 6.31500d0, 7.60457d0, 10.00000d0 /
45
46 data x_rosner / -69.900d0, -48.307d0, -21.850d0, -31.000d0, -21.200d0, &
47 -10.400d0, -21.940d0, -17.730d0, -26.602d0 /
48
49 data a_rosner / 11.70d0, 6.15d0, 0.00d0, 2.00d0, 0.00d0, &
50 -2.00d0, 0.00d0, -0.67d0, 0.50d0 /
51
52 data n_klimchuk / 7 /
53
54 data t_klimchuk / 3.00d0, 4.97d0, 5.67d0, 6.18d0, 6.55d0, 6.90d0, 7.63d0, 10.00d0 /
55
56 data x_klimchuk / -30.96257d0, -16.05208d0, -21.72125d0, -12.45223d0, &
57 -24.46092d0, -15.26043d0, -26.70774d0 /
58
59 data a_klimchuk / 2.00d0, -1.00d0, 0.00d0, -1.50d0, 0.33d0, -1.00d0, 0.50d0 /
60
61 data n_spex_dm_rough / 7 /
62
63 data t_spex_dm_rough / 1.000d0, 1.572d0, 3.992d0, 4.165d0, 5.221d0, 5.751d0, 7.295d0, 8.160d0 /
64
65 data x_spex_dm_rough / -34.286d0, -28.282d0, -108.273d0, -26.662d0, -9.729d0, -17.550d0, -24.767d0 /
66
67 data a_spex_dm_rough / 4.560d0, 0.740d0, 20.777d0, 1.182d0, -2.061d0, -0.701d0, 0.288d0 /
68
69 data n_spex_dm_fine / 14 /
70
71 data t_spex_dm_fine / 1.000d0, 1.422d0, 2.806d0, 3.980d0, 4.177d0, 4.443d0, 4.832d0, 5.397d0, &
72 5.570d0, 5.890d0, 6.232d0, 6.505d0, 6.941d0, 7.385d0, 8.160d0 /
73
74 data x_spex_dm_fine / -35.314d0, -29.195d0, -26.912d0, -108.273d0, -18.971d0, -32.195d0, -21.217d0, &
75 -0.247d0, -15.415d0, -19.275d0, -9.387d0, -22.476d0, -17.437d0, -25.026d0 /
76
77 data a_spex_dm_fine / 5.452d0, 1.150d0, 0.337d0, 20.777d0, -0.602d0, 2.374d0, 0.102d0, &
78 -3.784d0, -1.061d0, -0.406d0, -1.992d0, 0.020d0, -0.706d0, 0.321d0 /
79
80 ! Interpolatable tables
81
82 double precision :: t_dm(1:71), t_mb(1:51), t_mlcosmol(1:71) &
83 , t_mlwc(1:71), t_mlsolar1(1:71), t_spex(1:110) &
84 , t_jccorona(1:45), t_cl_ism(1:151), t_cl_solar(1:151)&
85 , t_dm_2(1:76), t_dere(1:101), t_colgan(1:55)
86
87 double precision :: l_dm(1:71), l_mb(1:51), l_mlcosmol(1:71) &
88 , l_mlwc(1:71), l_mlsolar1(1:71), l_spex(1:110) &
89 , l_jccorona(1:45), l_cl_ism(1:151), l_cl_solar(1:151) &
90 , l_dm_2(1:76), l_dere_corona(1:101), l_dere_photo(1:101)&
91 , l_colgan(1:55)
92
93 double precision :: lowfip_frac(1:101)
94
95 double precision :: nenh_spex(1:110)
96
97 integer :: n_dm , n_mb , n_mlcosmol &
98 , n_mlwc , n_mlsolar1, n_spex &
99 , n_jccorona, n_cl_ism , n_cl_solar &
100 , n_dm_2 , n_dere , n_colgan
101
102 data n_jccorona / 45 /
103
104 data t_jccorona / 4.00000d0, 4.14230d0, 4.21995d0, 4.29761d0, 4.37528d0, &
105 4.45294d0, 4.53061d0, 4.60827d0, 4.68593d0, 4.76359d0, &
106 4.79705d0, 4.83049d0, 4.86394d0, 4.89739d0, 4.93084d0, &
107 4.96428d0, 4.99773d0, 5.03117d0, 5.06461d0, 5.17574d0, &
108 5.28684d0, 5.39796d0, 5.50907d0, 5.62018d0, 5.73129d0, &
109 5.84240d0, 5.95351d0, 6.06461d0, 6.17574d0, 6.28684d0, &
110 6.39796d0, 6.50907d0, 6.62018d0, 6.73129d0, 6.84240d0, &
111 6.95351d0, 7.06461d0, 7.17574d0, 7.28684d0, 7.39796d0, &
112 7.50907d0, 7.62018d0, 7.73129d0, 7.84240d0, 7.95351d0 /
113
114 data l_jccorona / -200.18883d0, -100.78630d0, -30.60384d0, -22.68481d0, -21.76445d0, &
115 -21.67936d0, -21.54218d0, -21.37958d0, -21.25172d0, -21.17584d0, &
116 -21.15783d0, -21.14491d0, -21.13527d0, -21.12837d0, -21.12485d0, &
117 -21.12439d0, -21.12642d0, -21.12802d0, -21.12548d0, -21.08965d0, &
118 -21.08812d0, -21.19542d0, -21.34582d0, -21.34839d0, -21.31701d0, &
119 -21.29072d0, -21.28900d0, -21.34104d0, -21.43122d0, -21.62448d0, &
120 -21.86694d0, -22.02897d0, -22.08051d0, -22.06057d0, -22.01973d0, &
121 -22.00000d0, -22.05161d0, -22.22175d0, -22.41452d0, -22.52581d0, &
122 -22.56914d0, -22.57486d0, -22.56151d0, -22.53969d0, -22.51490d0 /
123
124 data n_dm / 71 /
125
126 data t_dm / 2.0d0, 2.1d0, 2.2d0, 2.3d0, 2.4d0, &
127 2.5d0, 2.6d0, 2.7d0, 2.8d0, 2.9d0, &
128 3.0d0, 3.1d0, 3.2d0, 3.3d0, 3.4d0, &
129 3.5d0, 3.6d0, 3.7d0, 3.8d0, 3.9d0, &
130 4.0d0, 4.1d0, 4.2d0, 4.3d0, 4.4d0, &
131 4.5d0, 4.6d0, 4.7d0, 4.8d0, 4.9d0, &
132 5.0d0, 5.1d0, 5.2d0, 5.3d0, 5.4d0, &
133 5.5d0, 5.6d0, 5.7d0, 5.8d0, 5.9d0, &
134 6.0d0, 6.1d0, 6.2d0, 6.3d0, 6.4d0, &
135 6.5d0, 6.6d0, 6.7d0, 6.8d0, 6.9d0, &
136 7.0d0, 7.1d0, 7.2d0, 7.3d0, 7.4d0, &
137 7.5d0, 7.6d0, 7.7d0, 7.8d0, 7.9d0, &
138 8.0d0, 8.1d0, 8.2d0, 8.3d0, 8.4d0, &
139 8.5d0, 8.6d0, 8.7d0, 8.8d0, 8.9d0, &
140 9.0d0 /
141
142 data l_dm / -26.523d0, -26.398d0, -26.301d0, -26.222d0, -26.097d0 &
143 , -26.011d0, -25.936d0, -25.866d0, -25.807d0, -25.754d0 &
144 , -25.708d0, -25.667d0, -25.630d0, -25.595d0, -25.564d0 &
145 , -25.534d0, -25.506d0, -25.479d0, -25.453d0, -25.429d0 &
146 , -25.407d0, -23.019d0, -21.762d0, -21.742d0, -21.754d0 &
147 , -21.730d0, -21.523d0, -21.455d0, -21.314d0, -21.229d0 &
148 , -21.163d0, -21.126d0, -21.092d0, -21.060d0, -21.175d0 &
149 , -21.280d0, -21.390d0, -21.547d0, -21.762d0, -22.050d0 &
150 , -22.271d0, -22.521d0, -22.646d0, -22.660d0, -22.676d0 &
151 , -22.688d0, -22.690d0, -22.662d0, -22.635d0, -22.609d0 &
152 , -22.616d0, -22.646d0, -22.697d0, -22.740d0, -22.788d0 &
153 , -22.815d0, -22.785d0, -22.754d0, -22.728d0, -22.703d0 &
154 , -22.680d0, -22.630d0, -22.580d0, -22.530d0, -22.480d0 &
155 , -22.430d0, -22.380d0, -22.330d0, -22.280d0, -22.230d0 &
156 , -22.180d0 /
157
158 data n_mb / 51 /
159
160 data t_mb / 4.0d0, 4.1d0, 4.2d0, 4.3d0, 4.4d0, &
161 4.5d0, 4.6d0, 4.7d0, 4.8d0, 4.9d0, &
162 5.0d0, 5.1d0, 5.2d0, 5.3d0, 5.4d0, &
163 5.5d0, 5.6d0, 5.7d0, 5.8d0, 5.9d0, &
164 6.0d0, 6.1d0, 6.2d0, 6.3d0, 6.4d0, &
165 6.5d0, 6.6d0, 6.7d0, 6.8d0, 6.9d0, &
166 7.0d0, 7.1d0, 7.2d0, 7.3d0, 7.4d0, &
167 7.5d0, 7.6d0, 7.7d0, 7.8d0, 7.9d0, &
168 8.0d0, 8.1d0, 8.2d0, 8.3d0, 8.4d0, &
169 8.5d0, 8.6d0, 8.7d0, 8.8d0, 8.9d0, &
170 9.0d0 /
171
172 data l_mb / -23.133d0, -22.895d0, -22.548d0, -22.285d0, -22.099d0 &
173 , -21.970d0, -21.918d0, -21.826d0, -21.743d0, -21.638d0 &
174 , -21.552d0, -21.447d0, -21.431d0, -21.418d0, -21.461d0 &
175 , -21.570d0, -21.743d0, -21.832d0, -21.908d0, -21.981d0 &
176 , -22.000d0, -21.998d0, -21.992d0, -22.013d0, -22.095d0 &
177 , -22.262d0, -22.397d0, -22.445d0, -22.448d0, -22.446d0 &
178 , -22.448d0, -22.465d0, -22.575d0, -22.725d0, -22.749d0 &
179 , -22.768d0, -22.753d0, -22.717d0, -22.678d0, -22.637d0 &
180 , -22.603d0, -22.553d0, -22.503d0, -22.453d0, -22.403d0 &
181 , -22.353d0, -22.303d0, -22.253d0, -22.203d0, -22.153d0 &
182 , -22.103d0 /
183
184 data n_mlcosmol / 71 /
185
186 data t_mlcosmol / &
187 2.0d0, 2.1d0, 2.2d0, 2.3d0, 2.4d0, &
188 2.5d0, 2.6d0, 2.7d0, 2.8d0, 2.9d0, &
189 3.0d0, 3.1d0, 3.2d0, 3.3d0, 3.4d0, &
190 3.5d0, 3.6d0, 3.7d0, 3.8d0, 3.9d0, &
191 4.0d0, 4.1d0, 4.2d0, 4.3d0, 4.4d0, &
192 4.5d0, 4.6d0, 4.7d0, 4.8d0, 4.9d0, &
193 5.0d0, 5.1d0, 5.2d0, 5.3d0, 5.4d0, &
194 5.5d0, 5.6d0, 5.7d0, 5.8d0, 5.9d0, &
195 6.0d0, 6.1d0, 6.2d0, 6.3d0, 6.4d0, &
196 6.5d0, 6.6d0, 6.7d0, 6.8d0, 6.9d0, &
197 7.0d0, 7.1d0, 7.2d0, 7.3d0, 7.4d0, &
198 7.5d0, 7.6d0, 7.7d0, 7.8d0, 7.9d0, &
199 8.0d0, 8.1d0, 8.2d0, 8.3d0, 8.4d0, &
200 8.5d0, 8.6d0, 8.7d0, 8.8d0, 8.9d0, &
201 9.0d0 /
202
203 data l_mlcosmol / &
204 -99.000d0, -99.000d0, -99.000d0, -99.000d0, -99.000d0 &
205 , -99.000d0, -99.000d0, -99.000d0, -99.000d0, -99.000d0 &
206 , -99.000d0, -99.000d0, -99.000d0, -99.000d0, -99.000d0 &
207 , -99.000d0, -44.649d0, -38.362d0, -33.324d0, -29.292d0 &
208 , -26.063d0, -23.532d0, -22.192d0, -22.195d0, -22.454d0 &
209 , -22.676d0, -22.909d0, -22.925d0, -22.499d0, -22.276d0 &
210 , -22.440d0, -22.688d0, -22.917d0, -23.116d0, -23.274d0 &
211 , -23.394d0, -23.472d0, -23.516d0, -23.530d0, -23.525d0 &
212 , -23.506d0, -23.478d0, -23.444d0, -23.408d0, -23.368d0 &
213 , -23.328d0, -23.286d0, -23.244d0, -23.201d0, -23.157d0 &
214 , -23.114d0, -23.070d0, -23.026d0, -22.981d0, -22.937d0 &
215 , -22.893d0, -22.848d0, -22.803d0, -22.759d0, -22.714d0 &
216 , -22.669d0, -22.619d0, -22.569d0, -22.519d0, -22.469d0 &
217 , -22.419d0, -22.369d0, -22.319d0, -22.269d0, -22.190d0 &
218 , -22.169d0 /
219
220 data n_mlwc / 71 /
221
222 data t_mlwc / &
223 2.0d0, 2.1d0, 2.2d0, 2.3d0, 2.4d0, &
224 2.5d0, 2.6d0, 2.7d0, 2.8d0, 2.9d0, &
225 3.0d0, 3.1d0, 3.2d0, 3.3d0, 3.4d0, &
226 3.5d0, 3.6d0, 3.7d0, 3.8d0, 3.9d0, &
227 4.0d0, 4.1d0, 4.2d0, 4.3d0, 4.4d0, &
228 4.5d0, 4.6d0, 4.7d0, 4.8d0, 4.9d0, &
229 5.0d0, 5.1d0, 5.2d0, 5.3d0, 5.4d0, &
230 5.5d0, 5.6d0, 5.7d0, 5.8d0, 5.9d0, &
231 6.0d0, 6.1d0, 6.2d0, 6.3d0, 6.4d0, &
232 6.5d0, 6.6d0, 6.7d0, 6.8d0, 6.9d0, &
233 7.0d0, 7.1d0, 7.2d0, 7.3d0, 7.4d0, &
234 7.5d0, 7.6d0, 7.7d0, 7.8d0, 7.9d0, &
235 8.0d0, 8.1d0, 8.2d0, 8.3d0, 8.4d0, &
236 8.5d0, 8.6d0, 8.7d0, 8.8d0, 8.9d0, &
237 9.0d0 /
238
239 data l_mlwc / &
240 -21.192d0, -21.160d0, -21.150d0, -21.150d0, -21.166d0 &
241 , -21.191d0, -21.222d0, -21.264d0, -21.308d0, -21.357d0 &
242 , -21.408d0, -21.449d0, -21.494d0, -21.544d0, -21.587d0 &
243 , -21.638d0, -21.686d0, -21.736d0, -21.780d0, -21.800d0 &
244 , -21.744d0, -21.547d0, -21.208d0, -20.849d0, -20.345d0 &
245 , -19.771d0, -19.409d0, -19.105d0, -18.827d0, -18.555d0 &
246 , -18.460d0, -18.763d0, -19.168d0, -19.334d0, -19.400d0 &
247 , -19.701d0, -20.090d0, -20.288d0, -20.337d0, -20.301d0 &
248 , -20.233d0, -20.275d0, -20.363d0, -20.508d0, -20.675d0 &
249 , -20.856d0, -21.025d0, -21.159d0, -21.256d0, -21.320d0 &
250 , -21.354d0, -21.366d0, -21.361d0, -21.343d0, -21.317d0 &
251 , -21.285d0, -21.250d0, -21.212d0, -21.172d0, -21.131d0 &
252 , -21.089d0, -21.039d0, -20.989d0, -20.939d0, -20.889d0 &
253 , -20.839d0, -20.789d0, -20.739d0, -20.689d0, -20.639d0 &
254 , -20.589d0 /
255
256 data n_mlsolar1 / 71 /
257
258 data t_mlsolar1 / &
259 2.0d0, 2.1d0, 2.2d0, 2.3d0, 2.4d0, &
260 2.5d0, 2.6d0, 2.7d0, 2.8d0, 2.9d0, &
261 3.0d0, 3.1d0, 3.2d0, 3.3d0, 3.4d0, &
262 3.5d0, 3.6d0, 3.7d0, 3.8d0, 3.9d0, &
263 4.0d0, 4.1d0, 4.2d0, 4.3d0, 4.4d0, &
264 4.5d0, 4.6d0, 4.7d0, 4.8d0, 4.9d0, &
265 5.0d0, 5.1d0, 5.2d0, 5.3d0, 5.4d0, &
266 5.5d0, 5.6d0, 5.7d0, 5.8d0, 5.9d0, &
267 6.0d0, 6.1d0, 6.2d0, 6.3d0, 6.4d0, &
268 6.5d0, 6.6d0, 6.7d0, 6.8d0, 6.9d0, &
269 7.0d0, 7.1d0, 7.2d0, 7.3d0, 7.4d0, &
270 7.5d0, 7.6d0, 7.7d0, 7.8d0, 7.9d0, &
271 8.0d0, 8.1d0, 8.2d0, 8.3d0, 8.4d0, &
272 8.5d0, 8.6d0, 8.7d0, 8.8d0, 8.9d0, &
273 9.0d0 /
274
275 data l_mlsolar1 / &
276 -26.983d0, -26.951d0, -26.941d0, -26.940d0, -26.956d0 &
277 , -26.980d0, -27.011d0, -27.052d0, -27.097d0, -27.145d0 &
278 , -27.195d0, -27.235d0, -27.279d0, -27.327d0, -27.368d0 &
279 , -27.415d0, -27.456d0, -27.485d0, -27.468d0, -27.223d0 &
280 , -25.823d0, -23.501d0, -22.162d0, -22.084d0, -22.157d0 &
281 , -22.101d0, -21.974d0, -21.782d0, -21.542d0, -21.335d0 &
282 , -21.251d0, -21.275d0, -21.236d0, -21.173d0, -21.167d0 &
283 , -21.407d0, -21.670d0, -21.788d0, -21.879d0, -22.008d0 &
284 , -22.192d0, -22.912d0, -22.918d0, -22.887d0, -22.929d0 &
285 , -23.023d0, -23.094d0, -23.117d0, -23.108d0, -23.083d0 &
286 , -23.049d0, -23.011d0, -22.970d0, -22.928d0, -22.885d0 &
287 , -22.842d0, -22.798d0, -22.754d0, -22.709d0, -22.665d0 &
288 , -22.620d0, -22.570d0, -22.520d0, -22.470d0, -22.420d0 &
289 , -22.370d0, -22.320d0, -22.270d0, -22.220d0, -22.170d0 &
290 , -22.120d0 /
291
292 data n_spex / 110 /
293
294 data t_spex / &
295 3.80d0, 3.84d0, 3.88d0, 3.92d0, 3.96d0, &
296 4.00d0, 4.04d0, 4.08d0, 4.12d0, 4.16d0, &
297 4.20d0, 4.24d0, 4.28d0, 4.32d0, 4.36d0, &
298 4.40d0, 4.44d0, 4.48d0, 4.52d0, 4.56d0, &
299 4.60d0, 4.64d0, 4.68d0, 4.72d0, 4.76d0, &
300 4.80d0, 4.84d0, 4.88d0, 4.92d0, 4.96d0, &
301 5.00d0, 5.04d0, 5.08d0, 5.12d0, 5.16d0, &
302 5.20d0, 5.24d0, 5.28d0, 5.32d0, 5.36d0, &
303 5.40d0, 5.44d0, 5.48d0, 5.52d0, 5.56d0, &
304 5.60d0, 5.64d0, 5.68d0, 5.72d0, 5.76d0, &
305 5.80d0, 5.84d0, 5.88d0, 5.92d0, 5.96d0, &
306 6.00d0, 6.04d0, 6.08d0, 6.12d0, 6.16d0, &
307 6.20d0, 6.24d0, 6.28d0, 6.32d0, 6.36d0, &
308 6.40d0, 6.44d0, 6.48d0, 6.52d0, 6.56d0, &
309 6.60d0, 6.64d0, 6.68d0, 6.72d0, 6.76d0, &
310 6.80d0, 6.84d0, 6.88d0, 6.92d0, 6.96d0, &
311 7.00d0, 7.04d0, 7.08d0, 7.12d0, 7.16d0, &
312 7.20d0, 7.24d0, 7.28d0, 7.32d0, 7.36d0, &
313 7.40d0, 7.44d0, 7.48d0, 7.52d0, 7.56d0, &
314 7.60d0, 7.64d0, 7.68d0, 7.72d0, 7.76d0, &
315 7.80d0, 7.84d0, 7.88d0, 7.92d0, 7.96d0, &
316 8.00d0, 8.04d0, 8.08d0, 8.12d0, 8.16d0 /
317
318 data l_spex / &
319 -25.7331d0, -25.0383d0, -24.4059d0, -23.8288d0, -23.3027d0 &
320 , -22.8242d0, -22.3917d0, -22.0067d0, -21.6818d0, -21.4529d0 &
321 , -21.3246d0, -21.3459d0, -21.4305d0, -21.5293d0, -21.6138d0 &
322 , -21.6615d0, -21.6551d0, -21.5919d0, -21.5092d0, -21.4124d0 &
323 , -21.3085d0, -21.2047d0, -21.1067d0, -21.0194d0, -20.9413d0 &
324 , -20.8735d0, -20.8205d0, -20.7805d0, -20.7547d0, -20.7455d0 &
325 , -20.7565d0, -20.7820d0, -20.8008d0, -20.7994d0, -20.7847d0 &
326 , -20.7687d0, -20.7590d0, -20.7544d0, -20.7505d0, -20.7545d0 &
327 , -20.7888d0, -20.8832d0, -21.0450d0, -21.2286d0, -21.3737d0 &
328 , -21.4573d0, -21.4935d0, -21.5098d0, -21.5345d0, -21.5863d0 &
329 , -21.6548d0, -21.7108d0, -21.7424d0, -21.7576d0, -21.7696d0 &
330 , -21.7883d0, -21.8115d0, -21.8303d0, -21.8419d0, -21.8514d0 &
331 , -21.8690d0, -21.9057d0, -21.9690d0, -22.0554d0, -22.1488d0 &
332 , -22.2355d0, -22.3084d0, -22.3641d0, -22.4033d0, -22.4282d0 &
333 , -22.4408d0, -22.4443d0, -22.4411d0, -22.4334d0, -22.4242d0 &
334 , -22.4164d0, -22.4134d0, -22.4168d0, -22.4267d0, -22.4418d0 &
335 , -22.4603d0, -22.4830d0, -22.5112d0, -22.5449d0, -22.5819d0 &
336 , -22.6177d0, -22.6483d0, -22.6719d0, -22.6883d0, -22.6985d0 &
337 , -22.7032d0, -22.7037d0, -22.7008d0, -22.6950d0, -22.6869d0 &
338 , -22.6769d0, -22.6655d0, -22.6531d0, -22.6397d0, -22.6258d0 &
339 , -22.6111d0, -22.5964d0, -22.5816d0, -22.5668d0, -22.5519d0 &
340 , -22.5367d0, -22.5216d0, -22.5062d0, -22.4912d0, -22.4753d0 /
341
342 data nenh_spex / &
343 0.000013264d0, 0.000042428d0, 0.000088276d0, 0.00017967d0 &
344 , 0.00084362d0, 0.0034295d0, 0.013283d0, 0.042008d0 &
345 , 0.12138d0, 0.30481d0, 0.53386d0, 0.76622d0 &
346 , 0.89459d0, 0.95414d0, 0.98342d0 &
347 , 1.0046d0, 1.0291d0, 1.0547d0 &
348 , 1.0767d0, 1.0888d0 &
349 , 1.0945d0, 1.0972d0, 1.0988d0 &
350 , 1.1004d0, 1.1034d0 &
351 , 1.1102d0, 1.1233d0, 1.1433d0 &
352 , 1.1638d0, 1.1791d0 &
353 , 1.1885d0, 1.1937d0, 1.1966d0 &
354 , 1.1983d0, 1.1993d0 &
355 , 1.1999d0, 1.2004d0, 1.2008d0 &
356 , 1.2012d0, 1.2015d0 &
357 , 1.2020d0, 1.2025d0, 1.2030d0 &
358 , 1.2035d0, 1.2037d0 &
359 , 1.2039d0, 1.2040d0, 1.2041d0 &
360 , 1.2042d0, 1.2044d0 &
361 , 1.2045d0, 1.2046d0, 1.2047d0 &
362 , 1.2049d0, 1.2050d0 &
363 , 1.2051d0, 1.2053d0, 1.2055d0 &
364 , 1.2056d0, 1.2058d0 &
365 , 1.2060d0, 1.2062d0, 1.2065d0 &
366 , 1.2067d0, 1.2070d0 &
367 , 1.2072d0, 1.2075d0, 1.2077d0 &
368 , 1.2078d0, 1.2079d0 &
369 , 1.2080d0, 1.2081d0, 1.2082d0 &
370 , 1.2083d0, 1.2083d0 &
371 , 1.2084d0, 1.2084d0, 1.2085d0 &
372 , 1.2085d0, 1.2086d0 &
373 , 1.2086d0, 1.2087d0, 1.2087d0 &
374 , 1.2088d0, 1.2088d0 &
375 , 1.2089d0, 1.2089d0, 1.2089d0 &
376 , 1.2089d0, 1.2089d0 &
377 , 1.2090d0, 1.2090d0, 1.2090d0 &
378 , 1.2090d0, 1.2090d0 &
379 , 1.2090d0, 1.2090d0, 1.2090d0 &
380 , 1.2090d0, 1.2090d0 &
381 , 1.2090d0, 1.2090d0, 1.2090d0 &
382 , 1.2090d0, 1.2090d0 &
383 , 1.2090d0, 1.2090d0, 1.2090d0 &
384 , 1.2090d0, 1.2090d0 /
385 !
386 ! To be used together with the SPEX table for the SPEX_DM option
387 ! Assuming an ionization fraction of 10^-3
388 !
389 data n_dm_2 / 76 /
390
391 data t_dm_2 / 1.00d0, 1.04d0, 1.08d0, 1.12d0, 1.16d0, 1.20d0 &
392 , 1.24d0, 1.28d0, 1.32d0, 1.36d0, 1.40d0 &
393 , 1.44d0, 1.48d0, 1.52d0, 1.56d0, 1.60d0 &
394 , 1.64d0, 1.68d0, 1.72d0, 1.76d0, 1.80d0 &
395 , 1.84d0, 1.88d0, 1.92d0, 1.96d0, 2.00d0 &
396 , 2.04d0, 2.08d0, 2.12d0, 2.16d0, 2.20d0 &
397 , 2.24d0, 2.28d0, 2.32d0, 2.36d0, 2.40d0 &
398 , 2.44d0, 2.48d0, 2.52d0, 2.56d0, 2.60d0 &
399 , 2.64d0, 2.68d0, 2.72d0, 2.76d0, 2.80d0 &
400 , 2.84d0, 2.88d0, 2.92d0, 2.96d0, 3.00d0 &
401 , 3.04d0, 3.08d0, 3.12d0, 3.16d0, 3.20d0 &
402 , 3.24d0, 3.28d0, 3.32d0, 3.36d0, 3.40d0 &
403 , 3.44d0, 3.48d0, 3.52d0, 3.56d0, 3.60d0 &
404 , 3.64d0, 3.68d0, 3.72d0, 3.76d0, 3.80d0 &
405 , 3.84d0, 3.88d0, 3.92d0, 3.96d0, 4.00d0 /
406
407 data l_dm_2 / -30.0377d0, -29.7062d0, -29.4055d0, -29.1331d0, -28.8864d0, -28.6631d0 &
408 , -28.4614d0, -28.2791d0, -28.1146d0, -27.9662d0, -27.8330d0 &
409 , -27.7129d0, -27.6052d0, -27.5088d0, -27.4225d0, -27.3454d0 &
410 , -27.2767d0, -27.2153d0, -27.1605d0, -27.1111d0, -27.0664d0 &
411 , -27.0251d0, -26.9863d0, -26.9488d0, -26.9119d0, -26.8742d0 &
412 , -26.8353d0, -26.7948d0, -26.7523d0, -26.7080d0, -26.6619d0 &
413 , -26.6146d0, -26.5666d0, -26.5183d0, -26.4702d0, -26.4229d0 &
414 , -26.3765d0, -26.3317d0, -26.2886d0, -26.2473d0, -26.2078d0 &
415 , -26.1704d0, -26.1348d0, -26.1012d0, -26.0692d0, -26.0389d0 &
416 , -26.0101d0, -25.9825d0, -25.9566d0, -25.9318d0, -25.9083d0 &
417 , -25.8857d0, -25.8645d0, -25.8447d0, -25.8259d0, -25.8085d0 &
418 , -25.7926d0, -25.7778d0, -25.7642d0, -25.7520d0, -25.7409d0 &
419 , -25.7310d0, -25.7222d0, -25.7142d0, -25.7071d0, -25.7005d0 &
420 , -25.6942d0, -25.6878d0, -25.6811d0, -25.6733d0, -25.6641d0 &
421 , -25.6525d0, -25.6325d0, -25.6080d0, -25.5367d0, -25.4806d0 /
422
423 data n_cl_solar / 151 /
424
425 data t_cl_solar / &
426 1.000d0 , 1.050d0 , 1.100d0 , 1.150d0 , &
427 1.200d0 , 1.250d0 , 1.300d0 , 1.350d0 , &
428 1.400d0 , 1.450d0 , 1.500d0 , 1.550d0 , &
429 1.600d0 , 1.650d0 , 1.700d0 , 1.750d0 , &
430 1.800d0 , 1.850d0 , 1.900d0 , 1.950d0 , &
431 2.000d0 , 2.050d0 , 2.100d0 , 2.150d0 , &
432 2.200d0 , 2.250d0 , 2.300d0 , 2.350d0 , &
433 2.400d0 , 2.450d0 , 2.500d0 , 2.550d0 , &
434 2.600d0 , 2.650d0 , 2.700d0 , 2.750d0 , &
435 2.800d0 , 2.850d0 , 2.900d0 , 2.950d0 , &
436 3.000d0 , 3.050d0 , 3.100d0 , 3.150d0 , &
437 3.200d0 , 3.250d0 , 3.300d0 , 3.350d0 , &
438 3.400d0 , 3.450d0 , 3.500d0 , 3.550d0 , &
439 3.600d0 , 3.650d0 , 3.700d0 , 3.750d0 , &
440 3.800d0 , 3.850d0 , 3.900d0 , 3.950d0 , &
441 4.000d0 , 4.050d0 , 4.100d0 , 4.150d0 , &
442 4.200d0 , 4.250d0 , 4.300d0 , 4.350d0 , &
443 4.400d0 , 4.450d0 , 4.500d0 , 4.550d0 , &
444 4.600d0 , 4.650d0 , 4.700d0 , 4.750d0 , &
445 4.800d0 , 4.850d0 , 4.900d0 , 4.950d0 , &
446 5.000d0 , 5.050d0 , 5.100d0 , 5.150d0 , &
447 5.200d0 , 5.250d0 , 5.300d0 , 5.350d0 , &
448 5.400d0 , 5.450d0 , 5.500d0 , 5.550d0 , &
449 5.600d0 , 5.650d0 , 5.700d0 , 5.750d0 , &
450 5.800d0 , 5.850d0 , 5.900d0 , 5.950d0 , &
451 6.000d0 , 6.050d0 , 6.100d0 , 6.150d0 , &
452 6.200d0 , 6.250d0 , 6.300d0 , 6.350d0 , &
453 6.400d0 , 6.450d0 , 6.500d0 , 6.550d0 , &
454 6.600d0 , 6.650d0 , 6.700d0 , 6.750d0 , &
455 6.800d0 , 6.850d0 , 6.900d0 , 6.950d0 , &
456 7.000d0 , 7.050d0 , 7.100d0 , 7.150d0 , &
457 7.200d0 , 7.250d0 , 7.300d0 , 7.350d0 , &
458 7.400d0 , 7.450d0 , 7.500d0 , 7.550d0 , &
459 7.600d0 , 7.650d0 , 7.700d0 , 7.750d0 , &
460 7.800d0 , 7.850d0 , 7.900d0 , 7.950d0 , &
461 8.000d0 , 8.100d0 , 8.200d0 , 8.300d0 , &
462 8.400d0 , 8.500d0 , 8.600d0 , 8.700d0 , &
463 8.800d0 , 8.900d0 , 9.00d0 /
464
465 data l_cl_solar / &
466 -28.375d0 , -28.251d0 , -28.137d0 , -28.029d0 , &
467 -27.929d0 , -27.834d0 , -27.745d0 , -27.662d0 , &
468 -27.584d0 , -27.512d0 , -27.445d0 , -27.383d0 , &
469 -27.326d0 , -27.273d0 , -27.223d0 , -27.175d0 , &
470 -27.128d0 , -27.079d0 , -27.027d0 , -26.972d0 , &
471 -26.911d0 , -26.846d0 , -26.777d0 , -26.705d0 , &
472 -26.632d0 , -26.554d0 , -26.479d0 , -26.407d0 , &
473 -26.338d0 , -26.274d0 , -26.213d0 , -26.156d0 , &
474 -26.101d0 , -26.049d0 , -25.999d0 , -25.949d0 , &
475 -25.901d0 , -25.852d0 , -25.803d0 , -25.754d0 , &
476 -25.707d0 , -25.662d0 , -25.621d0 , -25.588d0 , &
477 -25.561d0 , -25.538d0 , -25.518d0 , -25.497d0 , &
478 -25.475d0 , -25.452d0 , -25.426d0 , -25.400d0 , &
479 -25.374d0 , -25.333d0 , -25.295d0 , -25.261d0 , &
480 -25.228d0 , -25.189d0 , -25.136d0 , -25.053d0 , &
481 -24.888d0 , -24.454d0 , -23.480d0 , -22.562d0 , &
482 -22.009d0 , -21.826d0 , -21.840d0 , -21.905d0 , &
483 -21.956d0 , -21.971d0 , -21.958d0 , -21.928d0 , &
484 -21.879d0 , -21.810d0 , -21.724d0 , -21.623d0 , &
485 -21.512d0 , -21.404d0 , -21.321d0 , -21.273d0 , &
486 -21.250d0 , -21.253d0 , -21.275d0 , -21.287d0 , &
487 -21.282d0 , -21.275d0 , -21.272d0 , -21.267d0 , &
488 -21.281d0 , -21.357d0 , -21.496d0 , -21.616d0 , &
489 -21.677d0 , -21.698d0 , -21.708d0 , -21.730d0 , &
490 -21.767d0 , -21.793d0 , -21.794d0 , -21.787d0 , &
491 -21.787d0 , -21.802d0 , -21.826d0 , -21.859d0 , &
492 -21.911d0 , -21.987d0 , -22.082d0 , -22.173d0 , &
493 -22.253d0 , -22.325d0 , -22.392d0 , -22.448d0 , &
494 -22.487d0 , -22.512d0 , -22.524d0 , -22.528d0 , &
495 -22.524d0 , -22.516d0 , -22.507d0 , -22.501d0 , &
496 -22.502d0 , -22.511d0 , -22.533d0 , -22.565d0 , &
497 -22.600d0 , -22.630d0 , -22.648d0 , -22.656d0 , &
498 -22.658d0 , -22.654d0 , -22.647d0 , -22.634d0 , &
499 -22.619d0 , -22.602d0 , -22.585d0 , -22.566d0 , &
500 -22.546d0 , -22.525d0 , -22.505d0 , -22.480d0 , &
501 -22.465d0 , -22.415d0 , -22.365d0 , -22.315d0 , &
502 -22.265d0 , -22.215d0 , -22.165d0 , -22.115d0 , &
503 -22.065d0 , -22.015d0 , -21.965d0 /
504
505 data n_cl_ism / 151 /
506
507 data t_cl_ism / &
508 1.000d0 , 1.050d0 , 1.100d0 , 1.150d0 , &
509 1.200d0 , 1.250d0 , 1.300d0 , 1.350d0 , &
510 1.400d0 , 1.450d0 , 1.500d0 , 1.550d0 , &
511 1.600d0 , 1.650d0 , 1.700d0 , 1.750d0 , &
512 1.800d0 , 1.850d0 , 1.900d0 , 1.950d0 , &
513 2.000d0 , 2.050d0 , 2.100d0 , 2.150d0 , &
514 2.200d0 , 2.250d0 , 2.300d0 , 2.350d0 , &
515 2.400d0 , 2.450d0 , 2.500d0 , 2.550d0 , &
516 2.600d0 , 2.650d0 , 2.700d0 , 2.750d0 , &
517 2.800d0 , 2.850d0 , 2.900d0 , 2.950d0 , &
518 3.000d0 , 3.050d0 , 3.100d0 , 3.150d0 , &
519 3.200d0 , 3.250d0 , 3.300d0 , 3.350d0 , &
520 3.400d0 , 3.450d0 , 3.500d0 , 3.550d0 , &
521 3.600d0 , 3.650d0 , 3.700d0 , 3.750d0 , &
522 3.800d0 , 3.850d0 , 3.900d0 , 3.950d0 , &
523 4.000d0 , 4.050d0 , 4.100d0 , 4.150d0 , &
524 4.200d0 , 4.250d0 , 4.300d0 , 4.350d0 , &
525 4.400d0 , 4.450d0 , 4.500d0 , 4.550d0 , &
526 4.600d0 , 4.650d0 , 4.700d0 , 4.750d0 , &
527 4.800d0 , 4.850d0 , 4.900d0 , 4.950d0 , &
528 5.000d0 , 5.050d0 , 5.100d0 , 5.150d0 , &
529 5.200d0 , 5.250d0 , 5.300d0 , 5.350d0 , &
530 5.400d0 , 5.450d0 , 5.500d0 , 5.550d0 , &
531 5.600d0 , 5.650d0 , 5.700d0 , 5.750d0 , &
532 5.800d0 , 5.850d0 , 5.900d0 , 5.950d0 , &
533 6.000d0 , 6.050d0 , 6.100d0 , 6.150d0 , &
534 6.200d0 , 6.250d0 , 6.300d0 , 6.350d0 , &
535 6.400d0 , 6.450d0 , 6.500d0 , 6.550d0 , &
536 6.600d0 , 6.650d0 , 6.700d0 , 6.750d0 , &
537 6.800d0 , 6.850d0 , 6.900d0 , 6.950d0 , &
538 7.000d0 , 7.050d0 , 7.100d0 , 7.150d0 , &
539 7.200d0 , 7.250d0 , 7.300d0 , 7.350d0 , &
540 7.400d0 , 7.450d0 , 7.500d0 , 7.550d0 , &
541 7.600d0 , 7.650d0 , 7.700d0 , 7.750d0 , &
542 7.800d0 , 7.850d0 , 7.900d0 , 7.950d0 , &
543 8.000d0 , 8.100d0 , 8.200d0 , 8.300d0 , &
544 8.400d0 , 8.500d0 , 8.600d0 , 8.700d0 , &
545 8.800d0 , 8.900d0 , 9.00d0 /
546
547 data l_cl_ism / &
548 -28.365d0 , -28.242d0 , -28.127d0 , -28.020d0 , &
549 -27.919d0 , -27.825d0 , -27.736d0 , -27.653d0 , &
550 -27.575d0 , -27.504d0 , -27.437d0 , -27.376d0 , &
551 -27.319d0 , -27.267d0 , -27.220d0 , -27.176d0 , &
552 -27.134d0 , -27.095d0 , -27.058d0 , -27.021d0 , &
553 -26.985d0 , -26.948d0 , -26.910d0 , -26.870d0 , &
554 -26.827d0 , -26.775d0 , -26.721d0 , -26.664d0 , &
555 -26.608d0 , -26.552d0 , -26.495d0 , -26.437d0 , &
556 -26.378d0 , -26.317d0 , -26.255d0 , -26.190d0 , &
557 -26.123d0 , -26.053d0 , -25.984d0 , -25.913d0 , &
558 -25.847d0 , -25.786d0 , -25.736d0 , -25.702d0 , &
559 -25.678d0 , -25.662d0 , -25.649d0 , -25.636d0 , &
560 -25.621d0 , -25.604d0 , -25.587d0 , -25.571d0 , &
561 -25.562d0 , -25.526d0 , -25.505d0 , -25.499d0 , &
562 -25.499d0 , -25.491d0 , -25.468d0 , -25.410d0 , &
563 -25.268d0 , -24.888d0 , -23.702d0 , -22.624d0 , &
564 -22.036d0 , -21.843d0 , -21.854d0 , -21.924d0 , &
565 -21.986d0 , -22.017d0 , -22.021d0 , -22.005d0 , &
566 -21.964d0 , -21.896d0 , -21.806d0 , -21.699d0 , &
567 -21.580d0 , -21.463d0 , -21.370d0 , -21.312d0 , &
568 -21.284d0 , -21.290d0 , -21.322d0 , -21.345d0 , &
569 -21.354d0 , -21.366d0 , -21.385d0 , -21.396d0 , &
570 -21.414d0 , -21.483d0 , -21.600d0 , -21.696d0 , &
571 -21.742d0 , -21.759d0 , -21.776d0 , -21.816d0 , &
572 -21.885d0 , -21.939d0 , -21.946d0 , -21.918d0 , &
573 -21.873d0 , -21.818d0 , -21.756d0 , -21.689d0 , &
574 -21.618d0 , -21.547d0 , -21.475d0 , -21.403d0 , &
575 -21.331d0 , -21.260d0 , -21.188d0 , -21.114d0 , &
576 -21.039d0 , -20.963d0 , -20.887d0 , -20.810d0 , &
577 -20.734d0 , -20.657d0 , -20.581d0 , -20.505d0 , &
578 -20.429d0 , -20.352d0 , -20.276d0 , -20.200d0 , &
579 -20.125d0 , -20.049d0 , -19.973d0 , -19.898d0 , &
580 -19.822d0 , -19.747d0 , -19.671d0 , -19.596d0 , &
581 -19.520d0 , -19.445d0 , -19.370d0 , -19.295d0 , &
582 -19.220d0 , -19.144d0 , -19.069d0 , -18.994d0 , &
583 -18.919d0 , -18.869d0 , -18.819d0 , -18.769d0 , &
584 -18.719d0 , -18.669d0 , -18.619d0 , -18.569d0 , &
585 -18.519d0 , -18.469d0 , -18.419d0 /
586
587 data n_dere / 101 /
588
589 data t_dere / 4.00d0, 4.05d0, 4.10d0, 4.15d0, 4.20d0, 4.25d0, 4.30d0, 4.35d0, &
590 4.40d0, 4.45d0, 4.50d0, 4.55d0, 4.60d0, 4.65d0, 4.70d0, 4.75d0, &
591 4.80d0, 4.85d0, 4.90d0, 4.95d0, 5.00d0, 5.05d0, 5.10d0, 5.15d0, &
592 5.20d0, 5.25d0, 5.30d0, 5.35d0, 5.40d0, 5.45d0, 5.50d0, 5.55d0, &
593 5.60d0, 5.65d0, 5.70d0, 5.75d0, 5.80d0, 5.85d0, 5.90d0, 5.95d0, &
594 6.00d0, 6.05d0, 6.10d0, 6.15d0, 6.20d0, 6.25d0, 6.30d0, 6.35d0, &
595 6.40d0, 6.45d0, 6.50d0, 6.55d0, 6.60d0, 6.65d0, 6.70d0, 6.75d0, &
596 6.80d0, 6.85d0, 6.90d0, 6.95d0, 7.00d0, 7.05d0, 7.10d0, 7.15d0, &
597 7.20d0, 7.25d0, 7.30d0, 7.35d0, 7.40d0, 7.45d0, 7.50d0, 7.55d0, &
598 7.60d0, 7.65d0, 7.70d0, 7.75d0, 7.80d0, 7.85d0, 7.90d0, 7.95d0, &
599 8.00d0, 8.05d0, 8.10d0, 8.15d0, 8.20d0, 8.25d0, 8.30d0, 8.35d0, &
600 8.40d0, 8.45d0, 8.50d0, 8.55d0, 8.60d0, 8.65d0, 8.70d0, 8.75d0, &
601 8.80d0, 8.85d0, 8.90d0, 8.95d0, 9.00d0 /
602
603 data l_dere_corona / &
604 -23.00744648d0, -22.55439580d0, -22.15614458d0, -21.83268267d0, -21.64589156d0, &
605 -21.61618463d0, -21.68402965d0, -21.79048499d0, -21.87614836d0, -21.91009489d0, &
606 -21.89962945d0, -21.86012091d0, -21.79588002d0, -21.71669877d0, -21.62342304d0, &
607 -21.52143350d0, -21.41793664d0, -21.33068312d0, -21.27736608d0, -21.25181197d0, &
608 -21.24184538d0, -21.25806092d0, -21.27901426d0, -21.27164622d0, -21.24412514d0, &
609 -21.21467016d0, -21.19586057d0, -21.18309616d0, -21.18708664d0, -21.24108811d0, &
610 -21.35163999d0, -21.45099674d0, -21.48678240d0, -21.47625353d0, -21.45222529d0, &
611 -21.43297363d0, -21.42596873d0, -21.42021640d0, -21.41005040d0, -21.40120949d0, &
612 -21.40450378d0, -21.42250820d0, -21.44977165d0, -21.47755577d0, -21.51144928d0, &
613 -21.55909092d0, -21.63451202d0, -21.73754891d0, -21.85078089d0, -21.95467702d0, &
614 -22.03526908d0, -22.08990945d0, -22.11804503d0, -22.12436006d0, -22.11633856d0, &
615 -22.10072681d0, -22.08301995d0, -22.06701918d0, -22.05650548d0, -22.05551733d0, &
616 -22.06803389d0, -22.10072681d0, -22.15926677d0, -22.24033216d0, -22.32882716d0, &
617 -22.40782324d0, -22.47108330d0, -22.51855737d0, -22.54975089d0, -22.57024772d0, &
618 -22.58004425d0, -22.58335949d0, -22.58169871d0, -22.57511836d0, -22.56543110d0, &
619 -22.55439580d0, -22.54060751d0, -22.52724355d0, -22.51427857d0, -22.49894074d0, &
620 -22.48545225d0, -22.47108330d0, -22.45593196d0, -22.44009337d0, -22.42365865d0, &
621 -22.40671393d0, -22.38933984d0, -22.37059040d0, -22.35163999d0, -22.33161408d0, &
622 -22.31158018d0, -22.29073004d0, -22.26921772d0, -22.24795155d0, -22.22621356d0, &
623 -22.20411998d0, -22.18111459d0, -22.15864053d0, -22.13608262d0, -22.11294562d0, &
624 -22.08937560d0 /
625
626 data l_dere_photo / &
627 -23.25649024d0, -22.74232143d0, -22.28988263d0, -21.93554201d0, -21.74232143d0, &
628 -21.72353820d0, -21.80410035d0, -21.91009489d0, -22.00000000d0, -22.04143612d0, &
629 -22.04575749d0, -22.02502801d0, -21.97469413d0, -21.89962945d0, -21.80410035d0, &
630 -21.69250396d0, -21.57511836d0, -21.46852108d0, -21.38827669d0, -21.33629907d0, &
631 -21.31069114d0, -21.31966449d0, -21.33724217d0, -21.32882716d0, -21.30189945d0, &
632 -21.27408837d0, -21.25649024d0, -21.24718357d0, -21.25649024d0, -21.32148162d0, &
633 -21.46092390d0, -21.61083392d0, -21.70774393d0, -21.74958000d0, -21.76955108d0, &
634 -21.79588002d0, -21.84466396d0, -21.88941029d0, -21.91009489d0, -21.91721463d0, &
635 -21.92811799d0, -21.94692156d0, -21.96657624d0, -21.99139983d0, -22.01412464d0, &
636 -22.04963515d0, -22.10402527d0, -22.17848647d0, -22.26201267d0, -22.34390180d0, &
637 -22.41680123d0, -22.47237010d0, -22.50723961d0, -22.52578374d0, -22.53313238d0, &
638 -22.53165267d0, -22.52578374d0, -22.51855737d0, -22.51286162d0, -22.51286162d0, &
639 -22.52143350d0, -22.54211810d0, -22.57839607d0, -22.62708800d0, -22.67366414d0, &
640 -22.70996539d0, -22.73282827d0, -22.74714697d0, -22.75202673d0, -22.74958000d0, &
641 -22.74232143d0, -22.73282827d0, -22.72124640d0, -22.70553377d0, -22.69036983d0, &
642 -22.67162040d0, -22.65364703d0, -22.63638802d0, -22.61618463d0, -22.59687948d0, &
643 -22.57675413d0, -22.55752023d0, -22.53610701d0, -22.51570016d0, -22.49485002d0, &
644 -22.47366072d0, -22.45222529d0, -22.42945706d0, -22.40782324d0, -22.38510278d0, &
645 -22.36251027d0, -22.33913452d0, -22.31605287d0, -22.29242982d0, -22.26921772d0, &
646 -22.24565166d0, -22.22184875d0, -22.19859629d0, -22.17457388d0, -22.15058059d0, &
647 -22.12609840d0 /
648
649 data lowfip_frac / &
650 0.22546617d0, 0.14134011d0, 0.08795079d0, 0.06053747d0, 0.05337165d0, &
651 0.05879824d0, 0.06897329d0, 0.083151d0 , 0.10490873d0, 0.13426904d0, &
652 0.16729657d0, 0.19477968d0, 0.20424476d0, 0.19795407d0, 0.1850906d0 , &
653 0.16770414d0, 0.14401552d0, 0.1100625d0 , 0.06667179d0, 0.03197211d0, &
654 0.01537663d0, 0.0101183d0 , 0.0098894d0 , 0.01151467d0, 0.01383194d0, &
655 0.01702564d0, 0.02182203d0, 0.02865559d0, 0.03893335d0, 0.058658d0 , &
656 0.10065221d0, 0.17148029d0, 0.25484917d0, 0.33214691d0, 0.40393216d0, &
657 0.48846549d0, 0.60029338d0, 0.71388649d0, 0.79315821d0, 0.83508279d0, &
658 0.85167525d0, 0.85334779d0, 0.84727076d0, 0.83593845d0, 0.81732922d0, &
659 0.78699436d0, 0.73859693d0, 0.67075862d0, 0.60109553d0, 0.55162438d0, &
660 0.52663236d0, 0.52171481d0, 0.53231339d0, 0.55326768d0, 0.57915489d0, &
661 0.6053247d0 , 0.62870688d0, 0.64678595d0, 0.65724692d0, 0.65814188d0, &
662 0.64734694d0, 0.62185287d0, 0.5779309d0 , 0.51564068d0, 0.44378961d0, &
663 0.37428962d0, 0.31473933d0, 0.26719147d0, 0.23069006d0, 0.20302917d0, &
664 0.18205528d0, 0.16592395d0, 0.15317906d0, 0.14282262d0, 0.13400571d0, &
665 0.12609272d0, 0.11858749d0, 0.11115605d0, 0.10365161d0, 0.0960813d0 , &
666 0.08856139d0, 0.0813042d0 , 0.0744821d0 , 0.06825494d0, 0.06266579d0, &
667 0.05773579d0, 0.05343019d0, 0.04968256d0, 0.04643515d0, 0.04360624d0, &
668 0.04114288d0, 0.03899849d0, 0.03712123d0, 0.03547786d0, 0.03404048d0, &
669 0.0327726d0 , 0.03166073d0, 0.03068298d0, 0.02982491d0, 0.02907099d0, &
670 0.02841026d0 /
671
672 data n_colgan / 55 /
673
674 data t_colgan / 4.06460772d0, 4.14229559d0, 4.21995109d0, 4.29760733d0, 4.37527944d0, 4.45293587d0, &
675 4.53060946d0, 4.60826923d0, 4.68592974d0, 4.76359269d0, 4.79704583d0, 4.83049243d0, &
676 4.86394114d0, 4.89738514d0, 4.93083701d0, 4.96428321d0, 4.99773141d0, 5.03116600d0, &
677 5.06460772d0, 5.17574368d0, 5.28683805d0, 5.39795738d0, 5.50906805d0, 5.62017771d0, &
678 5.73129054d0, 5.84240328d0, 5.95351325d0, 6.06460772d0, 6.17574368d0, 6.28683805d0, &
679 6.39795738d0, 6.50906805d0, 6.62017771d0, 6.73129054d0, 6.84240328d0, 6.95351325d0, &
680 7.06460772d0, 7.17574368d0, 7.28683805d0, 7.39795738d0, 7.50906805d0, 7.62017771d0, &
681 7.73129054d0, 7.84240328d0, 7.95351325d0, 8.06460772d0, 8.17574368d0, 8.28683805d0, &
682 8.39795738d0, 8.50906805d0, 8.62017771d0, 8.73129054d0, 8.84240328d0, 8.95351325d0, &
683 9.06460772d0 /
684
685 data l_colgan / -22.18883401d0, -21.78629635d0, -21.60383554d0, -21.68480662d0, -21.76444630d0, &
686 -21.67935529d0, -21.54217864d0, -21.37958284d0, -21.25171892d0, -21.17584161d0, &
687 -21.15783402d0, -21.14491111d0, -21.13526945d0, -21.12837453d0, -21.12485189d0, &
688 -21.12438898d0, -21.12641785d0, -21.12802448d0, -21.12547760d0, -21.08964778d0, &
689 -21.08812360d0, -21.19542445d0, -21.34582346d0, -21.34839251d0, -21.31700703d0, &
690 -21.29072156d0, -21.28900309d0, -21.34104468d0, -21.43122351d0, -21.62448270d0, &
691 -21.86694036d0, -22.02897478d0, -22.08050874d0, -22.06057061d0, -22.01973295d0, &
692 -22.00000434d0, -22.05161149d0, -22.22175466d0, -22.41451671d0, -22.52581288d0, &
693 -22.56913516d0, -22.57485721d0, -22.56150512d0, -22.53968863d0, -22.51490350d0, &
694 -22.48895932d0, -22.46071057d0, -22.42908363d0, -22.39358639d0, -22.35456791d0, &
695 -22.31261375d0, -22.26827428d0, -22.22203698d0, -22.17422996d0, -22.12514145d0 /
696
697 ! Composite cooling curve: three-zone blend of Dere_corona (CHIANTI),
698 ! Colgan (Los Alamos), SPEX (Schure+2009), and cloudy_solar (Cloudy).
699 ! T > 12 kK : weighted average (Dere=6, Colgan=4, SPEX=2)
700 ! 7-12 kK : SPEX CIE through the Lya transition
701 ! T < 7 kK : cloudy_solar (fine-structure, molecular)
702 ! Smooth tanh blends at the zone boundaries; no cliff, no discontinuity.
703 ! Same T grid as cloudy_solar (151 points, log T = 1.0 to 9.0).
704 double precision :: t_composite(1:151), l_composite(1:151)
705 integer :: n_composite
706
707 data n_composite / 151 /
708
709 data t_composite / &
710 1.000d0, 1.050d0, 1.100d0, 1.150d0, &
711 1.200d0, 1.250d0, 1.300d0, 1.350d0, &
712 1.400d0, 1.450d0, 1.500d0, 1.550d0, &
713 1.600d0, 1.650d0, 1.700d0, 1.750d0, &
714 1.800d0, 1.850d0, 1.900d0, 1.950d0, &
715 2.000d0, 2.050d0, 2.100d0, 2.150d0, &
716 2.200d0, 2.250d0, 2.300d0, 2.350d0, &
717 2.400d0, 2.450d0, 2.500d0, 2.550d0, &
718 2.600d0, 2.650d0, 2.700d0, 2.750d0, &
719 2.800d0, 2.850d0, 2.900d0, 2.950d0, &
720 3.000d0, 3.050d0, 3.100d0, 3.150d0, &
721 3.200d0, 3.250d0, 3.300d0, 3.350d0, &
722 3.400d0, 3.450d0, 3.500d0, 3.550d0, &
723 3.600d0, 3.650d0, 3.700d0, 3.750d0, &
724 3.800d0, 3.850d0, 3.900d0, 3.950d0, &
725 4.000d0, 4.050d0, 4.100d0, 4.150d0, &
726 4.200d0, 4.250d0, 4.300d0, 4.350d0, &
727 4.400d0, 4.450d0, 4.500d0, 4.550d0, &
728 4.600d0, 4.650d0, 4.700d0, 4.750d0, &
729 4.800d0, 4.850d0, 4.900d0, 4.950d0, &
730 5.000d0, 5.050d0, 5.100d0, 5.150d0, &
731 5.200d0, 5.250d0, 5.300d0, 5.350d0, &
732 5.400d0, 5.450d0, 5.500d0, 5.550d0, &
733 5.600d0, 5.650d0, 5.700d0, 5.750d0, &
734 5.800d0, 5.850d0, 5.900d0, 5.950d0, &
735 6.000d0, 6.050d0, 6.100d0, 6.150d0, &
736 6.200d0, 6.250d0, 6.300d0, 6.350d0, &
737 6.400d0, 6.450d0, 6.500d0, 6.550d0, &
738 6.600d0, 6.650d0, 6.700d0, 6.750d0, &
739 6.800d0, 6.850d0, 6.900d0, 6.950d0, &
740 7.000d0, 7.050d0, 7.100d0, 7.150d0, &
741 7.200d0, 7.250d0, 7.300d0, 7.350d0, &
742 7.400d0, 7.450d0, 7.500d0, 7.550d0, &
743 7.600d0, 7.650d0, 7.700d0, 7.750d0, &
744 7.800d0, 7.850d0, 7.900d0, 7.950d0, &
745 8.000d0, 8.100d0, 8.200d0, 8.300d0, &
746 8.400d0, 8.500d0, 8.600d0, 8.700d0, &
747 8.800d0, 8.900d0, 9.000d0 /
748
749 data l_composite / &
750 -28.375d0, -28.251d0, -28.137d0, -28.029d0, &
751 -27.929d0, -27.834d0, -27.745d0, -27.662d0, &
752 -27.584d0, -27.512d0, -27.445d0, -27.383d0, &
753 -27.326d0, -27.273d0, -27.223d0, -27.175d0, &
754 -27.128d0, -27.079d0, -27.027d0, -26.971d0, &
755 -26.911d0, -26.846d0, -26.777d0, -26.705d0, &
756 -26.631d0, -26.554d0, -26.479d0, -26.407d0, &
757 -26.338d0, -26.274d0, -26.213d0, -26.156d0, &
758 -26.101d0, -26.049d0, -25.999d0, -25.949d0, &
759 -25.901d0, -25.852d0, -25.803d0, -25.754d0, &
760 -25.707d0, -25.662d0, -25.621d0, -25.588d0, &
761 -25.561d0, -25.538d0, -25.518d0, -25.497d0, &
762 -25.475d0, -25.451d0, -25.426d0, -25.400d0, &
763 -25.371d0, -25.333d0, -25.295d0, -25.261d0, &
764 -25.228d0, -25.098d0, -24.607d0, -23.796d0, &
765 -23.008d0, -22.456d0, -22.066d0, -21.760d0, &
766 -21.604d0, -21.584d0, -21.651d0, -21.736d0, &
767 -21.810d0, -21.872d0, -21.894d0, -21.860d0, &
768 -21.795d0, -21.715d0, -21.622d0, -21.520d0, &
769 -21.419d0, -21.335d0, -21.280d0, -21.252d0, &
770 -21.246d0, -21.253d0, -21.275d0, -21.270d0, &
771 -21.244d0, -21.217d0, -21.196d0, -21.183d0, &
772 -21.193d0, -21.250d0, -21.350d0, -21.442d0, &
773 -21.482d0, -21.476d0, -21.453d0, -21.435d0, &
774 -21.426d0, -21.419d0, -21.410d0, -21.403d0, &
775 -21.407d0, -21.424d0, -21.449d0, -21.478d0, &
776 -21.512d0, -21.562d0, -21.638d0, -21.739d0, &
777 -21.850d0, -21.952d0, -22.033d0, -22.087d0, &
778 -22.116d0, -22.123d0, -22.116d0, -22.101d0, &
779 -22.083d0, -22.068d0, -22.058d0, -22.057d0, &
780 -22.070d0, -22.103d0, -22.162d0, -22.241d0, &
781 -22.328d0, -22.406d0, -22.470d0, -22.517d0, &
782 -22.549d0, -22.569d0, -22.579d0, -22.583d0, &
783 -22.581d0, -22.575d0, -22.565d0, -22.554d0, &
784 -22.541d0, -22.525d0, -22.505d0, -22.480d0, &
785 -22.465d0, -22.415d0, -22.365d0, -22.315d0, &
786 -22.265d0, -22.215d0, -22.165d0, -22.115d0, &
787 -22.065d0, -22.015d0, -21.965d0 /
788
789end module mod_radloss_tables
mod_radloss_tables
module containing all optically thin radiative cooling tables
Definition mod_radloss_tables.t:7
mod_radloss_tables::l_dere_corona
double precision, dimension(1:101) l_dere_corona
Definition mod_radloss_tables.t:87
mod_radloss_tables::t_mlsolar1
double precision, dimension(1:71) t_mlsolar1
Definition mod_radloss_tables.t:82
mod_radloss_tables::n_mlcosmol
integer n_mlcosmol
Definition mod_radloss_tables.t:97
mod_radloss_tables::l_cl_solar
double precision, dimension(1:151) l_cl_solar
Definition mod_radloss_tables.t:87
mod_radloss_tables::t_fm
double precision, dimension(1:5) t_fm
Definition mod_radloss_tables.t:14
mod_radloss_tables::a_spex_dm_fine
double precision, dimension(1:14) a_spex_dm_fine
Definition mod_radloss_tables.t:20
mod_radloss_tables::n_cl_ism
integer n_cl_ism
Definition mod_radloss_tables.t:97
mod_radloss_tables::a_rosner
double precision, dimension(1:9) a_rosner
Definition mod_radloss_tables.t:20
mod_radloss_tables::l_spex
double precision, dimension(1:110) l_spex
Definition mod_radloss_tables.t:87
mod_radloss_tables::l_mb
double precision, dimension(1:51) l_mb
Definition mod_radloss_tables.t:87
mod_radloss_tables::n_composite
integer n_composite
Definition mod_radloss_tables.t:705
mod_radloss_tables::t_rosner
double precision, dimension(1:10) t_rosner
Definition mod_radloss_tables.t:14
mod_radloss_tables::a_hildner
double precision, dimension(1:5) a_hildner
Definition mod_radloss_tables.t:20
mod_radloss_tables::x_rosner
double precision, dimension(1:9) x_rosner
Definition mod_radloss_tables.t:17
mod_radloss_tables::x_klimchuk
double precision, dimension(1:7) x_klimchuk
Definition mod_radloss_tables.t:17
mod_radloss_tables::l_composite
double precision, dimension(1:151) l_composite
Definition mod_radloss_tables.t:704
mod_radloss_tables::l_cl_ism
double precision, dimension(1:151) l_cl_ism
Definition mod_radloss_tables.t:87
mod_radloss_tables::t_spex_dm_rough
double precision, dimension(1:8) t_spex_dm_rough
Definition mod_radloss_tables.t:14
mod_radloss_tables::n_colgan
integer n_colgan
Definition mod_radloss_tables.t:97
mod_radloss_tables::nenh_spex
double precision, dimension(1:110) nenh_spex
Definition mod_radloss_tables.t:95
mod_radloss_tables::t_spex
double precision, dimension(1:110) t_spex
Definition mod_radloss_tables.t:82
mod_radloss_tables::l_dm_2
double precision, dimension(1:76) l_dm_2
Definition mod_radloss_tables.t:87
mod_radloss_tables::t_composite
double precision, dimension(1:151) t_composite
Definition mod_radloss_tables.t:704
mod_radloss_tables::t_spex_dm_fine
double precision, dimension(1:15) t_spex_dm_fine
Definition mod_radloss_tables.t:14
mod_radloss_tables::n_rosner
integer n_rosner
Definition mod_radloss_tables.t:23
mod_radloss_tables::n_spex_dm_fine
integer n_spex_dm_fine
Definition mod_radloss_tables.t:23
mod_radloss_tables::x_spex_dm_rough
double precision, dimension(1:7) x_spex_dm_rough
Definition mod_radloss_tables.t:17
mod_radloss_tables::x_spex_dm_fine
double precision, dimension(1:14) x_spex_dm_fine
Definition mod_radloss_tables.t:17
mod_radloss_tables::l_mlsolar1
double precision, dimension(1:71) l_mlsolar1
Definition mod_radloss_tables.t:87
mod_radloss_tables::n_mlsolar1
integer n_mlsolar1
Definition mod_radloss_tables.t:97
mod_radloss_tables::t_jccorona
double precision, dimension(1:45) t_jccorona
Definition mod_radloss_tables.t:82
mod_radloss_tables::n_mb
integer n_mb
Definition mod_radloss_tables.t:97
mod_radloss_tables::n_hildner
integer n_hildner
Definition mod_radloss_tables.t:23
mod_radloss_tables::x_hildner
double precision, dimension(1:5) x_hildner
Definition mod_radloss_tables.t:17
mod_radloss_tables::t_mlcosmol
double precision, dimension(1:71) t_mlcosmol
Definition mod_radloss_tables.t:82
mod_radloss_tables::t_cl_ism
double precision, dimension(1:151) t_cl_ism
Definition mod_radloss_tables.t:82
mod_radloss_tables::t_cl_solar
double precision, dimension(1:151) t_cl_solar
Definition mod_radloss_tables.t:82
mod_radloss_tables::t_mb
double precision, dimension(1:51) t_mb
Definition mod_radloss_tables.t:82
mod_radloss_tables::t_klimchuk
double precision, dimension(1:8) t_klimchuk
Definition mod_radloss_tables.t:14
mod_radloss_tables::t_colgan
double precision, dimension(1:55) t_colgan
Definition mod_radloss_tables.t:82
mod_radloss_tables::l_colgan
double precision, dimension(1:55) l_colgan
Definition mod_radloss_tables.t:87
mod_radloss_tables::a_fm
double precision, dimension(1:4) a_fm
Definition mod_radloss_tables.t:20
mod_radloss_tables::n_dm
integer n_dm
Definition mod_radloss_tables.t:97
mod_radloss_tables::n_klimchuk
integer n_klimchuk
Definition mod_radloss_tables.t:23
mod_radloss_tables::l_dere_photo
double precision, dimension(1:101) l_dere_photo
Definition mod_radloss_tables.t:87
mod_radloss_tables::n_spex_dm_rough
integer n_spex_dm_rough
Definition mod_radloss_tables.t:23
mod_radloss_tables::n_dere
integer n_dere
Definition mod_radloss_tables.t:97
mod_radloss_tables::l_jccorona
double precision, dimension(1:45) l_jccorona
Definition mod_radloss_tables.t:87
mod_radloss_tables::l_mlwc
double precision, dimension(1:71) l_mlwc
Definition mod_radloss_tables.t:87
mod_radloss_tables::n_cl_solar
integer n_cl_solar
Definition mod_radloss_tables.t:97
mod_radloss_tables::l_dm
double precision, dimension(1:71) l_dm
Definition mod_radloss_tables.t:87
mod_radloss_tables::t_mlwc
double precision, dimension(1:71) t_mlwc
Definition mod_radloss_tables.t:82
mod_radloss_tables::a_spex_dm_rough
double precision, dimension(1:7) a_spex_dm_rough
Definition mod_radloss_tables.t:20
mod_radloss_tables::n_spex
integer n_spex
Definition mod_radloss_tables.t:97
mod_radloss_tables::t_dm
double precision, dimension(1:71) t_dm
Definition mod_radloss_tables.t:82
mod_radloss_tables::n_fm
integer n_fm
Definition mod_radloss_tables.t:23
mod_radloss_tables::a_klimchuk
double precision, dimension(1:7) a_klimchuk
Definition mod_radloss_tables.t:20
mod_radloss_tables::l_mlcosmol
double precision, dimension(1:71) l_mlcosmol
Definition mod_radloss_tables.t:87
mod_radloss_tables::t_dm_2
double precision, dimension(1:76) t_dm_2
Definition mod_radloss_tables.t:82
mod_radloss_tables::n_jccorona
integer n_jccorona
Definition mod_radloss_tables.t:97
mod_radloss_tables::n_mlwc
integer n_mlwc
Definition mod_radloss_tables.t:97
mod_radloss_tables::t_hildner
double precision, dimension(1:6) t_hildner
Definition mod_radloss_tables.t:14
mod_radloss_tables::x_fm
double precision, dimension(1:4) x_fm
Definition mod_radloss_tables.t:17
mod_radloss_tables::n_dm_2
integer n_dm_2
Definition mod_radloss_tables.t:97
mod_radloss_tables::t_dere
double precision, dimension(1:101) t_dere
Definition mod_radloss_tables.t:82
mod_radloss_tables::lowfip_frac
double precision, dimension(1:101) lowfip_frac
Definition mod_radloss_tables.t:93
Generated by doxygen 1.15.0