You can interact with this notebook online: Launch notebook
Quickstart for TARDIS¶
Every simulation run requires atomic data and a configuration file.
Atomic Data¶
We recommend using the kurucz_cd23_chianti_H_He.h5 dataset.
[1]:
from tardis.io.atom_data import download_atom_data
[2]:
download_atom_data('kurucz_cd23_chianti_H_He')
Atomic Data kurucz_cd23_chianti_H_He already exists in /home/runner/Downloads/tardis-data/kurucz_cd23_chianti_H_He.h5. Will not download - override with force_download=True.
You can also obtain a copy of the atomic data from the tardis-regression-data repository.
Example Configuration File¶
The configuration file tardis_example.yml is used throughout this Quickstart.
[3]:
!wget -q -nc https://raw.githubusercontent.com/tardis-sn/tardis/master/docs/tardis_example.yml
[4]:
!cat tardis_example.yml
# Example YAML configuration for TARDIS
tardis_config_version: v1.0
supernova:
luminosity_requested: 9.44 log_lsun
time_explosion: 13 day
atom_data: kurucz_cd23_chianti_H_He.h5
model:
structure:
type: specific
velocity:
start: 1.1e4 km/s
stop: 20000 km/s
num: 20
density:
type: branch85_w7
abundances:
type: uniform
O: 0.19
Mg: 0.03
Si: 0.52
S: 0.19
Ar: 0.04
Ca: 0.03
plasma:
disable_electron_scattering: no
ionization: lte
excitation: lte
radiative_rates_type: dilute-blackbody
line_interaction_type: macroatom
montecarlo:
seed: 23111963
no_of_packets: 4.0e+4
iterations: 20
nthreads: 1
last_no_of_packets: 1.e+5
no_of_virtual_packets: 10
convergence_strategy:
type: damped
damping_constant: 1.0
threshold: 0.05
fraction: 0.8
hold_iterations: 3
t_inner:
damping_constant: 0.5
spectrum:
start: 500 angstrom
stop: 20000 angstrom
num: 10000
Running the Simulation¶
To run the simulation, import the run_tardis function and create the sim object.
[5]:
from tardis import run_tardis
Note:
Get more information about the progress bars, logging configuration, and convergence plots.
[6]:
sim = run_tardis("tardis_example.yml",
virtual_packet_logging=True,
show_convergence_plots=True,
export_convergence_plots=True,
log_level="INFO")
[tardis.io.model.parse_atom_data][INFO ]
Reading Atomic Data from kurucz_cd23_chianti_H_He.h5 (parse_atom_data.py:40)
[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:34)
[tardis.io.atom_data.base][INFO ]
Reading Atom Data with: UUID = 6f7b09e887a311e7a06b246e96350010 MD5 = 864f1753714343c41f99cb065710cace (base.py:258)
[tardis.io.atom_data.base][INFO ]
Non provided Atomic Data: synpp_refs, photoionization_data, yg_data, two_photon_data, linelist_atoms, linelist_molecules (base.py:262)
[tardis.io.model.parse_density_configuration][WARNING]
Number of density points larger than number of shells. Assuming inner point irrelevant (parse_density_configuration.py:114)
[tardis.model.matter.decay][INFO ]
Decaying abundances for 1123200.0 seconds (decay.py:101)
[tardis.simulation.base][INFO ]
Starting iteration 1 of 20 (base.py:444)
[py.warnings ][WARNING]
/home/runner/work/tardis/tardis/tardis/transport/montecarlo/montecarlo_main_loop.py:123: NumbaTypeSafetyWarning:
unsafe cast from uint64 to int64. Precision may be lost.
(warnings.py:112)
[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:662)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 9.93e+03 K | 1.01e+04 K | 0.4 | 0.507 |
| 5 | 9.85e+03 K | 1.02e+04 K | 0.211 | 0.197 |
| 10 | 9.78e+03 K | 1.02e+04 K | 0.143 | 0.117 |
| 15 | 9.71e+03 K | 9.88e+03 K | 0.105 | 0.0868 |
[tardis.simulation.base][INFO ]
Current t_inner = 9933.952 K
Expected t_inner for next iteration = 10703.327 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 2 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.071e+43 erg / s
Luminosity absorbed = 3.577e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.01e+04 K | 1.08e+04 K | 0.507 | 0.524 |
| 5 | 1.02e+04 K | 1.1e+04 K | 0.197 | 0.202 |
| 10 | 1.02e+04 K | 1.08e+04 K | 0.117 | 0.124 |
| 15 | 9.88e+03 K | 1.05e+04 K | 0.0868 | 0.0926 |
[tardis.simulation.base][INFO ]
Current t_inner = 10703.327 K
Expected t_inner for next iteration = 10673.709 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 3 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.074e+43 erg / s
Luminosity absorbed = 3.397e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 1/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.08e+04 K | 1.1e+04 K | 0.524 | 0.484 |
| 5 | 1.1e+04 K | 1.12e+04 K | 0.202 | 0.19 |
| 10 | 1.08e+04 K | 1.1e+04 K | 0.124 | 0.118 |
| 15 | 1.05e+04 K | 1.06e+04 K | 0.0926 | 0.0894 |
[tardis.simulation.base][INFO ]
Current t_inner = 10673.709 K
Expected t_inner for next iteration = 10637.548 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 4 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.058e+43 erg / s
Luminosity absorbed = 3.365e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 2/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 K | 1.1e+04 K | 0.484 | 0.471 |
| 5 | 1.12e+04 K | 1.13e+04 K | 0.19 | 0.181 |
| 10 | 1.1e+04 K | 1.1e+04 K | 0.118 | 0.113 |
| 15 | 1.06e+04 K | 1.07e+04 K | 0.0894 | 0.0851 |
[tardis.simulation.base][INFO ]
Current t_inner = 10637.548 K
Expected t_inner for next iteration = 10641.159 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 5 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.057e+43 erg / s
Luminosity absorbed = 3.391e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 3/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 K | 1.1e+04 K | 0.471 | 0.48 |
| 5 | 1.13e+04 K | 1.13e+04 K | 0.181 | 0.179 |
| 10 | 1.1e+04 K | 1.1e+04 K | 0.113 | 0.115 |
| 15 | 1.07e+04 K | 1.07e+04 K | 0.0851 | 0.0838 |
[tardis.simulation.base][INFO ]
Current t_inner = 10641.159 K
Expected t_inner for next iteration = 10647.218 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 6 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.060e+43 erg / s
Luminosity absorbed = 3.394e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 4/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 K | 1.11e+04 K | 0.48 | 0.466 |
| 5 | 1.13e+04 K | 1.13e+04 K | 0.179 | 0.182 |
| 10 | 1.1e+04 K | 1.1e+04 K | 0.115 | 0.115 |
| 15 | 1.07e+04 K | 1.07e+04 K | 0.0838 | 0.0853 |
[tardis.simulation.base][INFO ]
Current t_inner = 10647.218 K
Expected t_inner for next iteration = 10645.645 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 7 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.060e+43 erg / s
Luminosity absorbed = 3.389e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 5/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 K | 1.1e+04 K | 0.466 | 0.469 |
| 5 | 1.13e+04 K | 1.14e+04 K | 0.182 | 0.173 |
| 10 | 1.1e+04 K | 1.1e+04 K | 0.115 | 0.113 |
| 15 | 1.07e+04 K | 1.06e+04 K | 0.0853 | 0.0865 |
[tardis.simulation.base][INFO ]
Current t_inner = 10645.645 K
Expected t_inner for next iteration = 10644.007 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 8 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.062e+43 erg / s
Luminosity absorbed = 3.356e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 6/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 K | 1.11e+04 K | 0.469 | 0.47 |
| 5 | 1.14e+04 K | 1.13e+04 K | 0.173 | 0.178 |
| 10 | 1.1e+04 K | 1.11e+04 K | 0.113 | 0.112 |
| 15 | 1.06e+04 K | 1.07e+04 K | 0.0865 | 0.0836 |
[tardis.simulation.base][INFO ]
Current t_inner = 10644.007 K
Expected t_inner for next iteration = 10636.775 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 9 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.053e+43 erg / s
Luminosity absorbed = 3.409e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 7/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 K | 1.1e+04 K | 0.47 | 0.471 |
| 5 | 1.13e+04 K | 1.14e+04 K | 0.178 | 0.172 |
| 10 | 1.11e+04 K | 1.12e+04 K | 0.112 | 0.109 |
| 15 | 1.07e+04 K | 1.07e+04 K | 0.0836 | 0.0836 |
[tardis.simulation.base][INFO ]
Current t_inner = 10636.775 K
Expected t_inner for next iteration = 10653.396 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 10 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.071e+43 erg / s
Luminosity absorbed = 3.320e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 8/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 K | 1.1e+04 K | 0.471 | 0.477 |
| 5 | 1.14e+04 K | 1.14e+04 K | 0.172 | 0.178 |
| 10 | 1.12e+04 K | 1.11e+04 K | 0.109 | 0.114 |
| 15 | 1.07e+04 K | 1.06e+04 K | 0.0836 | 0.0872 |
[tardis.simulation.base][INFO ]
Current t_inner = 10653.396 K
Expected t_inner for next iteration = 10624.877 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 11 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.050e+43 erg / s
Luminosity absorbed = 3.378e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 9/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 K | 1.1e+04 K | 0.477 | 0.473 |
| 5 | 1.14e+04 K | 1.12e+04 K | 0.178 | 0.182 |
| 10 | 1.11e+04 K | 1.1e+04 K | 0.114 | 0.115 |
| 15 | 1.06e+04 K | 1.06e+04 K | 0.0872 | 0.0858 |
[tardis.simulation.base][INFO ]
Current t_inner = 10624.877 K
Expected t_inner for next iteration = 10649.271 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 12 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.057e+43 erg / s
Luminosity absorbed = 3.427e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 10/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 K | 1.11e+04 K | 0.473 | 0.465 |
| 5 | 1.12e+04 K | 1.13e+04 K | 0.182 | 0.18 |
| 10 | 1.1e+04 K | 1.11e+04 K | 0.115 | 0.111 |
| 15 | 1.06e+04 K | 1.07e+04 K | 0.0858 | 0.083 |
[tardis.simulation.base][INFO ]
Current t_inner = 10649.271 K
Expected t_inner for next iteration = 10653.802 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 13 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.062e+43 erg / s
Luminosity absorbed = 3.405e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 11/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 K | 1.11e+04 K | 0.465 | 0.466 |
| 5 | 1.13e+04 K | 1.13e+04 K | 0.18 | 0.178 |
| 10 | 1.11e+04 K | 1.11e+04 K | 0.111 | 0.113 |
| 15 | 1.07e+04 K | 1.07e+04 K | 0.083 | 0.086 |
[tardis.simulation.base][INFO ]
Current t_inner = 10653.802 K
Expected t_inner for next iteration = 10646.746 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 14 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.061e+43 erg / s
Luminosity absorbed = 3.377e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 12/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 K | 1.11e+04 K | 0.466 | 0.468 |
| 5 | 1.13e+04 K | 1.13e+04 K | 0.178 | 0.181 |
| 10 | 1.11e+04 K | 1.1e+04 K | 0.113 | 0.114 |
| 15 | 1.07e+04 K | 1.06e+04 K | 0.086 | 0.0867 |
[tardis.simulation.base][INFO ]
Current t_inner = 10646.746 K
Expected t_inner for next iteration = 10641.333 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 15 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.055e+43 erg / s
Luminosity absorbed = 3.411e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 13/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 K | 1.1e+04 K | 0.468 | 0.484 |
| 5 | 1.13e+04 K | 1.13e+04 K | 0.181 | 0.18 |
| 10 | 1.1e+04 K | 1.1e+04 K | 0.114 | 0.113 |
| 15 | 1.06e+04 K | 1.07e+04 K | 0.0867 | 0.0852 |
[tardis.simulation.base][INFO ]
Current t_inner = 10641.333 K
Expected t_inner for next iteration = 10652.281 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 16 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.058e+43 erg / s
Luminosity absorbed = 3.437e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 14/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 K | 1.1e+04 K | 0.484 | 0.472 |
| 5 | 1.13e+04 K | 1.14e+04 K | 0.18 | 0.173 |
| 10 | 1.1e+04 K | 1.11e+04 K | 0.113 | 0.112 |
| 15 | 1.07e+04 K | 1.07e+04 K | 0.0852 | 0.0832 |
[tardis.simulation.base][INFO ]
Current t_inner = 10652.281 K
Expected t_inner for next iteration = 10654.516 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 17 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.066e+43 erg / s
Luminosity absorbed = 3.375e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 15/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.1e+04 K | 1.11e+04 K | 0.472 | 0.465 |
| 5 | 1.14e+04 K | 1.14e+04 K | 0.173 | 0.171 |
| 10 | 1.11e+04 K | 1.12e+04 K | 0.112 | 0.108 |
| 15 | 1.07e+04 K | 1.08e+04 K | 0.0832 | 0.0813 |
[tardis.simulation.base][INFO ]
Current t_inner = 10654.516 K
Expected t_inner for next iteration = 10637.793 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 18 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.059e+43 erg / s
Luminosity absorbed = 3.354e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 16/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 K | 1.11e+04 K | 0.465 | 0.463 |
| 5 | 1.14e+04 K | 1.14e+04 K | 0.171 | 0.174 |
| 10 | 1.12e+04 K | 1.11e+04 K | 0.108 | 0.11 |
| 15 | 1.08e+04 K | 1.07e+04 K | 0.0813 | 0.0844 |
[tardis.simulation.base][INFO ]
Current t_inner = 10637.793 K
Expected t_inner for next iteration = 10638.424 K
(base.py:657)
[tardis.simulation.base][INFO ]
Starting iteration 19 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.054e+43 erg / s
Luminosity absorbed = 3.403e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
[tardis.simulation.base][INFO ]
Iteration converged 17/4 consecutive times. (base.py:260)
[tardis.simulation.base][INFO ]
Plasma stratification: (base.py:630)
| Shell No. | t_rad | next_t_rad | w | next_w |
|---|---|---|---|---|
| 0 | 1.11e+04 K | 1.11e+04 K | 0.463 | 0.465 |
| 5 | 1.14e+04 K | 1.14e+04 K | 0.174 | 0.173 |
| 10 | 1.11e+04 K | 1.11e+04 K | 0.11 | 0.109 |
| 15 | 1.07e+04 K | 1.07e+04 K | 0.0844 | 0.0839 |
[tardis.simulation.base][INFO ]
Current t_inner = 10638.424 K
Expected t_inner for next iteration = 10651.315 K
(base.py:657)
[tardis.simulation.base][INFO ]
Simulation finished in 19 iterations
Simulation took 81.54 s
(base.py:547)
[tardis.simulation.base][INFO ]
Starting iteration 20 of 20 (base.py:444)
[tardis.simulation.base][INFO ]
Luminosity emitted = 1.061e+43 erg / s
Luminosity absorbed = 3.400e+42 erg / s
Luminosity requested = 1.059e+43 erg / s
(base.py:662)
Plotting the Spectrum¶
Finally, plot the generated spectrum with matplotlib.
[7]:
import matplotlib.pyplot as plt
[8]:
spectrum = sim.spectrum_solver.spectrum_real_packets
spectrum_virtual = sim.spectrum_solver.spectrum_virtual_packets
spectrum_integrated = sim.spectrum_solver.spectrum_integrated
[py.warnings ][WARNING]
/home/runner/work/tardis/tardis/tardis/spectrum/formal_integral.py:398: UserWarning:
The number of interpolate_shells was not specified. The value was set to 80.
(warnings.py:112)
[9]:
%matplotlib inline
plt.figure(figsize=(10, 6.5))
spectrum.plot(label="Normal packets")
spectrum_virtual.plot(label="Virtual packets")
spectrum_integrated.plot(label='Formal integral')
plt.xlim(500, 9000)
plt.title("TARDIS example model spectrum")
plt.xlabel("Wavelength [$\AA$]")
plt.ylabel("Luminosity density [erg/s/$\AA$]")
plt.legend()
plt.show()
[py.warnings ][WARNING]
<>:10: SyntaxWarning:
invalid escape sequence '\A'
(warnings.py:112)
[py.warnings ][WARNING]
<>:11: SyntaxWarning:
invalid escape sequence '\A'
(warnings.py:112)
[py.warnings ][WARNING]
<>:10: SyntaxWarning:
invalid escape sequence '\A'
(warnings.py:112)
[py.warnings ][WARNING]
<>:11: SyntaxWarning:
invalid escape sequence '\A'
(warnings.py:112)
[py.warnings ][WARNING]
/tmp/ipykernel_5932/3606673508.py:10: SyntaxWarning:
invalid escape sequence '\A'
(warnings.py:112)
[py.warnings ][WARNING]
/tmp/ipykernel_5932/3606673508.py:11: SyntaxWarning:
invalid escape sequence '\A'
(warnings.py:112)
[ ]: