You can interact with this notebook online: Launch interactive version
Quickstart for TARDIS¶
Every simulation run requires the atomic database (for more info refer to atomic data) and a configuration file (more info at configuration). You can obtain a copy of the atomic database from the (https://github.com/tardis-sn/tardis-refdata) repository (atom_data subfolder). We recommended to use the kurucz_cd23_chianti_H_He.h5 dataset (which is auto-downloaded in the second cell already). The
configuration file for this quickstart is tardis_example.yml, which can be downloaded here), though this file is auto-downloaded in the third cell.
After the installation, start a Jupyter server executing jupyter notebook on the command line in a directory that contains this quickstart.
[1]:
from tardis import run_tardis
from tardis.io.atom_data.util import download_atom_data
/usr/share/miniconda3/envs/tardis/lib/python3.8/site-packages/setuptools_scm/git.py:105: UserWarning: "/home/runner/work/tardis/tardis" is shallow and may cause errors
warnings.warn(f'"{wd.path}" is shallow and may cause errors')
Downloading the atomic data¶
[2]:
# The data are automatically downloaded
download_atom_data('kurucz_cd23_chianti_H_He')
Configuration File /home/runner/.astropy/config/tardis_internal_config.yml does not exist - creating new one from default
CRITICAL:root:
********************************************************************************
TARDIS will download different kinds of data (e.g. atomic) to its data directory /home/runner/Downloads/tardis-data
TARDIS DATA DIRECTORY not specified in /home/runner/.astropy/config/tardis_internal_config.yml:
ASSUMING DEFAULT DATA DIRECTORY /home/runner/Downloads/tardis-data
YOU CAN CHANGE THIS AT ANY TIME IN /home/runner/.astropy/config/tardis_internal_config.yml
********************************************************************************
Downloading the example file¶
[3]:
!curl -O https://raw.githubusercontent.com/tardis-sn/tardis/master/docs/tardis_example.yml
% Total % Received % Xferd Average Speed Time Time Time Current
Dload Upload Total Spent Left Speed
100 980 100 980 0 0 7716 0 --:--:-- --:--:-- --:--:-- 7716
Running the simulation (long output)¶
Note
The progress of the simulation can be tracked using progress bars which are displayed when the notebook is run, but are not displayed in the documentation.
[4]:
sim = run_tardis('tardis_example.yml')
[tardis.plasma.standard_plasmas][INFO ]
Reading Atomic Data from kurucz_cd23_chianti_H_He.h5 (standard_plasmas.py:93)
INFO:tardis.plasma.standard_plasmas:
Reading Atomic Data from kurucz_cd23_chianti_H_He.h5
[tardis.io.atom_data.util][INFO ]
Atom Data kurucz_cd23_chianti_H_He.h5 not found in local path.
Exists in TARDIS Data repo /home/runner/Downloads/tardis-data/kurucz_cd23_chianti_H_He.h5 (util.py:33)
INFO:tardis.io.atom_data.util:
Atom Data kurucz_cd23_chianti_H_He.h5 not found in local path.
Exists in TARDIS Data repo /home/runner/Downloads/tardis-data/kurucz_cd23_chianti_H_He.h5
[tardis.io.atom_data.base][INFO ]
Reading Atom Data with: UUID = 6f7b09e887a311e7a06b246e96350010 MD5 = 864f1753714343c41f99cb065710cace (base.py:249)
INFO:tardis.io.atom_data.base:Reading Atom Data with: UUID = 6f7b09e887a311e7a06b246e96350010 MD5 = 864f1753714343c41f99cb065710cace
[tardis.io.atom_data.base][INFO ]
Non provided Atomic Data: synpp_refs, photoionization_data, yg_data, two_photon_data (base.py:253)
INFO:tardis.io.atom_data.base:Non provided Atomic Data: synpp_refs, photoionization_data, yg_data, two_photon_data
[py.warnings ][WARNING]
/home/runner/work/tardis/tardis/tardis/plasma/properties/radiative_properties.py:93: RuntimeWarning: invalid value encountered in true_divide
(g_lower * n_upper) / (g_upper * n_lower)
(warnings.py:109)
WARNING:py.warnings:/home/runner/work/tardis/tardis/tardis/plasma/properties/radiative_properties.py:93: RuntimeWarning: invalid value encountered in true_divide
(g_lower * n_upper) / (g_upper * n_lower)
[py.warnings ][WARNING]
/home/runner/work/tardis/tardis/tardis/plasma/properties/radiative_properties.py:93: RuntimeWarning: invalid value encountered in true_divide
(g_lower * n_upper) / (g_upper * n_lower)
(warnings.py:109)
WARNING:py.warnings:/home/runner/work/tardis/tardis/tardis/plasma/properties/radiative_properties.py:93: RuntimeWarning: invalid value encountered in true_divide
(g_lower * n_upper) / (g_upper * n_lower)
[tardis.simulation.base][INFO ]
Starting iteration 1 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 1 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 7.942e+42 erg / s
Luminosity absorbed = 2.659e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 7.942e+42 erg / s
Luminosity absorbed = 2.659e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 9.93e+03 | 1.01e+04 | 0.4 | 0.507 |
| 5 | 9.85e+03 | 1.02e+04 | 0.211 | 0.197 |
| 10 | 9.78e+03 | 1.01e+04 | 0.143 | 0.117 |
| 15 | 9.71e+03 | 9.87e+03 | 0.105 | 0.0869 |
[tardis.simulation.base][INFO ]
Current t_inner = 9933.952 K
Expected t_inner for next iteration = 10703.212 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 9933.952 K
Expected t_inner for next iteration = 10703.212 K
[py.warnings ][WARNING]
/home/runner/work/tardis/tardis/tardis/plasma/properties/radiative_properties.py:93: RuntimeWarning: invalid value encountered in true_divide
(g_lower * n_upper) / (g_upper * n_lower)
(warnings.py:109)
WARNING:py.warnings:/home/runner/work/tardis/tardis/tardis/plasma/properties/radiative_properties.py:93: RuntimeWarning: invalid value encountered in true_divide
(g_lower * n_upper) / (g_upper * n_lower)
[tardis.simulation.base][INFO ]
Starting iteration 2 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 2 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.071e+43 erg / s
Luminosity absorbed = 3.576e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.071e+43 erg / s
Luminosity absorbed = 3.576e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.01e+04 | 1.08e+04 | 0.507 | 0.525 |
| 5 | 1.02e+04 | 1.1e+04 | 0.197 | 0.203 |
| 10 | 1.01e+04 | 1.08e+04 | 0.117 | 0.125 |
| 15 | 9.87e+03 | 1.05e+04 | 0.0869 | 0.0933 |
[tardis.simulation.base][INFO ]
Current t_inner = 10703.212 K
Expected t_inner for next iteration = 10673.712 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10703.212 K
Expected t_inner for next iteration = 10673.712 K
[tardis.simulation.base][INFO ]
Starting iteration 3 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 3 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.074e+43 erg / s
Luminosity absorbed = 3.391e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.074e+43 erg / s
Luminosity absorbed = 3.391e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.08e+04 | 1.1e+04 | 0.525 | 0.483 |
| 5 | 1.1e+04 | 1.12e+04 | 0.203 | 0.189 |
| 10 | 1.08e+04 | 1.1e+04 | 0.125 | 0.118 |
| 15 | 1.05e+04 | 1.06e+04 | 0.0933 | 0.0895 |
[tardis.simulation.base][INFO ]
Current t_inner = 10673.712 K
Expected t_inner for next iteration = 10635.953 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10673.712 K
Expected t_inner for next iteration = 10635.953 K
[tardis.simulation.base][INFO ]
Starting iteration 4 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 4 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.058e+43 erg / s
Luminosity absorbed = 3.352e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.058e+43 erg / s
Luminosity absorbed = 3.352e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 | 1.1e+04 | 0.483 | 0.469 |
| 5 | 1.12e+04 | 1.12e+04 | 0.189 | 0.182 |
| 10 | 1.1e+04 | 1.1e+04 | 0.118 | 0.113 |
| 15 | 1.06e+04 | 1.07e+04 | 0.0895 | 0.0861 |
[tardis.simulation.base][INFO ]
Current t_inner = 10635.953 K
Expected t_inner for next iteration = 10638.407 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10635.953 K
Expected t_inner for next iteration = 10638.407 K
[tardis.simulation.base][INFO ]
Starting iteration 5 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 5 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.055e+43 erg / s
Luminosity absorbed = 3.399e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.055e+43 erg / s
Luminosity absorbed = 3.399e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Iteration converged 1/4 consecutive times. (base.py:259)
INFO:tardis.simulation.base:Iteration converged 1/4 consecutive times.
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 | 1.1e+04 | 0.469 | 0.479 |
| 5 | 1.12e+04 | 1.13e+04 | 0.182 | 0.178 |
| 10 | 1.1e+04 | 1.1e+04 | 0.113 | 0.113 |
| 15 | 1.07e+04 | 1.07e+04 | 0.0861 | 0.0839 |
[tardis.simulation.base][INFO ]
Current t_inner = 10638.407 K
Expected t_inner for next iteration = 10650.202 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10638.407 K
Expected t_inner for next iteration = 10650.202 K
[tardis.simulation.base][INFO ]
Starting iteration 6 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 6 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.061e+43 erg / s
Luminosity absorbed = 3.398e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.061e+43 erg / s
Luminosity absorbed = 3.398e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Iteration converged 2/4 consecutive times. (base.py:259)
INFO:tardis.simulation.base:Iteration converged 2/4 consecutive times.
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 | 1.1e+04 | 0.479 | 0.47 |
| 5 | 1.13e+04 | 1.12e+04 | 0.178 | 0.185 |
| 10 | 1.1e+04 | 1.11e+04 | 0.113 | 0.112 |
| 15 | 1.07e+04 | 1.07e+04 | 0.0839 | 0.0856 |
[tardis.simulation.base][INFO ]
Current t_inner = 10650.202 K
Expected t_inner for next iteration = 10645.955 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10650.202 K
Expected t_inner for next iteration = 10645.955 K
[tardis.simulation.base][INFO ]
Starting iteration 7 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 7 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.061e+43 erg / s
Luminosity absorbed = 3.382e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.061e+43 erg / s
Luminosity absorbed = 3.382e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Iteration converged 3/4 consecutive times. (base.py:259)
INFO:tardis.simulation.base:Iteration converged 3/4 consecutive times.
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 | 1.1e+04 | 0.47 | 0.47 |
| 5 | 1.12e+04 | 1.13e+04 | 0.185 | 0.178 |
| 10 | 1.11e+04 | 1.11e+04 | 0.112 | 0.112 |
| 15 | 1.07e+04 | 1.07e+04 | 0.0856 | 0.086 |
[tardis.simulation.base][INFO ]
Current t_inner = 10645.955 K
Expected t_inner for next iteration = 10642.050 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10645.955 K
Expected t_inner for next iteration = 10642.050 K
[tardis.simulation.base][INFO ]
Starting iteration 8 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 8 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.062e+43 erg / s
Luminosity absorbed = 3.350e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.062e+43 erg / s
Luminosity absorbed = 3.350e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Iteration converged 4/4 consecutive times. (base.py:259)
INFO:tardis.simulation.base:Iteration converged 4/4 consecutive times.
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 | 1.11e+04 | 0.47 | 0.472 |
| 5 | 1.13e+04 | 1.14e+04 | 0.178 | 0.175 |
| 10 | 1.11e+04 | 1.11e+04 | 0.112 | 0.111 |
| 15 | 1.07e+04 | 1.07e+04 | 0.086 | 0.084 |
[tardis.simulation.base][INFO ]
Current t_inner = 10642.050 K
Expected t_inner for next iteration = 10636.106 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10642.050 K
Expected t_inner for next iteration = 10636.106 K
[tardis.simulation.base][INFO ]
Starting iteration 9 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 9 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.052e+43 erg / s
Luminosity absorbed = 3.411e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.052e+43 erg / s
Luminosity absorbed = 3.411e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Iteration converged 5/4 consecutive times. (base.py:259)
INFO:tardis.simulation.base:Iteration converged 5/4 consecutive times.
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 | 1.11e+04 | 0.472 | 0.469 |
| 5 | 1.14e+04 | 1.15e+04 | 0.175 | 0.17 |
| 10 | 1.11e+04 | 1.11e+04 | 0.111 | 0.109 |
| 15 | 1.07e+04 | 1.08e+04 | 0.084 | 0.0822 |
[tardis.simulation.base][INFO ]
Current t_inner = 10636.106 K
Expected t_inner for next iteration = 10654.313 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10636.106 K
Expected t_inner for next iteration = 10654.313 K
[tardis.simulation.base][INFO ]
Starting iteration 10 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 10 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.070e+43 erg / s
Luminosity absorbed = 3.335e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.070e+43 erg / s
Luminosity absorbed = 3.335e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 | 1.1e+04 | 0.469 | 0.475 |
| 5 | 1.15e+04 | 1.14e+04 | 0.17 | 0.177 |
| 10 | 1.11e+04 | 1.11e+04 | 0.109 | 0.112 |
| 15 | 1.08e+04 | 1.06e+04 | 0.0822 | 0.0878 |
[tardis.simulation.base][INFO ]
Current t_inner = 10654.313 K
Expected t_inner for next iteration = 10628.190 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10654.313 K
Expected t_inner for next iteration = 10628.190 K
[tardis.simulation.base][INFO ]
Starting iteration 11 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 11 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.053e+43 erg / s
Luminosity absorbed = 3.363e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.053e+43 erg / s
Luminosity absorbed = 3.363e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Iteration converged 1/4 consecutive times. (base.py:259)
INFO:tardis.simulation.base:Iteration converged 1/4 consecutive times.
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 | 1.1e+04 | 0.475 | 0.472 |
| 5 | 1.14e+04 | 1.12e+04 | 0.177 | 0.184 |
| 10 | 1.11e+04 | 1.1e+04 | 0.112 | 0.114 |
| 15 | 1.06e+04 | 1.06e+04 | 0.0878 | 0.0859 |
[tardis.simulation.base][INFO ]
Current t_inner = 10628.190 K
Expected t_inner for next iteration = 10644.054 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10628.190 K
Expected t_inner for next iteration = 10644.054 K
[tardis.simulation.base][INFO ]
Starting iteration 12 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 12 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.056e+43 erg / s
Luminosity absorbed = 3.420e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.056e+43 erg / s
Luminosity absorbed = 3.420e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 | 1.11e+04 | 0.472 | 0.467 |
| 5 | 1.12e+04 | 1.13e+04 | 0.184 | 0.176 |
| 10 | 1.1e+04 | 1.11e+04 | 0.114 | 0.11 |
| 15 | 1.06e+04 | 1.08e+04 | 0.0859 | 0.0821 |
[tardis.simulation.base][INFO ]
Current t_inner = 10644.054 K
Expected t_inner for next iteration = 10653.543 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10644.054 K
Expected t_inner for next iteration = 10653.543 K
[tardis.simulation.base][INFO ]
Starting iteration 13 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 13 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.062e+43 erg / s
Luminosity absorbed = 3.406e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.062e+43 erg / s
Luminosity absorbed = 3.406e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Iteration converged 1/4 consecutive times. (base.py:259)
INFO:tardis.simulation.base:Iteration converged 1/4 consecutive times.
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 | 1.11e+04 | 0.467 | 0.466 |
| 5 | 1.13e+04 | 1.13e+04 | 0.176 | 0.18 |
| 10 | 1.11e+04 | 1.11e+04 | 0.11 | 0.111 |
| 15 | 1.08e+04 | 1.08e+04 | 0.0821 | 0.0841 |
[tardis.simulation.base][INFO ]
Current t_inner = 10653.543 K
Expected t_inner for next iteration = 10647.277 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10653.543 K
Expected t_inner for next iteration = 10647.277 K
[tardis.simulation.base][INFO ]
Starting iteration 14 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 14 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.063e+43 erg / s
Luminosity absorbed = 3.369e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.063e+43 erg / s
Luminosity absorbed = 3.369e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Iteration converged 2/4 consecutive times. (base.py:259)
INFO:tardis.simulation.base:Iteration converged 2/4 consecutive times.
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 | 1.11e+04 | 0.466 | 0.469 |
| 5 | 1.13e+04 | 1.13e+04 | 0.18 | 0.182 |
| 10 | 1.11e+04 | 1.1e+04 | 0.111 | 0.113 |
| 15 | 1.08e+04 | 1.07e+04 | 0.0841 | 0.0854 |
[tardis.simulation.base][INFO ]
Current t_inner = 10647.277 K
Expected t_inner for next iteration = 10638.875 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10647.277 K
Expected t_inner for next iteration = 10638.875 K
[tardis.simulation.base][INFO ]
Starting iteration 15 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 15 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.053e+43 erg / s
Luminosity absorbed = 3.417e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.053e+43 erg / s
Luminosity absorbed = 3.417e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Iteration converged 3/4 consecutive times. (base.py:259)
INFO:tardis.simulation.base:Iteration converged 3/4 consecutive times.
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 | 1.1e+04 | 0.469 | 0.484 |
| 5 | 1.13e+04 | 1.13e+04 | 0.182 | 0.181 |
| 10 | 1.1e+04 | 1.1e+04 | 0.113 | 0.113 |
| 15 | 1.07e+04 | 1.07e+04 | 0.0854 | 0.0858 |
[tardis.simulation.base][INFO ]
Current t_inner = 10638.875 K
Expected t_inner for next iteration = 10655.125 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10638.875 K
Expected t_inner for next iteration = 10655.125 K
[tardis.simulation.base][INFO ]
Starting iteration 16 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 16 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.059e+43 erg / s
Luminosity absorbed = 3.445e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.059e+43 erg / s
Luminosity absorbed = 3.445e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Iteration converged 4/4 consecutive times. (base.py:259)
INFO:tardis.simulation.base:Iteration converged 4/4 consecutive times.
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 | 1.1e+04 | 0.484 | 0.472 |
| 5 | 1.13e+04 | 1.13e+04 | 0.181 | 0.177 |
| 10 | 1.1e+04 | 1.1e+04 | 0.113 | 0.113 |
| 15 | 1.07e+04 | 1.06e+04 | 0.0858 | 0.0858 |
[tardis.simulation.base][INFO ]
Current t_inner = 10655.125 K
Expected t_inner for next iteration = 10655.561 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10655.125 K
Expected t_inner for next iteration = 10655.561 K
[tardis.simulation.base][INFO ]
Starting iteration 17 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 17 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.067e+43 erg / s
Luminosity absorbed = 3.372e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.067e+43 erg / s
Luminosity absorbed = 3.372e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 | 1.11e+04 | 0.472 | 0.468 |
| 5 | 1.13e+04 | 1.14e+04 | 0.177 | 0.175 |
| 10 | 1.1e+04 | 1.11e+04 | 0.113 | 0.11 |
| 15 | 1.06e+04 | 1.08e+04 | 0.0858 | 0.0816 |
[tardis.simulation.base][INFO ]
Current t_inner = 10655.561 K
Expected t_inner for next iteration = 10636.536 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10655.561 K
Expected t_inner for next iteration = 10636.536 K
[tardis.simulation.base][INFO ]
Starting iteration 18 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 18 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.057e+43 erg / s
Luminosity absorbed = 3.365e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.057e+43 erg / s
Luminosity absorbed = 3.365e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Iteration converged 1/4 consecutive times. (base.py:259)
INFO:tardis.simulation.base:Iteration converged 1/4 consecutive times.
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 | 1.11e+04 | 0.468 | 0.464 |
| 5 | 1.14e+04 | 1.13e+04 | 0.175 | 0.177 |
| 10 | 1.11e+04 | 1.1e+04 | 0.11 | 0.113 |
| 15 | 1.08e+04 | 1.07e+04 | 0.0816 | 0.0848 |
[tardis.simulation.base][INFO ]
Current t_inner = 10636.536 K
Expected t_inner for next iteration = 10641.692 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10636.536 K
Expected t_inner for next iteration = 10641.692 K
[tardis.simulation.base][INFO ]
Starting iteration 19 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 19 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.056e+43 erg / s
Luminosity absorbed = 3.405e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.056e+43 erg / s
Luminosity absorbed = 3.405e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Iteration converged 2/4 consecutive times. (base.py:259)
INFO:tardis.simulation.base:Iteration converged 2/4 consecutive times.
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:515)
INFO:tardis.simulation.base:
Plasma stratification:
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 | 1.11e+04 | 0.464 | 0.466 |
| 5 | 1.13e+04 | 1.13e+04 | 0.177 | 0.177 |
| 10 | 1.1e+04 | 1.11e+04 | 0.113 | 0.111 |
| 15 | 1.07e+04 | 1.07e+04 | 0.0848 | 0.0853 |
[tardis.simulation.base][INFO ]
Current t_inner = 10641.692 K
Expected t_inner for next iteration = 10650.463 K
(base.py:542)
INFO:tardis.simulation.base:
Current t_inner = 10641.692 K
Expected t_inner for next iteration = 10650.463 K
[tardis.simulation.base][INFO ]
Starting iteration 20 of 20 (base.py:378)
INFO:tardis.simulation.base:
Starting iteration 20 of 20
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.061e+43 erg / s
Luminosity absorbed = 3.401e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:547)
INFO:tardis.simulation.base:
Luminosity emitted = 1.061e+43 erg / s
Luminosity absorbed = 3.401e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
[tardis.simulation.base][INFO ]
Simulation finished in 20 iterations
Simulation took 52.19 s
(base.py:468)
INFO:tardis.simulation.base:
Simulation finished in 20 iterations
Simulation took 52.19 s
Plotting the Spectrum¶
[5]:
%pylab inline
spectrum = sim.runner.spectrum
spectrum_virtual = sim.runner.spectrum_virtual
spectrum_integrated = sim.runner.spectrum_integrated
figure(figsize=(10,6))
spectrum.plot(label='normal packets')
spectrum_virtual.plot(label='virtual packets')
spectrum_integrated.plot(label='formal integral')
xlabel('Wavelength [$\AA$]')
ylabel('Luminosity Density [erg/s/$\AA$]')
legend()
xlim(500, 9000)
[py.warnings ][WARNING]
/usr/share/miniconda3/envs/tardis/lib/python3.8/site-packages/traitlets/traitlets.py:3258: FutureWarning: --rc={'figure.dpi': 96} for dict-traits is deprecated in traitlets 5.0. You can pass --rc <key=value> ... multiple times to add items to a dict.
warn(
(warnings.py:109)
WARNING:py.warnings:/usr/share/miniconda3/envs/tardis/lib/python3.8/site-packages/traitlets/traitlets.py:3258: FutureWarning: --rc={'figure.dpi': 96} for dict-traits is deprecated in traitlets 5.0. You can pass --rc <key=value> ... multiple times to add items to a dict.
warn(
%pylab is deprecated, use %matplotlib inline and import the required libraries.
Populating the interactive namespace from numpy and matplotlib
[py.warnings ][WARNING]
/home/runner/work/tardis/tardis/tardis/montecarlo/montecarlo_numba/formal_integral.py:356: UserWarning: The number of interpolate_shells was not specified. The value was set to 80.
warnings.warn(
(warnings.py:109)
WARNING:py.warnings:/home/runner/work/tardis/tardis/tardis/montecarlo/montecarlo_numba/formal_integral.py:356: UserWarning: The number of interpolate_shells was not specified. The value was set to 80.
warnings.warn(
[5]:
(500.0, 9000.0)
