This notebook is available at https://github.com/tardis-sn/tardis/tree/master/docs/research/code_comparison/toy_models/reading blondin toymodel.ipynb
Compare Blondin Toy Models¶
[1]:
import re
import base64
import pandas as pd
import numpy as np
import yaml
from tardis import run_tardis
from astropy import units as u
from IPython.display import HTML
from tardis.io.parsers import read_blondin_toymodel, convert_blondin_toymodel
/Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/tqdm/autonotebook/__init__.py:14: TqdmExperimentalWarning: Using `tqdm.autonotebook.tqdm` in notebook mode. Use `tqdm.tqdm` instead to force console mode (e.g. in jupyter console)
" (e.g. in jupyter console)", TqdmExperimentalWarning)
[2]:
def create_df_download_link( df, title="Download CSV file", filename="data.csv", type='zip'):
if type != 'zip':
csv = df.to_csv()
b64 = base64.b64encode(csv.encode())
payload = b64.decode()
else:
payload = df.decode()
html = '<a download="{filename}" href="data:text/csv;base64,{payload}" target="_blank">{title}</a>'
html = html.format(payload=payload,title=title,filename=filename)
return HTML(html)
def create_download_link(download_data, title="Download CSV file", filename="data.csv", type='zip'):
if type != 'zip':
csv = df.to_csv()
b64 = base64.b64encode(csv.encode())
payload = b64.decode()
else:
payload = df.decode()
html = '<a download="{filename}" href="data:text/csv;base64,{payload}" target="_blank">{title}</a>'
html = html.format(payload=payload,title=title,filename=filename)
return HTML(html)
pattern_remove_bracket = re.compile('\[.+\]')
t0_pattern = re.compile('tend = (.+)\n')
Comparison of snia_toy01.dat¶
[3]:
convert_blondin_toymodel('snia_toy01.dat', 'snia_toy01.test.csvy', 9000, 30000, conversion_t_electron_rad=1/0.9)
[4]:
!code snia_toy01.test.csvy
[4]:
blondin_dict['v_inner_boundary'] = '9000 km/s'
blondin_dict['v_outer_boundary'] = '20000 km/s'
[5]:
csvy_file = '---\n{0}\n---\n{1}'.format(yaml.dump(blondin_dict, default_flow_style=False), blondin_csv.to_csv(index=False))
[6]:
with open('blondin_compare_01.csvy', 'w') as fh:
fh.write(csvy_file)
Download the model and setup
Assumptions and explanations
We assume \(T_\textrm{electron} = 0.9 T_\textrm{radiative}\)
[7]:
sim = run_tardis('blondin_model_compare_01.yml')
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/importlib/_bootstrap.py:219: QAWarning: pyne.data is not yet QA compliant.
return f(*args, **kwds)
(warnings.py:99)
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/importlib/_bootstrap.py:219: QAWarning: pyne.material is not yet QA compliant.
return f(*args, **kwds)
(warnings.py:99)
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.
AstropyDeprecationWarning)
(warnings.py:99)
[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/wkerzend/projects/tardis/tardis-data/kurucz_cd23_chianti_H_He.h5 (util.py:29)
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/IPython/core/interactiveshell.py:3267: PerformanceWarning: indexing past lexsort depth may impact performance.
exec(code_obj, self.user_global_ns, self.user_ns)
(warnings.py:99)
[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)
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.
AstropyDeprecationWarning)
(warnings.py:99)
[py.warnings ][WARNING] /Users/wkerzend/python/tardis/tardis/plasma/properties/atomic.py:195: FutureWarning: .labels was deprecated in version 0.24.0. Use .codes instead.
zeta_data['atomic_number'] = zeta_data.index.labels[0] + 1
(warnings.py:99)
[py.warnings ][WARNING] /Users/wkerzend/python/tardis/tardis/plasma/properties/atomic.py:196: FutureWarning: .labels was deprecated in version 0.24.0. Use .codes instead.
zeta_data['ion_number'] = zeta_data.index.labels[1] + 1
(warnings.py:99)
[tardis.plasma.properties.atomic][WARNING] Zeta_data missing - replaced with 1s. Missing ions: [(14, 15), (16, 17), (20, 21), (26, 27), (27, 28), (28, 29)] (atomic.py:215)
[py.warnings ][WARNING] /Users/wkerzend/python/tardis/tardis/plasma/properties/ion_population.py:63: FutureWarning:
Passing list-likes to .loc or [] with any missing label will raise
KeyError in the future, you can use .reindex() as an alternative.
See the documentation here:
https://pandas.pydata.org/pandas-docs/stable/indexing.html#deprecate-loc-reindex-listlike
partition_function.index].dropna())
(warnings.py:99)
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/equivalencies.py:90: RuntimeWarning: divide by zero encountered in double_scalars
(si.m, si.Hz, lambda x: _si.c.value / x),
(warnings.py:99)
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.
AstropyDeprecationWarning)
(warnings.py:99)
[tardis.simulation.base][INFO ] Starting iteration 1/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.30150e+43 erg / s Luminosity absorbed = 7.43372e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 5000.0 7831.059753 0.462810 0.987065
5 5000.0 7106.913190 0.379149 1.304910
10 5000.0 6779.027752 0.336294 1.473708
15 5000.0 6591.160845 0.304883 1.556553
20 5000.0 6469.413444 0.279753 1.549488
25 5000.0 6369.832186 0.258765 1.548078
30 5000.0 6285.178795 0.240774 1.531938
35 5000.0 6216.656067 0.225078 1.509501
40 5000.0 6155.241379 0.211205 1.474743
45 5000.0 6105.918970 0.198821 1.414270
50 5000.0 6076.371342 0.187679 1.348556
55 5000.0 6030.241880 0.177588 1.311088
60 5000.0 6005.973722 0.168399 1.256169
65 5000.0 6002.125995 0.159992 1.169288
70 5000.0 5976.207283 0.152269 1.122986
75 5000.0 5951.666260 0.145148 1.077897
80 5000.0 5918.254016 0.138563 1.053437
85 5000.0 5918.208740 0.132453 0.989942
90 5000.0 5867.100159 0.126772 0.990185
95 5000.0 5859.007030 0.121475 0.929519
100 5000.0 5843.254166 0.116527 0.894981
105 5000.0 5841.131195 0.111894 0.849290
110 5000.0 5821.292215 0.107550 0.811782
115 5000.0 5820.272035 0.103468 0.769661
120 5000.0 5811.350827 0.099626 0.736181
125 5000.0 5793.908583 0.096006 0.716442
130 5000.0 5795.478604 0.092589 0.681608
135 5000.0 5773.849040 0.089359 0.667536
140 5000.0 5775.403713 0.086303 0.634047
145 5000.0 5736.539129 0.083408 0.630179
150 5000.0 5732.357290 0.080662 0.601033
155 5000.0 5743.364284 0.078054 0.570580
160 5000.0 5733.145769 0.075575 0.551961
165 5000.0 5718.331258 0.073217 0.529948
170 5000.0 5699.632837 0.070971 0.519324
175 5000.0 5695.349689 0.068830 0.499699
180 5000.0 5688.534677 0.066787 0.478921
185 5000.0 5683.055043 0.064837 0.463320
190 5000.0 5665.611101 0.062973 0.450357
195 5000.0 5654.224456 0.061191 0.437244
200 5000.0 5651.225192 0.059485 0.421169
205 5000.0 5626.660446 0.057852 0.414763
210 5000.0 5595.199895 0.056286 0.409042
215 5000.0 5577.506895 0.054785 0.400196
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[py.warnings ][WARNING] /Users/wkerzend/python/tardis/tardis/plasma/properties/ion_population.py:63: FutureWarning:
Passing list-likes to .loc or [] with any missing label will raise
KeyError in the future, you can use .reindex() as an alternative.
See the documentation here:
https://pandas.pydata.org/pandas-docs/stable/indexing.html#deprecate-loc-reindex-listlike
partition_function.index].dropna())
(warnings.py:99)
[tardis.simulation.base][INFO ] Starting iteration 2/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.38526e+43 erg / s Luminosity absorbed = 6.66321e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 7831.059753 9943.504253 0.987065 0.941908
5 7106.913190 9168.067423 1.304910 1.150356
10 6779.027752 8678.999816 1.473708 1.286543
15 6591.160845 8377.317578 1.556553 1.346125
20 6469.413444 8151.787425 1.549488 1.358795
25 6369.832186 7985.142325 1.548078 1.353925
30 6285.178795 7878.847609 1.531938 1.319217
35 6216.656067 7741.589421 1.509501 1.307884
40 6155.241379 7655.316512 1.474743 1.270084
45 6105.918970 7556.986602 1.414270 1.234470
50 6076.371342 7522.785118 1.348556 1.170937
55 6030.241880 7447.391749 1.311088 1.133516
60 6005.973722 7408.925630 1.256169 1.089917
65 6002.125995 7364.361995 1.169288 1.029933
70 5976.207283 7336.100563 1.122986 0.990663
75 5951.666260 7306.592489 1.077897 0.946869
80 5918.254016 7238.893988 1.053437 0.931127
85 5918.208740 7226.363390 0.989942 0.877828
90 5867.100159 7197.898366 0.990185 0.857879
95 5859.007030 7165.997446 0.929519 0.812933
100 5843.254166 7159.070470 0.894981 0.775213
105 5841.131195 7143.741089 0.849290 0.738824
110 5821.292215 7119.704519 0.811782 0.707784
115 5820.272035 7112.481717 0.769661 0.671211
120 5811.350827 7060.505651 0.736181 0.652498
125 5793.908583 7066.008522 0.716442 0.624559
130 5795.478604 7029.029121 0.681608 0.605209
135 5773.849040 7021.659449 0.667536 0.582533
140 5775.403713 7010.212505 0.634047 0.555628
145 5736.539129 6990.562369 0.630179 0.542516
150 5732.357290 6984.884673 0.601033 0.518933
155 5743.364284 6980.471223 0.570580 0.494955
160 5733.145769 6938.386837 0.551961 0.484650
165 5718.331258 6928.177798 0.529948 0.465554
170 5699.632837 6904.323015 0.519324 0.451990
175 5695.349689 6905.684428 0.499699 0.434141
180 5688.534677 6893.408541 0.478921 0.418797
185 5683.055043 6881.358878 0.463320 0.404739
190 5665.611101 6848.434616 0.450357 0.394703
195 5654.224456 6831.751052 0.437244 0.382167
200 5651.225192 6821.557616 0.421169 0.367724
205 5626.660446 6798.252052 0.414763 0.360225
210 5595.199895 6760.804314 0.409042 0.354866
215 5577.506895 6740.013534 0.400196 0.344957
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 3/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.32445e+43 erg / s Luminosity absorbed = 7.18316e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 9943.504253 11184.865494 0.941908 0.874243
5 9168.067423 10590.890516 1.150356 0.967498
10 8678.999816 10184.538479 1.286543 1.029334
15 8377.317578 9900.090529 1.346125 1.054300
20 8151.787425 9664.894816 1.358795 1.054968
25 7985.142325 9510.226327 1.353925 1.037032
30 7878.847609 9394.106182 1.319217 1.004565
35 7741.589421 9271.496145 1.307884 0.982576
40 7655.316512 9161.047735 1.270084 0.949923
45 7556.986602 9072.338911 1.234470 0.915208
50 7522.785118 9010.672618 1.170937 0.872802
55 7447.391749 8901.383373 1.133516 0.845834
60 7408.925630 8834.036103 1.089917 0.815912
65 7364.361995 8769.307119 1.029933 0.772520
70 7336.100563 8711.903127 0.990663 0.746307
75 7306.592489 8640.885761 0.946869 0.718271
80 7238.893988 8542.892115 0.931127 0.704924
85 7226.363390 8516.723270 0.877828 0.666870
90 7197.898366 8502.592440 0.857879 0.642794
95 7165.997446 8447.746462 0.812933 0.611320
100 7159.070470 8412.121792 0.775213 0.587311
105 7143.741089 8381.138885 0.738824 0.562031
110 7119.704519 8359.007818 0.707784 0.537153
115 7112.481717 8322.937242 0.671211 0.513156
120 7060.505651 8249.430763 0.652498 0.500380
125 7066.008522 8275.891240 0.624559 0.473512
130 7029.029121 8204.492926 0.605209 0.461398
135 7021.659449 8180.726449 0.582533 0.445042
140 7010.212505 8169.127793 0.555628 0.424694
145 6990.562369 8158.252968 0.542516 0.411711
150 6984.884673 8112.858348 0.518933 0.397343
155 6980.471223 8098.797522 0.494955 0.379947
160 6938.386837 8055.533873 0.484650 0.370261
165 6928.177798 8055.402146 0.465554 0.354906
170 6904.323015 8019.347064 0.451990 0.345449
175 6905.684428 8006.258863 0.434141 0.331780
180 6893.408541 8002.951619 0.418797 0.319606
185 6881.358878 7970.659310 0.404739 0.309280
190 6848.434616 7932.276565 0.394703 0.301617
195 6831.751052 7915.200790 0.382167 0.291916
200 6821.557616 7898.840956 0.367724 0.280425
205 6798.252052 7883.486836 0.360225 0.273069
210 6760.804314 7840.128677 0.354866 0.268436
215 6740.013534 7815.699388 0.344957 0.260891
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 4/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.30498e+43 erg / s Luminosity absorbed = 7.38203e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 11184.865494 11911.711939 0.874243 0.815342
5 10590.890516 11423.170586 0.967498 0.852032
10 10184.538479 11134.470932 1.029334 0.863305
15 9900.090529 10887.349224 1.054300 0.863750
20 9664.894816 10740.461976 1.054968 0.832290
25 9510.226327 10562.641071 1.037032 0.813050
30 9394.106182 10481.388048 1.004565 0.774264
35 9271.496145 10372.742105 0.982576 0.747951
40 9161.047735 10274.337236 0.949923 0.716987
45 9072.338911 10202.946659 0.915208 0.684430
50 9010.672618 10120.383319 0.872802 0.652331
55 8901.383373 10039.433120 0.845834 0.624391
60 8834.036103 9985.105525 0.815912 0.597331
65 8769.307119 9909.687407 0.772520 0.567351
70 8711.903127 9853.157362 0.746307 0.545114
75 8640.885761 9803.411756 0.718271 0.521662
80 8542.892115 9703.094119 0.704924 0.508055
85 8516.723270 9664.369661 0.666870 0.483056
90 8502.592440 9644.808583 0.642794 0.462651
95 8447.746462 9595.659098 0.611320 0.439338
100 8412.121792 9572.851164 0.587311 0.421222
105 8381.138885 9519.270385 0.562031 0.403866
110 8359.007818 9520.334232 0.537153 0.383152
115 8322.937242 9482.738472 0.513156 0.366423
120 8249.430763 9426.321672 0.500380 0.356164
125 8275.891240 9391.728992 0.473512 0.340621
130 8204.492926 9367.333274 0.461398 0.327122
135 8180.726449 9303.898069 0.445042 0.318214
140 8169.127793 9283.891027 0.424694 0.304527
145 8158.252968 9270.457341 0.411711 0.293142
150 8112.858348 9230.938730 0.397343 0.282674
155 8098.797522 9190.155128 0.379947 0.272658
160 8055.533873 9192.939613 0.370261 0.262679
165 8055.402146 9166.640157 0.354906 0.253095
170 8019.347064 9116.658009 0.345449 0.246487
175 8006.258863 9096.261076 0.331780 0.237289
180 8002.951619 9095.824747 0.319606 0.228009
185 7970.659310 9083.284785 0.309280 0.219454
190 7932.276565 9051.094660 0.301617 0.212873
195 7915.200790 9020.044932 0.291916 0.207175
200 7898.840956 8998.843933 0.280425 0.199036
205 7883.486836 8982.946944 0.273069 0.193372
210 7840.128677 8952.622269 0.268436 0.188983
215 7815.699388 8935.684680 0.260891 0.182965
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 5/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.28254e+43 erg / s Luminosity absorbed = 7.59247e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 11911.711939 12274.410751 0.815342 0.789790
5 11423.170586 11871.455425 0.852032 0.792676
10 11134.470932 11649.695187 0.863305 0.773668
15 10887.349224 11471.210556 0.863750 0.748405
20 10740.461976 11360.332943 0.832290 0.712461
25 10562.641071 11218.786535 0.813050 0.682594
30 10481.388048 11129.781624 0.774264 0.647152
35 10372.742105 11054.982577 0.747951 0.615252
40 10274.337236 10978.674047 0.716987 0.581735
45 10202.946659 10896.716952 0.684430 0.554543
50 10120.383319 10820.747796 0.652331 0.526362
55 10039.433120 10788.603828 0.624391 0.494951
60 9985.105525 10721.956293 0.597331 0.471665
65 9909.687407 10655.755485 0.567351 0.446374
70 9853.157362 10634.209059 0.545114 0.423620
75 9803.411756 10590.884258 0.521662 0.402309
80 9703.094119 10508.509946 0.508055 0.388143
85 9664.369661 10508.882030 0.483056 0.365331
90 9644.808583 10483.115608 0.462651 0.347868
95 9595.659098 10430.871132 0.439338 0.329624
100 9572.851164 10415.185286 0.421222 0.314476
105 9519.270385 10368.229649 0.403866 0.300532
110 9520.334232 10371.842110 0.383152 0.284290
115 9482.738472 10344.461437 0.366423 0.270853
120 9426.321672 10330.632364 0.356164 0.260494
125 9391.728992 10280.580995 0.340621 0.250351
130 9367.333274 10249.596303 0.327122 0.239675
135 9303.898069 10236.391898 0.318214 0.229879
140 9283.891027 10200.349641 0.304527 0.220010
145 9270.457341 10180.401007 0.293142 0.212525
150 9230.938730 10140.384246 0.282674 0.204939
155 9190.155128 10095.159566 0.272658 0.197206
160 9192.939613 10116.202875 0.262679 0.188372
165 9166.640157 10064.971970 0.253095 0.182719
170 9116.658009 10059.156300 0.246487 0.176322
175 9096.261076 10026.674915 0.237289 0.170347
180 9095.824747 9986.590882 0.228009 0.165054
185 9083.284785 9991.204987 0.219454 0.157743
190 9051.094660 9954.915111 0.212873 0.153163
195 9020.044932 9946.636259 0.207175 0.148020
200 8998.843933 9904.779162 0.199036 0.142722
205 8982.946944 9878.309679 0.193372 0.138820
210 8952.622269 9865.404965 0.188983 0.134775
215 8935.684680 9838.418760 0.182965 0.130692
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 6/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.28283e+43 erg / s Luminosity absorbed = 7.59368e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 12274.410751 12621.985366 0.789790 0.740187
5 11871.455425 12288.282401 0.792676 0.714197
10 11649.695187 12107.496695 0.773668 0.683174
15 11471.210556 11966.421542 0.748405 0.651469
20 11360.332943 11869.120463 0.712461 0.610504
25 11218.786535 11782.584290 0.682594 0.576377
30 11129.781624 11661.517833 0.647152 0.549776
35 11054.982577 11596.680082 0.615252 0.515964
40 10978.674047 11486.100778 0.581735 0.492717
45 10896.716952 11411.634100 0.554543 0.466075
50 10820.747796 11320.929647 0.526362 0.443230
55 10788.603828 11304.819988 0.494951 0.411733
60 10721.956293 11223.787471 0.471665 0.391653
65 10655.755485 11187.391580 0.446374 0.366972
70 10634.209059 11145.379508 0.423620 0.349026
75 10590.884258 11132.155630 0.402309 0.327809
80 10508.509946 11080.451217 0.388143 0.312311
85 10508.882030 11053.402190 0.365331 0.294910
90 10483.115608 11022.793935 0.347868 0.280215
95 10430.871132 11018.638959 0.329624 0.263074
100 10415.185286 10969.698101 0.314476 0.251057
105 10368.229649 10996.575238 0.300532 0.235733
110 10371.842110 10946.225411 0.284290 0.225373
115 10344.461437 10915.429238 0.270853 0.215369
120 10330.632364 10911.516887 0.260494 0.204801
125 10280.580995 10876.787083 0.250351 0.195768
130 10249.596303 10846.880881 0.239675 0.187419
135 10236.391898 10827.878807 0.229879 0.179927
140 10200.349641 10843.147042 0.220010 0.169929
145 10180.401007 10836.452058 0.212525 0.163079
150 10140.384246 10779.508126 0.204939 0.157783
155 10095.159566 10792.246053 0.197206 0.150208
160 10116.202875 10777.101313 0.188372 0.144407
165 10064.971970 10764.279483 0.182719 0.138515
170 10059.156300 10755.323524 0.176322 0.133355
175 10026.674915 10716.093969 0.170347 0.129014
180 9986.590882 10681.240209 0.165054 0.124735
185 9991.204987 10710.349167 0.157743 0.118591
190 9954.915111 10647.113250 0.153163 0.115792
195 9946.636259 10621.085728 0.148020 0.112042
200 9904.779162 10611.632892 0.142722 0.107582
205 9878.309679 10575.876336 0.138820 0.104696
210 9865.404965 10534.592000 0.134775 0.102013
215 9838.418760 10485.222357 0.130692 0.099331
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 7/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.28515e+43 erg / s Luminosity absorbed = 7.57676e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 12621.985366 12759.629863 0.740187 0.721720
5 12288.282401 12486.425455 0.714197 0.679950
10 12107.496695 12369.251952 0.683174 0.636363
15 11966.421542 12252.854693 0.651469 0.595850
20 11869.120463 12167.826060 0.610504 0.556063
25 11782.584290 12071.481834 0.576377 0.522488
30 11661.517833 11955.870966 0.549776 0.496104
35 11596.680082 11917.514987 0.515964 0.459262
40 11486.100778 11870.688289 0.492717 0.427744
45 11411.634100 11790.779683 0.466075 0.405292
50 11320.929647 11703.842476 0.443230 0.384510
55 11304.819988 11687.592974 0.411733 0.356790
60 11223.787471 11618.853302 0.391653 0.337614
65 11187.391580 11585.501944 0.366972 0.314979
70 11145.379508 11526.322819 0.349026 0.299224
75 11132.155630 11523.858835 0.327809 0.281351
80 11080.451217 11448.221615 0.312311 0.269039
85 11053.402190 11462.009705 0.294910 0.249884
90 11022.793935 11416.159361 0.280215 0.238194
95 11018.638959 11406.804183 0.263074 0.222978
100 10969.698101 11376.225310 0.251057 0.211479
105 10996.575238 11433.350181 0.235733 0.196343
110 10946.225411 11372.910863 0.225373 0.189059
115 10915.429238 11332.677175 0.215369 0.180570
120 10911.516887 11285.604610 0.204801 0.172781
125 10876.787083 11299.429914 0.195768 0.162189
130 10846.880881 11263.153673 0.187419 0.155143
135 10827.878807 11204.982403 0.179927 0.150862
140 10843.147042 11258.798865 0.169929 0.141258
145 10836.452058 11231.762696 0.163079 0.135271
150 10779.508126 11207.707717 0.157783 0.129487
155 10792.246053 11191.745451 0.150208 0.124238
160 10777.101313 11181.509283 0.144407 0.119067
165 10764.279483 11162.391880 0.138515 0.114391
170 10755.323524 11158.541462 0.133355 0.109545
175 10716.093969 11151.075575 0.129014 0.105748
180 10681.240209 11118.631655 0.124735 0.102085
185 10710.349167 11103.407863 0.118591 0.097843
190 10647.113250 11072.785328 0.115792 0.094608
195 10621.085728 11057.565986 0.112042 0.091442
200 10611.632892 11035.369976 0.107582 0.088212
205 10575.876336 11002.115565 0.104696 0.085741
210 10534.592000 10955.610971 0.102013 0.083497
215 10485.222357 10917.579162 0.099331 0.081022
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 8/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.27712e+43 erg / s Luminosity absorbed = 7.65528e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 12759.629863 12860.359189 0.721720 0.708992
5 12486.425455 12638.631249 0.679950 0.651342
10 12369.251952 12561.844673 0.636363 0.598991
15 12252.854693 12457.507444 0.595850 0.554747
20 12167.826060 12405.099383 0.556063 0.514695
25 12071.481834 12309.733919 0.522488 0.482479
30 11955.870966 12245.489227 0.496104 0.450183
35 11917.514987 12183.541188 0.459262 0.418285
40 11870.688289 12158.321534 0.427744 0.386364
45 11790.779683 12097.953291 0.405292 0.362136
50 11703.842476 12043.102357 0.384510 0.340574
55 11687.592974 11992.235264 0.356790 0.317048
60 11618.853302 11931.000829 0.337614 0.297179
65 11585.501944 11868.817978 0.314979 0.281821
70 11526.322819 11845.559084 0.299224 0.264268
75 11523.858835 11839.159985 0.281351 0.247294
80 11448.221615 11766.869353 0.269039 0.236873
85 11462.009705 11781.216970 0.249884 0.219222
90 11416.159361 11763.657217 0.238194 0.206710
95 11406.804183 11736.050598 0.222978 0.194684
100 11376.225310 11709.075282 0.211479 0.184521
105 11433.350181 11764.053125 0.196343 0.171069
110 11372.910863 11693.485667 0.189059 0.164673
115 11332.677175 11662.729622 0.180570 0.157151
120 11285.604610 11626.621706 0.172781 0.149390
125 11299.429914 11591.005256 0.162189 0.141987
130 11263.153673 11621.706183 0.155143 0.133331
135 11204.982403 11562.094749 0.150862 0.129172
140 11258.798865 11593.366562 0.141258 0.121613
145 11231.762696 11543.895750 0.135271 0.117417
150 11207.707717 11510.519401 0.129487 0.112442
155 11191.745451 11521.475871 0.124238 0.107162
160 11181.509283 11508.760099 0.119067 0.102595
165 11162.391880 11476.537118 0.114391 0.098371
170 11158.541462 11464.294275 0.109545 0.094355
175 11151.075575 11445.541800 0.105748 0.091276
180 11118.631655 11422.237088 0.102085 0.087966
185 11103.407863 11397.864821 0.097843 0.084664
190 11072.785328 11382.914370 0.094608 0.081583
195 11057.565986 11349.784124 0.091442 0.078976
200 11035.369976 11322.806026 0.088212 0.076562
205 11002.115565 11263.555912 0.085741 0.074833
210 10955.610971 11221.660698 0.083497 0.072828
215 10917.579162 11200.040416 0.081022 0.070359
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 9/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.26945e+43 erg / s Luminosity absorbed = 7.72465e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 12860.359189 12957.621989 0.708992 0.691150
5 12638.631249 12803.884316 0.651342 0.624900
10 12561.844673 12730.751370 0.598991 0.569845
15 12457.507444 12604.055249 0.554747 0.528149
20 12405.099383 12583.906104 0.514695 0.485502
25 12309.733919 12491.390514 0.482479 0.453408
30 12245.489227 12424.403292 0.450183 0.423553
35 12183.541188 12380.236496 0.418285 0.388916
40 12158.321534 12363.358789 0.386364 0.359253
45 12097.953291 12328.785784 0.362136 0.335280
50 12043.102357 12251.735081 0.340574 0.314188
55 11992.235264 12261.737757 0.317048 0.288563
60 11931.000829 12152.175154 0.297179 0.274420
65 11868.817978 12117.845990 0.281821 0.257238
70 11845.559084 12093.081412 0.264268 0.241895
75 11839.159985 12093.052004 0.247294 0.224445
80 11766.869353 12037.640748 0.236873 0.214143
85 11781.216970 12049.271412 0.219222 0.198544
90 11763.657217 12015.545530 0.206710 0.188042
95 11736.050598 11993.213751 0.194684 0.176258
100 11709.075282 11930.668715 0.184521 0.168432
105 11764.053125 12037.937411 0.171069 0.153223
110 11693.485667 11988.734420 0.164673 0.146614
115 11662.729622 11937.603116 0.157151 0.140315
120 11626.621706 11887.030622 0.149390 0.133416
125 11591.005256 11836.102537 0.141987 0.127737
130 11621.706183 11833.294958 0.133331 0.120973
135 11562.094749 11777.150834 0.129172 0.116478
140 11593.366562 11825.884043 0.121613 0.109144
145 11543.895750 11766.196357 0.117417 0.105560
150 11510.519401 11740.623271 0.112442 0.101324
155 11521.475871 11732.541118 0.107162 0.096637
160 11508.760099 11708.165413 0.102595 0.092791
165 11476.537118 11672.422129 0.098371 0.089207
170 11464.294275 11660.830520 0.094355 0.085588
175 11445.541800 11633.201479 0.091276 0.082922
180 11422.237088 11615.548095 0.087966 0.079616
185 11397.864821 11566.309428 0.084664 0.077288
190 11382.914370 11575.749336 0.081583 0.073956
195 11349.784124 11543.883213 0.078976 0.071674
200 11322.806026 11495.598254 0.076562 0.069762
205 11263.555912 11490.479524 0.074833 0.067104
210 11221.660698 11428.432630 0.072828 0.065627
215 11200.040416 11408.594202 0.070359 0.063437
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 10/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.26942e+43 erg / s Luminosity absorbed = 7.72233e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 12957.621989 13015.569839 0.691150 0.682918
5 12803.884316 12904.875006 0.624900 0.607993
10 12730.751370 12842.565429 0.569845 0.550057
15 12604.055249 12753.343916 0.528149 0.505693
20 12583.906104 12730.133218 0.485502 0.462479
25 12491.390514 12664.083226 0.453408 0.427473
30 12424.403292 12588.842740 0.423553 0.399502
35 12380.236496 12576.400081 0.388916 0.365150
40 12363.358789 12571.096532 0.359253 0.335182
45 12328.785784 12567.267055 0.335280 0.308863
50 12251.735081 12490.762479 0.314188 0.289858
55 12261.737757 12473.100822 0.288563 0.268680
60 12152.175154 12403.844620 0.274420 0.252799
65 12117.845990 12379.017248 0.257238 0.235858
70 12093.081412 12384.414417 0.241895 0.219277
75 12093.052004 12359.661933 0.224445 0.205534
80 12037.640748 12301.567714 0.214143 0.194798
85 12049.271412 12282.445145 0.198544 0.182032
90 12015.545530 12245.263730 0.188042 0.172696
95 11993.213751 12240.211097 0.176258 0.161126
100 11930.668715 12166.777597 0.168432 0.154560
105 12037.937411 12230.213020 0.153223 0.141698
110 11988.734420 12176.872696 0.146614 0.135522
115 11937.603116 12173.728933 0.140315 0.128104
120 11887.030622 12135.497772 0.133416 0.121988
125 11836.102537 12064.394959 0.127737 0.117167
130 11833.294958 12010.587301 0.120973 0.112690
135 11777.150834 11986.931983 0.116478 0.107698
140 11825.884043 12025.244790 0.109144 0.100791
145 11766.196357 11995.910253 0.105560 0.096670
150 11740.623271 11947.593265 0.101324 0.093294
155 11732.541118 11957.542102 0.096637 0.088421
160 11708.165413 11910.670130 0.092791 0.085384
165 11672.422129 11888.528792 0.089207 0.081976
170 11660.830520 11879.429183 0.085588 0.078285
175 11633.201479 11864.677062 0.082922 0.075451
180 11615.548095 11802.638064 0.079616 0.073662
185 11566.309428 11793.285395 0.077288 0.070677
190 11575.749336 11763.174836 0.073956 0.068157
195 11543.883213 11739.187171 0.071674 0.065917
200 11495.598254 11661.501839 0.069762 0.064739
205 11490.479524 11650.971928 0.067104 0.062392
210 11428.432630 11591.428842 0.065627 0.060987
215 11408.594202 11546.502844 0.063437 0.059336
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 11/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.26262e+43 erg / s Luminosity absorbed = 7.78847e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13015.569839 13032.401779 0.682918 0.684003
5 12904.875006 12990.461599 0.607993 0.593997
10 12842.565429 12898.821890 0.550057 0.540517
15 12753.343916 12849.193235 0.505693 0.490550
20 12730.133218 12828.953617 0.462479 0.448185
25 12664.083226 12759.303136 0.427473 0.414128
30 12588.842740 12696.103173 0.399502 0.385663
35 12576.400081 12697.730273 0.365150 0.352612
40 12571.096532 12662.177207 0.335182 0.324425
45 12567.267055 12656.118492 0.308863 0.298253
50 12490.762479 12628.623376 0.289858 0.274518
55 12473.100822 12623.212204 0.268680 0.252985
60 12403.844620 12591.031591 0.252799 0.236744
65 12379.017248 12568.788110 0.235858 0.220900
70 12384.414417 12504.218648 0.219277 0.209178
75 12359.661933 12472.497115 0.205534 0.196505
80 12301.567714 12476.433585 0.194798 0.182159
85 12282.445145 12444.116662 0.182032 0.170683
90 12245.263730 12408.486842 0.172696 0.162078
95 12240.211097 12333.668316 0.161126 0.154184
100 12166.777597 12303.324181 0.154560 0.145725
105 12230.213020 12324.197495 0.141698 0.135636
110 12176.872696 12297.477032 0.135522 0.128069
115 12173.728933 12298.340501 0.128104 0.120772
120 12135.497772 12244.175855 0.121988 0.116145
125 12064.394959 12220.660279 0.117167 0.110258
130 12010.587301 12203.437058 0.112690 0.104950
135 11986.931983 12204.633362 0.107698 0.099736
140 12025.244790 12194.368362 0.100791 0.094746
145 11995.910253 12188.022851 0.096670 0.090104
150 11947.593265 12106.440886 0.093294 0.087625
155 11957.542102 12108.357437 0.088421 0.083173
160 11910.670130 12056.905371 0.085384 0.080364
165 11888.528792 12075.433483 0.081976 0.076282
170 11879.429183 12050.882022 0.078285 0.073261
175 11864.677062 12025.339383 0.075451 0.070721
180 11802.638064 11946.287508 0.073662 0.069445
185 11793.285395 11931.799902 0.070677 0.066790
190 11763.174836 11906.870608 0.068157 0.064230
195 11739.187171 11849.898176 0.065917 0.062563
200 11661.501839 11805.166140 0.064739 0.060991
205 11650.971928 11750.970412 0.062392 0.059559
210 11591.428842 11729.652424 0.060987 0.057636
215 11546.502844 11651.688056 0.059336 0.056416
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 12/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.25945e+43 erg / s Luminosity absorbed = 7.81529e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13032.401779 13060.756844 0.684003 0.676972
5 12990.461599 13053.655799 0.593997 0.583799
10 12898.821890 12940.461814 0.540517 0.534608
15 12849.193235 12911.505202 0.490550 0.483174
20 12828.953617 12900.206203 0.448185 0.439200
25 12759.303136 12863.788022 0.414128 0.400384
30 12696.103173 12810.013851 0.385663 0.371252
35 12697.730273 12825.414714 0.352612 0.338216
40 12662.177207 12801.179422 0.324425 0.310208
45 12656.118492 12744.748623 0.298253 0.290688
50 12628.623376 12756.940903 0.274518 0.265593
55 12623.212204 12719.938507 0.252985 0.246729
60 12591.031591 12728.329492 0.236744 0.226363
65 12568.788110 12688.708886 0.220900 0.211208
70 12504.218648 12660.204209 0.209178 0.197615
75 12472.497115 12599.352691 0.196505 0.186919
80 12476.433585 12585.793885 0.182159 0.174973
85 12444.116662 12589.414222 0.170683 0.162163
90 12408.486842 12515.279884 0.162078 0.155383
95 12333.668316 12486.897379 0.154184 0.145659
100 12303.324181 12430.500086 0.145725 0.138556
105 12324.197495 12450.715847 0.135636 0.129425
110 12297.477032 12426.647529 0.128069 0.122309
115 12298.340501 12435.989959 0.120772 0.114494
120 12244.175855 12407.884661 0.116145 0.109221
125 12220.660279 12380.981505 0.110258 0.103989
130 12203.437058 12351.668125 0.104950 0.099268
135 12204.633362 12304.519778 0.099736 0.095737
140 12194.368362 12324.560098 0.094746 0.090400
145 12188.022851 12265.886677 0.090104 0.087024
150 12106.440886 12249.589898 0.087625 0.083028
155 12108.357437 12237.977955 0.083173 0.079321
160 12056.905371 12207.745803 0.080364 0.076236
165 12075.433483 12201.219362 0.076282 0.072831
170 12050.882022 12167.349382 0.073261 0.070259
175 12025.339383 12095.417592 0.070721 0.068696
180 11946.287508 12054.818787 0.069445 0.066319
185 11931.799902 12019.446420 0.066790 0.064144
190 11906.870608 11990.990568 0.064230 0.061757
195 11849.898176 11961.647695 0.062563 0.059858
200 11805.166140 11924.727184 0.060991 0.058232
205 11750.970412 11858.754302 0.059559 0.056888
210 11729.652424 11829.679705 0.057636 0.055114
215 11651.688056 11767.828592 0.056416 0.053853
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 13/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.24878e+43 erg / s Luminosity absorbed = 7.90910e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13060.756844 13083.224882 0.676972 0.673697
5 13053.655799 13088.802662 0.583799 0.577711
10 12940.461814 13027.101067 0.534608 0.522890
15 12911.505202 13006.462644 0.483174 0.473946
20 12900.206203 12983.888230 0.439200 0.429670
25 12863.788022 12971.745863 0.400384 0.388831
30 12810.013851 12913.558012 0.371252 0.361218
35 12825.414714 12948.249963 0.338216 0.325705
40 12801.179422 12912.855169 0.310208 0.301236
45 12744.748623 12866.250179 0.290688 0.279866
50 12756.940903 12848.652789 0.265593 0.258571
55 12719.938507 12848.859656 0.246729 0.236994
60 12728.329492 12851.601535 0.226363 0.217969
65 12688.708886 12808.658408 0.211208 0.204525
70 12660.204209 12799.122574 0.197615 0.189788
75 12599.352691 12726.546201 0.186919 0.179897
80 12585.793885 12747.782242 0.174973 0.165916
85 12589.414222 12750.918368 0.162163 0.154663
90 12515.279884 12668.223097 0.155383 0.148464
95 12486.897379 12675.791676 0.145659 0.137905
100 12430.500086 12620.707010 0.138556 0.130616
105 12450.715847 12593.289701 0.129425 0.123022
110 12426.647529 12585.676771 0.122309 0.116424
115 12435.989959 12537.431630 0.114494 0.110995
120 12407.884661 12520.786382 0.109221 0.105135
125 12380.981505 12516.716438 0.103989 0.099087
130 12351.668125 12472.741826 0.099268 0.095006
135 12304.519778 12422.139907 0.095737 0.091867
140 12324.560098 12453.435118 0.090400 0.086369
145 12265.886677 12389.456714 0.087024 0.083330
150 12249.589898 12363.390556 0.083028 0.079554
155 12237.977955 12341.590207 0.079321 0.076214
160 12207.745803 12291.756528 0.076236 0.073703
165 12201.219362 12252.749433 0.072831 0.070982
170 12167.349382 12256.286840 0.070259 0.067874
175 12095.417592 12182.336237 0.068696 0.066291
180 12054.818787 12145.611052 0.066319 0.063949
185 12019.446420 12120.021352 0.064144 0.061681
190 11990.990568 12095.073710 0.061757 0.059340
195 11961.647695 12049.179195 0.059858 0.057897
200 11924.727184 11992.337972 0.058232 0.056356
205 11858.754302 11949.219275 0.056888 0.054787
210 11829.679705 11885.344010 0.055114 0.053671
215 11767.828592 11823.561736 0.053853 0.052411
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 14/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.24323e+43 erg / s Luminosity absorbed = 7.95660e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13083.224882 13110.758198 0.673697 0.670180
5 13088.802662 13113.028903 0.577711 0.574875
10 13027.101067 13093.088462 0.522890 0.515120
15 13006.462644 13078.630835 0.473946 0.464486
20 12983.888230 13073.991303 0.429670 0.419770
25 12971.745863 13110.381168 0.388831 0.375311
30 12913.558012 13055.574882 0.361218 0.344822
35 12948.249963 13041.575051 0.325705 0.318035
40 12912.855169 13022.005042 0.301236 0.292756
45 12866.250179 13003.812183 0.279866 0.269639
50 12848.652789 12963.032879 0.258571 0.250234
55 12848.859656 12982.389906 0.236994 0.229901
60 12851.601535 12975.350494 0.217969 0.210624
65 12808.658408 12965.368804 0.204525 0.196254
70 12799.122574 12890.584285 0.189788 0.185073
75 12726.546201 12893.430696 0.179897 0.171994
80 12747.782242 12901.177598 0.165916 0.159412
85 12750.918368 12894.708333 0.154663 0.148871
90 12668.223097 12834.988058 0.148464 0.140906
95 12675.791676 12833.456390 0.137905 0.131452
100 12620.707010 12807.914014 0.130616 0.123631
105 12593.289701 12725.023233 0.123022 0.118612
110 12585.676771 12715.728557 0.116424 0.112363
115 12537.431630 12681.404241 0.110995 0.106558
120 12520.786382 12659.267376 0.105135 0.100885
125 12516.716438 12647.301366 0.099087 0.095555
130 12472.741826 12601.531043 0.095006 0.091565
135 12422.139907 12567.359999 0.091867 0.087506
140 12453.435118 12518.962460 0.086369 0.084548
145 12389.456714 12478.839839 0.083330 0.080795
150 12363.390556 12494.070660 0.079554 0.076008
155 12341.590207 12463.099230 0.076214 0.073105
160 12291.756528 12383.175568 0.073703 0.071259
165 12252.749433 12349.074839 0.070982 0.068442
170 12256.286840 12352.713395 0.067874 0.065349
175 12182.336237 12260.653843 0.066291 0.064265
180 12145.611052 12228.542720 0.063949 0.062020
185 12120.021352 12199.166752 0.061681 0.059829
190 12095.073710 12177.586624 0.059340 0.057465
195 12049.179195 12107.631819 0.057897 0.056255
200 11992.337972 12076.794692 0.056356 0.054539
205 11949.219275 12033.176341 0.054787 0.053045
210 11885.344010 11961.453925 0.053671 0.052030
215 11823.561736 11886.901099 0.052411 0.051071
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 15/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.24191e+43 erg / s Luminosity absorbed = 7.96929e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13110.758198 13157.348776 0.670180 0.660872
5 13113.028903 13154.537628 0.574875 0.567941
10 13093.088462 13122.466176 0.515120 0.508848
15 13078.630835 13143.615683 0.464486 0.455821
20 13073.991303 13184.013843 0.419770 0.406711
25 13110.381168 13204.092069 0.375311 0.367136
30 13055.574882 13126.833694 0.344822 0.339156
35 13041.575051 13129.174789 0.318035 0.310268
40 13022.005042 13109.493361 0.292756 0.285613
45 13003.812183 13130.769795 0.269639 0.259172
50 12963.032879 13078.921737 0.250234 0.241394
55 12982.389906 13088.297094 0.229901 0.222381
60 12975.350494 13098.649549 0.210624 0.203764
65 12965.368804 13049.630202 0.196254 0.189580
70 12890.584285 12990.355074 0.185073 0.179681
75 12893.430696 12968.256977 0.171994 0.167990
80 12901.177598 12970.293312 0.159412 0.156160
85 12894.708333 12968.271549 0.148871 0.145543
90 12834.988058 12945.709957 0.140906 0.136494
95 12833.456390 12931.106463 0.131452 0.127930
100 12807.914014 12919.726451 0.123631 0.120351
105 12725.023233 12876.617692 0.118612 0.113582
110 12715.728557 12815.972255 0.112363 0.108565
115 12681.404241 12832.576329 0.106558 0.102167
120 12659.267376 12757.169548 0.100885 0.097927
125 12647.301366 12772.943537 0.095555 0.091839
130 12601.531043 12694.520268 0.091565 0.088882
135 12567.359999 12652.718277 0.087506 0.084953
140 12518.962460 12658.308413 0.084548 0.080415
145 12478.839839 12630.562968 0.080795 0.077078
150 12494.070660 12597.378869 0.076008 0.073622
155 12463.099230 12559.021645 0.073105 0.070842
160 12383.175568 12498.714833 0.071259 0.068965
165 12349.074839 12440.298871 0.068442 0.066375
170 12352.713395 12430.528073 0.065349 0.063553
175 12260.653843 12372.900884 0.064265 0.062067
180 12228.542720 12348.207758 0.062020 0.059442
185 12199.166752 12294.218663 0.059829 0.057943
190 12177.586624 12281.713874 0.057465 0.055475
195 12107.631819 12194.852959 0.056255 0.054501
200 12076.794692 12140.252777 0.054539 0.053392
205 12033.176341 12102.119966 0.053045 0.051561
210 11961.453925 12043.807386 0.052030 0.050461
215 11886.901099 11969.967457 0.051071 0.049585
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 16/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.23023e+43 erg / s Luminosity absorbed = 8.07641e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13157.348776 13208.359288 0.660872 0.651688
5 13154.537628 13167.132794 0.567941 0.565322
10 13122.466176 13185.668397 0.508848 0.502269
15 13143.615683 13211.047516 0.455821 0.448037
20 13184.013843 13236.736583 0.406711 0.401473
25 13204.092069 13206.804495 0.367136 0.366627
30 13126.833694 13210.054340 0.339156 0.333981
35 13129.174789 13223.266295 0.310268 0.303712
40 13109.493361 13144.587177 0.285613 0.283888
45 13130.769795 13206.257030 0.259172 0.253767
50 13078.921737 13138.179491 0.241394 0.237449
55 13088.297094 13162.707185 0.222381 0.217089
60 13098.649549 13139.213259 0.203764 0.202532
65 13049.630202 13075.253357 0.189580 0.188542
70 12990.355074 12998.806255 0.179681 0.178885
75 12968.256977 13007.626047 0.167990 0.165915
80 12970.293312 12989.924028 0.156160 0.154532
85 12968.271549 12950.778512 0.145543 0.144925
90 12945.709957 12971.021631 0.136494 0.135208
95 12931.106463 12931.545423 0.127930 0.126700
100 12919.726451 12967.279881 0.120351 0.118255
105 12876.617692 12910.093006 0.113582 0.112558
110 12815.972255 12836.599003 0.108565 0.107565
115 12832.576329 12807.968042 0.102167 0.102752
120 12757.169548 12793.413777 0.097927 0.096426
125 12772.943537 12790.716105 0.091839 0.091588
130 12694.520268 12762.147943 0.088882 0.086533
135 12652.718277 12707.570217 0.084953 0.083250
140 12658.308413 12723.808197 0.080415 0.078572
145 12630.562968 12631.971267 0.077078 0.076837
150 12597.378869 12613.669514 0.073622 0.072776
155 12559.021645 12571.136723 0.070842 0.070026
160 12498.714833 12582.788478 0.068965 0.066758
165 12440.298871 12515.372784 0.066375 0.064675
170 12430.528073 12463.372449 0.063553 0.062916
175 12372.900884 12414.922501 0.062067 0.060897
180 12348.207758 12374.558711 0.059442 0.058655
185 12294.218663 12325.073521 0.057943 0.056945
190 12281.713874 12271.369723 0.055475 0.055274
195 12194.852959 12215.455284 0.054501 0.054039
200 12140.252777 12184.530283 0.053392 0.052251
205 12102.119966 12150.063082 0.051561 0.050542
210 12043.807386 12106.155134 0.050461 0.049397
215 11969.967457 12017.087711 0.049585 0.048618
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 17/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.23593e+43 erg / s Luminosity absorbed = 8.01895e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13208.359288 13180.535471 0.651688 0.658012
5 13167.132794 13197.522520 0.565322 0.562610
10 13185.668397 13174.694866 0.502269 0.505905
15 13211.047516 13218.240350 0.448037 0.447772
20 13236.736583 13241.225892 0.401473 0.403604
25 13206.804495 13281.145141 0.366627 0.360381
30 13210.054340 13267.313536 0.333981 0.328384
35 13223.266295 13286.738073 0.303712 0.297458
40 13144.587177 13181.026751 0.283888 0.279818
45 13206.257030 13250.028942 0.253767 0.251366
50 13138.179491 13221.355944 0.237449 0.232221
55 13162.707185 13184.596993 0.217089 0.216556
60 13139.213259 13167.869456 0.202532 0.200682
65 13075.253357 13157.266030 0.188542 0.184969
70 12998.806255 13067.082659 0.178885 0.176158
75 13007.626047 13063.765451 0.165915 0.163820
80 12989.924028 13107.403255 0.154532 0.150006
85 12950.778512 13019.626254 0.144925 0.142593
90 12971.021631 13075.761908 0.135208 0.131608
95 12931.545423 12974.276546 0.126700 0.125865
100 12967.279881 12999.965074 0.118255 0.117449
105 12910.093006 12960.919781 0.112558 0.111632
110 12836.599003 12862.321828 0.107565 0.106997
115 12807.968042 12888.724834 0.102752 0.100015
120 12793.413777 12843.322883 0.096426 0.094920
125 12790.716105 12840.666968 0.091588 0.089666
130 12762.147943 12747.676705 0.086533 0.087161
135 12707.570217 12724.372823 0.083250 0.083115
140 12723.808197 12748.930414 0.078572 0.078337
145 12631.971267 12717.078441 0.076837 0.074678
150 12613.669514 12663.247909 0.072776 0.071613
155 12571.136723 12639.626393 0.070026 0.069149
160 12582.788478 12605.013851 0.066758 0.066224
165 12515.372784 12529.037411 0.064675 0.064437
170 12463.372449 12512.441055 0.062916 0.061855
175 12414.922501 12504.982434 0.060897 0.059360
180 12374.558711 12421.316875 0.058655 0.057855
185 12325.073521 12407.198976 0.056945 0.055458
190 12271.369723 12340.725468 0.055274 0.054148
195 12215.455284 12284.031240 0.054039 0.052758
200 12184.530283 12240.133374 0.052251 0.051248
205 12150.063082 12197.305230 0.050542 0.049849
210 12106.155134 12155.057432 0.049397 0.048580
215 12017.087711 12048.113577 0.048618 0.048016
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 18/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.23486e+43 erg / s Luminosity absorbed = 8.02956e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13180.535471 13166.826599 0.658012 0.661773
5 13197.522520 13161.225439 0.562610 0.570444
10 13174.694866 13164.033321 0.505905 0.508207
15 13218.240350 13201.849346 0.447772 0.450656
20 13241.225892 13207.232370 0.403604 0.408133
25 13281.145141 13265.773013 0.360381 0.360941
30 13267.313536 13256.622272 0.328384 0.330341
35 13286.738073 13303.315390 0.297458 0.296541
40 13181.026751 13239.259991 0.279818 0.276023
45 13250.028942 13284.146284 0.251366 0.251606
50 13221.355944 13281.863495 0.232221 0.229871
55 13184.596993 13291.777175 0.216556 0.210540
60 13167.869456 13273.989481 0.200682 0.195987
65 13157.266030 13234.724167 0.184969 0.181474
70 13067.082659 13073.022860 0.176158 0.176104
75 13063.765451 13104.206567 0.163820 0.161397
80 13107.403255 13144.009619 0.150006 0.148310
85 13019.626254 13095.391976 0.142593 0.139576
90 13075.761908 13127.347738 0.131608 0.129448
95 12974.276546 13047.484719 0.125865 0.123275
100 12999.965074 13054.715390 0.117449 0.115835
105 12960.919781 13030.075993 0.111632 0.109099
110 12862.321828 12933.735701 0.106997 0.104811
115 12888.724834 12932.599283 0.100015 0.099248
120 12843.322883 12887.477861 0.094920 0.094004
125 12840.666968 12856.250561 0.089666 0.089628
130 12747.676705 12868.962924 0.087161 0.083908
135 12724.372823 12788.877268 0.083115 0.081364
140 12748.930414 12795.467458 0.078337 0.077004
145 12717.078441 12763.950779 0.074678 0.073739
150 12663.247909 12718.648984 0.071613 0.070597
155 12639.626393 12715.537044 0.069149 0.067710
160 12605.013851 12650.932274 0.066224 0.065329
165 12529.037411 12607.924877 0.064437 0.062762
170 12512.441055 12608.704017 0.061855 0.059987
175 12504.982434 12586.366176 0.059360 0.057683
180 12421.316875 12528.864095 0.057855 0.055953
185 12407.198976 12477.716927 0.055458 0.054166
190 12340.725468 12414.012399 0.054148 0.052969
195 12284.031240 12346.489093 0.052758 0.051603
200 12240.133374 12285.180423 0.051248 0.050515
205 12197.305230 12227.179556 0.049849 0.049302
210 12155.057432 12166.263653 0.048580 0.048302
215 12048.113577 12071.857550 0.048016 0.047522
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 19/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.22794e+43 erg / s Luminosity absorbed = 8.08940e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13166.826599 13195.519766 0.661773 0.655923
5 13161.225439 13210.741764 0.570444 0.561298
10 13164.033321 13235.889833 0.508207 0.496477
15 13201.849346 13257.976739 0.450656 0.442108
20 13207.232370 13244.393460 0.408133 0.404152
25 13265.773013 13296.492302 0.360941 0.360900
30 13256.622272 13297.435443 0.330341 0.328098
35 13303.315390 13345.323924 0.296541 0.294771
40 13239.259991 13291.955653 0.276023 0.273309
45 13284.146284 13330.280434 0.251606 0.248861
50 13281.863495 13308.757907 0.229871 0.230862
55 13291.777175 13304.285435 0.210540 0.211845
60 13273.989481 13296.137731 0.195987 0.196026
65 13234.724167 13270.557153 0.181474 0.180536
70 13073.022860 13177.400597 0.176104 0.171785
75 13104.206567 13183.943549 0.161397 0.158635
80 13144.009619 13178.261898 0.148310 0.148387
85 13095.391976 13152.438819 0.139576 0.138476
90 13127.347738 13215.251923 0.129448 0.126976
95 13047.484719 13142.709904 0.123275 0.121220
100 13054.715390 13117.686046 0.115835 0.114414
105 13030.075993 13084.549148 0.109099 0.107408
110 12933.735701 13016.108830 0.104811 0.103248
115 12932.599283 13034.016323 0.099248 0.095773
120 12887.477861 12996.130306 0.094004 0.091197
125 12856.250561 12941.186178 0.089628 0.087218
130 12868.962924 12925.433879 0.083908 0.082603
135 12788.877268 12880.728193 0.081364 0.079023
140 12795.467458 12847.699135 0.077004 0.075342
145 12763.950779 12794.111032 0.073739 0.072687
150 12718.648984 12736.936646 0.070597 0.069979
155 12715.537044 12710.772364 0.067710 0.067228
160 12650.932274 12649.594065 0.065329 0.065043
165 12607.924877 12622.622036 0.062762 0.062659
170 12608.704017 12611.817207 0.059987 0.059724
175 12586.366176 12585.348534 0.057683 0.057595
180 12528.864095 12550.454375 0.055953 0.055510
185 12477.716927 12478.754739 0.054166 0.053976
190 12414.012399 12444.539262 0.052969 0.052245
195 12346.489093 12373.740602 0.051603 0.050962
200 12285.180423 12313.701285 0.050515 0.049849
205 12227.179556 12241.245001 0.049302 0.048828
210 12166.263653 12183.394427 0.048302 0.047714
215 12071.857550 12099.018816 0.047522 0.046976
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 20/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.22565e+43 erg / s Luminosity absorbed = 8.10622e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13195.519766 13218.392110 0.655923 0.649349
5 13210.741764 13203.590883 0.561298 0.562997
10 13235.889833 13261.694418 0.496477 0.494806
15 13257.976739 13298.064876 0.442108 0.438235
20 13244.393460 13250.974631 0.404152 0.404452
25 13296.492302 13305.202082 0.360900 0.360661
30 13297.435443 13329.187270 0.328098 0.326259
35 13345.323924 13351.564942 0.294771 0.295542
40 13291.955653 13350.505065 0.273309 0.269488
45 13330.280434 13381.278649 0.248861 0.244352
50 13308.757907 13350.401391 0.230862 0.227866
55 13304.285435 13322.777854 0.211845 0.212514
60 13296.137731 13338.479356 0.196026 0.194321
65 13270.557153 13301.899148 0.180536 0.180111
70 13177.400597 13281.385654 0.171785 0.167801
75 13183.943549 13294.307652 0.158635 0.154838
80 13178.261898 13267.966077 0.148387 0.145283
85 13152.438819 13231.259161 0.138476 0.135669
90 13215.251923 13243.603339 0.126976 0.126946
95 13142.709904 13245.685361 0.121220 0.117504
100 13117.686046 13177.508230 0.114414 0.112469
105 13084.549148 13166.889212 0.107408 0.105195
110 13016.108830 13166.758929 0.103248 0.099258
115 13034.016323 13147.285482 0.095773 0.092737
120 12996.130306 13132.027948 0.091197 0.087745
125 12941.186178 13024.016775 0.087218 0.085256
130 12925.433879 13018.325987 0.082603 0.080448
135 12880.728193 13030.706325 0.079023 0.075968
140 12847.699135 12937.788695 0.075342 0.073766
145 12794.111032 12911.995175 0.072687 0.069983
150 12736.936646 12823.855253 0.069979 0.068040
155 12710.772364 12801.548304 0.067228 0.065332
160 12649.594065 12716.134743 0.065043 0.063471
165 12622.622036 12684.472936 0.062659 0.061499
170 12611.817207 12681.171176 0.059724 0.058254
175 12585.348534 12622.297246 0.057595 0.056786
180 12550.454375 12571.964450 0.055510 0.054959
185 12478.754739 12526.205063 0.053976 0.053273
190 12444.539262 12472.393760 0.052245 0.051593
195 12373.740602 12392.121290 0.050962 0.050498
200 12313.701285 12324.192334 0.049849 0.049509
205 12241.245001 12255.227451 0.048828 0.048525
210 12183.394427 12234.356154 0.047714 0.046820
215 12099.018816 12116.932757 0.046976 0.046558
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 21/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.22962e+43 erg / s Luminosity absorbed = 8.07205e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13218.392110 13193.465240 0.649349 0.653651
5 13203.590883 13227.367950 0.562997 0.559593
10 13261.694418 13274.159820 0.494806 0.492442
15 13298.064876 13293.745040 0.438235 0.439388
20 13250.974631 13288.270014 0.404452 0.399678
25 13305.202082 13292.192089 0.360661 0.361950
30 13329.187270 13317.502835 0.326259 0.325734
35 13351.564942 13401.993759 0.295542 0.290505
40 13350.505065 13359.932654 0.269488 0.270203
45 13381.278649 13407.014769 0.244352 0.243462
50 13350.401391 13409.046073 0.227866 0.224079
55 13322.777854 13340.715423 0.212514 0.210605
60 13338.479356 13381.415152 0.194321 0.191791
65 13301.899148 13377.188342 0.180111 0.176521
70 13281.385654 13334.189209 0.167801 0.165216
75 13294.307652 13348.055264 0.154838 0.152572
80 13267.966077 13322.318785 0.145283 0.142626
85 13231.259161 13272.891041 0.135669 0.133746
90 13243.603339 13280.211749 0.126946 0.125982
95 13245.685361 13267.189239 0.117504 0.116621
100 13177.508230 13215.630844 0.112469 0.111315
105 13166.889212 13175.010755 0.105195 0.105045
110 13166.758929 13200.228786 0.099258 0.097975
115 13147.285482 13162.211276 0.092737 0.092293
120 13132.027948 13106.849919 0.087745 0.088225
125 13024.016775 13055.657025 0.085256 0.084607
130 13018.325987 13012.306316 0.080448 0.080817
135 13030.706325 12997.425810 0.075968 0.076277
140 12937.788695 12945.009959 0.073766 0.073460
145 12911.995175 12918.379242 0.069983 0.070102
150 12823.855253 12811.234197 0.068040 0.068201
155 12801.548304 12773.399696 0.065332 0.066079
160 12716.134743 12751.645015 0.063471 0.063025
165 12684.472936 12724.598687 0.061499 0.060578
170 12681.171176 12676.535769 0.058254 0.058497
175 12622.297246 12617.758328 0.056786 0.056922
180 12571.964450 12566.081512 0.054959 0.055095
185 12526.205063 12519.386548 0.053273 0.053572
190 12472.393760 12501.909282 0.051593 0.051469
195 12392.121290 12437.036377 0.050498 0.049955
200 12324.192334 12352.874273 0.049509 0.049206
205 12255.227451 12282.392519 0.048525 0.048163
210 12234.356154 12217.553356 0.046820 0.047145
215 12116.932757 12106.194987 0.046558 0.046731
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 22/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.22206e+43 erg / s Luminosity absorbed = 8.13794e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13193.465240 13195.466506 0.653651 0.655681
5 13227.367950 13234.149938 0.559593 0.558807
10 13274.159820 13250.119206 0.492442 0.497110
15 13293.745040 13283.570386 0.439388 0.443546
20 13288.270014 13285.196098 0.399678 0.401518
25 13292.192089 13297.903037 0.361950 0.363788
30 13317.502835 13304.903433 0.325734 0.330382
35 13401.993759 13348.630577 0.290505 0.296154
40 13359.932654 13373.208881 0.270203 0.268559
45 13407.014769 13409.696334 0.243462 0.243838
50 13409.046073 13446.583648 0.224079 0.221260
55 13340.715423 13363.747733 0.210605 0.209475
60 13381.415152 13445.718362 0.191791 0.188191
65 13377.188342 13451.866500 0.176521 0.173362
70 13334.189209 13420.709807 0.165216 0.161404
75 13348.055264 13404.576966 0.152572 0.150490
80 13322.318785 13376.205532 0.142626 0.141152
85 13272.891041 13367.313511 0.133746 0.130660
90 13280.211749 13362.190573 0.125982 0.122816
95 13267.189239 13315.019271 0.116621 0.114927
100 13215.630844 13284.482648 0.111315 0.108049
105 13175.010755 13228.775249 0.105045 0.103409
110 13200.228786 13250.373347 0.097975 0.096366
115 13162.211276 13163.949582 0.092293 0.092277
120 13106.849919 13195.199414 0.088225 0.086233
125 13055.657025 13168.769528 0.084607 0.081787
130 13012.306316 13085.969357 0.080817 0.079046
135 12997.425810 13075.101627 0.076277 0.074932
140 12945.009959 12966.424621 0.073460 0.073250
145 12918.379242 12951.468346 0.070102 0.069209
150 12811.234197 12890.777445 0.068201 0.066757
155 12773.399696 12846.206270 0.066079 0.064555
160 12751.645015 12803.339348 0.063025 0.062149
165 12724.598687 12797.147056 0.060578 0.059464
170 12676.535769 12797.437041 0.058497 0.056845
175 12617.758328 12708.051794 0.056922 0.055482
180 12566.081512 12637.715416 0.055095 0.054127
185 12519.386548 12574.299005 0.053572 0.052499
190 12501.909282 12531.863048 0.051469 0.050803
195 12437.036377 12496.462195 0.049955 0.049082
200 12352.874273 12416.783340 0.049206 0.048305
205 12282.392519 12349.739513 0.048163 0.047137
210 12217.553356 12254.084604 0.047145 0.046408
215 12106.194987 12160.167476 0.046731 0.045780
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 23/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.22210e+43 erg / s Luminosity absorbed = 8.13819e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13195.466506 13216.402157 0.655681 0.652541
5 13234.149938 13237.225427 0.558807 0.559506
10 13250.119206 13252.209818 0.497110 0.497424
15 13283.570386 13243.043291 0.443546 0.449203
20 13285.196098 13279.986369 0.401518 0.401513
25 13297.903037 13307.974117 0.363788 0.360489
30 13304.903433 13330.259306 0.330382 0.326376
35 13348.630577 13361.962084 0.296154 0.295078
40 13373.208881 13394.308987 0.268559 0.266887
45 13409.696334 13436.815696 0.243838 0.242688
50 13446.583648 13461.659991 0.221260 0.221027
55 13363.747733 13443.756797 0.209475 0.205664
60 13445.718362 13439.851226 0.188191 0.188685
65 13451.866500 13447.748126 0.173362 0.174022
70 13420.709807 13423.536597 0.161404 0.161184
75 13404.576966 13423.527026 0.150490 0.149976
80 13376.205532 13431.673399 0.141152 0.138622
85 13367.313511 13367.261282 0.130660 0.130878
90 13362.190573 13392.504047 0.122816 0.121643
95 13315.019271 13379.838946 0.114927 0.113123
100 13284.482648 13305.320906 0.108049 0.107698
105 13228.775249 13299.509838 0.103409 0.101494
110 13250.373347 13303.155862 0.096366 0.095391
115 13163.949582 13239.401607 0.092277 0.090675
120 13195.199414 13225.853168 0.086233 0.085758
125 13168.769528 13182.779625 0.081787 0.082242
130 13085.969357 13158.856980 0.079046 0.077699
135 13075.101627 13080.616663 0.074932 0.075268
140 12966.424621 13027.168954 0.073250 0.072122
145 12951.468346 13000.819102 0.069209 0.068449
150 12890.777445 12985.503070 0.066757 0.065123
155 12846.206270 12930.672063 0.064555 0.062974
160 12803.339348 12854.302745 0.062149 0.061227
165 12797.147056 12825.373038 0.059464 0.058748
170 12797.437041 12819.415527 0.056845 0.056058
175 12708.051794 12708.497312 0.055482 0.055067
180 12637.715416 12639.246258 0.054127 0.054019
185 12574.299005 12631.362226 0.052499 0.051403
190 12531.863048 12559.773216 0.050803 0.050265
195 12496.462195 12612.585828 0.049082 0.047577
200 12416.783340 12489.714343 0.048305 0.047266
205 12349.739513 12414.869027 0.047137 0.046149
210 12254.084604 12333.119968 0.046408 0.045325
215 12160.167476 12218.394834 0.045780 0.044907
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 24/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.21424e+43 erg / s Luminosity absorbed = 8.21243e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13216.402157 13214.412341 0.652541 0.654355
5 13237.225427 13266.546018 0.559506 0.557950
10 13252.209818 13266.925438 0.497424 0.495958
15 13243.043291 13321.265822 0.449203 0.439403
20 13279.986369 13300.999426 0.401513 0.399679
25 13307.974117 13340.522690 0.360489 0.359001
30 13330.259306 13360.096557 0.326376 0.324063
35 13361.962084 13425.244338 0.295078 0.290336
40 13394.308987 13406.012752 0.266887 0.267246
45 13436.815696 13470.738156 0.242688 0.239683
50 13461.659991 13514.566779 0.221027 0.216757
55 13443.756797 13446.809708 0.205664 0.204399
60 13439.851226 13495.051511 0.188685 0.185571
65 13447.748126 13478.049586 0.174022 0.171728
70 13423.536597 13478.243438 0.161184 0.158297
75 13423.527026 13447.494233 0.149976 0.148463
80 13431.673399 13436.926051 0.138622 0.138342
85 13367.261282 13405.984926 0.130878 0.129554
90 13392.504047 13411.477844 0.121643 0.120318
95 13379.838946 13397.359911 0.113123 0.112874
100 13305.320906 13360.714798 0.107698 0.106382
105 13299.509838 13331.840059 0.101494 0.100378
110 13303.155862 13283.695143 0.095391 0.095469
115 13239.401607 13269.496703 0.090675 0.089693
120 13225.853168 13216.630334 0.085758 0.085353
125 13182.779625 13188.100904 0.082242 0.081391
130 13158.856980 13171.671286 0.077699 0.077299
135 13080.616663 13114.333697 0.075268 0.074280
140 13027.168954 13087.991186 0.072122 0.070821
145 13000.819102 13030.686619 0.068449 0.067675
150 12985.503070 12997.630280 0.065123 0.064682
155 12930.672063 12953.381884 0.062974 0.062689
160 12854.302745 12869.969481 0.061227 0.060758
165 12825.373038 12850.813423 0.058748 0.058063
170 12819.415527 12822.577348 0.056058 0.055673
175 12708.497312 12719.584918 0.055067 0.054788
180 12639.246258 12641.476771 0.054019 0.053491
185 12631.362226 12625.848643 0.051403 0.051273
190 12559.773216 12633.826899 0.050265 0.049295
195 12612.585828 12568.263736 0.047577 0.047966
200 12489.714343 12497.684613 0.047266 0.046947
205 12414.869027 12403.984719 0.046149 0.046160
210 12333.119968 12336.380238 0.045325 0.045221
215 12218.394834 12223.584261 0.044907 0.044754
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 25/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.21658e+43 erg / s Luminosity absorbed = 8.19057e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13214.412341 13226.637262 0.654355 0.649960
5 13266.546018 13273.317748 0.557950 0.557229
10 13266.925438 13308.672690 0.495958 0.489640
15 13321.265822 13362.955586 0.439403 0.433103
20 13300.999426 13295.175381 0.399679 0.399624
25 13340.522690 13344.623968 0.359001 0.358202
30 13360.096557 13391.442684 0.324063 0.322013
35 13425.244338 13466.503250 0.290336 0.288319
40 13406.012752 13436.991080 0.267246 0.265779
45 13470.738156 13513.082253 0.239683 0.237983
50 13514.566779 13499.586715 0.216757 0.218713
55 13446.809708 13471.751080 0.204399 0.203875
60 13495.051511 13474.365905 0.185571 0.187673
65 13478.049586 13495.089584 0.171728 0.171978
70 13478.243438 13506.466553 0.158297 0.158035
75 13447.494233 13472.618979 0.148463 0.147689
80 13436.926051 13426.092569 0.138342 0.138483
85 13405.984926 13435.279967 0.129554 0.127994
90 13411.477844 13398.549637 0.120318 0.120507
95 13397.359911 13392.899837 0.112874 0.112153
100 13360.714798 13355.280323 0.106382 0.106066
105 13331.840059 13350.876967 0.100378 0.100158
110 13283.695143 13262.987994 0.095469 0.096234
115 13269.496703 13269.261900 0.089693 0.089592
120 13216.630334 13220.271865 0.085353 0.085028
125 13188.100904 13227.413382 0.081391 0.080367
130 13171.671286 13183.623645 0.077299 0.076439
135 13114.333697 13154.819805 0.074280 0.072986
140 13087.991186 13073.616423 0.070821 0.070857
145 13030.686619 13043.634237 0.067675 0.067414
150 12997.630280 12994.570180 0.064682 0.064596
155 12953.381884 12946.036078 0.062689 0.062456
160 12869.969481 12899.812989 0.060758 0.060078
165 12850.813423 12860.432622 0.058063 0.057899
170 12822.577348 12758.171162 0.055673 0.056694
175 12719.584918 12736.508057 0.054788 0.054428
180 12641.476771 12693.757431 0.053491 0.052712
185 12625.848643 12627.684428 0.051273 0.051319
190 12633.826899 12605.649726 0.049295 0.049639
195 12568.263736 12519.805078 0.047966 0.048512
200 12497.684613 12441.375388 0.046947 0.047643
205 12403.984719 12390.758682 0.046160 0.046193
210 12336.380238 12341.490806 0.045221 0.045034
215 12223.584261 12269.951528 0.044754 0.044087
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 26/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.22024e+43 erg / s Luminosity absorbed = 8.15492e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13226.637262 13224.038802 0.649960 0.649671
5 13273.317748 13254.407087 0.557229 0.561481
10 13308.672690 13284.603566 0.489640 0.492310
15 13362.955586 13302.545105 0.433103 0.439342
20 13295.175381 13292.262849 0.399624 0.398956
25 13344.623968 13342.035745 0.358202 0.357113
30 13391.442684 13395.651625 0.322013 0.320898
35 13466.503250 13464.557456 0.288319 0.289676
40 13436.991080 13440.515426 0.265779 0.266129
45 13513.082253 13515.022982 0.237983 0.237762
50 13499.586715 13510.184201 0.218713 0.218257
55 13471.751080 13504.086537 0.203875 0.201539
60 13474.365905 13468.713278 0.187673 0.187599
65 13495.089584 13499.187263 0.171978 0.172129
70 13506.466553 13503.375732 0.158035 0.158278
75 13472.618979 13471.049950 0.147689 0.147953
80 13426.092569 13469.814316 0.138483 0.137510
85 13435.279967 13466.753912 0.127994 0.127882
90 13398.549637 13414.734938 0.120507 0.121051
95 13392.899837 13424.523000 0.112153 0.112245
100 13355.280323 13380.318364 0.106066 0.106457
105 13350.876967 13352.257978 0.100158 0.100303
110 13262.987994 13320.346400 0.096234 0.094864
115 13269.261900 13299.997262 0.089592 0.089149
120 13220.271865 13232.512535 0.085028 0.084875
125 13227.413382 13242.894290 0.080367 0.080056
130 13183.623645 13202.446665 0.076439 0.076082
135 13154.819805 13156.812152 0.072986 0.073040
140 13073.616423 13076.700422 0.070857 0.070527
145 13043.634237 13053.722896 0.067414 0.067399
150 12994.570180 13006.987684 0.064596 0.064253
155 12946.036078 12967.210224 0.062456 0.062075
160 12899.812989 12945.463999 0.060078 0.059531
165 12860.432622 12854.793486 0.057899 0.057858
170 12758.171162 12796.353980 0.056694 0.056173
175 12736.508057 12792.868597 0.054428 0.053488
180 12693.757431 12696.129047 0.052712 0.052492
185 12627.684428 12656.106175 0.051319 0.050890
190 12605.649726 12598.256389 0.049639 0.049612
195 12519.805078 12518.316319 0.048512 0.048575
200 12441.375388 12456.817563 0.047643 0.047361
205 12390.758682 12358.693871 0.046193 0.046701
210 12341.490806 12305.940437 0.045034 0.045632
215 12269.951528 12229.092749 0.044087 0.044648
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 27/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.21669e+43 erg / s Luminosity absorbed = 8.18670e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13224.038802 13227.901174 0.649671 0.649570
5 13254.407087 13253.437801 0.561481 0.562809
10 13284.603566 13315.822899 0.492310 0.487354
15 13302.545105 13308.504973 0.439342 0.441817
20 13292.262849 13313.086028 0.398956 0.398092
25 13342.035745 13362.257074 0.357113 0.355488
30 13395.651625 13385.926898 0.320898 0.322538
35 13464.557456 13455.152434 0.289676 0.290891
40 13440.515426 13468.663938 0.266129 0.263483
45 13515.022982 13519.866805 0.237762 0.238079
50 13510.184201 13506.898673 0.218257 0.219233
55 13504.086537 13549.458452 0.201539 0.199181
60 13468.713278 13509.018362 0.187599 0.185527
65 13499.187263 13486.072442 0.172129 0.172556
70 13503.375732 13525.678046 0.158278 0.157214
75 13471.049950 13482.036613 0.147953 0.146788
80 13469.814316 13503.408657 0.137510 0.135660
85 13466.753912 13518.750962 0.127882 0.125845
90 13414.734938 13448.183311 0.121051 0.118639
95 13424.523000 13466.878706 0.112245 0.110304
100 13380.318364 13393.365097 0.106457 0.105151
105 13352.257978 13320.206391 0.100303 0.100640
110 13320.346400 13296.148156 0.094864 0.095177
115 13299.997262 13277.749438 0.089149 0.089847
120 13232.512535 13255.594266 0.084875 0.084687
125 13242.894290 13216.769722 0.080056 0.080475
130 13202.446665 13134.280941 0.076082 0.077466
135 13156.812152 13143.270696 0.073040 0.073214
140 13076.700422 13118.062526 0.070527 0.069721
145 13053.722896 13129.778037 0.067399 0.066201
150 13006.987684 13016.758140 0.064253 0.064372
155 12967.210224 12952.608913 0.062075 0.062374
160 12945.463999 12972.195659 0.059531 0.058840
165 12854.793486 12901.782640 0.057858 0.057205
170 12796.353980 12836.110004 0.056173 0.055623
175 12792.868597 12813.578279 0.053488 0.053283
180 12696.129047 12746.183324 0.052492 0.052010
185 12656.106175 12642.400168 0.050890 0.051181
190 12598.256389 12623.518838 0.049612 0.049208
195 12518.316319 12545.983504 0.048575 0.048173
200 12456.817563 12486.360144 0.047361 0.046895
205 12358.693871 12395.242697 0.046701 0.046184
210 12305.940437 12343.814221 0.045632 0.045038
215 12229.092749 12244.749981 0.044648 0.044401
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 28/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.21499e+43 erg / s Luminosity absorbed = 8.20205e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13227.901174 13233.150013 0.649570 0.650204
5 13253.437801 13246.153944 0.562809 0.561930
10 13315.822899 13321.909830 0.487354 0.489863
15 13308.504973 13364.961838 0.441817 0.435686
20 13313.086028 13391.543417 0.398092 0.390502
25 13362.257074 13380.970341 0.355488 0.354792
30 13385.926898 13424.798756 0.322538 0.319017
35 13455.152434 13461.928210 0.290891 0.291085
40 13468.663938 13503.499504 0.263483 0.261981
45 13519.866805 13495.866176 0.238079 0.240783
50 13506.898673 13555.871313 0.219233 0.216347
55 13549.458452 13571.589585 0.199181 0.197516
60 13509.018362 13533.420573 0.185527 0.184918
65 13486.072442 13516.614447 0.172556 0.171709
70 13525.678046 13524.108128 0.157214 0.158169
75 13482.036613 13504.251745 0.146788 0.146163
80 13503.408657 13508.314541 0.135660 0.135793
85 13518.750962 13538.331114 0.125845 0.125334
90 13448.183311 13460.652047 0.118639 0.118766
95 13466.878706 13463.743682 0.110304 0.111095
100 13393.365097 13382.045279 0.105151 0.105867
105 13320.206391 13378.572343 0.100640 0.099143
110 13296.148156 13379.622903 0.095177 0.092730
115 13277.749438 13274.010882 0.089847 0.089717
120 13255.594266 13289.780395 0.084687 0.083850
125 13216.769722 13209.527751 0.080475 0.080657
130 13134.280941 13174.722631 0.077466 0.076697
135 13143.270696 13135.588395 0.073214 0.073239
140 13118.062526 13089.038929 0.069721 0.070019
145 13129.778037 13107.025066 0.066201 0.066344
150 13016.758140 13017.046389 0.064372 0.064109
155 12952.608913 12909.922348 0.062374 0.062857
160 12972.195659 12948.664891 0.058840 0.059189
165 12901.782640 12903.756393 0.057205 0.057100
170 12836.110004 12845.378254 0.055623 0.055307
175 12813.578279 12803.645892 0.053283 0.053442
180 12746.183324 12723.611287 0.052010 0.052088
185 12642.400168 12634.318515 0.051181 0.051017
190 12623.518838 12601.606254 0.049208 0.049256
195 12545.983504 12565.203662 0.048173 0.047612
200 12486.360144 12494.861637 0.046895 0.046639
205 12395.242697 12397.928412 0.046184 0.045991
210 12343.814221 12336.477565 0.045038 0.045007
215 12244.749981 12249.970878 0.044401 0.044239
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 29/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.21179e+43 erg / s Luminosity absorbed = 8.22822e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 13233.150013 13208.310381 0.650204 0.657074
5 13246.153944 13242.926346 0.561930 0.561964
10 13321.909830 13330.763628 0.489863 0.489069
15 13364.961838 13368.288133 0.435686 0.433913
20 13391.543417 13384.044346 0.390502 0.389619
25 13380.970341 13398.903276 0.354792 0.351207
30 13424.798756 13448.522939 0.319017 0.316476
35 13461.928210 13467.082069 0.291085 0.290285
40 13503.499504 13508.888789 0.261981 0.263058
45 13495.866176 13525.001658 0.240783 0.239508
50 13555.871313 13562.425891 0.216347 0.217229
55 13571.589585 13577.718177 0.197516 0.197810
60 13533.420573 13597.504461 0.184918 0.181660
65 13516.614447 13601.489270 0.171709 0.168048
70 13524.108128 13568.698227 0.158169 0.156417
75 13504.251745 13546.834728 0.146163 0.145109
80 13508.314541 13555.479137 0.135793 0.133817
85 13538.331114 13571.340270 0.125334 0.124472
90 13460.652047 13472.928583 0.118766 0.119160
95 13463.743682 13457.100412 0.111095 0.111430
100 13382.045279 13428.556594 0.105867 0.105074
105 13378.572343 13366.311341 0.099143 0.100007
110 13379.622903 13379.869001 0.092730 0.093092
115 13274.010882 13275.411938 0.089717 0.090079
120 13289.780395 13246.539043 0.083850 0.085249
125 13209.527751 13249.121435 0.080657 0.079777
130 13174.722631 13135.663120 0.076697 0.077560
135 13135.588395 13110.646608 0.073239 0.073567
140 13089.038929 13087.479523 0.070019 0.070224
145 13107.025066 13044.884829 0.066344 0.067281
150 13017.046389 13037.484772 0.064109 0.063956
155 12909.922348 12941.347490 0.062857 0.062526
160 12948.664891 12942.316969 0.059189 0.059498
165 12903.756393 12933.641487 0.057100 0.056895
170 12845.378254 12859.107981 0.055307 0.055298
175 12803.645892 12800.474636 0.053442 0.053562
180 12723.611287 12750.762210 0.052088 0.052045
185 12634.318515 12681.428533 0.051017 0.050617
190 12601.606254 12600.917221 0.049256 0.049396
195 12565.203662 12543.499820 0.047612 0.047945
200 12494.861637 12490.808295 0.046639 0.046748
205 12397.928412 12415.000297 0.045991 0.045719
210 12336.477565 12336.243932 0.045007 0.044929
215 12249.970878 12253.441375 0.044239 0.044120
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 30/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 1.21579e+43 erg / s Luminosity absorbed = 8.19099e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Simulation finished in 30 iterations and took 302.03 s (base.py:306)
[8]:
spectrum = pd.DataFrame(data = sim.runner.spectrum_integrated.frequency.value, columns=['frequency'])
spectrum['luminosity'] = sim.runner.spectrum_integrated.luminosity_density_nu.value
create_df_download_link(spectrum, title='Download Spectrum', filename='spectrum.csv', type='csv')
[8]:
[9]:
spec_artis_toy = pd.read_csv('spec_artis_toy.csv')
[10]:
%pylab notebook
plot(sim.runner.spectrum_integrated.frequency, sim.runner.spectrum_integrated.luminosity_density_nu, label='TARDIS')
#plot(sim.runner.spectrum_integrated.wavelength, sim.runner.spectrum_integrated.luminosity_density_lambda)
plot(spec_artis_toy['nu'], spec_artis_toy.iloc[:,53] * np.pi * 4 * u.Mpc.to(u.cm)**2, label='ARTIS')
legend()
xlim(0, 1e15)
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/matplotlib/style/core.py:167: UserWarning: In /Users/wkerzend/.matplotlib/stylelib/talk_default.mplstyle:
The text.latex.unicode rcparam was deprecated in Matplotlib 2.2 and will be removed in 3.1.
styles = read_style_directory(stylelib_path)
(warnings.py:99)
Populating the interactive namespace from numpy and matplotlib
[10]:
(0, 1000000000000000.0)
Comparison of snia_toy06.dat¶
[11]:
blondin_dict, blondin_csv = read_blondin_toymodel('snia_toy06.dat')
[12]:
blondin_dict['v_inner_boundary'] = '9000 km/s'
blondin_dict['v_outer_boundary'] = '20000 km/s'
[13]:
csvy_file = '---\n{0}\n---\n{1}'.format(yaml.dump(blondin_dict, default_flow_style=False), blondin_csv.to_csv(index=False))
Download the model and setup
[14]:
with open('blondin_compare_06.csvy', 'w') as fh:
fh.write(csvy_file)
[15]:
sim = run_tardis('blondin_model_compare_06.yml')
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.
AstropyDeprecationWarning)
(warnings.py:99)
[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/wkerzend/projects/tardis/tardis-data/kurucz_cd23_chianti_H_He.h5 (util.py:29)
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/IPython/core/interactiveshell.py:3267: PerformanceWarning: indexing past lexsort depth may impact performance.
exec(code_obj, self.user_global_ns, self.user_ns)
(warnings.py:99)
[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)
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.
AstropyDeprecationWarning)
(warnings.py:99)
[py.warnings ][WARNING] /Users/wkerzend/python/tardis/tardis/plasma/properties/atomic.py:195: FutureWarning: .labels was deprecated in version 0.24.0. Use .codes instead.
zeta_data['atomic_number'] = zeta_data.index.labels[0] + 1
(warnings.py:99)
[py.warnings ][WARNING] /Users/wkerzend/python/tardis/tardis/plasma/properties/atomic.py:196: FutureWarning: .labels was deprecated in version 0.24.0. Use .codes instead.
zeta_data['ion_number'] = zeta_data.index.labels[1] + 1
(warnings.py:99)
[tardis.plasma.properties.atomic][WARNING] Zeta_data missing - replaced with 1s. Missing ions: [(14, 15), (16, 17), (20, 21), (26, 27), (27, 28), (28, 29)] (atomic.py:215)
[py.warnings ][WARNING] /Users/wkerzend/python/tardis/tardis/plasma/properties/ion_population.py:63: FutureWarning:
Passing list-likes to .loc or [] with any missing label will raise
KeyError in the future, you can use .reindex() as an alternative.
See the documentation here:
https://pandas.pydata.org/pandas-docs/stable/indexing.html#deprecate-loc-reindex-listlike
partition_function.index].dropna())
(warnings.py:99)
[py.warnings ][WARNING] /Users/wkerzend/python/tardis/tardis/plasma/properties/ion_population.py:134: UserWarning: t_rads outside of zeta factor interpolation zeta_min=2000.00 zeta_max=40000.00 - replacing with 1s
t_rad))
(warnings.py:99)
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/equivalencies.py:90: RuntimeWarning: divide by zero encountered in double_scalars
(si.m, si.Hz, lambda x: _si.c.value / x),
(warnings.py:99)
[py.warnings ][WARNING] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.
AstropyDeprecationWarning)
(warnings.py:99)
[tardis.simulation.base][INFO ] Starting iteration 1/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 8.58316e+42 erg / s Luminosity absorbed = 1.14927e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 48232.0 33078.660848 0.462810 0.339407
5 45376.0 33359.208475 0.379149 0.240036
10 42491.0 32247.702823 0.336294 0.209552
15 39588.0 30940.509627 0.304883 0.189092
20 36670.0 29419.035049 0.279753 0.175165
25 33741.0 27862.998412 0.258765 0.162926
30 30797.0 26056.304566 0.240774 0.155731
35 27831.0 24233.374024 0.225078 0.150757
40 24827.0 22259.666051 0.211205 0.148148
45 21754.0 20052.519837 0.198821 0.149908
50 18561.0 17572.332905 0.187679 0.159651
55 15146.0 14878.504699 0.177588 0.178974
60 11272.0 11894.797719 0.168399 0.219075
65 6063.9 8017.697008 0.159992 0.321883
70 5000.0 6857.681232 0.152269 0.409785
75 5000.0 6706.225807 0.145148 0.420703
80 5000.0 6624.827933 0.138563 0.412427
85 5000.0 6559.646701 0.132453 0.401986
90 5000.0 6498.094778 0.126772 0.394748
95 5000.0 6486.173446 0.121475 0.375484
100 5000.0 6442.241085 0.116527 0.367131
105 5000.0 6396.521785 0.111894 0.359269
110 5000.0 6369.954386 0.107550 0.348460
115 5000.0 6365.069785 0.103468 0.330502
120 5000.0 6335.033359 0.099626 0.324085
125 5000.0 6335.240062 0.096006 0.307071
130 5000.0 6312.003849 0.092589 0.298308
135 5000.0 6315.460113 0.089359 0.283135
140 5000.0 6287.968101 0.086303 0.273750
145 5000.0 6283.131863 0.083408 0.262965
150 5000.0 6288.219630 0.080662 0.249380
155 5000.0 6273.416628 0.078054 0.237795
160 5000.0 6252.886407 0.075575 0.233189
165 5000.0 6239.557304 0.073217 0.224598
170 5000.0 6217.990468 0.070971 0.218918
175 5000.0 6204.862349 0.068830 0.210945
180 5000.0 6190.728150 0.066787 0.204948
185 5000.0 6194.861311 0.064837 0.195274
190 5000.0 6186.214267 0.062973 0.187556
195 5000.0 6163.168514 0.061191 0.183025
200 5000.0 6151.945647 0.059485 0.177255
205 5000.0 6155.450367 0.057852 0.169708
210 5000.0 6145.333602 0.056286 0.163614
215 5000.0 6112.801838 0.054785 0.160652
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 2/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 8.22032e+42 erg / s Luminosity absorbed = 1.18799e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 33078.660848 25646.127577 0.339407 0.272341
5 33359.208475 27779.011309 0.240036 0.166440
10 32247.702823 27647.402177 0.209552 0.142136
15 30940.509627 26956.243849 0.189092 0.130006
20 29419.035049 26110.749534 0.175165 0.121196
25 27862.998412 25184.797764 0.162926 0.114157
30 26056.304566 23903.310867 0.155731 0.111695
35 24233.374024 22464.949449 0.150757 0.111841
40 22259.666051 20807.749039 0.148148 0.114950
45 20052.519837 19138.312801 0.149908 0.120132
50 17572.332905 17049.758069 0.159651 0.135918
55 14878.504699 14744.241654 0.178974 0.161523
60 11894.797719 12368.873229 0.219075 0.204673
65 8017.697008 9747.118793 0.321883 0.282912
70 6857.681232 8504.295946 0.409785 0.360680
75 6706.225807 8261.138040 0.420703 0.372982
80 6624.827933 8131.977329 0.412427 0.367405
85 6559.646701 7986.892879 0.401986 0.364399
90 6498.094778 7900.213479 0.394748 0.359227
95 6486.173446 7866.326390 0.375484 0.341438
100 6442.241085 7801.310057 0.367131 0.334603
105 6396.521785 7734.044408 0.359269 0.328341
110 6369.954386 7687.907215 0.348460 0.320147
115 6365.069785 7693.182800 0.330502 0.301387
120 6335.033359 7655.515919 0.324085 0.293835
125 6335.240062 7620.713571 0.307071 0.282380
130 6312.003849 7580.429413 0.298308 0.274228
135 6315.460113 7566.309064 0.283135 0.261058
140 6287.968101 7550.789323 0.273750 0.250787
145 6283.131863 7535.937992 0.262965 0.240000
150 6288.219630 7512.312121 0.249380 0.230931
155 6273.416628 7478.890013 0.237795 0.223031
160 6252.886407 7475.047982 0.233189 0.215373
165 6239.557304 7468.079400 0.224598 0.206387
170 6217.990468 7451.004860 0.218918 0.198988
175 6204.862349 7409.557929 0.210945 0.194654
180 6190.728150 7404.275923 0.204948 0.187906
185 6194.861311 7390.741925 0.195274 0.180247
190 6186.214267 7383.089094 0.187556 0.172665
195 6163.168514 7341.277893 0.183025 0.169132
200 6151.945647 7320.617789 0.177255 0.164262
205 6155.450367 7305.229567 0.169708 0.157482
210 6145.333602 7287.420458 0.163614 0.153307
215 6112.801838 7260.224191 0.160652 0.148549
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 3/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 7.53288e+42 erg / s Luminosity absorbed = 1.25536e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 25646.127577 21687.341639 0.272341 0.253627
5 27779.011309 24469.580180 0.166440 0.140872
10 27647.402177 24909.976327 0.142136 0.116365
15 26956.243849 24656.427206 0.130006 0.105583
20 26110.749534 24165.648308 0.121196 0.098350
25 25184.797764 23409.557036 0.114157 0.094399
30 23903.310867 22273.546772 0.111695 0.094505
35 22464.949449 21015.842734 0.111841 0.097335
40 20807.749039 19598.445762 0.114950 0.102435
45 19138.312801 18194.697753 0.120132 0.108754
50 17049.758069 16474.690030 0.135918 0.123592
55 14744.241654 14426.658678 0.161523 0.151056
60 12368.873229 12409.246992 0.204673 0.193485
65 9747.118793 10566.675904 0.282912 0.250930
70 8504.295946 9641.544377 0.360680 0.298794
75 8261.138040 9425.314448 0.372982 0.303862
80 8131.977329 9237.406931 0.367405 0.302762
85 7986.892879 9076.762770 0.364399 0.302978
90 7900.213479 8945.885086 0.359227 0.299366
95 7866.326390 8894.168206 0.341438 0.286980
100 7801.310057 8808.146138 0.334603 0.279830
105 7734.044408 8709.008609 0.328341 0.275653
110 7687.907215 8641.546213 0.320147 0.269139
115 7693.182800 8594.743373 0.301387 0.257106
120 7655.515919 8547.204869 0.293835 0.251014
125 7620.713571 8497.141809 0.282380 0.243241
130 7580.429413 8470.974386 0.274228 0.234358
135 7566.309064 8436.190944 0.261058 0.223792
140 7550.789323 8426.221081 0.250787 0.213680
145 7535.937992 8395.605145 0.240000 0.205239
150 7512.312121 8367.423110 0.230931 0.197109
155 7478.890013 8314.965938 0.223031 0.190926
160 7475.047982 8303.325524 0.215373 0.185388
165 7468.079400 8294.385045 0.206387 0.176284
170 7451.004860 8277.336042 0.198988 0.169414
175 7409.557929 8222.330703 0.194654 0.165971
180 7404.275923 8204.174904 0.187906 0.160077
185 7390.741925 8183.696969 0.180247 0.154208
190 7383.089094 8159.654987 0.172665 0.148128
195 7341.277893 8129.675443 0.169132 0.144477
200 7320.617789 8108.159374 0.164262 0.140187
205 7305.229567 8087.266758 0.157482 0.135504
210 7287.420458 8055.881172 0.153307 0.130894
215 7260.224191 8039.969032 0.148549 0.126358
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 4/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.94811e+42 erg / s Luminosity absorbed = 1.31524e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 21687.341639 19337.291241 0.253627 0.270153
5 24469.580180 22592.595297 0.140872 0.136920
10 24909.976327 23285.393456 0.116365 0.109207
15 24656.427206 23122.516303 0.105583 0.098554
20 24165.648308 22709.273345 0.098350 0.092178
25 23409.557036 21980.507524 0.094399 0.089448
30 22273.546772 20907.299125 0.094505 0.090564
35 21015.842734 19878.350356 0.097335 0.092653
40 19598.445762 18595.512520 0.102435 0.098847
45 18194.697753 17339.057651 0.108754 0.106149
50 16474.690030 15794.453426 0.123592 0.120273
55 14426.658678 14033.865213 0.151056 0.144822
60 12409.246992 12258.337451 0.193485 0.187633
65 10566.675904 10849.973249 0.250930 0.233686
70 9641.544377 10196.773086 0.298794 0.260404
75 9425.314448 10011.550967 0.303862 0.258982
80 9237.406931 9888.790308 0.302762 0.252716
85 9076.762770 9760.412746 0.302978 0.249461
90 8945.885086 9596.708364 0.299366 0.247828
95 8894.168206 9549.777257 0.286980 0.235729
100 8808.146138 9432.561960 0.279830 0.232594
105 8709.008609 9347.553652 0.275653 0.226183
110 8641.546213 9278.843901 0.269139 0.220521
115 8594.743373 9253.390705 0.257106 0.209843
120 8547.204869 9200.145781 0.251014 0.204575
125 8497.141809 9119.450917 0.243241 0.199674
130 8470.974386 9068.333536 0.234358 0.193760
135 8436.190944 9000.651683 0.223792 0.188192
140 8426.221081 9003.427215 0.213680 0.177094
145 8395.605145 8990.573400 0.205239 0.168436
150 8367.423110 8941.331178 0.197109 0.162812
155 8314.965938 8903.292918 0.190926 0.158086
160 8303.325524 8872.487593 0.185388 0.153814
165 8294.385045 8853.965174 0.176284 0.146998
170 8277.336042 8831.227420 0.169414 0.140738
175 8222.330703 8751.598029 0.165971 0.139141
180 8204.174904 8757.994139 0.160077 0.133106
185 8183.696969 8748.758957 0.154208 0.127799
190 8159.654987 8739.461589 0.148128 0.122879
195 8129.675443 8705.022057 0.144477 0.119449
200 8108.159374 8682.105330 0.140187 0.115308
205 8087.266758 8647.041572 0.135504 0.111847
210 8055.881172 8611.942072 0.130894 0.108252
215 8039.969032 8603.779801 0.126358 0.103887
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 5/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.89650e+42 erg / s Luminosity absorbed = 1.32067e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 19337.291241 17785.274576 0.270153 0.306363
5 22592.595297 21098.722865 0.136920 0.150394
10 23285.393456 21889.872853 0.109207 0.117945
15 23122.516303 21814.999257 0.098554 0.105379
20 22709.273345 21484.857487 0.092178 0.097576
25 21980.507524 20868.094624 0.089448 0.094251
30 20907.299125 19961.904347 0.090564 0.095136
35 19878.350356 19040.763859 0.092653 0.097343
40 18595.512520 17894.437326 0.098847 0.103005
45 17339.057651 16807.693319 0.106149 0.107974
50 15794.453426 15362.328757 0.120273 0.122100
55 14033.865213 13619.150411 0.144822 0.150211
60 12258.337451 12017.493716 0.187633 0.193773
65 10849.973249 10910.161548 0.233686 0.231460
70 10196.773086 10460.227307 0.260404 0.243385
75 10011.550967 10319.638202 0.258982 0.235796
80 9888.790308 10171.471779 0.252716 0.232786
85 9760.412746 10095.034363 0.249461 0.223191
90 9596.708364 10001.062608 0.247828 0.215620
95 9549.777257 9957.963976 0.235729 0.204434
100 9432.561960 9835.949600 0.232594 0.202949
105 9347.553652 9790.333119 0.226183 0.193053
110 9278.843901 9685.887970 0.220521 0.189031
115 9253.390705 9668.390274 0.209843 0.179259
120 9200.145781 9625.488665 0.204575 0.172673
125 9119.450917 9544.269710 0.199674 0.169047
130 9068.333536 9482.576963 0.193760 0.163112
135 9000.651683 9434.093136 0.188192 0.158375
140 9003.427215 9427.040736 0.177094 0.149212
145 8990.573400 9393.195721 0.168436 0.143074
150 8941.331178 9363.940241 0.162812 0.137745
155 8903.292918 9308.334814 0.158086 0.134036
160 8872.487593 9284.914172 0.153814 0.129720
165 8853.965174 9262.204725 0.146998 0.124334
170 8831.227420 9257.308024 0.140738 0.118703
175 8751.598029 9181.075956 0.139141 0.117517
180 8757.994139 9137.539664 0.133106 0.113970
185 8748.758957 9117.715299 0.127799 0.109614
190 8739.461589 9099.968913 0.122879 0.105135
195 8705.022057 9125.450166 0.119449 0.100419
200 8682.105330 9088.203900 0.115308 0.097487
205 8647.041572 9035.134612 0.111847 0.095245
210 8611.942072 9014.882318 0.108252 0.091627
215 8603.779801 8988.733605 0.103887 0.088560
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 6/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.81705e+42 erg / s Luminosity absorbed = 1.32768e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 17785.274576 16849.548861 0.306363 0.339130
5 21098.722865 19982.766676 0.150394 0.168907
10 21889.872853 20814.585317 0.117945 0.132588
15 21814.999257 20836.878477 0.105379 0.116916
20 21484.857487 20562.035214 0.097576 0.108511
25 20868.094624 20073.230073 0.094251 0.103236
30 19961.904347 19318.326762 0.095136 0.102658
35 19040.763859 18473.243508 0.097343 0.103772
40 17894.437326 17473.192259 0.103005 0.108271
45 16807.693319 16351.507427 0.107974 0.114545
50 15362.328757 14936.801780 0.122100 0.128970
55 13619.150411 13344.023456 0.150211 0.156169
60 12017.493716 11852.938239 0.193773 0.199002
65 10910.161548 10902.808440 0.231460 0.231027
70 10460.227307 10539.550594 0.243385 0.237316
75 10319.638202 10409.216001 0.235796 0.230407
80 10171.471779 10345.447196 0.232786 0.220309
85 10095.034363 10265.687128 0.223191 0.210362
90 10001.062608 10200.853852 0.215620 0.200839
95 9957.963976 10128.975173 0.204434 0.191325
100 9835.949600 10083.193955 0.202949 0.183343
105 9790.333119 10065.440414 0.193053 0.173337
110 9685.887970 9946.066164 0.189031 0.169955
115 9668.390274 9941.685970 0.179259 0.160889
120 9625.488665 9855.119991 0.172673 0.156500
125 9544.269710 9767.913949 0.169047 0.153230
130 9482.576963 9741.390095 0.163112 0.145382
135 9434.093136 9699.488937 0.158375 0.141292
140 9427.040736 9699.638927 0.149212 0.133026
145 9393.195721 9647.447500 0.143074 0.127877
150 9363.940241 9658.721246 0.137745 0.121130
155 9308.334814 9557.671246 0.134036 0.119500
160 9284.914172 9548.568868 0.129720 0.114677
165 9262.204725 9518.814190 0.124334 0.110231
170 9257.308024 9536.709699 0.118703 0.104798
175 9181.075956 9472.954238 0.117517 0.102829
180 9137.539664 9430.024609 0.113970 0.099642
185 9117.715299 9448.801554 0.109614 0.094832
190 9099.968913 9382.195031 0.105135 0.092690
195 9125.450166 9369.232768 0.100419 0.089226
200 9088.203900 9326.009435 0.097487 0.086697
205 9035.134612 9283.305714 0.095245 0.084703
210 9014.882318 9286.349681 0.091627 0.080969
215 8988.733605 9238.697233 0.088560 0.078968
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 7/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.92481e+42 erg / s Luminosity absorbed = 1.31825e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 16849.548861 16209.755588 0.339130 0.370417
5 19982.766676 19077.723696 0.168907 0.191210
10 20814.585317 19911.636284 0.132588 0.149223
15 20836.878477 19980.515105 0.116916 0.131862
20 20562.035214 19730.173414 0.108511 0.122104
25 20073.230073 19351.203251 0.103236 0.114934
30 19318.326762 18708.163597 0.102658 0.112862
35 18473.243508 17970.785739 0.103772 0.112251
40 17473.192259 17014.462048 0.108271 0.115784
45 16351.507427 15913.141457 0.114545 0.123177
50 14936.801780 14504.845122 0.128970 0.139496
55 13344.023456 12991.456904 0.156169 0.170236
60 11852.938239 11647.196885 0.199002 0.211901
65 10902.808440 10883.735491 0.231027 0.233586
70 10539.550594 10569.881836 0.237316 0.237468
75 10409.216001 10440.537574 0.230407 0.230231
80 10345.447196 10380.171949 0.220309 0.217920
85 10265.687128 10330.778351 0.210362 0.207626
90 10200.853852 10261.329728 0.200839 0.199073
95 10128.975173 10197.887886 0.191325 0.189033
100 10083.193955 10136.908527 0.183343 0.183036
105 10065.440414 10178.032441 0.173337 0.167191
110 9946.066164 10053.893696 0.169955 0.165158
115 9941.685970 10048.474963 0.160889 0.154498
120 9855.119991 9977.360250 0.156500 0.149005
125 9767.913949 9930.208870 0.153230 0.143528
130 9741.390095 9915.603925 0.145382 0.136017
135 9699.488937 9866.200366 0.141292 0.131511
140 9699.638927 9814.951031 0.133026 0.126961
145 9647.447500 9798.071238 0.127877 0.120533
150 9658.721246 9808.797163 0.121130 0.113663
155 9557.671246 9793.162005 0.119500 0.108977
160 9548.568868 9725.343672 0.114677 0.106289
165 9518.814190 9698.623879 0.110231 0.102227
170 9536.709699 9689.497056 0.104798 0.097970
175 9472.954238 9639.823072 0.102829 0.095198
180 9430.024609 9580.700640 0.099642 0.093124
185 9448.801554 9567.030309 0.094832 0.089470
190 9382.195031 9515.358467 0.092690 0.087009
195 9369.232768 9542.712152 0.089226 0.082585
200 9326.009435 9489.230843 0.086697 0.080734
205 9283.305714 9459.468892 0.084703 0.078167
210 9286.349681 9417.719379 0.080969 0.075863
215 9238.697233 9392.331097 0.078968 0.073609
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 8/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.81528e+42 erg / s Luminosity absorbed = 1.32994e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 16209.755588 15749.985522 0.370417 0.402278
5 19077.723696 18407.128630 0.191210 0.213512
10 19911.636284 19251.095237 0.149223 0.165039
15 19980.515105 19350.147374 0.131862 0.144273
20 19730.173414 19223.160804 0.122104 0.131226
25 19351.203251 18903.086783 0.114934 0.122602
30 18708.163597 18263.389823 0.112862 0.120931
35 17970.785739 17583.304766 0.112251 0.118426
40 17014.462048 16653.952428 0.115784 0.121839
45 15913.141457 15489.395480 0.123177 0.131262
50 14504.845122 14098.052167 0.139496 0.150432
55 12991.456904 12695.327152 0.170236 0.181755
60 11647.196885 11464.075634 0.211901 0.221590
65 10883.735491 10765.874609 0.233586 0.241479
70 10569.881836 10471.508089 0.237468 0.244550
75 10440.537574 10381.748027 0.230231 0.233189
80 10380.171949 10333.422176 0.217920 0.221514
85 10330.778351 10316.315191 0.207626 0.206450
90 10261.329728 10263.249639 0.199073 0.197000
95 10197.887886 10202.942850 0.189033 0.187554
100 10136.908527 10114.185235 0.183036 0.182414
105 10178.032441 10116.105048 0.167191 0.169095
110 10053.893696 10016.312332 0.165158 0.165532
115 10048.474963 10003.578670 0.154498 0.156078
120 9977.360250 10000.442227 0.149005 0.145977
125 9930.208870 9936.034025 0.143528 0.140875
130 9915.603925 9958.866790 0.136017 0.132945
135 9866.200366 9870.055986 0.131511 0.129994
140 9814.951031 9831.256270 0.126961 0.123960
145 9798.071238 9849.863589 0.120533 0.116870
150 9808.797163 9852.694348 0.113663 0.110614
155 9793.162005 9819.575708 0.108977 0.106679
160 9725.343672 9758.739520 0.106289 0.103636
165 9698.623879 9728.138440 0.102227 0.099698
170 9689.497056 9715.576792 0.097970 0.095781
175 9639.823072 9666.766972 0.095198 0.093382
180 9580.700640 9630.084835 0.093124 0.090494
185 9567.030309 9583.778725 0.089470 0.087729
190 9515.358467 9551.138946 0.087009 0.085357
195 9542.712152 9573.820919 0.082585 0.080667
200 9489.230843 9506.186949 0.080734 0.079204
205 9459.468892 9469.368983 0.078167 0.076713
210 9417.719379 9433.468467 0.075863 0.074497
215 9392.331097 9411.402869 0.073609 0.072252
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 9/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.74528e+42 erg / s Luminosity absorbed = 1.33788e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 15749.985522 15458.227551 0.402278 0.427376
5 18407.128630 17798.747665 0.213512 0.235997
10 19251.095237 18657.315752 0.165039 0.181372
15 19350.147374 18805.169772 0.144273 0.158406
20 19223.160804 18723.570753 0.131226 0.142359
25 18903.086783 18431.317315 0.122602 0.132699
30 18263.389823 17838.583491 0.120931 0.129882
35 17583.304766 17223.765177 0.118426 0.126423
40 16653.952428 16303.461774 0.121839 0.129604
45 15489.395480 15147.982602 0.131262 0.138048
50 14098.052167 13738.955797 0.150432 0.160394
55 12695.327152 12406.424011 0.181755 0.193471
60 11464.075634 11270.721996 0.221590 0.233103
65 10765.874609 10574.233335 0.241479 0.256338
70 10471.508089 10298.222974 0.244550 0.259000
75 10381.748027 10214.316462 0.233189 0.248513
80 10333.422176 10202.758039 0.221514 0.232142
85 10316.315191 10166.178983 0.206450 0.216863
90 10263.249639 10115.022778 0.197000 0.205195
95 10202.942850 10080.125264 0.187554 0.192643
100 10114.185235 10048.036098 0.182414 0.184170
105 10116.105048 10046.911911 0.169095 0.171504
110 10016.312332 9981.713697 0.165532 0.164493
115 10003.578670 9927.658845 0.156078 0.157250
120 10000.442227 9945.885957 0.145977 0.146737
125 9936.034025 9855.570556 0.140875 0.143805
130 9958.866790 9879.558216 0.132945 0.135379
135 9870.055986 9839.980238 0.129994 0.128972
140 9831.256270 9791.478622 0.123960 0.123739
145 9849.863589 9780.099754 0.116870 0.117938
150 9852.694348 9778.970512 0.110614 0.111651
155 9819.575708 9751.095673 0.106679 0.107421
160 9758.739520 9730.283796 0.103636 0.103129
165 9728.138440 9673.844772 0.099698 0.100064
170 9715.576792 9701.542658 0.095781 0.094503
175 9666.766972 9631.239591 0.093382 0.093347
180 9630.084835 9596.608164 0.090494 0.090629
185 9583.778725 9555.849242 0.087729 0.087554
190 9551.138946 9540.583353 0.085357 0.084597
195 9573.820919 9546.132704 0.080667 0.080753
200 9506.186949 9492.322767 0.079204 0.078968
205 9469.368983 9428.838067 0.076713 0.077329
210 9433.468467 9399.859205 0.074497 0.074898
215 9411.402869 9390.115487 0.072252 0.072139
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 10/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.62875e+42 erg / s Luminosity absorbed = 1.35048e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 15458.227551 15259.085998 0.427376 0.444473
5 17798.747665 17355.691821 0.235997 0.255545
10 18657.315752 18157.173276 0.181372 0.196750
15 18805.169772 18349.330494 0.158406 0.169240
20 18723.570753 18277.884283 0.142359 0.151949
25 18431.317315 17992.846563 0.132699 0.141698
30 17838.583491 17444.630868 0.129882 0.136951
35 17223.765177 16745.087671 0.126423 0.135394
40 16303.461774 15888.514666 0.129604 0.138124
45 15147.982602 14787.685134 0.138048 0.146045
50 13738.955797 13402.515209 0.160394 0.171287
55 12406.424011 12139.191744 0.193471 0.204516
60 11270.721996 11086.341547 0.233103 0.242949
65 10574.233335 10451.528280 0.256338 0.266372
70 10298.222974 10217.549030 0.259000 0.263445
75 10214.316462 10121.607904 0.248513 0.252799
80 10202.758039 10056.072995 0.232142 0.242837
85 10166.178983 10027.549221 0.216863 0.224978
90 10115.022778 10006.260486 0.205195 0.211536
95 10080.125264 9970.058892 0.192643 0.199547
100 10048.036098 9952.284172 0.184170 0.188162
105 10046.911911 9922.490726 0.171504 0.177315
110 9981.713697 9851.368434 0.164493 0.170877
115 9927.658845 9855.665720 0.157250 0.160089
120 9945.885957 9777.411625 0.146737 0.156330
125 9855.570556 9749.340898 0.143805 0.149945
130 9879.558216 9749.231328 0.135379 0.139874
135 9839.980238 9751.632383 0.128972 0.131955
140 9791.478622 9741.260648 0.123739 0.124980
145 9780.099754 9708.857324 0.117938 0.119626
150 9778.970512 9654.790111 0.111651 0.116149
155 9751.095673 9655.682112 0.107421 0.110371
160 9730.283796 9627.661123 0.103129 0.106483
165 9673.844772 9619.013629 0.100064 0.101149
170 9701.542658 9598.251328 0.094503 0.096994
175 9631.239591 9542.422787 0.093347 0.095185
180 9596.608164 9503.474428 0.090629 0.092382
185 9555.849242 9503.126726 0.087554 0.088157
190 9540.583353 9466.285373 0.084597 0.085764
195 9546.132704 9474.500261 0.080753 0.081921
200 9492.322767 9429.140151 0.078968 0.079748
205 9428.838067 9377.256444 0.077329 0.077953
210 9399.859205 9321.736943 0.074898 0.076356
215 9390.115487 9308.095919 0.072139 0.073361
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 11/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.63952e+42 erg / s Luminosity absorbed = 1.34988e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 15259.085998 15069.724222 0.444473 0.464318
5 17355.691821 17013.949284 0.255545 0.272917
10 18157.173276 17767.853217 0.196750 0.211868
15 18349.330494 17986.842558 0.169240 0.181022
20 18277.884283 17925.349868 0.151949 0.162523
25 17992.846563 17673.353533 0.141698 0.150667
30 17444.630868 17128.254440 0.136951 0.145715
35 16745.087671 16443.113346 0.135394 0.142570
40 15888.514666 15546.604684 0.138124 0.146858
45 14787.685134 14412.527255 0.146045 0.156775
50 13402.515209 13084.050786 0.171287 0.184305
55 12139.191744 11903.396592 0.204516 0.217070
60 11086.341547 10971.633564 0.242949 0.249812
65 10451.528280 10371.162612 0.266372 0.269907
70 10217.549030 10144.176823 0.263445 0.266762
75 10121.607904 10026.528323 0.252799 0.258907
80 10056.072995 9967.778334 0.242837 0.246971
85 10027.549221 9947.688810 0.224978 0.230995
90 10006.260486 9922.877208 0.211536 0.217542
95 9970.058892 9902.150412 0.199547 0.202696
100 9952.284172 9878.071641 0.188162 0.191287
105 9922.490726 9807.207571 0.177315 0.182578
110 9851.368434 9776.057031 0.170877 0.174491
115 9855.665720 9793.114726 0.160089 0.161623
120 9777.411625 9720.241714 0.156330 0.157286
125 9749.340898 9700.181396 0.149945 0.149985
130 9749.231328 9695.255664 0.139874 0.140861
135 9751.632383 9651.303851 0.131955 0.136572
140 9741.260648 9636.967179 0.124980 0.129646
145 9708.857324 9635.475972 0.119626 0.121945
150 9654.790111 9571.647264 0.116149 0.118896
155 9655.682112 9586.792983 0.110371 0.112436
160 9627.661123 9538.190414 0.106483 0.109002
165 9619.013629 9517.650604 0.101149 0.104434
170 9598.251328 9532.764735 0.096994 0.099100
175 9542.422787 9441.089708 0.095185 0.098744
180 9503.474428 9430.644418 0.092382 0.094151
185 9503.126726 9412.404881 0.088157 0.090928
190 9466.285373 9369.906921 0.085764 0.088580
195 9474.500261 9389.026063 0.081921 0.084101
200 9429.140151 9341.755304 0.079748 0.082064
205 9377.256444 9306.842129 0.077953 0.079302
210 9321.736943 9271.674811 0.076356 0.076932
215 9308.095919 9265.772634 0.073361 0.074147
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 12/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.58967e+42 erg / s Luminosity absorbed = 1.35563e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 15069.724222 14917.931458 0.464318 0.481757
5 17013.949284 16681.387543 0.272917 0.292093
10 17767.853217 17391.570403 0.211868 0.226847
15 17986.842558 17612.334432 0.181022 0.193661
20 17925.349868 17584.544355 0.162523 0.172830
25 17673.353533 17308.500502 0.150667 0.161009
30 17128.254440 16848.358222 0.145715 0.152924
35 16443.113346 16088.743313 0.142570 0.151807
40 15546.604684 15092.105047 0.146858 0.158782
45 14412.527255 13989.676053 0.156775 0.171054
50 13084.050786 12767.134507 0.184305 0.197526
55 11903.396592 11665.941714 0.217070 0.230936
60 10971.633564 10763.915769 0.249812 0.268287
65 10371.162612 10208.482170 0.269907 0.284396
70 10144.176823 10000.818128 0.266762 0.278384
75 10026.528323 9953.990206 0.258907 0.262874
80 9967.778334 9892.024533 0.246971 0.250412
85 9947.688810 9879.684118 0.230995 0.235809
90 9922.877208 9832.746752 0.217542 0.222649
95 9902.150412 9789.877020 0.202696 0.210288
100 9878.071641 9725.587076 0.191287 0.201976
105 9807.207571 9711.994277 0.182578 0.188699
110 9776.057031 9672.974010 0.174491 0.180883
115 9793.114726 9719.612578 0.161623 0.166094
120 9720.241714 9608.051693 0.157286 0.164237
125 9700.181396 9566.806427 0.149985 0.157185
130 9695.255664 9562.261448 0.140861 0.146888
135 9651.303851 9538.230091 0.136572 0.141509
140 9636.967179 9571.976011 0.129646 0.131045
145 9635.475972 9524.539575 0.121945 0.127051
150 9571.647264 9507.344602 0.118896 0.121224
155 9586.792983 9434.355287 0.112436 0.118623
160 9538.190414 9420.685142 0.109002 0.112938
165 9517.650604 9440.729541 0.104434 0.107137
170 9532.764735 9416.143475 0.099100 0.104008
175 9441.089708 9346.265534 0.098744 0.102112
180 9430.644418 9390.130265 0.094151 0.095653
185 9412.404881 9320.179683 0.090928 0.093619
190 9369.906921 9291.694311 0.088580 0.090622
195 9389.026063 9299.466455 0.084101 0.086813
200 9341.755304 9246.966845 0.082064 0.084704
205 9306.842129 9214.024411 0.079302 0.082041
210 9271.674811 9170.332366 0.076932 0.079649
215 9265.772634 9162.214771 0.074147 0.077050
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 13/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.56386e+42 erg / s Luminosity absorbed = 1.35896e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14917.931458 14816.291445 0.481757 0.495850
5 16681.387543 16445.160515 0.292093 0.306152
10 17391.570403 17149.320789 0.226847 0.236870
15 17612.334432 17371.480447 0.193661 0.202580
20 17584.544355 17309.773020 0.172830 0.181312
25 17308.500502 17014.985923 0.161009 0.169232
30 16848.358222 16530.310429 0.152924 0.159637
35 16088.743313 15829.826213 0.151807 0.157504
40 15092.105047 14808.324702 0.158782 0.164809
45 13989.676053 13718.544582 0.171054 0.179536
50 12767.134507 12550.334846 0.197526 0.206525
55 11665.941714 11504.542052 0.230936 0.240235
60 10763.915769 10588.447548 0.268287 0.281583
65 10208.482170 10050.750010 0.284396 0.297889
70 10000.818128 9892.426668 0.278384 0.288697
75 9953.990206 9829.570089 0.262874 0.274457
80 9892.024533 9814.637274 0.250412 0.255496
85 9879.684118 9784.637874 0.235809 0.242381
90 9832.746752 9737.840160 0.222649 0.228370
95 9789.877020 9708.488504 0.210288 0.214617
100 9725.587076 9674.176935 0.201976 0.203487
105 9711.994277 9628.213177 0.188699 0.194663
110 9672.974010 9613.777761 0.180883 0.182497
115 9719.612578 9625.509249 0.166094 0.170664
120 9608.051693 9559.893129 0.164237 0.166325
125 9566.806427 9553.080449 0.157185 0.156412
130 9562.261448 9494.968899 0.146888 0.150355
135 9538.230091 9495.810034 0.141509 0.143296
140 9571.976011 9543.012317 0.131045 0.131393
145 9524.539575 9481.775076 0.127051 0.128468
150 9507.344602 9431.919378 0.121224 0.124537
155 9434.355287 9388.814400 0.118623 0.119913
160 9420.685142 9373.393445 0.112938 0.114866
165 9440.729541 9377.920730 0.107137 0.109589
170 9416.143475 9342.015533 0.104008 0.106559
175 9346.265534 9321.429425 0.102112 0.102646
180 9390.130265 9300.726449 0.095653 0.098686
185 9320.179683 9291.181373 0.093619 0.094208
190 9291.694311 9254.708746 0.090622 0.091497
195 9299.466455 9235.683723 0.086813 0.088512
200 9246.966845 9193.500046 0.084704 0.085947
205 9214.024411 9152.647188 0.082041 0.083762
210 9170.332366 9098.416126 0.079649 0.081872
215 9162.214771 9086.843405 0.077050 0.079206
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 14/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.51880e+42 erg / s Luminosity absorbed = 1.36449e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14816.291445 14741.524854 0.495850 0.505536
5 16445.160515 16302.427661 0.306152 0.315487
10 17149.320789 16921.612407 0.236870 0.247221
15 17371.480447 17079.673386 0.202580 0.212766
20 17309.773020 17010.805933 0.181312 0.189605
25 17014.985923 16681.660854 0.169232 0.176962
30 16530.310429 16188.860361 0.159637 0.167142
35 15829.826213 15452.165154 0.157504 0.166547
40 14808.324702 14505.502335 0.164809 0.173630
45 13718.544582 13400.103464 0.179536 0.191037
50 12550.334846 12278.064504 0.206525 0.218663
55 11504.542052 11302.809564 0.240235 0.251424
60 10588.447548 10464.116701 0.281583 0.293720
65 10050.750010 9919.047392 0.297889 0.311649
70 9892.426668 9754.714369 0.288697 0.304236
75 9829.570089 9662.426041 0.274457 0.292435
80 9814.637274 9721.976230 0.255496 0.263300
85 9784.637874 9655.293326 0.242381 0.251554
90 9737.840160 9573.766158 0.228370 0.242013
95 9708.488504 9556.790991 0.214617 0.227626
100 9674.176935 9543.057395 0.203487 0.214705
105 9628.213177 9490.271762 0.194663 0.204473
110 9613.777761 9491.404647 0.182497 0.190737
115 9625.509249 9471.396094 0.170664 0.180957
120 9559.893129 9461.898387 0.166325 0.170591
125 9553.080449 9447.671563 0.156412 0.161765
130 9494.968899 9406.334636 0.150355 0.155673
135 9495.810034 9360.098680 0.143296 0.150447
140 9543.012317 9388.773609 0.131393 0.140886
145 9481.775076 9367.112154 0.128468 0.134344
150 9431.919378 9347.950436 0.124537 0.128462
155 9388.814400 9303.913877 0.119913 0.124519
160 9373.393445 9246.031233 0.114866 0.121280
165 9377.920730 9265.111188 0.109589 0.114830
170 9342.015533 9233.722837 0.106559 0.110679
175 9321.429425 9202.246731 0.102646 0.107131
180 9300.726449 9218.089627 0.098686 0.101477
185 9291.181373 9208.648590 0.094208 0.096988
190 9254.708746 9187.962618 0.091497 0.093813
195 9235.683723 9162.184582 0.088512 0.090789
200 9193.500046 9154.961055 0.085947 0.086963
205 9152.647188 9100.935891 0.083762 0.085243
210 9098.416126 9049.060820 0.081872 0.083494
215 9086.843405 9007.634087 0.079206 0.081732
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 15/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.46813e+42 erg / s Luminosity absorbed = 1.37003e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14741.524854 14700.797294 0.505536 0.508817
5 16302.427661 16129.029953 0.315487 0.325829
10 16921.612407 16739.045746 0.247221 0.255196
15 17079.673386 16877.600406 0.212766 0.219726
20 17010.805933 16753.163119 0.189605 0.199363
25 16681.660854 16440.254561 0.176962 0.183494
30 16188.860361 15944.187087 0.167142 0.175580
35 15452.165154 15198.689814 0.166547 0.174046
40 14505.502335 14284.787110 0.173630 0.179963
45 13400.103464 13136.785727 0.191037 0.201713
50 12278.064504 12122.287873 0.218663 0.225248
55 11302.809564 11201.141803 0.251424 0.257274
60 10464.116701 10422.125928 0.293720 0.293363
65 9919.047392 9858.883913 0.311649 0.314089
70 9754.714369 9668.876272 0.304236 0.307468
75 9662.426041 9638.238290 0.292435 0.287820
80 9721.976230 9619.594622 0.263300 0.267199
85 9655.293326 9612.433962 0.251554 0.251954
90 9573.766158 9518.071270 0.242013 0.243992
95 9556.790991 9482.387456 0.227626 0.231885
100 9543.057395 9480.895313 0.214705 0.217102
105 9490.271762 9417.286049 0.204473 0.208232
110 9491.404647 9390.263510 0.190737 0.197420
115 9471.396094 9390.444305 0.180957 0.184848
120 9461.898387 9366.964720 0.170591 0.175256
125 9447.671563 9364.346951 0.161765 0.165437
130 9406.334636 9359.909760 0.155673 0.157065
135 9360.098680 9331.935249 0.150447 0.151009
140 9388.773609 9295.094988 0.140886 0.145477
145 9367.112154 9286.546131 0.134344 0.138151
150 9347.950436 9250.332848 0.128462 0.133313
155 9303.913877 9224.966666 0.124519 0.127288
160 9246.031233 9192.896201 0.121280 0.122356
165 9265.111188 9169.143728 0.114830 0.117800
170 9233.722837 9133.827174 0.110679 0.114357
175 9202.246731 9151.399328 0.107131 0.107990
180 9218.089627 9159.106703 0.101477 0.103143
185 9208.648590 9129.238399 0.096988 0.099730
190 9187.962618 9104.711875 0.093813 0.097049
195 9162.184582 9088.501384 0.090789 0.093356
200 9154.961055 9052.867418 0.086963 0.090758
205 9100.935891 9006.932103 0.085243 0.088057
210 9049.060820 8943.939575 0.083494 0.087197
215 9007.634087 8929.434956 0.081732 0.084054
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 16/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.55666e+42 erg / s Luminosity absorbed = 1.36172e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14700.797294 14675.336961 0.508817 0.510940
5 16129.029953 15994.651940 0.325829 0.333397
10 16739.045746 16570.419311 0.255196 0.264239
15 16877.600406 16688.537285 0.219726 0.227988
20 16753.163119 16527.258321 0.199363 0.207736
25 16440.254561 16196.252686 0.183494 0.191596
30 15944.187087 15704.334770 0.175580 0.184110
35 15198.689814 14943.780055 0.174046 0.182320
40 14284.787110 13949.133956 0.179963 0.191907
45 13136.785727 12893.597371 0.201713 0.212320
50 12122.287873 11912.477314 0.225248 0.240345
55 11201.141803 11083.883451 0.257274 0.269699
60 10422.125928 10344.582215 0.293363 0.304216
65 9858.883913 9763.974272 0.314089 0.330538
70 9668.876272 9604.801670 0.307468 0.317302
75 9638.238290 9559.687164 0.287820 0.298306
80 9619.594622 9527.009613 0.267199 0.279763
85 9612.433962 9506.391059 0.251954 0.263241
90 9518.071270 9504.460716 0.243992 0.245333
95 9482.387456 9406.362307 0.231885 0.239755
100 9480.895313 9408.214083 0.217102 0.222648
105 9417.286049 9364.636404 0.208232 0.212189
110 9390.263510 9307.327448 0.197420 0.204706
115 9390.444305 9279.999842 0.184848 0.193863
120 9366.964720 9238.041769 0.175256 0.185912
125 9364.346951 9257.451708 0.165437 0.175557
130 9359.909760 9260.204316 0.157065 0.164358
135 9331.935249 9213.576964 0.151009 0.159131
140 9295.094988 9199.068252 0.145477 0.151302
145 9286.546131 9205.022557 0.138151 0.144217
150 9250.332848 9204.560955 0.133313 0.136431
155 9224.966666 9169.216870 0.127288 0.130946
160 9192.896201 9146.924436 0.122356 0.125373
165 9169.143728 9108.426283 0.117800 0.121309
170 9133.827174 9081.036369 0.114357 0.118747
175 9151.399328 9057.063766 0.107990 0.113409
180 9159.106703 9050.563627 0.103143 0.108158
185 9129.238399 9042.510978 0.099730 0.103755
190 9104.711875 9018.719150 0.097049 0.100633
195 9088.501384 9028.408547 0.093356 0.095858
200 9052.867418 8959.354153 0.090758 0.095162
205 9006.932103 8922.859135 0.088057 0.092022
210 8943.939575 8894.757395 0.087197 0.089289
215 8929.434956 8845.408250 0.084054 0.087283
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 17/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.49691e+42 erg / s Luminosity absorbed = 1.36833e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14675.336961 14626.156540 0.510940 0.519047
5 15994.651940 15798.406058 0.333397 0.347915
10 16570.419311 16382.176183 0.264239 0.274867
15 16688.537285 16504.579428 0.227988 0.236788
20 16527.258321 16330.206061 0.207736 0.214864
25 16196.252686 15991.828329 0.191596 0.199845
30 15704.334770 15434.233842 0.184110 0.192475
35 14943.780055 14736.494432 0.182320 0.188556
40 13949.133956 13699.679329 0.191907 0.201966
45 12893.597371 12644.427534 0.212320 0.225883
50 11912.477314 11721.726841 0.240345 0.255238
55 11083.883451 10896.082659 0.269699 0.284146
60 10344.582215 10227.732980 0.304216 0.313576
65 9763.974272 9685.186908 0.330538 0.336527
70 9604.801670 9502.090724 0.317302 0.327945
75 9559.687164 9489.897966 0.298306 0.304092
80 9527.009613 9441.992597 0.279763 0.290588
85 9506.391059 9401.961200 0.263241 0.274427
90 9504.460716 9363.768277 0.245333 0.259647
95 9406.362307 9351.826940 0.239755 0.243144
100 9408.214083 9318.497010 0.222648 0.229994
105 9364.636404 9242.492713 0.212189 0.221486
110 9307.327448 9196.472404 0.204706 0.211307
115 9279.999842 9213.199956 0.193863 0.196662
120 9238.041769 9194.100721 0.185912 0.188332
125 9257.451708 9190.854195 0.175557 0.178530
130 9260.204316 9172.882103 0.164358 0.169729
135 9213.576964 9124.013241 0.159131 0.163895
140 9199.068252 9104.156431 0.151302 0.156559
145 9205.022557 9089.286509 0.144217 0.150299
150 9204.560955 9119.597239 0.136431 0.139739
155 9169.216870 9084.318714 0.130946 0.134658
160 9146.924436 9087.474700 0.125373 0.127661
165 9108.426283 9043.318543 0.121309 0.124503
170 9081.036369 9023.675631 0.118747 0.120826
175 9057.063766 8989.289310 0.113409 0.116466
180 9050.563627 8945.305962 0.108158 0.113446
185 9042.510978 8946.064687 0.103755 0.108459
190 9018.719150 8935.589065 0.100633 0.103492
195 9028.408547 8930.965039 0.095858 0.100165
200 8959.354153 8856.477621 0.095162 0.099208
205 8922.859135 8827.972351 0.092022 0.095769
210 8894.757395 8807.030814 0.089289 0.092573
215 8845.408250 8772.677924 0.087283 0.089805
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 18/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.51856e+42 erg / s Luminosity absorbed = 1.36659e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14626.156540 14563.147157 0.519047 0.527477
5 15798.406058 15672.389708 0.347915 0.359316
10 16382.176183 16220.017575 0.274867 0.283670
15 16504.579428 16333.737367 0.236788 0.244483
20 16330.206061 16161.729309 0.214864 0.221440
25 15991.828329 15728.408075 0.199845 0.209432
30 15434.233842 15188.263972 0.192475 0.200810
35 14736.494432 14495.726168 0.188556 0.197478
40 13699.679329 13462.900081 0.201966 0.212977
45 12644.427534 12427.498300 0.225883 0.240283
50 11721.726841 11591.690666 0.255238 0.263974
55 10896.082659 10772.459489 0.284146 0.295610
60 10227.732980 10121.320815 0.313576 0.324126
65 9685.186908 9584.404138 0.336527 0.349799
70 9502.090724 9393.681202 0.327945 0.342690
75 9489.897966 9346.974373 0.304092 0.322710
80 9441.992597 9286.769360 0.290588 0.311528
85 9401.961200 9303.547786 0.274427 0.284364
90 9363.768277 9266.201840 0.259647 0.268510
95 9351.826940 9254.889332 0.243144 0.252807
100 9318.497010 9220.120143 0.229994 0.238726
105 9242.492713 9207.942238 0.221486 0.224330
110 9196.472404 9192.809108 0.211307 0.213342
115 9213.199956 9204.304596 0.196662 0.197957
120 9194.100721 9150.351587 0.188332 0.191046
125 9190.854195 9119.837117 0.178530 0.184071
130 9172.882103 9133.071948 0.169729 0.172065
135 9124.013241 9073.685301 0.163895 0.167074
140 9104.156431 9055.870731 0.156559 0.160187
145 9089.286509 9028.264139 0.150299 0.154956
150 9119.597239 9025.478387 0.139739 0.144432
155 9084.318714 9030.253642 0.134658 0.137974
160 9087.474700 9018.519081 0.127661 0.132430
165 9043.318543 8956.523647 0.124503 0.130171
170 9023.675631 8964.270422 0.120826 0.123201
175 8989.289310 8922.482934 0.116466 0.119316
180 8945.305962 8909.570644 0.113446 0.115127
185 8946.064687 8917.350045 0.108459 0.109979
190 8935.589065 8875.604508 0.103492 0.106627
195 8930.965039 8851.674478 0.100165 0.103824
200 8856.477621 8814.122706 0.099208 0.100959
205 8827.972351 8779.145997 0.095769 0.097804
210 8807.030814 8753.530725 0.092573 0.094794
215 8772.677924 8733.551320 0.089805 0.091637
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 19/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.47707e+42 erg / s Luminosity absorbed = 1.37123e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14563.147157 14558.443269 0.527477 0.526808
5 15672.389708 15591.220418 0.359316 0.364858
10 16220.017575 16105.831741 0.283670 0.289515
15 16333.737367 16182.627782 0.244483 0.250794
20 16161.729309 16019.028642 0.221440 0.225371
25 15728.408075 15546.172496 0.209432 0.213277
30 15188.263972 14986.708483 0.200810 0.205112
35 14495.726168 14206.597030 0.197478 0.206946
40 13462.900081 13214.438927 0.212977 0.224297
45 12427.498300 12258.061088 0.240283 0.248411
50 11591.690666 11425.143136 0.263974 0.275795
55 10772.459489 10651.207876 0.295610 0.305247
60 10121.320815 10044.565974 0.324126 0.331662
65 9584.404138 9482.481764 0.349799 0.361361
70 9393.681202 9343.869393 0.342690 0.349708
75 9346.974373 9254.723777 0.322710 0.334508
80 9286.769360 9218.220310 0.311528 0.316625
85 9303.547786 9205.216824 0.284364 0.294723
90 9266.201840 9152.494346 0.268510 0.281090
95 9254.889332 9182.662852 0.252807 0.259026
100 9220.120143 9153.846783 0.238726 0.244104
105 9207.942238 9149.291329 0.224330 0.229578
110 9192.809108 9151.623093 0.213342 0.217216
115 9204.304596 9135.507252 0.197957 0.203908
120 9150.351587 9084.836347 0.191046 0.194848
125 9119.837117 9078.196409 0.184071 0.186688
130 9133.071948 9058.951103 0.172065 0.178687
135 9073.685301 9026.263992 0.167074 0.171056
140 9055.870731 9012.748184 0.160187 0.163922
145 9028.264139 8954.207909 0.154956 0.159354
150 9025.478387 8928.402657 0.144432 0.151011
155 9030.253642 8963.528102 0.137974 0.142079
160 9018.519081 8938.075415 0.132430 0.136851
165 8956.523647 8888.234477 0.130171 0.133802
170 8964.270422 8879.311161 0.123201 0.127947
175 8922.482934 8852.385238 0.119316 0.122870
180 8909.570644 8835.395724 0.115127 0.118942
185 8917.350045 8801.279301 0.109979 0.116151
190 8875.604508 8786.595443 0.106627 0.110797
195 8851.674478 8774.440587 0.103824 0.106874
200 8814.122706 8750.446628 0.100959 0.103055
205 8779.145997 8714.174646 0.097804 0.100483
210 8753.530725 8676.651315 0.094794 0.097982
215 8733.551320 8647.921184 0.091637 0.095233
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 20/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.54309e+42 erg / s Luminosity absorbed = 1.36510e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14558.443269 14527.510345 0.526808 0.532630
5 15591.220418 15483.937166 0.364858 0.372757
10 16105.831741 15969.646152 0.289515 0.297564
15 16182.627782 16056.775058 0.250794 0.256721
20 16019.028642 15865.180054 0.225371 0.232714
25 15546.172496 15387.050754 0.213277 0.220639
30 14986.708483 14774.849411 0.205112 0.214992
35 14206.597030 13953.260106 0.206946 0.217126
40 13214.438927 12972.799163 0.224297 0.239591
45 12258.061088 12083.399950 0.248411 0.262562
50 11425.143136 11263.594252 0.275795 0.292612
55 10651.207876 10576.864392 0.305247 0.314595
60 10044.565974 9987.435062 0.331662 0.337785
65 9482.481764 9398.676439 0.361361 0.373720
70 9343.869393 9239.473273 0.349708 0.362801
75 9254.723777 9199.715933 0.334508 0.341280
80 9218.220310 9166.418629 0.316625 0.321497
85 9205.216824 9173.953482 0.294723 0.297312
90 9152.494346 9124.813448 0.281090 0.284189
95 9182.662852 9131.809028 0.259026 0.262673
100 9153.846783 9094.894714 0.244104 0.250903
105 9149.291329 9087.344461 0.229578 0.236913
110 9151.623093 9041.887966 0.217216 0.228464
115 9135.507252 9013.462778 0.203908 0.215648
120 9084.836347 8996.927961 0.194848 0.205739
125 9078.196409 9011.093827 0.186688 0.193318
130 9058.951103 9011.613656 0.178687 0.181391
135 9026.263992 8946.943223 0.171056 0.177542
140 9012.748184 8941.099654 0.163922 0.168864
145 8954.207909 8880.697126 0.159354 0.164240
150 8928.402657 8869.726647 0.151011 0.155424
155 8963.528102 8900.349932 0.142079 0.146134
160 8938.075415 8850.320461 0.136851 0.142646
165 8888.234477 8838.440444 0.133802 0.137430
170 8879.311161 8829.209084 0.127947 0.131442
175 8852.385238 8791.232594 0.122870 0.127573
180 8835.395724 8766.217387 0.118942 0.123450
185 8801.279301 8750.839919 0.116151 0.118872
190 8786.595443 8726.856451 0.110797 0.114460
195 8774.440587 8700.005387 0.106874 0.111333
200 8750.446628 8697.315947 0.103055 0.106434
205 8714.174646 8692.351416 0.100483 0.102481
210 8676.651315 8653.201242 0.097982 0.099687
215 8647.921184 8602.118014 0.095233 0.097555
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 21/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.53652e+42 erg / s Luminosity absorbed = 1.36571e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14527.510345 14473.750936 0.532630 0.539420
5 15483.937166 15394.891469 0.372757 0.382897
10 15969.646152 15830.219839 0.297564 0.308088
15 16056.775058 15921.169104 0.256721 0.264342
20 15865.180054 15695.876561 0.232714 0.239606
25 15387.050754 15232.148714 0.220639 0.227645
30 14774.849411 14581.300817 0.214992 0.222243
35 13953.260106 13712.675877 0.217126 0.229782
40 12972.799163 12790.709179 0.239591 0.248271
45 12083.399950 11957.930823 0.262562 0.270906
50 11263.594252 11208.295764 0.292612 0.296500
55 10576.864392 10525.823724 0.314595 0.317217
60 9987.435062 9912.190537 0.337785 0.347295
65 9398.676439 9370.834115 0.373720 0.375902
70 9239.473273 9188.390217 0.362801 0.366786
75 9199.715933 9150.887903 0.341280 0.346870
80 9166.418629 9116.253822 0.321497 0.330431
85 9173.953482 9074.482396 0.297312 0.312219
90 9124.813448 9083.906896 0.284189 0.289921
95 9131.809028 9071.406685 0.262673 0.271828
100 9094.894714 9034.050780 0.250903 0.259545
105 9087.344461 9020.845402 0.236913 0.243679
110 9041.887966 9008.943817 0.228464 0.229227
115 9013.462778 8993.272549 0.215648 0.218263
120 8996.927961 8935.633387 0.205739 0.212411
125 9011.093827 8948.980261 0.193318 0.200199
130 9011.613656 8943.734949 0.181391 0.187033
135 8946.943223 8881.660936 0.177542 0.183569
140 8941.099654 8892.798021 0.168864 0.173944
145 8880.697126 8844.964402 0.164240 0.167491
150 8869.726647 8835.188012 0.155424 0.159309
155 8900.349932 8801.211139 0.146134 0.154711
160 8850.320461 8808.910540 0.142646 0.147143
165 8838.440444 8809.512980 0.137430 0.139741
170 8829.209084 8782.298215 0.131442 0.134675
175 8791.232594 8718.629476 0.127573 0.132374
180 8766.217387 8730.921502 0.123450 0.126125
185 8750.839919 8681.938497 0.118872 0.122719
190 8726.856451 8680.800852 0.114460 0.117310
195 8700.005387 8644.028642 0.111333 0.114752
200 8697.315947 8646.380875 0.106434 0.109396
205 8692.351416 8620.944955 0.102481 0.106071
210 8653.201242 8615.237945 0.099687 0.101676
215 8602.118014 8574.353124 0.097555 0.098878
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 22/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.47268e+42 erg / s Luminosity absorbed = 1.37300e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14473.750936 14466.699060 0.539420 0.541765
5 15394.891469 15379.257104 0.382897 0.382662
10 15830.219839 15730.556762 0.308088 0.312481
15 15921.169104 15839.987581 0.264342 0.266450
20 15695.876561 15549.993314 0.239606 0.244379
25 15232.148714 15034.693551 0.227645 0.234763
30 14581.300817 14386.349902 0.222243 0.229501
35 13712.675877 13537.848480 0.229782 0.238753
40 12790.709179 12672.745873 0.248271 0.254498
45 11957.930823 11883.354388 0.270906 0.274462
50 11208.295764 11143.884144 0.296500 0.299847
55 10525.823724 10478.176725 0.317217 0.322281
60 9912.190537 9864.967730 0.347295 0.349796
65 9370.834115 9316.705476 0.375902 0.379320
70 9188.390217 9181.792936 0.366786 0.366714
75 9150.887903 9108.642539 0.346870 0.348132
80 9116.253822 9084.295103 0.330431 0.328849
85 9074.482396 9111.357500 0.312219 0.302365
90 9083.906896 9069.522485 0.289921 0.286715
95 9071.406685 9022.555097 0.271828 0.273829
100 9034.050780 9029.769617 0.259545 0.257064
105 9020.845402 8976.672363 0.243679 0.245541
110 9008.943817 9012.549872 0.229227 0.226280
115 8993.272549 8971.977997 0.218263 0.216976
120 8935.633387 8904.790154 0.212411 0.212322
125 8948.980261 8919.544981 0.200199 0.200177
130 8943.734949 8907.158270 0.187033 0.188602
135 8881.660936 8885.055501 0.183569 0.180874
140 8892.798021 8833.700643 0.173944 0.176308
145 8844.964402 8801.127863 0.167491 0.168953
150 8835.188012 8824.797408 0.159309 0.158740
155 8801.211139 8801.151948 0.154711 0.153259
160 8808.910540 8790.722508 0.147143 0.147153
165 8809.512980 8766.134921 0.139741 0.140835
170 8782.298215 8733.853914 0.134675 0.136649
175 8718.629476 8678.905305 0.132374 0.133883
180 8730.921502 8713.240812 0.126125 0.126614
185 8681.938497 8690.867087 0.122719 0.121146
190 8680.800852 8674.415597 0.117310 0.117001
195 8644.028642 8625.598670 0.114752 0.114779
200 8646.380875 8611.261354 0.109396 0.110876
205 8620.944955 8581.897782 0.106071 0.107376
210 8615.237945 8566.974887 0.101676 0.103806
215 8574.353124 8523.003297 0.098878 0.100800
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 23/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.51623e+42 erg / s Luminosity absorbed = 1.36856e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14466.699060 14470.495196 0.541765 0.540916
5 15379.257104 15344.013647 0.382662 0.384238
10 15730.556762 15650.287302 0.312481 0.316687
15 15839.987581 15649.137611 0.266450 0.275412
20 15549.993314 15355.037039 0.244379 0.255026
25 15034.693551 14855.790242 0.234763 0.242837
30 14386.349902 14198.679220 0.229501 0.239303
35 13537.848480 13337.263084 0.238753 0.253282
40 12672.745873 12471.883194 0.254498 0.272017
45 11883.354388 11709.781442 0.274462 0.290670
50 11143.884144 11062.479762 0.299847 0.308937
55 10478.176725 10407.487680 0.322281 0.334168
60 9864.967730 9774.616033 0.349796 0.364399
65 9316.705476 9271.982358 0.379320 0.385156
70 9181.792936 9121.059366 0.366714 0.378376
75 9108.642539 9056.972704 0.348132 0.356722
80 9084.295103 9016.594635 0.328849 0.339635
85 9111.357500 9006.654975 0.302365 0.318891
90 9069.522485 8972.763135 0.286715 0.300559
95 9022.555097 8968.768294 0.273829 0.281449
100 9029.769617 8942.190881 0.257064 0.267309
105 8976.672363 8900.738938 0.245541 0.254386
110 9012.549872 8896.915088 0.226280 0.239292
115 8971.977997 8900.791271 0.216976 0.227092
120 8904.790154 8848.896836 0.212322 0.218636
125 8919.544981 8843.978465 0.200177 0.207149
130 8907.158270 8824.093831 0.188602 0.197145
135 8885.055501 8809.102746 0.180874 0.188480
140 8833.700643 8771.761343 0.176308 0.182293
145 8801.127863 8783.647599 0.168953 0.171703
150 8824.797408 8736.470439 0.158740 0.166432
155 8801.151948 8718.525135 0.153259 0.161367
160 8790.722508 8694.231608 0.147153 0.154423
165 8766.134921 8649.088026 0.140835 0.150506
170 8733.853914 8621.804914 0.136649 0.144998
175 8678.905305 8629.947937 0.133883 0.137221
180 8713.240812 8618.940400 0.126614 0.132390
185 8690.867087 8632.160122 0.121146 0.125620
190 8674.415597 8589.174521 0.117001 0.121948
195 8625.598670 8586.796554 0.114779 0.117890
200 8611.261354 8550.096150 0.110876 0.114517
205 8581.897782 8546.014530 0.107376 0.109796
210 8566.974887 8487.939548 0.103806 0.108179
215 8523.003297 8461.893246 0.100800 0.104419
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 24/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.52062e+42 erg / s Luminosity absorbed = 1.36911e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14470.495196 14439.383342 0.540916 0.545080
5 15344.013647 15287.824637 0.384238 0.388787
10 15650.287302 15554.430528 0.316687 0.322520
15 15649.137611 15548.412490 0.275412 0.280082
20 15355.037039 15278.361758 0.255026 0.256587
25 14855.790242 14722.099510 0.242837 0.247162
30 14198.679220 14030.505526 0.239303 0.246144
35 13337.263084 13202.971353 0.253282 0.258033
40 12471.883194 12379.235050 0.272017 0.277166
45 11709.781442 11653.492906 0.290670 0.294545
50 11062.479762 10961.781850 0.308937 0.317106
55 10407.487680 10314.959040 0.334168 0.341358
60 9774.616033 9788.937173 0.364399 0.360526
65 9271.982358 9252.044247 0.385156 0.388297
70 9121.059366 9094.933115 0.378376 0.378186
75 9056.972704 9061.701483 0.356722 0.353253
80 9016.594635 8995.518124 0.339635 0.340515
85 9006.654975 8959.328296 0.318891 0.323211
90 8972.763135 8920.517810 0.300559 0.305046
95 8968.768294 8912.435094 0.281449 0.289535
100 8942.190881 8894.124908 0.267309 0.273591
105 8900.738938 8873.382355 0.254386 0.257205
110 8896.915088 8868.773717 0.239292 0.241261
115 8900.791271 8854.993081 0.227092 0.230346
120 8848.896836 8811.078110 0.218636 0.220893
125 8843.978465 8829.644050 0.207149 0.207952
130 8824.093831 8788.689000 0.197145 0.201215
135 8809.102746 8787.612627 0.188480 0.189961
140 8771.761343 8740.167292 0.182293 0.185001
145 8783.647599 8732.804547 0.171703 0.175840
150 8736.470439 8712.980426 0.166432 0.167455
155 8718.525135 8694.556614 0.161367 0.162622
160 8694.231608 8688.460087 0.154423 0.154070
165 8649.088026 8652.726333 0.150506 0.149510
170 8621.804914 8664.425423 0.144998 0.141927
175 8629.947937 8645.616623 0.137221 0.135993
180 8618.940400 8617.566209 0.132390 0.132218
185 8632.160122 8586.371277 0.125620 0.127863
190 8589.174521 8583.876718 0.121948 0.122524
195 8586.796554 8565.133414 0.117890 0.118756
200 8550.096150 8587.038961 0.114517 0.112478
205 8546.014530 8525.571532 0.109796 0.110728
210 8487.939548 8476.450531 0.108179 0.108114
215 8461.893246 8441.936810 0.104419 0.104832
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 25/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.60156e+42 erg / s Luminosity absorbed = 1.36095e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14439.383342 14437.116780 0.545080 0.545538
5 15287.824637 15208.388456 0.388787 0.394694
10 15554.430528 15465.034944 0.322520 0.329390
15 15548.412490 15412.457674 0.280082 0.287618
20 15278.361758 15145.831208 0.256587 0.263885
25 14722.099510 14553.799657 0.247162 0.255353
30 14030.505526 13823.788145 0.246144 0.258975
35 13202.971353 12984.359020 0.258033 0.272115
40 12379.235050 12237.744937 0.277166 0.287929
45 11653.492906 11476.331616 0.294545 0.312857
50 10961.781850 10851.867539 0.317106 0.331583
55 10314.959040 10213.094685 0.341358 0.358849
60 9788.937173 9627.334087 0.360526 0.387457
65 9252.044247 9122.601885 0.388297 0.412794
70 9094.933115 8965.806200 0.378186 0.403706
75 9061.701483 8919.477778 0.353253 0.378108
80 8995.518124 8851.061232 0.340515 0.364432
85 8959.328296 8852.046845 0.323211 0.344019
90 8920.517810 8831.845485 0.305046 0.321067
95 8912.435094 8812.422594 0.289535 0.305359
100 8894.124908 8759.712819 0.273591 0.294303
105 8873.382355 8787.423778 0.257205 0.269151
110 8868.773717 8805.995679 0.241261 0.251308
115 8854.993081 8753.994314 0.230346 0.243327
120 8811.078110 8722.670406 0.220893 0.231825
125 8829.644050 8704.253189 0.207952 0.222245
130 8788.689000 8686.428399 0.201215 0.212143
135 8787.612627 8668.273534 0.189961 0.202056
140 8740.167292 8623.070495 0.185001 0.196536
145 8732.804547 8659.909430 0.175840 0.183356
150 8712.980426 8652.405947 0.167455 0.175147
155 8694.556614 8611.353196 0.162622 0.168822
160 8688.460087 8592.789426 0.154070 0.161985
165 8652.726333 8570.267904 0.149510 0.156621
170 8664.425423 8570.576999 0.141927 0.150070
175 8645.616623 8538.715386 0.135993 0.145110
180 8617.566209 8539.488067 0.132218 0.138459
185 8586.371277 8483.200715 0.127863 0.136243
190 8583.876718 8473.866237 0.122524 0.130449
195 8565.133414 8459.491337 0.118756 0.125984
200 8587.038961 8471.330226 0.112478 0.120085
205 8525.571532 8442.341282 0.110728 0.116065
210 8476.450531 8383.081663 0.108114 0.114197
215 8441.936810 8349.425816 0.104832 0.110679
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 26/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.57549e+42 erg / s Luminosity absorbed = 1.36399e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14437.116780 14447.800661 0.545538 0.543522
5 15208.388456 15187.337939 0.394694 0.396558
10 15465.034944 15372.013233 0.329390 0.332810
15 15412.457674 15255.096011 0.287618 0.296214
20 15145.831208 14975.570690 0.263885 0.271200
25 14553.799657 14386.021279 0.255353 0.261695
30 13823.788145 13656.578942 0.258975 0.266606
35 12984.359020 12866.983707 0.272115 0.279457
40 12237.744937 12177.931229 0.287929 0.292067
45 11476.331616 11429.490010 0.312857 0.316759
50 10851.867539 10845.027168 0.331583 0.331514
55 10213.094685 10235.032432 0.358849 0.356731
60 9627.334087 9669.110203 0.387457 0.383669
65 9122.601885 9097.954569 0.412794 0.419061
70 8965.806200 8968.475708 0.403706 0.404309
75 8919.477778 8930.583222 0.378108 0.376530
80 8851.061232 8883.311408 0.364432 0.361427
85 8852.046845 8853.670421 0.344019 0.341969
90 8831.845485 8838.048898 0.321067 0.321419
95 8812.422594 8827.250729 0.305359 0.304083
100 8759.712819 8801.457588 0.294303 0.288052
105 8787.423778 8767.655170 0.269151 0.270964
110 8805.995679 8762.291342 0.251308 0.256165
115 8753.994314 8770.631556 0.243327 0.240710
120 8722.670406 8740.191499 0.231825 0.230293
125 8704.253189 8715.880329 0.222245 0.221697
130 8686.428399 8680.240453 0.212143 0.212627
135 8668.273534 8682.604150 0.202056 0.200232
140 8623.070495 8640.082127 0.196536 0.195253
145 8659.909430 8650.251766 0.183356 0.183872
150 8652.405947 8649.814169 0.175147 0.175618
155 8611.353196 8619.046960 0.168822 0.168567
160 8592.789426 8610.523603 0.161985 0.162408
165 8570.267904 8559.554560 0.156621 0.157239
170 8570.576999 8558.020831 0.150070 0.151176
175 8538.715386 8513.431640 0.145110 0.146494
180 8539.488067 8545.040169 0.138459 0.138197
185 8483.200715 8491.199260 0.136243 0.135345
190 8473.866237 8474.738077 0.130449 0.129895
195 8459.491337 8430.986868 0.125984 0.127507
200 8471.330226 8430.450850 0.120085 0.122018
205 8442.341282 8419.764772 0.116065 0.117632
210 8383.081663 8366.360454 0.114197 0.114850
215 8349.425816 8331.657948 0.110679 0.111728
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 27/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.51806e+42 erg / s Luminosity absorbed = 1.36987e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14447.800661 14417.219086 0.543522 0.549349
5 15187.337939 15181.195510 0.396558 0.398061
10 15372.013233 15328.307051 0.332810 0.334343
15 15255.096011 15203.511026 0.296214 0.298500
20 14975.570690 14808.409736 0.271200 0.279158
25 14386.021279 14261.242828 0.261695 0.269496
30 13656.578942 13521.940483 0.266606 0.273082
35 12866.983707 12680.830794 0.279457 0.290405
40 12177.931229 12012.770662 0.292067 0.301843
45 11429.490010 11314.197665 0.316759 0.325619
50 10845.027168 10718.769264 0.331514 0.344833
55 10235.032432 10175.844754 0.356731 0.362672
60 9669.110203 9659.545119 0.383669 0.384407
65 9097.954569 9133.583181 0.419061 0.411875
70 8968.475708 8983.550656 0.404309 0.400423
75 8930.583222 8916.883150 0.376530 0.376356
80 8883.311408 8881.311910 0.361427 0.359073
85 8853.670421 8848.822085 0.341969 0.338215
90 8838.048898 8828.974564 0.321419 0.319146
95 8827.250729 8807.691955 0.304083 0.302713
100 8801.457588 8785.221475 0.288052 0.285162
105 8767.655170 8751.387856 0.270964 0.270968
110 8762.291342 8780.883817 0.256165 0.254049
115 8770.631556 8751.480700 0.240710 0.242601
120 8740.191499 8734.457617 0.230293 0.230096
125 8715.880329 8704.113589 0.221697 0.221491
130 8680.240453 8653.751674 0.212627 0.213909
135 8682.604150 8661.234925 0.200232 0.201672
140 8640.082127 8658.835032 0.195253 0.192969
145 8650.251766 8641.881040 0.183872 0.184367
150 8649.814169 8613.747196 0.175618 0.177634
155 8619.046960 8555.853044 0.168567 0.172421
160 8610.523603 8579.102688 0.162408 0.164319
165 8559.554560 8542.446129 0.157239 0.157480
170 8558.020831 8533.533199 0.151176 0.151698
175 8513.431640 8496.442832 0.146494 0.146762
180 8545.040169 8524.482824 0.138197 0.138349
185 8491.199260 8487.251017 0.135345 0.135365
190 8474.738077 8460.332814 0.129895 0.130966
195 8430.986868 8434.379977 0.127507 0.126942
200 8430.450850 8412.147879 0.122018 0.122168
205 8419.764772 8426.552341 0.117632 0.116774
210 8366.360454 8383.211923 0.114850 0.113791
215 8331.657948 8351.844797 0.111728 0.110146
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 28/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.63368e+42 erg / s Luminosity absorbed = 1.35831e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14417.219086 14429.531917 0.549349 0.545585
5 15181.195510 15120.350364 0.398061 0.404133
10 15328.307051 15281.732868 0.334343 0.339321
15 15203.511026 15117.848316 0.298500 0.304393
20 14808.409736 14694.996090 0.279158 0.287399
25 14261.242828 14183.269980 0.269496 0.276461
30 13521.940483 13407.675917 0.273082 0.281302
35 12680.830794 12617.861166 0.290405 0.297246
40 12012.770662 11969.509130 0.301843 0.306747
45 11314.197665 11231.733197 0.325619 0.336649
50 10718.769264 10650.678996 0.344833 0.355468
55 10175.844754 10105.622482 0.362672 0.373025
60 9659.545119 9622.146412 0.384407 0.390830
65 9133.583181 9111.276842 0.411875 0.417521
70 8983.550656 8968.479051 0.400423 0.402578
75 8916.883150 8927.469964 0.376356 0.378112
80 8881.311910 8901.508883 0.359073 0.357221
85 8848.822085 8868.433596 0.338215 0.337688
90 8828.974564 8817.757551 0.319146 0.324247
95 8807.691955 8798.092101 0.302713 0.307739
100 8785.221475 8763.887474 0.285162 0.290757
105 8751.387856 8739.862756 0.270968 0.276310
110 8780.883817 8768.800462 0.254049 0.257413
115 8751.480700 8735.955476 0.242601 0.246462
120 8734.457617 8756.530513 0.230096 0.231126
125 8704.113589 8703.467418 0.221491 0.224571
130 8653.751674 8646.364373 0.213909 0.217122
135 8661.234925 8653.142702 0.201672 0.203528
140 8658.835032 8678.613405 0.192969 0.192395
145 8641.881040 8626.451663 0.184367 0.186953
150 8613.747196 8597.106111 0.177634 0.179226
155 8555.853044 8584.136268 0.172421 0.171616
160 8579.102688 8589.260883 0.164319 0.164299
165 8542.446129 8541.224977 0.157480 0.158702
170 8533.533199 8522.338318 0.151698 0.152727
175 8496.442832 8491.754843 0.146762 0.148037
180 8524.482824 8474.412917 0.138349 0.141995
185 8487.251017 8460.237138 0.135365 0.137380
190 8460.332814 8415.212522 0.130966 0.134130
195 8434.379977 8399.982738 0.126942 0.129201
200 8412.147879 8372.301416 0.122168 0.125072
205 8426.552341 8361.337581 0.116774 0.120960
210 8383.211923 8326.692458 0.113791 0.117308
215 8351.844797 8289.240860 0.110146 0.114643
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 29/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.63865e+42 erg / s Luminosity absorbed = 1.35826e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Plasma stratification:
t_rad next_t_rad w next_w
Shell
0 14429.531917 14393.560704 0.545585 0.549791
5 15120.350364 15047.460876 0.404133 0.407487
10 15281.732868 15152.376157 0.339321 0.348021
15 15117.848316 15011.073356 0.304393 0.311534
20 14694.996090 14594.083857 0.287399 0.293825
25 14183.269980 14050.725200 0.276461 0.284887
30 13407.675917 13300.100420 0.281302 0.290669
35 12617.861166 12506.921436 0.297246 0.306635
40 11969.509130 11810.035443 0.306747 0.325268
45 11231.733197 11126.963002 0.336649 0.347387
50 10650.678996 10537.289550 0.355468 0.369531
55 10105.622482 10022.685606 0.373025 0.383605
60 9622.146412 9557.446390 0.390830 0.401964
65 9111.276842 9056.805821 0.417521 0.429370
70 8968.479051 8921.840347 0.402578 0.411349
75 8927.469964 8907.031457 0.378112 0.384014
80 8901.508883 8831.590758 0.357221 0.371482
85 8868.433596 8815.898171 0.337688 0.349352
90 8817.757551 8777.689182 0.324247 0.332044
95 8798.092101 8777.981323 0.307739 0.310802
100 8763.887474 8735.735990 0.290757 0.294033
105 8739.862756 8736.781935 0.276310 0.276744
110 8768.800462 8703.852852 0.257413 0.266497
115 8735.955476 8704.749811 0.246462 0.250630
120 8756.530513 8682.672809 0.231126 0.238432
125 8703.467418 8617.135647 0.224571 0.234778
130 8646.364373 8602.578401 0.217122 0.222175
135 8653.142702 8626.725758 0.203528 0.207815
140 8678.613405 8598.950945 0.192395 0.202267
145 8626.451663 8594.447574 0.186953 0.190719
150 8597.106111 8561.702948 0.179226 0.184284
155 8584.136268 8535.674486 0.171616 0.176376
160 8589.260883 8539.167829 0.164299 0.169219
165 8541.224977 8531.513673 0.158702 0.160695
170 8522.338318 8484.097260 0.152727 0.158442
175 8491.754843 8474.481043 0.148037 0.150545
180 8474.412917 8438.276493 0.141995 0.145788
185 8460.237138 8435.079927 0.137380 0.140271
190 8415.212522 8395.762702 0.134130 0.137036
195 8399.982738 8366.507497 0.129201 0.132203
200 8372.301416 8351.094536 0.125072 0.127218
205 8361.337581 8341.067731 0.120960 0.122846
210 8326.692458 8311.116368 0.117308 0.118902
215 8289.240860 8279.957190 0.114643 0.115706
(base.py:348)
[tardis.simulation.base][INFO ] t_inner 13000.000 K -- next t_inner 13000.000 K (base.py:350)
[tardis.simulation.base][INFO ] Starting iteration 30/30 (base.py:266)
[tardis.simulation.base][INFO ] Luminosity emitted = 6.60081e+42 erg / s Luminosity absorbed = 1.36197e+43 erg / s Luminosity requested = 1.05928e+43 erg / s (base.py:357)
[tardis.simulation.base][INFO ] Simulation finished in 30 iterations and took 337.79 s (base.py:306)
[16]:
spectrum = pd.DataFrame(data = sim.runner.spectrum_integrated.frequency.value, columns=['frequency'])
spectrum['luminosity'] = sim.runner.spectrum_integrated.luminosity_density_nu.value
create_df_download_link(spectrum, title='Download Spectrum', filename='spectrum.csv', type='csv')
[16]:
[23]:
sim.plasma.ion_number_density / sim.plasma.number_density
[23]:
| 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | ... | 210 | 211 | 212 | 213 | 214 | 215 | 216 | 217 | 218 | 219 | ||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| atomic_number | ion_number | |||||||||||||||||||||
| 14 | 0 | 3.371573e-15 | 2.836444e-15 | 2.419784e-15 | 2.136677e-15 | 1.999989e-15 | 2.054821e-15 | 2.110358e-15 | 2.153097e-15 | 2.178870e-15 | 2.396495e-15 | ... | 3.099967e-09 | 3.105191e-09 | 3.119057e-09 | 3.211761e-09 | 3.355139e-09 | 3.226573e-09 | 3.442658e-09 | 3.458125e-09 | 3.453447e-09 | 3.465355e-09 |
| 1 | 2.022453e-06 | 1.705701e-06 | 1.462930e-06 | 1.291823e-06 | 1.208121e-06 | 1.225786e-06 | 1.248424e-06 | 1.263770e-06 | 1.269590e-06 | 1.375301e-06 | ... | 7.935683e-02 | 7.986164e-02 | 8.051595e-02 | 8.248215e-02 | 8.528176e-02 | 8.360084e-02 | 8.762037e-02 | 8.833458e-02 | 8.872181e-02 | 8.938419e-02 | |
| 2 | 3.887190e-01 | 3.393144e-01 | 3.008827e-01 | 2.707304e-01 | 2.546273e-01 | 2.514531e-01 | 2.505741e-01 | 2.486651e-01 | 2.454469e-01 | 2.545289e-01 | ... | 9.206432e-01 | 9.201384e-01 | 9.194840e-01 | 9.175178e-01 | 9.147182e-01 | 9.163992e-01 | 9.123796e-01 | 9.116654e-01 | 9.112782e-01 | 9.106158e-01 | |
| 3 | 6.112531e-01 | 6.606446e-01 | 6.990609e-01 | 7.291956e-01 | 7.452862e-01 | 7.484558e-01 | 7.493327e-01 | 7.512385e-01 | 7.544522e-01 | 7.453764e-01 | ... | 3.653376e-09 | 3.585998e-09 | 3.481458e-09 | 3.266309e-09 | 2.989207e-09 | 3.137543e-09 | 2.781066e-09 | 2.722738e-09 | 2.685372e-09 | 2.637792e-09 | |
| 4 | 2.587464e-05 | 3.930435e-05 | 5.492756e-05 | 7.275498e-05 | 8.533674e-05 | 8.983376e-05 | 9.201249e-05 | 9.513098e-05 | 9.956385e-05 | 9.331045e-05 | ... | 3.671725e-25 | 3.512867e-25 | 3.268472e-25 | 2.814316e-25 | 2.286085e-25 | 2.557571e-25 | 1.930661e-25 | 1.837662e-25 | 1.778037e-25 | 1.706548e-25 | |
| 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 | |
| 16 | 0 | 2.016802e-12 | 1.658852e-12 | 1.393517e-12 | 1.205457e-12 | 1.110193e-12 | 1.112130e-12 | 1.121180e-12 | 1.122548e-12 | 1.114228e-12 | 1.198344e-12 | ... | 1.249653e-06 | 1.246017e-06 | 1.245134e-06 | 1.258252e-06 | 1.280630e-06 | 1.251831e-06 | 1.286763e-06 | 1.284397e-06 | 1.278783e-06 | 1.275507e-06 |
| 1 | 2.394294e-04 | 2.040317e-04 | 1.771086e-04 | 1.569568e-04 | 1.464626e-04 | 1.456358e-04 | 1.459572e-04 | 1.454903e-04 | 1.440266e-04 | 1.517962e-04 | ... | 9.957734e-01 | 9.958178e-01 | 9.958872e-01 | 9.960498e-01 | 9.962669e-01 | 9.961392e-01 | 9.964290e-01 | 9.964735e-01 | 9.964994e-01 | 9.965363e-01 | |
| 2 | 8.434871e-01 | 8.111234e-01 | 7.807144e-01 | 7.526197e-01 | 7.357674e-01 | 7.320337e-01 | 7.308866e-01 | 7.286039e-01 | 7.248339e-01 | 7.346920e-01 | ... | 4.225326e-03 | 4.180994e-03 | 4.111542e-03 | 3.948930e-03 | 3.731861e-03 | 3.859506e-03 | 3.569744e-03 | 3.525171e-03 | 3.499367e-03 | 3.462394e-03 | |
| 3 | 1.562730e-01 | 1.886718e-01 | 2.191073e-01 | 2.472216e-01 | 2.640839e-01 | 2.678184e-01 | 2.689651e-01 | 2.712481e-01 | 2.750194e-01 | 2.651538e-01 | ... | 8.812221e-13 | 8.555966e-13 | 8.161771e-13 | 7.350285e-13 | 6.352517e-13 | 6.895675e-13 | 5.650228e-13 | 5.461861e-13 | 5.345989e-13 | 5.195702e-13 | |
| 4 | 4.416652e-07 | 7.660508e-07 | 1.195779e-06 | 1.740700e-06 | 2.154429e-06 | 2.302852e-06 | 2.373892e-06 | 2.478209e-06 | 2.630008e-06 | 2.404086e-06 | ... | 0.000000e+00 | 0.000000e+00 | 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 | 1.740392e-20 | 5.708941e-20 | 1.493104e-19 | 3.421705e-19 | 5.532661e-19 | 6.679461e-19 | 7.371750e-19 | 8.355242e-19 | 9.808088e-19 | 8.453664e-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 | |
| 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 | |
| ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
| 27 | 27 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN |
| 28 | 0 | 3.867508e-15 | 3.345109e-15 | 2.918472e-15 | 2.607352e-15 | 2.446026e-15 | 2.476290e-15 | 2.512104e-15 | 2.533291e-15 | 2.535729e-15 | 2.725434e-15 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN |
| 1 | 4.768144e-06 | 4.199161e-06 | 3.730841e-06 | 3.373334e-06 | 3.185213e-06 | 3.198249e-06 | 3.225992e-06 | 3.237521e-06 | 3.228287e-06 | 3.416276e-06 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 2 | 7.608692e-01 | 7.150802e-01 | 6.735910e-01 | 6.363741e-01 | 6.145266e-01 | 6.092583e-01 | 6.073952e-01 | 6.041015e-01 | 5.989427e-01 | 6.108105e-01 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 3 | 2.391261e-01 | 2.849156e-01 | 3.264053e-01 | 3.636225e-01 | 3.854702e-01 | 3.907384e-01 | 3.926015e-01 | 3.958952e-01 | 4.010540e-01 | 3.891860e-01 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 4 | 5.988713e-09 | 1.116299e-08 | 1.840678e-08 | 2.814457e-08 | 3.588567e-08 | 3.917661e-08 | 4.094012e-08 | 4.337653e-08 | 4.679478e-08 | 4.298233e-08 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 5 | 1.638349e-24 | 6.005792e-24 | 1.712472e-23 | 4.242307e-23 | 7.190971e-23 | 8.956895e-23 | 1.008317e-22 | 1.168034e-22 | 1.405184e-22 | 1.216072e-22 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 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 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 21 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 22 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 23 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 24 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 25 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 26 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 27 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | |
| 28 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | 0.000000e+00 | ... | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN | NaN |
137 rows × 220 columns
[ ]: