You can interact with this notebook online: Launch notebook

Using Custom Callbacks When Running TARDISΒΆ

The function run_tardis allows users to provide a set of callbacks to the simulation. These callbacks are functions that will be run at the end of each iteration, and can do a variety of things, such as printing information about the simulation, storing data to a table, or even changing simulation parameters between iterations. This tutorial will show three examples of callbacks and how they can be used in TARDIS. One important thing to note is that the first argument of the callback must be the Simulation object being run.

Our first callback example will compute the (volume-weighted) average radiative temperature in the supernova ejecta (outside of the photosphere) and will print its value:

[1]:

def average_temp(sim):
    t_rads = sim.simulation_state.t_radiative
    volumes = sim.simulation_state.volume
    avg = sum(t_rads*volumes) / sum(volumes)
    print(f"Average temperature for iteration {sim.iterations_executed}: {avg}")

Now we give the callback to run_tardis. run_tardis offers the keyword argument simulation_callbacks which takes a list of lists containing the callback as well as any optional arguments you wish to include with your callback. For this example our function requires no extra arguments and we only have a single callback, so we give run_tardis a 2D list containing the callback as its only element:

[2]:

# We filter out warnings in this notebook
import warnings
warnings.filterwarnings('ignore')

from tardis import run_tardis
from tardis.io.atom_data import download_atom_data

# We download the atomic data needed to run the simulation
download_atom_data('kurucz_cd23_chianti_H_He')

# We run the simulation with our callback
sim = run_tardis('tardis_example.yml',
                 simulation_callbacks=[[average_temp]])

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.

[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)

[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)
Average temperature for iteration 1: 10002.345335990905 K
[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)
Average temperature for iteration 2: 10636.306724631839 K
[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)
Average temperature for iteration 3: 10781.451660994117 K
[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)
Average temperature for iteration 4: 10847.470538493748 K
[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)
Average temperature for iteration 5: 10852.893010414327 K
[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)
Average temperature for iteration 6: 10855.104475524438 K
[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)
Average temperature for iteration 7: 10872.983731448345 K
[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)
Average temperature for iteration 8: 10900.252453404333 K
[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)
Average temperature for iteration 9: 10928.447457469569 K
[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)
Average temperature for iteration 10: 10887.139946956373 K
[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)
Average temperature for iteration 11: 10808.256295066209 K
[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)
Average temperature for iteration 12: 10916.035502125074 K
[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)
Average temperature for iteration 13: 10890.679856236407 K
[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)
Average temperature for iteration 14: 10839.310459658938 K
[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)
Average temperature for iteration 15: 10835.970184787593 K
[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)
Average temperature for iteration 16: 10898.815372772322 K
[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)
Average temperature for iteration 17: 10993.425063230141 K
[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)
Average temperature for iteration 18: 10923.229603253854 K
[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)
Average temperature for iteration 19: 10915.386068997337 K
[tardis.simulation.base][INFO   ]

        Simulation finished in 19 iterations
        Simulation took 79.11 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)
Average temperature for iteration 20: 10915.386068997337 K

Running Callbacks with Extra ArgumentsΒΆ

The callbacks provided to run_tardis can also take extra arguments. As an example, we’ll make a callback that appends the number of monte carlo packets emitted by the supernova for each iteration to a list so we can plot the number of emitted packets for each iteration. We will also specify that we want this information for all but the last iteration, as more packets are used in the last iteration than are used in the other iterations. The callback will take a list we want to append to as an argument. We’ll send both this new callback and our original average_temp callback to run_tardis as an example of using multiple callbacks at once:

[3]:

def append_num_emitted_to_list(sim, lst):
    if sim.iterations_executed < sim.iterations:
        num_emitted_packets = len(sim.transport.transport_state.emitted_packet_nu)
        lst.append(num_emitted_packets)

In order to add our new callback, we just create another entry in our list of callbacks. Since append_num_emitted_to_list takes an extra argument, we will provide that argument in the inner list containing the callback:

[4]:

# Initialize a list to store the number of emitted packets
num_emitted_list = []

# Make our list of callbacks
callbacks = [[average_temp],
             [append_num_emitted_to_list, num_emitted_list]]

# Run the simulation with both of our callbacks
sim = run_tardis('tardis_example.yml',
                 simulation_callbacks=callbacks)

[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)

[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)
Average temperature for iteration 1: 10002.345335990905 K
[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)
Average temperature for iteration 2: 10636.306724631839 K
[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)
Average temperature for iteration 3: 10781.451660994117 K
[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)
Average temperature for iteration 4: 10847.470538493748 K
[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)
Average temperature for iteration 5: 10852.893010414327 K
[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)
Average temperature for iteration 6: 10855.104475524438 K
[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)
Average temperature for iteration 7: 10872.983731448345 K
[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)
Average temperature for iteration 8: 10900.252453404333 K
[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)
Average temperature for iteration 9: 10928.447457469569 K
[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)
Average temperature for iteration 10: 10887.139946956373 K
[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)
Average temperature for iteration 11: 10808.256295066209 K
[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)
Average temperature for iteration 12: 10916.035502125074 K
[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)
Average temperature for iteration 13: 10890.679856236407 K
[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)
Average temperature for iteration 14: 10839.310459658938 K
[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)
Average temperature for iteration 15: 10835.970184787593 K
[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)
Average temperature for iteration 16: 10898.815372772322 K
[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)
Average temperature for iteration 17: 10993.425063230141 K
[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)
Average temperature for iteration 18: 10923.229603253854 K
[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)
Average temperature for iteration 19: 10915.386068997337 K
[tardis.simulation.base][INFO   ]

        Simulation finished in 19 iterations
        Simulation took 16.58 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)
Average temperature for iteration 20: 10915.386068997337 K

Now we can look at how many packets are emitted after each iteration:

[5]:

import matplotlib.pyplot as plt

# Generate a list of each iteration number for the x-axis
iterations = list(range(1, len(num_emitted_list)+1))

# Plot the number of emitted packets
plt.plot(iterations, num_emitted_list)
plt.xlabel("Iteration")
plt.ylabel("Number of emitted packets")

[5]:

Text(0, 0.5, 'Number of emitted packets')
../../_images/io_optional_tutorial_callback_example_10_1.svg

Using Callbacks to Add New FunctionalityΒΆ

Callbacks can also add new functionality to the code. For example, we introduce one final callback inc_packets that will increase the number of packets in the following iteration by a number \(N\) (which is an argument to the callback, in our example we shall use \(N=1000\)):

[6]:

def inc_packets(sim, N):
    sim.no_of_packets += N

[7]:

# Initialize a new list to store the number of emitted packets
num_emitted_list_new = []

# Make our new list of callbacks
callbacks = [[average_temp],
             [append_num_emitted_to_list, num_emitted_list_new],
             [inc_packets, 1000]]

# Run the simulation with all three of our callbacks
sim = run_tardis('tardis_example.yml',
                 simulation_callbacks=callbacks)

[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)

[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)
Average temperature for iteration 1: 10002.345335990905 K
[tardis.simulation.base][INFO   ]

        Starting iteration 2 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.075e+43 erg / s
        Luminosity absorbed  = 3.538e+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.09e+04 K 0.507 0.515
5 1.02e+04 K 1.1e+04 K 0.197 0.203
10 1.02e+04 K 1.07e+04 K 0.117 0.13
15 9.88e+03 K 1.05e+04 K 0.0868 0.0916 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10703.327 K
        Expected t_inner for next iteration = 10663.259 K
 (base.py:657)
Average temperature for iteration 2: 10613.908754846558 K
[tardis.simulation.base][INFO   ]

        Starting iteration 3 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.070e+43 erg / s
        Luminosity absorbed  = 3.381e+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.09e+04 K 1.1e+04 K 0.515 0.486
5 1.1e+04 K 1.13e+04 K 0.203 0.18
10 1.07e+04 K 1.1e+04 K 0.13 0.116
15 1.05e+04 K 1.07e+04 K 0.0916 0.0851 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10663.259 K
        Expected t_inner for next iteration = 10637.291 K
 (base.py:657)
Average temperature for iteration 3: 10850.348012607232 K
[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.363e+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.1e+04 K 1.1e+04 K 0.486 0.473
5 1.13e+04 K 1.13e+04 K 0.18 0.178
10 1.1e+04 K 1.11e+04 K 0.116 0.111
15 1.07e+04 K 1.07e+04 K 0.0851 0.0841 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10637.291 K
        Expected t_inner for next iteration = 10640.804 K
 (base.py:657)
Average temperature for iteration 4: 10887.571829540853 K
[tardis.simulation.base][INFO   ]

        Starting iteration 5 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.056e+43 erg / s
        Luminosity absorbed  = 3.402e+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.473 0.473
5 1.13e+04 K 1.12e+04 K 0.178 0.184
10 1.11e+04 K 1.1e+04 K 0.111 0.114
15 1.07e+04 K 1.07e+04 K 0.0841 0.0843 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10640.804 K
        Expected t_inner for next iteration = 10650.106 K
 (base.py:657)
Average temperature for iteration 5: 10830.30318618678 K
[tardis.simulation.base][INFO   ]

        Starting iteration 6 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.062e+43 erg / s
        Luminosity absorbed  = 3.388e+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.473 0.474
5 1.12e+04 K 1.13e+04 K 0.184 0.18
10 1.1e+04 K 1.11e+04 K 0.114 0.112
15 1.07e+04 K 1.07e+04 K 0.0843 0.0859 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10650.106 K
        Expected t_inner for next iteration = 10643.291 K
 (base.py:657)
Average temperature for iteration 6: 10858.234351239738 K
[tardis.simulation.base][INFO   ]

        Starting iteration 7 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.059e+43 erg / s
        Luminosity absorbed  = 3.388e+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.1e+04 K 0.474 0.474
5 1.13e+04 K 1.12e+04 K 0.18 0.185
10 1.11e+04 K 1.1e+04 K 0.112 0.115
15 1.07e+04 K 1.05e+04 K 0.0859 0.0891 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10643.291 K
        Expected t_inner for next iteration = 10644.932 K
 (base.py:657)
Average temperature for iteration 7: 10792.481079535482 K
[tardis.simulation.base][INFO   ]

        Starting iteration 8 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.058e+43 erg / s
        Luminosity absorbed  = 3.396e+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.1e+04 K 1.1e+04 K 0.474 0.481
5 1.12e+04 K 1.14e+04 K 0.185 0.177
10 1.1e+04 K 1.09e+04 K 0.115 0.118
15 1.05e+04 K 1.07e+04 K 0.0891 0.0849 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10644.932 K
        Expected t_inner for next iteration = 10647.629 K
 (base.py:657)
Average temperature for iteration 8: 10838.693927807315 K
[tardis.simulation.base][INFO   ]

        Starting iteration 9 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.060e+43 erg / s
        Luminosity absorbed  = 3.398e+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.481 0.463
5 1.14e+04 K 1.13e+04 K 0.177 0.178
10 1.09e+04 K 1.12e+04 K 0.118 0.109
15 1.07e+04 K 1.08e+04 K 0.0849 0.0824 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10647.629 K
        Expected t_inner for next iteration = 10646.761 K
 (base.py:657)
Average temperature for iteration 9: 10947.325467890487 K
[tardis.simulation.base][INFO   ]

        Starting iteration 10 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.060e+43 erg / s
        Luminosity absorbed  = 3.390e+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.11e+04 K 0.463 0.472
5 1.13e+04 K 1.14e+04 K 0.178 0.174
10 1.12e+04 K 1.11e+04 K 0.109 0.112
15 1.08e+04 K 1.07e+04 K 0.0824 0.0853 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10646.761 K
        Expected t_inner for next iteration = 10645.487 K
 (base.py:657)
Average temperature for iteration 10: 10904.884829767207 K
[tardis.simulation.base][INFO   ]

        Starting iteration 11 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.059e+43 erg / s
        Luminosity absorbed  = 3.392e+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.11e+04 K 1.1e+04 K 0.472 0.473
5 1.14e+04 K 1.14e+04 K 0.174 0.173
10 1.11e+04 K 1.11e+04 K 0.112 0.112
15 1.07e+04 K 1.07e+04 K 0.0853 0.0843 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10645.487 K
        Expected t_inner for next iteration = 10646.183 K
 (base.py:657)
Average temperature for iteration 11: 10909.70977691611 K
[tardis.simulation.base][INFO   ]

        Starting iteration 12 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.059e+43 erg / s
        Luminosity absorbed  = 3.393e+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.11e+04 K 0.473 0.466
5 1.14e+04 K 1.15e+04 K 0.173 0.17
10 1.11e+04 K 1.11e+04 K 0.112 0.111
15 1.07e+04 K 1.07e+04 K 0.0843 0.0846 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10646.183 K
        Expected t_inner for next iteration = 10645.990 K
 (base.py:657)
Average temperature for iteration 12: 10936.034767331808 K
[tardis.simulation.base][INFO   ]

        Starting iteration 13 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.065e+43 erg / s
        Luminosity absorbed  = 3.336e+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.11e+04 K 1.1e+04 K 0.466 0.473
5 1.15e+04 K 1.13e+04 K 0.17 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.0846 0.085 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10645.990 K
        Expected t_inner for next iteration = 10631.183 K
 (base.py:657)
Average temperature for iteration 13: 10914.993081485425 K
[tardis.simulation.base][INFO   ]

        Starting iteration 14 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.054e+43 erg / s
        Luminosity absorbed  = 3.367e+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.1e+04 K 1.1e+04 K 0.473 0.467
5 1.13e+04 K 1.13e+04 K 0.178 0.176
10 1.11e+04 K 1.11e+04 K 0.113 0.111
15 1.07e+04 K 1.08e+04 K 0.085 0.0824 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10631.183 K
        Expected t_inner for next iteration = 10644.219 K
 (base.py:657)
Average temperature for iteration 14: 10926.52320767812 K
[tardis.simulation.base][INFO   ]

        Starting iteration 15 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.057e+43 erg / s
        Luminosity absorbed  = 3.400e+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.1e+04 K 1.1e+04 K 0.467 0.479
5 1.13e+04 K 1.14e+04 K 0.176 0.172
10 1.11e+04 K 1.11e+04 K 0.111 0.111
15 1.08e+04 K 1.08e+04 K 0.0824 0.0834 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10644.219 K
        Expected t_inner for next iteration = 10649.204 K
 (base.py:657)
Average temperature for iteration 15: 10919.420818412638 K
[tardis.simulation.base][INFO   ]

        Starting iteration 16 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.064e+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 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.1e+04 K 1.1e+04 K 0.479 0.47
5 1.14e+04 K 1.13e+04 K 0.172 0.177
10 1.11e+04 K 1.11e+04 K 0.111 0.113
15 1.08e+04 K 1.07e+04 K 0.0834 0.0838 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10649.204 K
        Expected t_inner for next iteration = 10637.338 K
 (base.py:657)
Average temperature for iteration 16: 10898.993427208812 K
[tardis.simulation.base][INFO   ]

        Starting iteration 17 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.055e+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 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.47 0.47
5 1.13e+04 K 1.14e+04 K 0.177 0.171
10 1.11e+04 K 1.11e+04 K 0.113 0.112
15 1.07e+04 K 1.07e+04 K 0.0838 0.0858 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10637.338 K
        Expected t_inner for next iteration = 10649.272 K
 (base.py:657)
Average temperature for iteration 17: 10868.841791849643 K
[tardis.simulation.base][INFO   ]

        Starting iteration 18 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.063e+43 erg / s
        Luminosity absorbed  = 3.373e+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.47 0.467
5 1.14e+04 K 1.14e+04 K 0.171 0.173
10 1.11e+04 K 1.11e+04 K 0.112 0.112
15 1.07e+04 K 1.06e+04 K 0.0858 0.0866 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10649.272 K
        Expected t_inner for next iteration = 10638.973 K
 (base.py:657)
Average temperature for iteration 18: 10868.404619085652 K
[tardis.simulation.base][INFO   ]

        Starting iteration 19 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.058e+43 erg / s
        Luminosity absorbed  = 3.371e+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.467 0.465
5 1.14e+04 K 1.14e+04 K 0.173 0.176
10 1.11e+04 K 1.11e+04 K 0.112 0.109
15 1.06e+04 K 1.08e+04 K 0.0866 0.0822 
[tardis.simulation.base][INFO   ]

        Current t_inner = 10638.973 K
        Expected t_inner for next iteration = 10643.052 K
 (base.py:657)
Average temperature for iteration 19: 10935.815023903202 K
[tardis.simulation.base][INFO   ]

        Simulation finished in 19 iterations
        Simulation took 19.27 s
 (base.py:547)
[tardis.simulation.base][INFO   ]

        Starting iteration 20 of 20 (base.py:444)
[tardis.simulation.base][INFO   ]

        Luminosity emitted   = 1.058e+43 erg / s
        Luminosity absorbed  = 3.390e+42 erg / s
        Luminosity requested = 1.059e+43 erg / s
 (base.py:662)
Average temperature for iteration 20: 10935.815023903202 K

Now, let’s see how this affected our plot for packets emitted in each iteration:

[8]:

plt.plot(iterations, num_emitted_list_new)
plt.xlabel("Iteration")
plt.ylabel("Number of emitted packets")

[8]:

Text(0, 0.5, 'Number of emitted packets')
../../_images/io_optional_tutorial_callback_example_15_1.svg

As expected, the number of packets emitted will keep on increasing as 1000 more packets are run each iteration.