In order to compute the synthetic spectrum, TARDIS must either be told or must calculate many physical properties of the model. To understand and test the code it can be important to look at these values. One easy way to do this is to run TARDIS in an interactive mode and then inspect the model properties.
In [1]:
# Download the atomic data
from tardis.io.atom_data.util import download_atom_data
download_atom_data('kurucz_cd23_chianti_H_He')
# Download the example configuration file
!curl -O https://raw.githubusercontent.com/tardis-sn/tardis/master/docs/models/examples/tardis_example.yml
/Users/sashmish/anaconda3/envs/tardis/lib/python3.6/importlib/_bootstrap.py:219: QAWarning: pyne.data is not yet QA compliant.
return f(*args, **kwds)
[tardis.io.atom_data.atom_web_download][INFO ] Downloading atomic data from https://media.githubusercontent.com/media/tardis-sn/tardis-refdata/master/atom_data/kurucz_cd23_chianti_H_He.h5 to /Users/sashmish/Downloads/tardis-data/kurucz_cd23_chianti_H_He.h5 (atom_web_download.py:47)
% Total % Received % Xferd Average Speed Time Time Time Current
Dload Upload Total Spent Left Speed
100 980 100 980 0 0 2222 0 --:--:-- --:--:-- --:--:-- 2222
In [2]:
from tardis import run_tardis
simulation = run_tardis('tardis_example.yml')
[tardis.plasma.standard_plasmas][INFO ] Reading Atomic Data from kurucz_cd23_chianti_H_He.h5 (standard_plasmas.py:74)
[tardis.io.atom_data.util][INFO ] Atom Data kurucz_cd23_chianti_H_He.h5 not found in local path. Exists in TARDIS Data repo /Users/sashmish/Downloads/tardis-data/kurucz_cd23_chianti_H_He.h5 (util.py:29)
[tardis.io.atom_data.base][INFO ] Read Atom Data with UUID=6f7b09e887a311e7a06b246e96350010 and MD5=864f1753714343c41f99cb065710cace. (base.py:184)
[tardis.io.atom_data.base][INFO ] Non provided atomic data: synpp_refs, photoionization_data (base.py:187)
[tardis.simulation.base][INFO ] Starting iteration 1/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 7.93730e+42 erg / s Luminosity absorbed = 2.66400e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 9926.501965 10168.422346 0.400392 0.501001
5 9852.611678 10367.510165 0.211205 0.186556
10 9779.813302 10222.882075 0.142695 0.113985
15 9708.082813 9928.649533 0.104556 0.085217
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 9933.952 K -- next t_inner 11475.997 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 2/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.41371e+43 erg / s Luminosity absorbed = 4.73752e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 10168.422346 11497.903617 0.501001 0.545136
5 10367.510165 11621.417600 0.186556 0.212874
10 10222.882075 11335.692704 0.113985 0.135635
15 9928.649533 11018.841465 0.085217 0.101203
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 11475.997 K -- next t_inner 9933.797 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 3/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 8.04734e+42 erg / s Luminosity absorbed = 2.55640e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 11497.903617 10477.653343 0.545136 0.439296
5 11621.417600 10856.518852 0.212874 0.159058
10 11335.692704 10635.474976 0.135635 0.100253
15 11018.841465 10290.254692 0.101203 0.075529
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 9933.797 K -- next t_inner 11397.084 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 4/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.38199e+43 erg / s Luminosity absorbed = 4.54590e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 10477.653343 11608.211738 0.439296 0.506501
5 10856.518852 11859.259666 0.159058 0.192802
10 10635.474976 11479.604579 0.100253 0.125779
15 10290.254692 11149.192411 0.075529 0.094052
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 11397.084 K -- next t_inner 9978.058 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 5/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 8.16431e+42 erg / s Luminosity absorbed = 2.62454e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 11608.211738 10556.913414 0.506501 0.440229
5 11859.259666 11021.144022 0.192802 0.153136
10 11479.604579 10713.104080 0.125779 0.099268
15 11149.192411 10399.079825 0.094052 0.073482
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 9978.058 K -- next t_inner 11365.561 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 6/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.36136e+43 erg / s Luminosity absorbed = 4.54590e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 10556.913414 11562.006340 0.440229 0.509603
5 11021.144022 11799.980145 0.153136 0.194437
10 10713.104080 11482.742597 0.099268 0.124721
15 10399.079825 11124.799749 0.073482 0.093597
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 11365.561 K -- next t_inner 10025.564 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 7/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 8.33278e+42 erg / s Luminosity absorbed = 2.66283e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 11562.006340 10655.549804 0.509603 0.430187
5 11799.980145 11063.665621 0.194437 0.155127
10 11482.742597 10926.058943 0.124721 0.094642
15 11124.799749 10477.041314 0.093597 0.072578
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 10025.564 K -- next t_inner 11303.643 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 8/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.34650e+43 erg / s Luminosity absorbed = 4.30670e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 10655.549804 11572.119736 0.430187 0.496558
5 11063.665621 11696.304249 0.155127 0.198228
10 10926.058943 11482.609682 0.094642 0.124382
15 10477.041314 11098.080451 0.072578 0.093537
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 11303.643 K -- next t_inner 10025.810 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 9/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 8.36983e+42 erg / s Luminosity absorbed = 2.63049e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 11572.119736 10590.775937 0.496558 0.439574
5 11696.304249 11060.842307 0.198228 0.154917
10 11482.609682 10761.665022 0.124382 0.099840
15 11098.080451 10436.629671 0.093537 0.074676
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 10025.810 K -- next t_inner 11278.874 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 10/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.32624e+43 erg / s Luminosity absorbed = 4.34986e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 10590.775937 11489.769556 0.439574 0.511403
5 11060.842307 11791.517605 0.154917 0.191127
10 10761.665022 11527.732208 0.099840 0.120106
15 10436.629671 11111.114390 0.074676 0.091713
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 11278.874 K -- next t_inner 10079.967 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 11/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 8.51915e+42 erg / s Luminosity absorbed = 2.72009e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 11489.769556 10664.234475 0.511403 0.437272
5 11791.517605 11091.518839 0.191127 0.156638
10 11527.732208 10885.771059 0.120106 0.096973
15 11111.114390 10484.177251 0.091713 0.074261
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 10079.967 K -- next t_inner 11239.980 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 12/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.30886e+43 erg / s Luminosity absorbed = 4.28617e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 10664.234475 11510.842068 0.437272 0.495944
5 11091.518839 11716.009425 0.156638 0.189856
10 10885.771059 11479.636914 0.096973 0.120067
15 10484.177251 11117.963837 0.074261 0.090628
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 11239.980 K -- next t_inner 10111.668 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 13/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 8.65047e+42 erg / s Luminosity absorbed = 2.72874e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 11510.842068 10636.025493 0.495944 0.446635
5 11716.009425 11175.711706 0.189856 0.154916
10 11479.636914 10916.188023 0.120067 0.097843
15 11117.963837 10431.720831 0.090628 0.077032
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 10111.668 K -- next t_inner 11189.420 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 14/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.28712e+43 erg / s Luminosity absorbed = 4.19242e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 10636.025493 11423.118214 0.446635 0.503706
5 11175.711706 11682.224028 0.154916 0.189577
10 10916.188023 11389.994295 0.097843 0.121114
15 10431.720831 11029.450833 0.077032 0.091350
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 11189.420 K -- next t_inner 10150.866 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 15/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 8.76501e+42 erg / s Luminosity absorbed = 2.79433e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 11423.118214 10634.831735 0.503706 0.452948
5 11682.224028 11112.411355 0.189577 0.159316
10 11389.994295 10831.887277 0.121114 0.101805
15 11029.450833 10396.372864 0.091350 0.078776
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 10150.866 K -- next t_inner 11159.164 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 16/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.27462e+43 erg / s Luminosity absorbed = 4.13373e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 10634.831735 11415.925493 0.452948 0.501153
5 11112.411355 11688.667754 0.159316 0.188430
10 10831.887277 11397.519913 0.101805 0.120134
15 10396.372864 11045.200458 0.078776 0.090180
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 11159.164 K -- next t_inner 10172.913 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 17/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 8.86785e+42 erg / s Luminosity absorbed = 2.79094e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 11415.925493 10678.510651 0.501153 0.453127
5 11688.667754 11152.822099 0.188430 0.159726
10 11397.519913 10907.112682 0.120134 0.099725
15 11045.200458 10508.816973 0.090180 0.076667
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 10172.913 K -- next t_inner 11118.359 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 18/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.24749e+43 erg / s Luminosity absorbed = 4.15445e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 10678.510651 11482.446420 0.453127 0.481917
5 11152.822099 11693.549165 0.159726 0.184456
10 10907.112682 11472.013555 0.099725 0.115208
15 10508.816973 10997.309982 0.076667 0.089730
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 11118.359 K -- next t_inner 10245.337 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 19/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 9.09476e+42 erg / s Luminosity absorbed = 2.89713e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 11482.446420 10834.729769 0.481917 0.438322
5 11693.549165 11346.133670 0.184456 0.154204
10 11472.013555 10996.421420 0.115208 0.099677
15 10997.309982 10625.479432 0.089730 0.075531
(base.py:350)
[tardis.simulation.base][INFO ] t_inner 10245.337 K -- next t_inner 11056.949 K (base.py:352)
[tardis.simulation.base][INFO ] Starting iteration 20/20 (base.py:268)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.22623e+43 erg / s Luminosity absorbed = 4.00603e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:359)
[tardis.simulation.base][INFO ] Simulation finished in 20 iterations and took 23.27 s (base.py:308)
If all goes well, the simulation should run as usual. Afterwards, the
information from the simulation will all exist in Simulation and
can be examined.
Some examples for useful/interesting quantities are
given below (but much more information is available: contact us via
tardis-sn-users if you need
further help).
For example, two of our important quantities are the parameters of the
radiation field model, $T_{\rm rad}$ and $W$. These exist as numpy.ndarray
Thus simulation.plasma.t_rad will give you a list of the $T_{\rm rad}$-values for the model zones in cgs units.
In [3]:
simulation.plasma.t_rad
Out[3]:
array([10834.7297692 , 11024.52012206, 11134.85912792, 11299.31893687,
11304.06901637, 11346.13367044, 11273.8729839 , 11201.3282228 ,
11080.35051086, 11053.1815164 , 10996.42142 , 10906.50207322,
10852.60134857, 10768.23991178, 10701.07288864, 10625.47943241,
10578.25968193, 10462.71630881, 10328.2464022 , 10198.53837217])
Similarly, the $W$-values can be accessed using simulation.plasma.w
In [4]:
simulation.plasma.w
Out[4]:
array([0.43832195, 0.32677702, 0.26128919, 0.20770296, 0.17725587,
0.15420374, 0.1387194 , 0.1264399 , 0.11830376, 0.10744254,
0.09967675, 0.09360326, 0.08756804, 0.08325674, 0.07930315,
0.07553141, 0.07138798, 0.06935054, 0.06784402, 0.06650938])
Several important quantities that were setup when the model was defined by the configuration file are located in the model section of the simulation. For example, the inner and outer velocity boundaries of the zones in the model is given by simulation.model.v_inner.cgs and simulation.model.v_outer.cgs respectively. These exist as Astropy Quantities.
In [5]:
simulation.model.v_inner.cgs
Out[5]:
$[1.1 \times 10^{9},~1.145 \times 10^{9},~1.19 \times 10^{9},~1.235 \times 10^{9},~1.28 \times 10^{9},~1.325 \times 10^{9},~1.37 \times 10^{9},~1.415 \times 10^{9},~1.46 \times 10^{9},~1.505 \times 10^{9},~1.55 \times 10^{9},~1.595 \times 10^{9},~1.64 \times 10^{9},~1.685 \times 10^{9},~1.73 \times 10^{9},~1.775 \times 10^{9},~1.82 \times 10^{9},~1.865 \times 10^{9},~1.91 \times 10^{9},~1.955 \times 10^{9}] \; \mathrm{\frac{cm}{s}}$
In [6]:
simulation.model.v_outer.cgs
Out[6]:
$[1.145 \times 10^{9},~1.19 \times 10^{9},~1.235 \times 10^{9},~1.28 \times 10^{9},~1.325 \times 10^{9},~1.37 \times 10^{9},~1.415 \times 10^{9},~1.46 \times 10^{9},~1.505 \times 10^{9},~1.55 \times 10^{9},~1.595 \times 10^{9},~1.64 \times 10^{9},~1.685 \times 10^{9},~1.73 \times 10^{9},~1.775 \times 10^{9},~1.82 \times 10^{9},~1.865 \times 10^{9},~1.91 \times 10^{9},~1.955 \times 10^{9},~2 \times 10^{9}] \; \mathrm{\frac{cm}{s}}$
The average density in the zones is given by simulation.model.density.cgs. These also exist as Astropy Quantities.
In [7]:
simulation.model.density.cgs
Out[7]:
$[7.5428036 \times 10^{-14},~5.728475 \times 10^{-14},~4.3960742 \times 10^{-14},~3.4062874 \times 10^{-14},~2.6631346 \times 10^{-14},~2.0995965 \times 10^{-14},~1.6682872 \times 10^{-14},~1.3353105 \times 10^{-14},~1.0761538 \times 10^{-14},~8.7290848 \times 10^{-15},~7.1236516 \times 10^{-15},~5.8469209 \times 10^{-15},~4.8250928 \times 10^{-15},~4.0023242 \times 10^{-15},~3.3360386 \times 10^{-15},~2.7935404 \times 10^{-15},~2.3495504 \times 10^{-15},~1.9843968 \times 10^{-15},~1.6826769 \times 10^{-15},~1.4322598 \times 10^{-15}] \; \mathrm{\frac{g}{cm^{3}}}$
Many other interesting quantities are stored in the plasma.
For example the calculated ion populations and level populations is given by simulation.plasma.ion_number_density and simulation.plasma.level_number_density respectively.
In [8]:
simulation.plasma.ion_number_density
Out[8]:
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
atomic_number
ion_number
8
0
1.380280e+03
6.045221e+02
3.046902e+02
1.456868e+02
8.855874e+01
5.201972e+01
3.629170e+01
2.573196e+01
1.983458e+01
1.357656e+01
9.817542e+00
7.543514e+00
5.565623e+00
4.345404e+00
3.343570e+00
2.633578e+00
2.005298e+00
1.715762e+00
1.531726e+00
1.374292e+00
1
5.393290e+08
4.094843e+08
3.141101e+08
2.431228e+08
1.899693e+08
1.495973e+08
1.188702e+08
9.515165e+07
7.673358e+07
6.221905e+07
5.077488e+07
4.169038e+07
3.440446e+07
2.854724e+07
2.379873e+07
1.993347e+07
1.676551e+07
1.416739e+07
1.201946e+07
1.023437e+07
2
9.780989e+04
1.904997e+05
2.775463e+05
4.797901e+05
4.863297e+05
5.565803e+05
4.383339e+05
3.438730e+05
2.281928e+05
2.075235e+05
1.703380e+05
1.242273e+05
1.024814e+05
7.562351e+04
5.915039e+04
4.470640e+04
3.743945e+04
2.414008e+04
1.430706e+04
8.524308e+03
3
3.252281e-09
2.343311e-08
7.976996e-08
4.162940e-07
5.521097e-07
9.904266e-07
6.786694e-07
4.570699e-07
1.993626e-07
1.932109e-07
1.432370e-07
7.806358e-08
5.796411e-08
3.221132e-08
2.066853e-08
1.209496e-08
9.155287e-09
3.541742e-09
1.101099e-09
3.457566e-10
4
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
5
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
6
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
7
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
8
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
12
0
1.346333e-06
3.950660e-07
1.426881e-07
4.794894e-08
2.277172e-08
1.031969e-08
5.976987e-09
3.547825e-09
2.376965e-09
1.345053e-09
8.243636e-10
5.520734e-10
3.490056e-10
2.397691e-10
1.614181e-10
1.125574e-10
7.472232e-11
5.898938e-11
4.971474e-11
4.231294e-11
1
6.084662e+02
2.608497e+02
1.299483e+02
6.113372e+01
3.716500e+01
2.176241e+01
1.529509e+01
1.092661e+01
8.527359e+00
5.856457e+00
4.261848e+00
3.307046e+00
2.455226e+00
1.935492e+00
1.501131e+00
1.193128e+00
9.138621e-01
7.931379e-01
7.204117e-01
6.577398e-01
2
5.606550e+07
4.257984e+07
3.267615e+07
2.531908e+07
1.979520e+07
1.560640e+07
1.240046e+07
9.925432e+06
7.999107e+06
6.488375e+06
5.295049e+06
4.346048e+06
3.586518e+06
2.974950e+06
2.479696e+06
2.076454e+06
1.746434e+06
1.475013e+06
1.250743e+06
1.064607e+06
3
4.120107e-17
1.847504e-16
4.315753e-16
1.483199e-15
1.533951e-15
2.087760e-15
1.222248e-15
7.092049e-16
2.819936e-16
2.283522e-16
1.467030e-16
7.216755e-17
4.688294e-17
2.368231e-17
1.364184e-17
7.273552e-18
4.886724e-18
1.821574e-18
5.619339e-19
1.755039e-19
4
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
5
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
6
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
7
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
8
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
9
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
10
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
11
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
12
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
14
0
7.121410e-04
1.921013e-04
6.599549e-05
2.054382e-05
9.708417e-06
4.296360e-06
2.567117e-06
1.572651e-06
1.113705e-06
6.364852e-07
3.997533e-07
2.788974e-07
1.804863e-07
1.288509e-07
8.939595e-08
6.450650e-08
4.371505e-08
3.639895e-08
3.262155e-08
2.943795e-08
1
1.057234e+05
4.409044e+04
2.159442e+04
9.886849e+03
5.989365e+03
3.468313e+03
2.460782e+03
1.775245e+03
1.412253e+03
9.716312e+02
7.123876e+02
5.606352e+02
4.193053e+02
3.350009e+02
2.624996e+02
2.111578e+02
1.628349e+02
1.441288e+02
1.339320e+02
1.249635e+02
2
8.401073e+08
6.371392e+08
4.879292e+08
3.760672e+08
2.931690e+08
2.298363e+08
1.826186e+08
1.461958e+08
1.181696e+08
9.566890e+07
7.805477e+07
6.417624e+07
5.295412e+07
4.399297e+07
3.669628e+07
3.076475e+07
2.587415e+07
2.191312e+07
1.863251e+07
1.589078e+07
3
8.192553e+05
1.549018e+06
2.217086e+06
3.728327e+06
3.767905e+06
4.267904e+06
3.395076e+06
2.691213e+06
1.821504e+06
1.660538e+06
1.374144e+06
1.017116e+06
8.458608e+05
6.331180e+05
5.007190e+05
3.833675e+05
3.234797e+05
2.130184e+05
1.294170e+05
7.897841e+04
4
7.678176e-04
4.490784e-03
1.357295e-02
5.947951e-02
7.832168e-02
1.340207e-01
9.884805e-02
7.175832e-02
3.562515e-02
3.548489e-02
2.795368e-02
1.682397e-02
1.325599e-02
8.104905e-03
5.615334e-03
3.587664e-03
2.868875e-03
1.275024e-03
4.675402e-04
1.727755e-04
5
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
6
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
7
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
18
10
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
11
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
12
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
13
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
14
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
15
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
16
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
17
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
18
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
20
0
1.332681e-08
4.248991e-09
1.606786e-09
5.764976e-10
2.742923e-10
1.263310e-10
7.115362e-11
4.104046e-11
2.617119e-11
1.464219e-11
8.760639e-12
5.642345e-12
3.482507e-12
2.302501e-12
1.502345e-12
1.010506e-12
6.556032e-13
4.885474e-13
3.839339e-13
3.045743e-13
1
2.335159e+01
1.077167e+01
5.593754e+00
2.795791e+00
1.702578e+00
1.012249e+00
6.930320e-01
4.820902e-01
3.596454e-01
2.444799e-01
1.741170e-01
1.305166e-01
9.488691e-02
7.235695e-02
5.463538e-02
4.212014e-02
3.164596e-02
2.618354e-02
2.246849e-02
1.939451e-02
2
3.400155e+07
2.582290e+07
1.981669e+07
1.535491e+07
1.200492e+07
9.464593e+06
7.520331e+06
6.019333e+06
4.851103e+06
3.934910e+06
3.211211e+06
2.635684e+06
2.175063e+06
1.804174e+06
1.503825e+06
1.259276e+06
1.059134e+06
8.945295e+05
7.585197e+05
6.456363e+05
3
7.014243e-04
1.836190e-03
3.163463e-03
6.982618e-03
7.131144e-03
8.685252e-03
6.156902e-03
4.339911e-03
2.400688e-03
2.095729e-03
1.577122e-03
1.000142e-03
7.580963e-04
4.890214e-04
3.431818e-04
2.291922e-04
1.775305e-04
9.424970e-05
4.431810e-05
2.100575e-05
4
0.000000e+00
0.000000e+00
0.000000e+00
4.670974e-20
6.275559e-20
1.256827e-19
7.138208e-20
3.972085e-20
1.252727e-20
1.128297e-20
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
5
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
6
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
7
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
8
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
9
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
10
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
11
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
12
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
13
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
14
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
15
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
16
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
17
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
18
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
19
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
20
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
94 rows × 20 columns
In [9]:
simulation.plasma.level_number_density
Out[9]:
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
atomic_number
ion_number
level_number
8
0
0
724.299083
316.384041
159.218359
75.955544
46.168086
27.103405
18.927817
13.433996
10.372597
7.102624
5.140143
3.954488
2.919823
2.282355e+00
1.757802e+00
1.385998e+00
1.056041e+00
9.050186e-01
8.094569e-01
7.275697e-01
1
425.541419
185.949751
93.597245
44.664108
27.148431
15.938929
11.129610
7.898195
6.096966
4.174679
3.020879
2.323669
1.715519
1.340759e+00
1.032475e+00
8.139671e-01
6.201311e-01
5.313209e-01
4.750836e-01
4.269029e-01
2
140.558746
61.429907
30.923305
14.758344
8.970679
5.266885
3.677480
2.609598
2.014270
1.379169
0.997947
0.767567
0.566654
4.428354e-01
3.409936e-01
2.688096e-01
2.047875e-01
1.754414e-01
1.568526e-01
1.409282e-01
3
88.063868
39.888569
20.489867
10.070845
6.126565
3.623688
2.498195
1.749989
1.321457
0.900295
0.644631
0.487520
0.356241
2.739129e-01
2.081708e-01
1.616668e-01
1.220039e-01
1.020938e-01
8.875591e-02
7.756551e-02
4
1.629618
0.769011
0.404288
0.205522
0.125148
0.074651
0.050720
0.035007
0.025778
0.017462
0.012354
0.009164
0.006618
4.994795e-03
3.739352e-03
2.854639e-03
2.131074e-03
1.735927e-03
1.461477e-03
1.237289e-03
5
0.040328
0.020852
0.011544
0.006327
0.003861
0.002347
0.001544
0.001031
0.000718
0.000480
0.000331
0.000235
0.000165
1.196208e-04
8.660153e-05
6.363173e-05
4.636896e-05
3.556994e-05
2.787792e-05
2.198773e-05
6
0.016188
0.008428
0.004684
0.002582
0.001576
0.000959
0.000629
0.000419
0.000291
0.000194
0.000133
0.000095
0.000066
4.789850e-05
3.458903e-05
2.534132e-05
1.843267e-05
1.407569e-05
1.097219e-05
8.607646e-06
7
0.004388
0.002336
0.001315
0.000738
0.000451
0.000276
0.000179
0.000119
0.000081
0.000054
0.000037
0.000026
0.000018
1.287856e-05
9.223661e-06
6.694378e-06
4.840483e-06
3.642144e-06
2.789568e-06
2.150620e-06
8
0.007311
0.003893
0.002191
0.001230
0.000751
0.000459
0.000299
0.000198
0.000135
0.000090
0.000061
0.000043
0.000030
2.145848e-05
1.536859e-05
1.115424e-05
8.065254e-06
6.068553e-06
4.647971e-06
3.583345e-06
9
0.010231
0.005448
0.003066
0.001722
0.001051
0.000643
0.000418
0.000277
0.000189
0.000126
0.000086
0.000060
0.000042
3.002712e-05
2.150540e-05
1.560817e-05
1.128571e-05
8.491681e-06
6.503828e-06
5.014082e-06
10
0.003362
0.001799
0.001015
0.000572
0.000349
0.000214
0.000139
0.000092
0.000063
0.000042
0.000028
0.000020
0.000014
9.852180e-06
7.044304e-06
5.102842e-06
3.685227e-06
2.764553e-06
2.109822e-06
1.620802e-06
11
0.005603
0.002997
0.001692
0.000953
0.000582
0.000356
0.000232
0.000153
0.000104
0.000069
0.000047
0.000033
0.000023
1.641907e-05
1.173962e-05
8.504093e-06
6.141578e-06
4.607234e-06
3.516097e-06
2.701124e-06
12
0.001121
0.000599
0.000338
0.000191
0.000116
0.000071
0.000046
0.000031
0.000021
0.000014
0.000009
0.000007
0.000005
3.283746e-06
2.347875e-06
1.700782e-06
1.228289e-06
9.214269e-07
7.032039e-07
5.402128e-07
13
0.002258
0.001227
0.000698
0.000399
0.000244
0.000150
0.000097
0.000063
0.000043
0.000028
0.000019
0.000013
0.000009
6.578356e-06
4.676593e-06
3.365578e-06
2.420553e-06
1.797254e-06
1.354899e-06
1.028308e-06
14
0.001226
0.000668
0.000380
0.000218
0.000133
0.000082
0.000053
0.000035
0.000023
0.000015
0.000011
0.000007
0.000005
3.571458e-06
2.537378e-06
1.824754e-06
1.311786e-06
9.729029e-07
7.324624e-07
5.551699e-07
15
0.003139
0.001713
0.000978
0.000560
0.000342
0.000210
0.000136
0.000089
0.000060
0.000040
0.000027
0.000019
0.000013
9.132513e-06
6.481784e-06
4.656043e-06
3.344735e-06
2.476215e-06
1.860264e-06
1.407004e-06
16
0.002442
0.001333
0.000761
0.000436
0.000266
0.000164
0.000106
0.000069
0.000047
0.000031
0.000021
0.000015
0.000010
7.102984e-06
5.041329e-06
3.621325e-06
2.601430e-06
1.925923e-06
1.446855e-06
1.094323e-06
17
0.001744
0.000952
0.000543
0.000311
0.000190
0.000117
0.000075
0.000049
0.000033
0.000022
0.000015
0.000010
0.000007
5.073477e-06
3.600890e-06
2.586618e-06
1.858133e-06
1.375636e-06
1.033450e-06
7.816457e-07
18
0.001046
0.000571
0.000326
0.000187
0.000114
0.000070
0.000045
0.000030
0.000020
0.000013
0.000009
0.000006
0.000004
3.044044e-06
2.160504e-06
1.551949e-06
1.114864e-06
8.253697e-07
6.200611e-07
4.689806e-07
19
0.000349
0.000190
0.000109
0.000062
0.000038
0.000023
0.000015
0.000010
0.000007
0.000004
0.000003
0.000002
0.000001
1.014675e-06
7.201632e-07
5.173129e-07
3.716190e-07
2.751214e-07
2.066856e-07
1.563258e-07
20
0.001037
0.000566
0.000323
0.000185
0.000113
0.000070
0.000045
0.000029
0.000020
0.000013
0.000009
0.000006
0.000004
3.016739e-06
2.141003e-06
1.537842e-06
1.104685e-06
8.177511e-07
6.142634e-07
4.645400e-07
21
0.001728
0.000943
0.000539
0.000309
0.000189
0.000116
0.000075
0.000049
0.000033
0.000022
0.000015
0.000010
0.000007
5.027852e-06
3.568305e-06
2.563046e-06
1.841125e-06
1.362905e-06
1.023762e-06
7.742258e-07
22
0.002420
0.001321
0.000754
0.000432
0.000264
0.000162
0.000105
0.000069
0.000046
0.000031
0.000021
0.000014
0.000010
7.038908e-06
4.995567e-06
3.588220e-06
2.577543e-06
1.908044e-06
1.433249e-06
1.083902e-06
23
0.000838
0.000459
0.000263
0.000151
0.000092
0.000057
0.000037
0.000024
0.000016
0.000011
0.000007
0.000005
0.000003
2.434600e-06
1.725531e-06
1.237515e-06
8.880885e-07
6.558311e-07
4.912221e-07
3.704348e-07
24
0.001397
0.000765
0.000438
0.000252
0.000154
0.000095
0.000061
0.000040
0.000027
0.000018
0.000012
0.000008
0.000006
4.057302e-06
2.875625e-06
2.062337e-06
1.480012e-06
1.092951e-06
8.186267e-07
6.173327e-07
25
0.001955
0.001071
0.000613
0.000352
0.000215
0.000132
0.000085
0.000056
0.000037
0.000025
0.000017
0.000012
0.000008
5.679315e-06
4.025227e-06
2.886805e-06
2.071680e-06
1.529879e-06
1.145886e-06
8.641199e-07
26
0.000775
0.000425
0.000243
0.000140
0.000086
0.000053
0.000034
0.000022
0.000015
0.000010
0.000007
0.000005
0.000003
2.250648e-06
1.594368e-06
1.142805e-06
8.198292e-07
6.048888e-07
4.525892e-07
3.409460e-07
27
0.001292
0.000709
0.000406
0.000233
0.000143
0.000088
0.000057
0.000037
0.000025
0.000016
0.000011
0.000008
0.000005
3.751039e-06
2.657251e-06
1.904654e-06
1.366367e-06
1.008137e-06
7.543069e-07
5.682369e-07
28
0.000258
0.000142
0.000081
0.000047
0.000029
0.000018
0.000011
0.000007
0.000005
0.000003
0.000002
0.000002
0.000001
7.501817e-07
5.314315e-07
3.809173e-07
2.732637e-07
2.016202e-07
1.508559e-07
1.136432e-07
29
0.001491
0.000821
0.000471
0.000272
0.000166
0.000102
0.000066
0.000043
0.000029
0.000019
0.000013
0.000009
0.000006
4.324013e-06
3.059414e-06
2.189866e-06
1.569593e-06
1.155555e-06
8.623594e-07
6.479636e-07
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
20
8
53
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
54
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
55
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
56
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
57
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
58
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
59
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
60
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
61
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
62
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
63
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
64
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
65
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
66
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
67
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
68
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
69
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
70
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
9
0
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
10
0
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
11
0
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
12
0
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
13
0
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
14
0
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
15
0
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
16
0
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
17
0
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
18
0
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
19
0
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
20
0
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
0.000000e+00
4435 rows × 20 columns
These are stored as Pandas DataFrames. An index can be supplied to obtain the population in a particular zone. E.g., for the ion populations of the innermost zone (index = 0), we will use
simulation.plasma.ion_number_density[0]
In [10]:
simulation.plasma.ion_number_density[0]
Out[10]:
atomic_number ion_number
8 0 1.380280e+03
1 5.393290e+08
2 9.780989e+04
3 3.252281e-09
4 0.000000e+00
5 0.000000e+00
6 0.000000e+00
7 0.000000e+00
8 0.000000e+00
12 0 1.346333e-06
1 6.084662e+02
2 5.606550e+07
3 4.120107e-17
4 0.000000e+00
5 0.000000e+00
6 0.000000e+00
7 0.000000e+00
8 0.000000e+00
9 0.000000e+00
10 0.000000e+00
11 0.000000e+00
12 0.000000e+00
14 0 7.121410e-04
1 1.057234e+05
2 8.401073e+08
3 8.192553e+05
4 7.678176e-04
5 0.000000e+00
6 0.000000e+00
7 0.000000e+00
...
18 10 0.000000e+00
11 0.000000e+00
12 0.000000e+00
13 0.000000e+00
14 0.000000e+00
15 0.000000e+00
16 0.000000e+00
17 0.000000e+00
18 0.000000e+00
20 0 1.332681e-08
1 2.335159e+01
2 3.400155e+07
3 7.014243e-04
4 0.000000e+00
5 0.000000e+00
6 0.000000e+00
7 0.000000e+00
8 0.000000e+00
9 0.000000e+00
10 0.000000e+00
11 0.000000e+00
12 0.000000e+00
13 0.000000e+00
14 0.000000e+00
15 0.000000e+00
16 0.000000e+00
17 0.000000e+00
18 0.000000e+00
19 0.000000e+00
20 0.000000e+00
Name: 0, Length: 94, dtype: float64
Ion populations for a particular ionization stage of a particular element can be accessed by specifying an appropriate tuple (𝑍,𝐶), which identifies the element (via atomic number 𝑍 ) and the charge (via the ion charge 𝐶 ). Thus, simulation.plasma.ion_number_density.loc[14,1] will identify the ion popuations for Si II (𝑍=14,𝐶=1) in all the zones.
In [11]:
simulation.plasma.ion_number_density.loc[14,1]
Out[11]:
0 105723.411638
1 44090.436208
2 21594.422001
3 9886.848793
4 5989.365036
5 3468.313242
6 2460.782422
7 1775.244868
8 1412.252955
9 971.631213
10 712.387551
11 560.635175
12 419.305282
13 335.000930
14 262.499557
15 211.157774
16 162.834905
17 144.128751
18 133.932017
19 124.963505
Name: (14, 1), dtype: float64
The above examples can be combined to obtain e.g. the Si II population in the innermost zone
can be obtained by
simulation.plasma.ion_number_density[0].loc[14,1]
In [12]:
simulation.plasma.ion_number_density[0].loc[14,1]
Out[12]:
105723.41163844905
The level populations are stored (and can be accessed) in a similar way - a third label can be used to pick out a particular atomic level. E.g., to pull out the population of the ground state (index 0) of Si II we can use simulation.plasma.level_number_density.loc[14,1,0]
In [13]:
simulation.plasma.level_number_density.loc[14,1,0]
Out[13]:
0 35851.432649
1 14931.643658
2 7307.391148
3 3341.597977
4 2024.238524
5 1171.823393
6 831.862975
7 600.439318
8 478.083068
9 328.985257
10 241.304587
11 190.021082
12 142.171548
13 113.651958
14 89.095253
15 71.705046
16 55.312585
17 48.994766
18 45.567026
19 42.549623
Name: (14, 1, 0), dtype: float64
Content source: kaushik94/tardis
Similar notebooks: