From 009b589db14e75e040874abbdc9f19bbea4271a4 Mon Sep 17 00:00:00 2001 From: Wolfgang Kerzendorf Date: Tue, 23 Jul 2019 21:09:21 -0400 Subject: [PATCH] add code comparison (#953) --- docs/research/code_comparison/index.rst | 4 + .../toy_models/reading blondin toymodel.ipynb | 2635 +++++++++++++++++ 2 files changed, 2639 insertions(+) create mode 100644 docs/research/code_comparison/toy_models/reading blondin toymodel.ipynb diff --git a/docs/research/code_comparison/index.rst b/docs/research/code_comparison/index.rst index 0eaeb23a7d6..da2ec52f850 100644 --- a/docs/research/code_comparison/index.rst +++ b/docs/research/code_comparison/index.rst @@ -6,3 +6,7 @@ codes to understand systematics uncertainties and different strategies. The first meeting is summarized `here `_. +.. toctree:: + + toy_models/reading blondin toymodel + diff --git a/docs/research/code_comparison/toy_models/reading blondin toymodel.ipynb b/docs/research/code_comparison/toy_models/reading blondin toymodel.ipynb new file mode 100644 index 00000000000..3769031248d --- /dev/null +++ b/docs/research/code_comparison/toy_models/reading blondin toymodel.ipynb @@ -0,0 +1,2635 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### Compare Blondin Toy Models ###" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/tqdm/autonotebook/__init__.py:14: TqdmExperimentalWarning: Using `tqdm.autonotebook.tqdm` in notebook mode. Use `tqdm.tqdm` instead to force console mode (e.g. in jupyter console)\n", + " \" (e.g. in jupyter console)\", TqdmExperimentalWarning)\n" + ] + } + ], + "source": [ + "import pandas as pd\n", + "import re\n", + "from tardis.util.base import parse_quantity\n", + "import numpy as np\n", + "import yaml\n", + "from tardis import run_tardis\n", + "from astropy import units as u" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [], + "source": [ + "\n", + "pattern_remove_bracket = re.compile('\\[.+\\]')\n", + "t0_pattern = re.compile('tend = (.+)\\n')\n", + "\n", + "\n", + "def read_blondin_toymodel(fname):\n", + " with open(fname, 'r') as fh:\n", + " for line in fh:\n", + " if line.startswith(\"#idx\"):\n", + " break\n", + " else:\n", + " raise ValueError('File {0} does not conform to Toy Model format as it does not contain #idx')\n", + " columns = [pattern_remove_bracket.sub('', item) for item in line[1:].split()]\n", + " \n", + " raw_blondin_csv = pd.read_csv(fname, delim_whitespace=True, comment='#', header=None, names=columns)\n", + " raw_blondin_csv.set_index('idx', inplace=True)\n", + " \n", + " blondin_csv = raw_blondin_csv.loc[:, ['vel', 'dens', 'temp', 'X_56Ni0', 'X_Ti', 'X_Ca', 'X_S', 'X_Si', 'X_O', 'X_C']]\n", + " rename_col_dict = {'vel':'velocity', 'dens':'density', 'temp':'t_electron'}\n", + " rename_col_dict.update({item:item[2:] for item in blondin_csv.columns[3:]})\n", + " rename_col_dict['X_56Ni0'] = 'Ni56'\n", + " blondin_csv.rename(columns=rename_col_dict, inplace=True)\n", + " blondin_csv.iloc[:, 3:] = blondin_csv.iloc[:,3:].divide(blondin_csv.iloc[:,3:].sum(axis=1), axis=0)\n", + " \n", + " \n", + " #changing velocities to outer boundary\n", + " new_velocities = 0.5 * (blondin_csv.velocity.iloc[:-1].values + blondin_csv.velocity.iloc[1:].values)\n", + " new_velocities = np.hstack((new_velocities, [2 * new_velocities[-1] - new_velocities[-2]]))\n", + " blondin_csv['velocity'] = new_velocities\n", + " \n", + " \n", + " \n", + " with open(fname, 'r') as fh:\n", + " t0_string = t0_pattern.findall(fh.read())[0]\n", + " \n", + " t0 = parse_quantity(t0_string.replace('DAYS', 'day'))\n", + " blondin_dict = {}\n", + " blondin_dict['model_density_time_0'] = str(t0)\n", + " blondin_dict['description'] = 'Converted {0} to csvy format'.format(fname)\n", + " blondin_dict['tardis_model_config_version'] = 'v1.0'\n", + " blondin_dict_fields = [dict(name='velocity', unit='km/s', desc='velocities of shell outer bounderies.')]\n", + " blondin_dict_fields.append(dict(name='density', unit='g/cm^3', desc='mean density of shell.'))\n", + " blondin_dict_fields.append(dict(name='t_electron', unit='K', desc='electron temperature.'))\n", + " \n", + " for abund in blondin_csv.columns[3:]:\n", + " blondin_dict_fields.append(dict(name=abund, desc='Fraction {0} abundance'.format(abund)))\n", + " blondin_dict['datatype'] = {'fields':blondin_dict_fields}\n", + " \n", + " return blondin_dict, blondin_csv" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [], + "source": [ + "blondin_dict, blondin_csv = read_blondin_toymodel('snia_toy01.dat')" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": {}, + "outputs": [], + "source": [ + "blondin_dict['v_inner_boundary'] = '9000 km/s'\n", + "blondin_dict['v_outer_boundary'] = '20000 km/s'\n", + "\n" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "outputs": [], + "source": [ + "csvy_file = '---\\n{0}\\n---\\n{1}'.format(yaml.dump(blondin_dict, default_flow_style=False), blondin_csv.to_csv(index=False))" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": {}, + "outputs": [], + "source": [ + "with open('blondin_compare.csvy', 'w') as fh:\n", + " fh.write(csvy_file)" + ] + }, + { + "cell_type": "code", + "execution_count": 7, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/importlib/_bootstrap.py:219: QAWarning: pyne.data is not yet QA compliant.\n", + " return f(*args, **kwds)\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n", + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/importlib/_bootstrap.py:219: QAWarning: pyne.material is not yet QA compliant.\n", + " return f(*args, **kwds)\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n", + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.\n", + " AstropyDeprecationWarning)\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n", + "[\u001b[1mtardis.plasma.standard_plasmas\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Reading Atomic Data from kurucz_cd23_chianti_H_He.h5 (\u001b[1mstandard_plasmas.py\u001b[0m:74)\n", + "[\u001b[1mtardis.io.atom_data.util\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Atom Data kurucz_cd23_chianti_H_He.h5 not found in local path. Exists in TARDIS Data repo /Users/wkerzend/projects/tardis/tardis-data/kurucz_cd23_chianti_H_He.h5 (\u001b[1mutil.py\u001b[0m:29)\n", + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/IPython/core/interactiveshell.py:3267: PerformanceWarning: indexing past lexsort depth may impact performance.\n", + " exec(code_obj, self.user_global_ns, self.user_ns)\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n", + "[\u001b[1mtardis.io.atom_data.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Read Atom Data with UUID=6f7b09e887a311e7a06b246e96350010 and MD5=864f1753714343c41f99cb065710cace. (\u001b[1mbase.py\u001b[0m:184)\n", + "[\u001b[1mtardis.io.atom_data.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Non provided atomic data: synpp_refs, photoionization_data (\u001b[1mbase.py\u001b[0m:187)\n", + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.\n", + " AstropyDeprecationWarning)\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n", + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/python/tardis/tardis/plasma/properties/atomic.py:195: FutureWarning: .labels was deprecated in version 0.24.0. Use .codes instead.\n", + " zeta_data['atomic_number'] = zeta_data.index.labels[0] + 1\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n", + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/python/tardis/tardis/plasma/properties/atomic.py:196: FutureWarning: .labels was deprecated in version 0.24.0. Use .codes instead.\n", + " zeta_data['ion_number'] = zeta_data.index.labels[1] + 1\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n", + "[\u001b[1mtardis.plasma.properties.atomic\u001b[0m][\u001b[1;33mWARNING\u001b[0m] Zeta_data missing - replaced with 1s. Missing ions: [(14, 15), (16, 17), (20, 21), (26, 27), (27, 28), (28, 29)] (\u001b[1matomic.py\u001b[0m:215)\n", + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/python/tardis/tardis/plasma/properties/ion_population.py:63: FutureWarning: \n", + "Passing list-likes to .loc or [] with any missing label will raise\n", + "KeyError in the future, you can use .reindex() as an alternative.\n", + "\n", + "See the documentation here:\n", + "https://pandas.pydata.org/pandas-docs/stable/indexing.html#deprecate-loc-reindex-listlike\n", + " partition_function.index].dropna())\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n", + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/equivalencies.py:90: RuntimeWarning: divide by zero encountered in double_scalars\n", + " (si.m, si.Hz, lambda x: _si.c.value / x),\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n", + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.\n", + " AstropyDeprecationWarning)\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 1/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.22587e+43 erg / s Luminosity absorbed = 8.12489e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 12998.916007 13061.816463 0.462810 0.567717\n", + "\t5 12988.086012 13056.730857 0.379149 0.480626\n", + "\t10 12977.274048 13025.637830 0.336294 0.428170\n", + "\t15 12966.480071 13033.048290 0.304883 0.382503\n", + "\t20 12955.704034 13065.468766 0.279753 0.343327\n", + "\t25 12944.945893 13026.575503 0.258765 0.318805\n", + "\t30 12934.205605 13031.454768 0.240774 0.290452\n", + "\t35 12923.483123 13012.104561 0.225078 0.269456\n", + "\t40 12912.778405 12984.321293 0.211205 0.251158\n", + "\t45 12902.091406 12960.901645 0.198821 0.232334\n", + "\t50 12891.422083 12941.563429 0.187679 0.217530\n", + "\t55 12880.770390 12921.758766 0.177588 0.203183\n", + "\t60 12870.136285 12942.591441 0.168399 0.186982\n", + "\t65 12859.519724 12916.751988 0.159992 0.176438\n", + "\t70 12848.920665 12981.212652 0.152269 0.162411\n", + "\t75 12838.339062 12932.407924 0.145148 0.153813\n", + "\t80 12827.774874 12916.331855 0.138563 0.145349\n", + "\t85 12817.228057 12917.001822 0.132453 0.136457\n", + "\t90 12806.698569 12887.128678 0.126772 0.129621\n", + "\t95 12796.186367 12888.946086 0.121475 0.122397\n", + "\t100 12785.691409 12881.697128 0.116527 0.116324\n", + "\t105 12775.213651 12895.681422 0.111894 0.109162\n", + "\t110 12764.753053 12851.913702 0.107550 0.105005\n", + "\t115 12754.309570 12795.235678 0.103468 0.101558\n", + "\t120 12743.883163 12778.835290 0.099626 0.097184\n", + "\t125 12733.473789 12747.043599 0.096006 0.093453\n", + "\t130 12723.081405 12743.322681 0.092589 0.089142\n", + "\t135 12712.705972 12704.356267 0.089359 0.085928\n", + "\t140 12702.347446 12731.984118 0.086303 0.081671\n", + "\t145 12692.005787 12690.395130 0.083408 0.079170\n", + "\t150 12681.680954 12675.569081 0.080662 0.075764\n", + "\t155 12671.372906 12673.798388 0.078054 0.072396\n", + "\t160 12661.081601 12688.321977 0.075575 0.069417\n", + "\t165 12650.807000 12602.489519 0.073217 0.068280\n", + "\t170 12640.549060 12575.649977 0.070971 0.066077\n", + "\t175 12630.307743 12552.103659 0.068830 0.064002\n", + "\t180 12620.083007 12524.123590 0.066787 0.062033\n", + "\t185 12609.874813 12470.157871 0.064837 0.060517\n", + "\t190 12599.683119 12452.067488 0.062973 0.058731\n", + "\t195 12589.507887 12415.450898 0.061191 0.057065\n", + "\t200 12579.349076 12377.598670 0.059485 0.055637\n", + "\t205 12569.206647 12369.427919 0.057852 0.053812\n", + "\t210 12559.080559 12308.221817 0.056286 0.052806\n", + "\t215 12548.970775 12274.329208 0.054785 0.051496\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/python/tardis/tardis/plasma/properties/ion_population.py:63: FutureWarning: \n", + "Passing list-likes to .loc or [] with any missing label will raise\n", + "KeyError in the future, you can use .reindex() as an alternative.\n", + "\n", + "See the documentation here:\n", + "https://pandas.pydata.org/pandas-docs/stable/indexing.html#deprecate-loc-reindex-listlike\n", + " partition_function.index].dropna())\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 2/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.24347e+43 erg / s Luminosity absorbed = 7.94642e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13061.816463 13117.652070 0.567717 0.613290\n", + "\t5 13056.730857 13128.710972 0.480626 0.520966\n", + "\t10 13025.637830 13087.946452 0.428170 0.469492\n", + "\t15 13033.048290 13117.221906 0.382503 0.416974\n", + "\t20 13065.468766 13114.807782 0.343327 0.377475\n", + "\t25 13026.575503 13096.095422 0.318805 0.347216\n", + "\t30 13031.454768 13094.720257 0.290452 0.315541\n", + "\t35 13012.104561 13060.231563 0.269456 0.292868\n", + "\t40 12984.321293 13059.642743 0.251158 0.267728\n", + "\t45 12960.901645 13036.502049 0.232334 0.248840\n", + "\t50 12941.563429 13088.971015 0.217530 0.225795\n", + "\t55 12921.758766 13057.048962 0.203183 0.211304\n", + "\t60 12942.591441 13012.342405 0.186982 0.197369\n", + "\t65 12916.751988 13011.652080 0.176438 0.182820\n", + "\t70 12981.212652 13053.319564 0.162411 0.169516\n", + "\t75 12932.407924 13074.224102 0.153813 0.156755\n", + "\t80 12916.331855 13029.565434 0.145349 0.147823\n", + "\t85 12917.001822 13042.421261 0.136457 0.137239\n", + "\t90 12887.128678 12990.734303 0.129621 0.130452\n", + "\t95 12888.946086 12973.415343 0.122397 0.122766\n", + "\t100 12881.697128 12959.651664 0.116324 0.116410\n", + "\t105 12895.681422 12960.317508 0.109162 0.109050\n", + "\t110 12851.913702 12929.769910 0.105005 0.103908\n", + "\t115 12795.235678 12913.972612 0.101558 0.098979\n", + "\t120 12778.835290 12886.386512 0.097184 0.094654\n", + "\t125 12747.043599 12843.837749 0.093453 0.090621\n", + "\t130 12743.322681 12816.084072 0.089142 0.086600\n", + "\t135 12704.356267 12790.203288 0.085928 0.082578\n", + "\t140 12731.984118 12796.350067 0.081671 0.078901\n", + "\t145 12690.395130 12791.502994 0.079170 0.075382\n", + "\t150 12675.569081 12704.857834 0.075764 0.073887\n", + "\t155 12673.798388 12655.990720 0.072396 0.070855\n", + "\t160 12688.321977 12610.883401 0.069417 0.068569\n", + "\t165 12602.489519 12558.138741 0.068280 0.066569\n", + "\t170 12575.649977 12553.564630 0.066077 0.063869\n", + "\t175 12552.103659 12520.095372 0.064002 0.061832\n", + "\t180 12524.123590 12457.646020 0.062033 0.060085\n", + "\t185 12470.157871 12434.132208 0.060517 0.058028\n", + "\t190 12452.067488 12380.739186 0.058731 0.056670\n", + "\t195 12415.450898 12353.512233 0.057065 0.054970\n", + "\t200 12377.598670 12315.990345 0.055637 0.053348\n", + "\t205 12369.427919 12289.918274 0.053812 0.051717\n", + "\t210 12308.221817 12224.340350 0.052806 0.050777\n", + "\t215 12274.329208 12164.080107 0.051496 0.049729\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 3/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.23675e+43 erg / s Luminosity absorbed = 8.01326e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13117.652070 13122.102732 0.613290 0.641024\n", + "\t5 13128.710972 13164.042522 0.520966 0.543481\n", + "\t10 13087.946452 13113.990367 0.469492 0.494977\n", + "\t15 13117.221906 13173.538351 0.416974 0.432339\n", + "\t20 13114.807782 13144.097172 0.377475 0.392818\n", + "\t25 13096.095422 13164.934622 0.347216 0.356523\n", + "\t30 13094.720257 13184.070310 0.315541 0.319199\n", + "\t35 13060.231563 13136.832241 0.292868 0.297554\n", + "\t40 13059.642743 13136.025386 0.267728 0.273157\n", + "\t45 13036.502049 13155.265509 0.248840 0.249419\n", + "\t50 13088.971015 13166.480491 0.225795 0.228277\n", + "\t55 13057.048962 13166.079767 0.211304 0.210512\n", + "\t60 13012.342405 13127.443828 0.197369 0.197161\n", + "\t65 13011.652080 13141.214262 0.182820 0.181772\n", + "\t70 13053.319564 13174.463561 0.169516 0.167142\n", + "\t75 13074.224102 13150.237468 0.156755 0.157056\n", + "\t80 13029.565434 13097.524561 0.147823 0.147829\n", + "\t85 13042.421261 13127.609151 0.137239 0.137357\n", + "\t90 12990.734303 13040.008757 0.130452 0.130839\n", + "\t95 12973.415343 13029.856575 0.122766 0.122425\n", + "\t100 12959.651664 12982.871415 0.116410 0.116097\n", + "\t105 12960.317508 12979.071109 0.109050 0.109282\n", + "\t110 12929.769910 12943.014800 0.103908 0.104322\n", + "\t115 12913.972612 12927.071313 0.098979 0.098592\n", + "\t120 12886.386512 12886.480175 0.094654 0.094474\n", + "\t125 12843.837749 12864.414625 0.090621 0.089431\n", + "\t130 12816.084072 12874.140628 0.086600 0.084479\n", + "\t135 12790.203288 12818.832488 0.082578 0.081002\n", + "\t140 12796.350067 12826.258359 0.078901 0.076974\n", + "\t145 12791.502994 12795.258432 0.075382 0.073890\n", + "\t150 12704.857834 12722.903444 0.073887 0.072095\n", + "\t155 12655.990720 12672.627773 0.070855 0.069298\n", + "\t160 12610.883401 12629.678970 0.068569 0.066939\n", + "\t165 12558.138741 12586.628880 0.066569 0.064793\n", + "\t170 12553.564630 12573.279154 0.063869 0.062097\n", + "\t175 12520.095372 12536.734639 0.061832 0.059908\n", + "\t180 12457.646020 12499.514570 0.060085 0.057760\n", + "\t185 12434.132208 12489.203566 0.058028 0.055605\n", + "\t190 12380.739186 12388.978106 0.056670 0.054919\n", + "\t195 12353.512233 12311.880538 0.054970 0.054015\n", + "\t200 12315.990345 12300.791731 0.053348 0.051990\n", + "\t205 12289.918274 12254.716475 0.051717 0.050505\n", + "\t210 12224.340350 12226.649994 0.050777 0.049039\n", + "\t215 12164.080107 12126.030478 0.049729 0.048560\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 4/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.23576e+43 erg / s Luminosity absorbed = 8.01378e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13122.102732 13166.764592 0.641024 0.647005\n", + "\t5 13164.042522 13160.763682 0.543481 0.557020\n", + "\t10 13113.990367 13167.628696 0.494977 0.498091\n", + "\t15 13173.538351 13190.751769 0.432339 0.444245\n", + "\t20 13144.097172 13184.024122 0.392818 0.401907\n", + "\t25 13164.934622 13227.151619 0.356523 0.359620\n", + "\t30 13184.070310 13180.743685 0.319199 0.331947\n", + "\t35 13136.832241 13204.580495 0.297554 0.300127\n", + "\t40 13136.025386 13194.812326 0.273157 0.275673\n", + "\t45 13155.265509 13230.326739 0.249419 0.248915\n", + "\t50 13166.480491 13234.880733 0.228277 0.229173\n", + "\t55 13166.079767 13233.024725 0.210512 0.211402\n", + "\t60 13127.443828 13173.721736 0.197161 0.199245\n", + "\t65 13141.214262 13236.160576 0.181772 0.180157\n", + "\t70 13174.463561 13220.204788 0.167142 0.166780\n", + "\t75 13150.237468 13206.871230 0.157056 0.156329\n", + "\t80 13097.524561 13191.695756 0.147829 0.146271\n", + "\t85 13127.609151 13132.044552 0.137357 0.138585\n", + "\t90 13040.008757 13107.559438 0.130839 0.129874\n", + "\t95 13029.856575 13096.184312 0.122425 0.121443\n", + "\t100 12982.871415 13051.565010 0.116097 0.115060\n", + "\t105 12979.071109 13026.202501 0.109282 0.108914\n", + "\t110 12943.014800 12971.433048 0.104322 0.103877\n", + "\t115 12927.071313 12951.898476 0.098592 0.098483\n", + "\t120 12886.480175 12938.304729 0.094474 0.093136\n", + "\t125 12864.414625 12899.761435 0.089431 0.088597\n", + "\t130 12874.140628 12882.151103 0.084479 0.084103\n", + "\t135 12818.832488 12882.960846 0.081002 0.079350\n", + "\t140 12826.258359 12869.549641 0.076974 0.075580\n", + "\t145 12795.258432 12842.281777 0.073890 0.072636\n", + "\t150 12722.903444 12749.157914 0.072095 0.070611\n", + "\t155 12672.627773 12703.241916 0.069298 0.067873\n", + "\t160 12629.678970 12670.297608 0.066939 0.065557\n", + "\t165 12586.628880 12636.895690 0.064793 0.063296\n", + "\t170 12573.279154 12597.918830 0.062097 0.061002\n", + "\t175 12536.734639 12550.742876 0.059908 0.058575\n", + "\t180 12499.514570 12528.304583 0.057760 0.056474\n", + "\t185 12489.203566 12475.332811 0.055605 0.054925\n", + "\t190 12388.978106 12397.080100 0.054919 0.054007\n", + "\t195 12311.880538 12372.673017 0.054015 0.052291\n", + "\t200 12300.791731 12355.614341 0.051990 0.050307\n", + "\t205 12254.716475 12286.492951 0.050505 0.049312\n", + "\t210 12226.649994 12227.262848 0.049039 0.048146\n", + "\t215 12126.030478 12110.655363 0.048560 0.047822\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 5/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.23135e+43 erg / s Luminosity absorbed = 8.05807e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13166.764592 13133.133359 0.647005 0.661105\n", + "\t5 13160.763682 13180.154270 0.557020 0.561426\n", + "\t10 13167.628696 13187.356843 0.498091 0.499942\n", + "\t15 13190.751769 13230.009983 0.444245 0.443721\n", + "\t20 13184.024122 13221.404988 0.401907 0.401814\n", + "\t25 13227.151619 13279.927924 0.359620 0.358538\n", + "\t30 13180.743685 13213.416642 0.331947 0.331098\n", + "\t35 13204.580495 13307.286771 0.300127 0.293396\n", + "\t40 13194.812326 13257.000678 0.275673 0.272724\n", + "\t45 13230.326739 13266.454328 0.248915 0.249218\n", + "\t50 13234.880733 13285.344929 0.229173 0.228049\n", + "\t55 13233.024725 13311.193287 0.211402 0.209753\n", + "\t60 13173.721736 13279.875959 0.199245 0.194716\n", + "\t65 13236.160576 13264.863530 0.180157 0.180622\n", + "\t70 13220.204788 13336.570629 0.166780 0.162591\n", + "\t75 13206.871230 13264.422259 0.156329 0.154915\n", + "\t80 13191.695756 13249.784141 0.146271 0.144610\n", + "\t85 13132.044552 13218.036151 0.138585 0.135840\n", + "\t90 13107.559438 13177.846374 0.129874 0.127847\n", + "\t95 13096.184312 13175.012522 0.121443 0.119499\n", + "\t100 13051.565010 13135.169311 0.115060 0.112385\n", + "\t105 13026.202501 13096.897799 0.108914 0.106725\n", + "\t110 12971.433048 13083.418933 0.103877 0.100297\n", + "\t115 12951.898476 13056.777956 0.098483 0.095490\n", + "\t120 12938.304729 13032.014535 0.093136 0.090889\n", + "\t125 12899.761435 13018.849241 0.088597 0.086436\n", + "\t130 12882.151103 12981.466541 0.084103 0.082128\n", + "\t135 12882.960846 12920.917764 0.079350 0.078529\n", + "\t140 12869.549641 12909.860356 0.075580 0.074552\n", + "\t145 12842.281777 12851.799676 0.072636 0.072046\n", + "\t150 12749.157914 12782.647581 0.070611 0.069667\n", + "\t155 12703.241916 12738.044835 0.067873 0.066928\n", + "\t160 12670.297608 12706.411554 0.065557 0.064376\n", + "\t165 12636.895690 12647.891447 0.063296 0.062578\n", + "\t170 12597.918830 12625.602574 0.061002 0.059937\n", + "\t175 12550.742876 12555.971329 0.058575 0.058448\n", + "\t180 12528.304583 12552.544848 0.056474 0.055863\n", + "\t185 12475.332811 12498.211140 0.054925 0.054189\n", + "\t190 12397.080100 12456.284077 0.054007 0.052761\n", + "\t195 12372.673017 12363.777196 0.052291 0.051798\n", + "\t200 12355.614341 12365.014946 0.050307 0.049562\n", + "\t205 12286.492951 12267.941719 0.049312 0.048978\n", + "\t210 12227.262848 12211.404200 0.048146 0.047722\n", + "\t215 12110.655363 12156.535212 0.047822 0.046684\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 6/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.22714e+43 erg / s Luminosity absorbed = 8.09114e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13133.133359 13196.025024 0.661105 0.650910\n", + "\t5 13180.154270 13255.698653 0.561426 0.556263\n", + "\t10 13187.356843 13239.518470 0.499942 0.494814\n", + "\t15 13230.009983 13288.479637 0.443721 0.438216\n", + "\t20 13221.404988 13326.424072 0.401814 0.392288\n", + "\t25 13279.927924 13327.703183 0.358538 0.357912\n", + "\t30 13213.416642 13290.593772 0.331098 0.326841\n", + "\t35 13307.286771 13379.166431 0.293396 0.290233\n", + "\t40 13257.000678 13332.043792 0.272724 0.269385\n", + "\t45 13266.454328 13342.544624 0.249218 0.245700\n", + "\t50 13285.344929 13393.368804 0.228049 0.223222\n", + "\t55 13311.193287 13405.465136 0.209753 0.205613\n", + "\t60 13279.875959 13337.444936 0.194716 0.192957\n", + "\t65 13264.863530 13322.573830 0.180622 0.179129\n", + "\t70 13336.570629 13376.807762 0.162591 0.162838\n", + "\t75 13264.422259 13341.931624 0.154915 0.152609\n", + "\t80 13249.784141 13323.359370 0.144610 0.142282\n", + "\t85 13218.036151 13320.624255 0.135840 0.132022\n", + "\t90 13177.846374 13276.335581 0.127847 0.124791\n", + "\t95 13175.012522 13291.680639 0.119499 0.116002\n", + "\t100 13135.169311 13231.017853 0.112385 0.110371\n", + "\t105 13096.897799 13210.946440 0.106725 0.104111\n", + "\t110 13083.418933 13207.592577 0.100297 0.097619\n", + "\t115 13056.777956 13119.069100 0.095490 0.094113\n", + "\t120 13032.014535 13109.529071 0.090889 0.088884\n", + "\t125 13018.849241 13095.270885 0.086436 0.084279\n", + "\t130 12981.466541 13033.700377 0.082128 0.080103\n", + "\t135 12920.917764 12957.809294 0.078529 0.077480\n", + "\t140 12909.860356 12931.802369 0.074552 0.073862\n", + "\t145 12851.799676 12925.100614 0.072046 0.070224\n", + "\t150 12782.647581 12855.814961 0.069667 0.068059\n", + "\t155 12738.044835 12839.979057 0.066928 0.064786\n", + "\t160 12706.411554 12827.797549 0.064376 0.061897\n", + "\t165 12647.891447 12755.000233 0.062578 0.060413\n", + "\t170 12625.602574 12693.096249 0.059937 0.058603\n", + "\t175 12555.971329 12601.717780 0.058448 0.057507\n", + "\t180 12552.544848 12615.109130 0.055863 0.054529\n", + "\t185 12498.211140 12586.653068 0.054189 0.052663\n", + "\t190 12456.284077 12537.507236 0.052761 0.051090\n", + "\t195 12363.777196 12473.427562 0.051798 0.049878\n", + "\t200 12365.014946 12477.999050 0.049562 0.047867\n", + "\t205 12267.941719 12362.115770 0.048978 0.047396\n", + "\t210 12211.404200 12253.089156 0.047722 0.046878\n", + "\t215 12156.535212 12180.868013 0.046684 0.045927\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 7/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.22118e+43 erg / s Luminosity absorbed = 8.15225e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13196.025024 13210.192507 0.650910 0.651204\n", + "\t5 13255.698653 13279.272431 0.556263 0.553787\n", + "\t10 13239.518470 13325.730308 0.494814 0.485617\n", + "\t15 13288.479637 13305.186251 0.438216 0.438444\n", + "\t20 13326.424072 13343.413094 0.392288 0.392317\n", + "\t25 13327.703183 13329.023277 0.357912 0.358184\n", + "\t30 13290.593772 13327.705187 0.326841 0.323893\n", + "\t35 13379.166431 13366.676710 0.290233 0.291474\n", + "\t40 13332.043792 13347.637776 0.269385 0.270268\n", + "\t45 13342.544624 13391.584255 0.245700 0.243847\n", + "\t50 13393.368804 13371.818600 0.223222 0.224589\n", + "\t55 13405.465136 13429.111170 0.205613 0.205312\n", + "\t60 13337.444936 13385.068591 0.192957 0.190720\n", + "\t65 13322.573830 13377.164416 0.179129 0.177522\n", + "\t70 13376.807762 13389.265265 0.162838 0.162714\n", + "\t75 13341.931624 13335.296594 0.152609 0.152417\n", + "\t80 13323.359370 13317.898827 0.142282 0.142520\n", + "\t85 13320.624255 13328.550643 0.132022 0.131243\n", + "\t90 13276.335581 13310.694565 0.124791 0.123615\n", + "\t95 13291.680639 13291.914034 0.116002 0.115864\n", + "\t100 13231.017853 13270.160350 0.110371 0.108953\n", + "\t105 13210.946440 13206.626648 0.104111 0.103936\n", + "\t110 13207.592577 13204.813252 0.097619 0.096775\n", + "\t115 13119.069100 13133.099872 0.094113 0.092629\n", + "\t120 13109.529071 13152.191134 0.088884 0.087390\n", + "\t125 13095.270885 13119.431535 0.084279 0.082767\n", + "\t130 13033.700377 13054.659703 0.080103 0.079317\n", + "\t135 12957.809294 13030.381168 0.077480 0.075530\n", + "\t140 12931.802369 12994.145013 0.073862 0.072263\n", + "\t145 12925.100614 12958.166533 0.070224 0.069155\n", + "\t150 12855.814961 12883.286594 0.068059 0.066754\n", + "\t155 12839.979057 12863.007784 0.064786 0.064032\n", + "\t160 12827.797549 12827.137770 0.061897 0.061571\n", + "\t165 12755.000233 12799.473800 0.060413 0.059162\n", + "\t170 12693.096249 12717.620831 0.058603 0.057730\n", + "\t175 12601.717780 12641.557351 0.057507 0.056527\n", + "\t180 12615.109130 12613.094183 0.054529 0.054178\n", + "\t185 12586.653068 12581.699392 0.052663 0.052259\n", + "\t190 12537.507236 12556.326916 0.051090 0.050396\n", + "\t195 12473.427562 12463.328882 0.049878 0.049704\n", + "\t200 12477.999050 12421.543590 0.047867 0.048312\n", + "\t205 12362.115770 12339.759839 0.047396 0.047358\n", + "\t210 12253.089156 12243.798783 0.046878 0.046663\n", + "\t215 12180.868013 12163.326591 0.045927 0.045862\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 8/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.22375e+43 erg / s Luminosity absorbed = 8.12095e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13210.192507 13194.127608 0.651204 0.656434\n", + "\t5 13279.272431 13228.388565 0.553787 0.559574\n", + "\t10 13325.730308 13289.414493 0.485617 0.489004\n", + "\t15 13305.186251 13346.761126 0.438444 0.433787\n", + "\t20 13343.413094 13363.742432 0.392317 0.390643\n", + "\t25 13329.023277 13385.269116 0.358184 0.351372\n", + "\t30 13327.705187 13359.128491 0.323893 0.321469\n", + "\t35 13366.676710 13398.262188 0.291474 0.291795\n", + "\t40 13347.637776 13415.204376 0.270268 0.266293\n", + "\t45 13391.584255 13428.721895 0.243847 0.242985\n", + "\t50 13371.818600 13413.695497 0.224589 0.222309\n", + "\t55 13429.111170 13406.001338 0.205312 0.206976\n", + "\t60 13385.068591 13406.277223 0.190720 0.191089\n", + "\t65 13377.164416 13447.307416 0.177522 0.174859\n", + "\t70 13389.265265 13427.486749 0.162714 0.162477\n", + "\t75 13335.296594 13400.407701 0.152417 0.150796\n", + "\t80 13317.898827 13381.691224 0.142520 0.140777\n", + "\t85 13328.550643 13371.963378 0.131243 0.131026\n", + "\t90 13310.694565 13375.266319 0.123615 0.121836\n", + "\t95 13291.914034 13356.498239 0.115864 0.114540\n", + "\t100 13270.160350 13323.180618 0.108953 0.107497\n", + "\t105 13206.626648 13246.701015 0.103936 0.102790\n", + "\t110 13204.813252 13208.360067 0.096775 0.097150\n", + "\t115 13133.099872 13148.966352 0.092629 0.092266\n", + "\t120 13152.191134 13174.872225 0.087390 0.086782\n", + "\t125 13119.431535 13142.241261 0.082767 0.082443\n", + "\t130 13054.659703 13072.396994 0.079317 0.079388\n", + "\t135 13030.381168 13051.399307 0.075530 0.075497\n", + "\t140 12994.145013 13026.259101 0.072263 0.071752\n", + "\t145 12958.166533 12996.797161 0.069155 0.068470\n", + "\t150 12883.286594 12925.202949 0.066754 0.065936\n", + "\t155 12863.007784 12864.354520 0.064032 0.064026\n", + "\t160 12827.137770 12859.747285 0.061571 0.061331\n", + "\t165 12799.473800 12793.706958 0.059162 0.059341\n", + "\t170 12717.620831 12733.667948 0.057730 0.057447\n", + "\t175 12641.557351 12669.138594 0.056527 0.055941\n", + "\t180 12613.094183 12647.126180 0.054178 0.053803\n", + "\t185 12581.699392 12589.173650 0.052259 0.052224\n", + "\t190 12556.326916 12550.587983 0.050396 0.050632\n", + "\t195 12463.328882 12450.144160 0.049704 0.049861\n", + "\t200 12421.543590 12405.865846 0.048312 0.048493\n", + "\t205 12339.759839 12351.709722 0.047358 0.047307\n", + "\t210 12243.798783 12253.994038 0.046663 0.046529\n", + "\t215 12163.326591 12193.782185 0.045862 0.045482\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 9/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21908e+43 erg / s Luminosity absorbed = 8.16872e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13194.127608 13207.712484 0.656434 0.655960\n", + "\t5 13228.388565 13224.770429 0.559574 0.562303\n", + "\t10 13289.414493 13275.361500 0.489004 0.494601\n", + "\t15 13346.761126 13317.607311 0.433787 0.435487\n", + "\t20 13363.742432 13352.395556 0.390643 0.391052\n", + "\t25 13385.269116 13386.828235 0.351372 0.350598\n", + "\t30 13359.128491 13410.026678 0.321469 0.315679\n", + "\t35 13398.262188 13389.009476 0.291795 0.290711\n", + "\t40 13415.204376 13418.441455 0.266293 0.263041\n", + "\t45 13428.721895 13378.162118 0.242985 0.245205\n", + "\t50 13413.695497 13443.730262 0.222309 0.220039\n", + "\t55 13406.001338 13424.303697 0.206976 0.205223\n", + "\t60 13406.277223 13419.555276 0.191089 0.188897\n", + "\t65 13447.307416 13460.875005 0.174859 0.172920\n", + "\t70 13427.486749 13450.337756 0.162477 0.160245\n", + "\t75 13400.407701 13448.630047 0.150796 0.148399\n", + "\t80 13381.691224 13438.668160 0.140777 0.138444\n", + "\t85 13371.963378 13419.643621 0.131026 0.130018\n", + "\t90 13375.266319 13410.569390 0.121836 0.121143\n", + "\t95 13356.498239 13411.775288 0.114540 0.112150\n", + "\t100 13323.180618 13345.602702 0.107497 0.106578\n", + "\t105 13246.701015 13309.641634 0.102790 0.100505\n", + "\t110 13208.360067 13267.237025 0.097150 0.095508\n", + "\t115 13148.966352 13245.989409 0.092266 0.089846\n", + "\t120 13174.872225 13192.393305 0.086782 0.086280\n", + "\t125 13142.241261 13206.993275 0.082443 0.081211\n", + "\t130 13072.396994 13109.780991 0.079388 0.078693\n", + "\t135 13051.399307 13105.515227 0.075497 0.074463\n", + "\t140 13026.259101 13080.013791 0.071752 0.070828\n", + "\t145 12996.797161 13034.310835 0.068470 0.067873\n", + "\t150 12925.202949 12954.520922 0.065936 0.065855\n", + "\t155 12864.354520 12892.530978 0.064026 0.063870\n", + "\t160 12859.747285 12871.554047 0.061331 0.061258\n", + "\t165 12793.706958 12829.332539 0.059341 0.058871\n", + "\t170 12733.667948 12764.372956 0.057447 0.057038\n", + "\t175 12669.138594 12758.719713 0.055941 0.054577\n", + "\t180 12647.126180 12643.886903 0.053803 0.053704\n", + "\t185 12589.173650 12584.799996 0.052224 0.052303\n", + "\t190 12550.587983 12539.090757 0.050632 0.050683\n", + "\t195 12450.144160 12482.487445 0.049861 0.049185\n", + "\t200 12405.865846 12385.358724 0.048493 0.048675\n", + "\t205 12351.709722 12362.404518 0.047307 0.046954\n", + "\t210 12253.994038 12303.329057 0.046529 0.045802\n", + "\t215 12193.782185 12217.514244 0.045482 0.045041\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 10/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21621e+43 erg / s Luminosity absorbed = 8.18951e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13207.712484 13250.044474 0.655960 0.648306\n", + "\t5 13224.770429 13268.260885 0.562303 0.554688\n", + "\t10 13275.361500 13343.630262 0.494601 0.485191\n", + "\t15 13317.607311 13372.127852 0.435487 0.433539\n", + "\t20 13352.395556 13394.404128 0.391052 0.387859\n", + "\t25 13386.828235 13423.447190 0.350598 0.348268\n", + "\t30 13410.026678 13435.790768 0.315679 0.315451\n", + "\t35 13389.009476 13459.456710 0.290711 0.284964\n", + "\t40 13418.441455 13459.066834 0.263041 0.262896\n", + "\t45 13378.162118 13429.157522 0.245205 0.244167\n", + "\t50 13443.730262 13447.760292 0.220039 0.221943\n", + "\t55 13424.303697 13475.944820 0.205223 0.203639\n", + "\t60 13419.555276 13456.185170 0.188897 0.187987\n", + "\t65 13460.875005 13505.405683 0.172920 0.171399\n", + "\t70 13450.337756 13471.586824 0.160245 0.159227\n", + "\t75 13448.630047 13509.654307 0.148399 0.146613\n", + "\t80 13438.668160 13485.367813 0.138444 0.136661\n", + "\t85 13419.643621 13459.486003 0.130018 0.128842\n", + "\t90 13410.569390 13413.151984 0.121143 0.120975\n", + "\t95 13411.775288 13372.719193 0.112150 0.113967\n", + "\t100 13345.602702 13347.397283 0.106578 0.107160\n", + "\t105 13309.641634 13357.042867 0.100505 0.099694\n", + "\t110 13267.237025 13292.283955 0.095508 0.095070\n", + "\t115 13245.989409 13254.313335 0.089846 0.089742\n", + "\t120 13192.393305 13203.352329 0.086280 0.086462\n", + "\t125 13206.993275 13200.112218 0.081211 0.081222\n", + "\t130 13109.780991 13116.657467 0.078693 0.078214\n", + "\t135 13105.515227 13120.957986 0.074463 0.073713\n", + "\t140 13080.013791 13117.413156 0.070828 0.070116\n", + "\t145 13034.310835 13042.554199 0.067873 0.067352\n", + "\t150 12954.520922 12983.088194 0.065855 0.065417\n", + "\t155 12892.530978 12953.555946 0.063870 0.062668\n", + "\t160 12871.554047 12916.451771 0.061258 0.060032\n", + "\t165 12829.332539 12858.080133 0.058871 0.058206\n", + "\t170 12764.372956 12812.042920 0.057038 0.056128\n", + "\t175 12758.719713 12747.638758 0.054577 0.054594\n", + "\t180 12643.886903 12709.889498 0.053704 0.052654\n", + "\t185 12584.799996 12659.454233 0.052303 0.050867\n", + "\t190 12539.090757 12603.737771 0.050683 0.049698\n", + "\t195 12482.487445 12542.866664 0.049185 0.048387\n", + "\t200 12385.358724 12457.714659 0.048675 0.047626\n", + "\t205 12362.404518 12404.158500 0.046954 0.046263\n", + "\t210 12303.329057 12314.481431 0.045802 0.045468\n", + "\t215 12217.514244 12219.633343 0.045041 0.044806\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 11/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.22099e+43 erg / s Luminosity absorbed = 8.14722e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13250.044474 13230.303793 0.648306 0.653347\n", + "\t5 13268.260885 13278.921453 0.554688 0.557141\n", + "\t10 13343.630262 13297.867380 0.485191 0.494971\n", + "\t15 13372.127852 13338.977010 0.433539 0.437654\n", + "\t20 13394.404128 13388.041689 0.387859 0.389896\n", + "\t25 13423.447190 13399.607318 0.348268 0.352616\n", + "\t30 13435.790768 13402.658723 0.315451 0.318820\n", + "\t35 13459.456710 13450.995960 0.284964 0.286188\n", + "\t40 13459.066834 13448.503156 0.262896 0.263254\n", + "\t45 13429.157522 13465.392891 0.244167 0.240151\n", + "\t50 13447.760292 13467.648960 0.221943 0.221032\n", + "\t55 13475.944820 13475.701645 0.203639 0.202742\n", + "\t60 13456.185170 13493.350980 0.187987 0.185672\n", + "\t65 13505.405683 13518.310710 0.171399 0.170532\n", + "\t70 13471.586824 13494.025081 0.159227 0.158113\n", + "\t75 13509.654307 13539.078602 0.146613 0.144461\n", + "\t80 13485.367813 13548.423948 0.136661 0.133634\n", + "\t85 13459.486003 13467.039561 0.128842 0.127736\n", + "\t90 13413.151984 13451.638423 0.120975 0.119543\n", + "\t95 13372.719193 13457.248316 0.113967 0.111862\n", + "\t100 13347.397283 13412.863481 0.107160 0.105685\n", + "\t105 13357.042867 13398.455904 0.099694 0.098452\n", + "\t110 13292.283955 13302.921440 0.095070 0.094659\n", + "\t115 13254.313335 13282.193580 0.089742 0.089329\n", + "\t120 13203.352329 13219.401175 0.086462 0.086073\n", + "\t125 13200.112218 13201.420216 0.081222 0.081001\n", + "\t130 13116.657467 13139.211214 0.078214 0.077949\n", + "\t135 13120.957986 13164.116084 0.073713 0.073320\n", + "\t140 13117.413156 13143.426144 0.070116 0.069725\n", + "\t145 13042.554199 13072.451666 0.067352 0.067062\n", + "\t150 12983.088194 13012.224631 0.065417 0.064803\n", + "\t155 12953.555946 12948.682525 0.062668 0.062814\n", + "\t160 12916.451771 12932.493342 0.060032 0.059969\n", + "\t165 12858.080133 12837.412003 0.058206 0.058573\n", + "\t170 12812.042920 12822.790631 0.056128 0.056038\n", + "\t175 12747.638758 12754.419790 0.054594 0.054566\n", + "\t180 12709.889498 12690.188142 0.052654 0.052906\n", + "\t185 12659.454233 12657.419216 0.050867 0.051031\n", + "\t190 12603.737771 12567.241073 0.049698 0.050083\n", + "\t195 12542.866664 12506.494295 0.048387 0.048923\n", + "\t200 12457.714659 12458.719247 0.047626 0.047384\n", + "\t205 12404.158500 12402.111064 0.046263 0.046297\n", + "\t210 12314.481431 12309.834715 0.045468 0.045563\n", + "\t215 12219.633343 12225.965187 0.044806 0.044749\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 12/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21227e+43 erg / s Luminosity absorbed = 8.22264e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13230.303793 13208.849994 0.653347 0.659598\n", + "\t5 13278.921453 13272.158210 0.557141 0.561212\n", + "\t10 13297.867380 13342.748204 0.494971 0.488881\n", + "\t15 13338.977010 13390.263252 0.437654 0.432805\n", + "\t20 13388.041689 13424.191029 0.389896 0.388292\n", + "\t25 13399.607318 13458.184819 0.352616 0.348057\n", + "\t30 13402.658723 13470.170887 0.318820 0.314842\n", + "\t35 13450.995960 13486.075663 0.286188 0.284269\n", + "\t40 13448.503156 13522.925178 0.263254 0.259422\n", + "\t45 13465.392891 13497.508334 0.240151 0.239422\n", + "\t50 13467.648960 13489.284340 0.221032 0.220600\n", + "\t55 13475.701645 13512.115225 0.202742 0.202209\n", + "\t60 13493.350980 13496.883357 0.185672 0.186368\n", + "\t65 13518.310710 13519.225487 0.170532 0.172357\n", + "\t70 13494.025081 13459.411751 0.158113 0.160356\n", + "\t75 13539.078602 13486.183066 0.144461 0.147243\n", + "\t80 13548.423948 13486.241355 0.133634 0.136797\n", + "\t85 13467.039561 13451.472664 0.127736 0.128289\n", + "\t90 13451.638423 13456.663239 0.119543 0.119819\n", + "\t95 13457.248316 13468.237878 0.111862 0.111512\n", + "\t100 13412.863481 13417.948513 0.105685 0.105811\n", + "\t105 13398.455904 13361.202313 0.098452 0.099821\n", + "\t110 13302.921440 13316.056762 0.094659 0.094989\n", + "\t115 13282.193580 13243.700355 0.089329 0.090345\n", + "\t120 13219.401175 13243.163337 0.086073 0.085463\n", + "\t125 13201.420216 13250.057352 0.081001 0.080080\n", + "\t130 13139.211214 13164.311347 0.077949 0.077706\n", + "\t135 13164.116084 13184.287993 0.073320 0.072553\n", + "\t140 13143.426144 13154.157361 0.069725 0.069428\n", + "\t145 13072.451666 13079.781830 0.067062 0.067008\n", + "\t150 13012.224631 13030.156864 0.064803 0.064555\n", + "\t155 12948.682525 12944.264398 0.062814 0.062676\n", + "\t160 12932.493342 12949.815017 0.059969 0.059988\n", + "\t165 12837.412003 12897.310883 0.058573 0.057579\n", + "\t170 12822.790631 12922.159742 0.056038 0.054708\n", + "\t175 12754.419790 12821.949007 0.054566 0.053511\n", + "\t180 12690.188142 12730.870651 0.052906 0.052254\n", + "\t185 12657.419216 12722.247624 0.051031 0.050018\n", + "\t190 12567.241073 12642.740779 0.050083 0.048842\n", + "\t195 12506.494295 12560.850177 0.048923 0.048038\n", + "\t200 12458.719247 12511.262912 0.047384 0.046644\n", + "\t205 12402.111064 12421.083202 0.046297 0.045961\n", + "\t210 12309.834715 12344.680043 0.045563 0.044970\n", + "\t215 12225.965187 12226.871851 0.044749 0.044510\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 13/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21689e+43 erg / s Luminosity absorbed = 8.18137e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13208.849994 13221.314409 0.659598 0.656564\n", + "\t5 13272.158210 13266.903886 0.561212 0.559145\n", + "\t10 13342.748204 13343.312042 0.488881 0.485847\n", + "\t15 13390.263252 13396.321205 0.432805 0.432258\n", + "\t20 13424.191029 13431.814773 0.388292 0.387402\n", + "\t25 13458.184819 13514.656129 0.348057 0.342639\n", + "\t30 13470.170887 13503.249569 0.314842 0.311340\n", + "\t35 13486.075663 13497.806703 0.284269 0.285753\n", + "\t40 13522.925178 13536.387470 0.259422 0.259071\n", + "\t45 13497.508334 13526.018045 0.239422 0.237719\n", + "\t50 13489.284340 13479.479891 0.220600 0.222196\n", + "\t55 13512.115225 13531.789695 0.202209 0.199924\n", + "\t60 13496.883357 13524.251543 0.186368 0.185102\n", + "\t65 13519.225487 13531.626249 0.172357 0.170451\n", + "\t70 13459.411751 13491.592574 0.160356 0.159119\n", + "\t75 13486.183066 13491.954060 0.147243 0.148060\n", + "\t80 13486.241355 13493.132689 0.136797 0.136837\n", + "\t85 13451.472664 13485.073739 0.128289 0.127372\n", + "\t90 13456.663239 13466.187776 0.119819 0.119361\n", + "\t95 13468.237878 13451.422770 0.111512 0.111297\n", + "\t100 13417.948513 13454.427565 0.105811 0.103878\n", + "\t105 13361.202313 13386.326056 0.099821 0.099301\n", + "\t110 13316.056762 13366.800827 0.094989 0.093696\n", + "\t115 13243.700355 13306.582006 0.090345 0.088835\n", + "\t120 13243.163337 13245.016932 0.085463 0.084972\n", + "\t125 13250.057352 13255.717688 0.080080 0.080067\n", + "\t130 13164.311347 13246.575711 0.077706 0.075900\n", + "\t135 13184.287993 13181.247980 0.072553 0.072767\n", + "\t140 13154.157361 13122.812683 0.069428 0.069772\n", + "\t145 13079.781830 13037.707680 0.067008 0.068075\n", + "\t150 13030.156864 13026.386344 0.064555 0.064898\n", + "\t155 12944.264398 12925.145144 0.062676 0.063013\n", + "\t160 12949.815017 12885.376665 0.059988 0.060850\n", + "\t165 12897.310883 12861.447381 0.057579 0.058018\n", + "\t170 12922.159742 12833.429801 0.054708 0.056175\n", + "\t175 12821.949007 12783.510362 0.053511 0.054147\n", + "\t180 12730.870651 12717.340806 0.052254 0.052455\n", + "\t185 12722.247624 12680.884053 0.050018 0.050573\n", + "\t190 12642.740779 12654.520515 0.048842 0.048821\n", + "\t195 12560.850177 12556.882955 0.048038 0.048179\n", + "\t200 12511.262912 12509.126284 0.046644 0.046719\n", + "\t205 12421.083202 12446.538483 0.045961 0.045501\n", + "\t210 12344.680043 12344.694905 0.044970 0.044970\n", + "\t215 12226.871851 12231.908348 0.044510 0.044500\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 14/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21370e+43 erg / s Luminosity absorbed = 8.21618e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13221.314409 13248.085322 0.656564 0.652078\n", + "\t5 13266.903886 13271.262850 0.559145 0.557217\n", + "\t10 13343.312042 13323.100432 0.485847 0.488675\n", + "\t15 13396.321205 13364.000693 0.432258 0.434756\n", + "\t20 13431.814773 13415.622861 0.387402 0.388410\n", + "\t25 13514.656129 13519.302099 0.342639 0.342131\n", + "\t30 13503.249569 13522.659188 0.311340 0.309502\n", + "\t35 13497.806703 13549.670505 0.285753 0.280875\n", + "\t40 13536.387470 13564.411362 0.259071 0.258257\n", + "\t45 13526.018045 13547.677730 0.237719 0.236120\n", + "\t50 13479.479891 13531.242465 0.222196 0.217546\n", + "\t55 13531.789695 13562.478514 0.199924 0.199222\n", + "\t60 13524.251543 13560.618390 0.185102 0.183767\n", + "\t65 13531.626249 13595.350076 0.170451 0.168717\n", + "\t70 13491.592574 13583.685267 0.159119 0.155018\n", + "\t75 13491.954060 13568.688277 0.148060 0.144394\n", + "\t80 13493.132689 13511.447091 0.136837 0.136209\n", + "\t85 13485.073739 13524.307991 0.127372 0.126028\n", + "\t90 13466.187776 13518.242204 0.119361 0.116933\n", + "\t95 13451.422770 13473.696003 0.111297 0.110181\n", + "\t100 13454.427565 13470.970672 0.103878 0.103747\n", + "\t105 13386.326056 13439.625097 0.099301 0.097789\n", + "\t110 13366.800827 13341.589885 0.093696 0.094411\n", + "\t115 13306.582006 13284.506102 0.088835 0.088979\n", + "\t120 13245.016932 13212.934919 0.084972 0.085421\n", + "\t125 13255.717688 13236.311034 0.080067 0.079694\n", + "\t130 13246.575711 13209.300239 0.075900 0.076463\n", + "\t135 13181.247980 13147.404552 0.072767 0.072989\n", + "\t140 13122.812683 13123.191470 0.069772 0.069362\n", + "\t145 13037.707680 13062.442351 0.068075 0.067216\n", + "\t150 13026.386344 13002.501066 0.064898 0.064495\n", + "\t155 12925.145144 12941.430006 0.063013 0.062345\n", + "\t160 12885.376665 12885.248161 0.060850 0.060475\n", + "\t165 12861.447381 12856.926441 0.058018 0.057861\n", + "\t170 12833.429801 12831.914851 0.056175 0.055705\n", + "\t175 12783.510362 12783.093570 0.054147 0.053784\n", + "\t180 12717.340806 12690.757185 0.052455 0.052567\n", + "\t185 12680.884053 12641.084383 0.050573 0.050898\n", + "\t190 12654.520515 12572.460696 0.048821 0.049764\n", + "\t195 12556.882955 12516.398175 0.048179 0.048572\n", + "\t200 12509.126284 12454.357431 0.046719 0.047301\n", + "\t205 12446.538483 12367.350498 0.045501 0.046470\n", + "\t210 12344.694905 12318.940270 0.044970 0.045142\n", + "\t215 12231.908348 12221.586509 0.044500 0.044518\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 15/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.20596e+43 erg / s Luminosity absorbed = 8.28887e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13248.085322 13269.950479 0.652078 0.645755\n", + "\t5 13271.262850 13301.816652 0.557217 0.549777\n", + "\t10 13323.100432 13356.938920 0.488675 0.484341\n", + "\t15 13364.000693 13387.126609 0.434756 0.432017\n", + "\t20 13415.622861 13474.324367 0.388410 0.380816\n", + "\t25 13519.302099 13530.987543 0.342131 0.342549\n", + "\t30 13522.659188 13531.828247 0.309502 0.310338\n", + "\t35 13549.670505 13545.706924 0.280875 0.281759\n", + "\t40 13564.411362 13546.619910 0.258257 0.258598\n", + "\t45 13547.677730 13498.710008 0.236120 0.240528\n", + "\t50 13531.242465 13568.778885 0.217546 0.215732\n", + "\t55 13562.478514 13606.052514 0.199222 0.195722\n", + "\t60 13560.618390 13624.818819 0.183767 0.178910\n", + "\t65 13595.350076 13619.048574 0.168717 0.166233\n", + "\t70 13583.685267 13592.981823 0.155018 0.154182\n", + "\t75 13568.688277 13620.893813 0.144394 0.141796\n", + "\t80 13511.447091 13597.052224 0.136209 0.133395\n", + "\t85 13524.307991 13585.058724 0.126028 0.123785\n", + "\t90 13518.242204 13591.535826 0.116933 0.114551\n", + "\t95 13473.696003 13520.772286 0.110181 0.108313\n", + "\t100 13470.970672 13457.777873 0.103747 0.103595\n", + "\t105 13439.625097 13454.673879 0.097789 0.097445\n", + "\t110 13341.589885 13422.703934 0.094411 0.091966\n", + "\t115 13284.506102 13389.027818 0.088979 0.086189\n", + "\t120 13212.934919 13309.391317 0.085421 0.083230\n", + "\t125 13236.311034 13240.145769 0.079694 0.079493\n", + "\t130 13209.300239 13204.256727 0.076463 0.076387\n", + "\t135 13147.404552 13181.230898 0.072989 0.072152\n", + "\t140 13123.191470 13140.558809 0.069362 0.068672\n", + "\t145 13062.442351 13100.506478 0.067216 0.066137\n", + "\t150 13002.501066 13031.203765 0.064495 0.063556\n", + "\t155 12941.430006 13001.184124 0.062345 0.060949\n", + "\t160 12885.248161 12942.508886 0.060475 0.059074\n", + "\t165 12856.926441 12912.024518 0.057861 0.056502\n", + "\t170 12831.914851 12891.299428 0.055705 0.054352\n", + "\t175 12783.093570 12829.756525 0.053784 0.052885\n", + "\t180 12690.757185 12724.080499 0.052567 0.051865\n", + "\t185 12641.084383 12693.543905 0.050898 0.050000\n", + "\t190 12572.460696 12633.066604 0.049764 0.048720\n", + "\t195 12516.398175 12597.395784 0.048572 0.047184\n", + "\t200 12454.357431 12524.995233 0.047301 0.046111\n", + "\t205 12367.350498 12411.617488 0.046470 0.045838\n", + "\t210 12318.940270 12367.171116 0.045142 0.044366\n", + "\t215 12221.586509 12267.188024 0.044518 0.043741\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 16/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21019e+43 erg / s Luminosity absorbed = 8.24203e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13269.950479 13294.332464 0.645755 0.639894\n", + "\t5 13301.816652 13334.225375 0.549777 0.545794\n", + "\t10 13356.938920 13341.731754 0.484341 0.489215\n", + "\t15 13387.126609 13392.884162 0.432017 0.433385\n", + "\t20 13474.324367 13449.916507 0.380816 0.384548\n", + "\t25 13530.987543 13504.961340 0.342549 0.345200\n", + "\t30 13531.828247 13515.546456 0.310338 0.311512\n", + "\t35 13545.706924 13519.588682 0.281759 0.284215\n", + "\t40 13546.619910 13540.852305 0.258598 0.259600\n", + "\t45 13498.710008 13551.504707 0.240528 0.237025\n", + "\t50 13568.778885 13595.189072 0.215732 0.214145\n", + "\t55 13606.052514 13640.400297 0.195722 0.194052\n", + "\t60 13624.818819 13650.180777 0.178910 0.178014\n", + "\t65 13619.048574 13682.597008 0.166233 0.163665\n", + "\t70 13592.981823 13628.687357 0.154182 0.153787\n", + "\t75 13620.893813 13686.415370 0.141796 0.140718\n", + "\t80 13597.052224 13667.638767 0.133395 0.131648\n", + "\t85 13585.058724 13614.508146 0.123785 0.123306\n", + "\t90 13591.535826 13592.293834 0.114551 0.115055\n", + "\t95 13520.772286 13591.458834 0.108313 0.106462\n", + "\t100 13457.777873 13516.525837 0.103595 0.102804\n", + "\t105 13454.673879 13493.432984 0.097445 0.097303\n", + "\t110 13422.703934 13506.262863 0.091966 0.089921\n", + "\t115 13389.027818 13462.951218 0.086189 0.084904\n", + "\t120 13309.391317 13401.175756 0.083230 0.081100\n", + "\t125 13240.145769 13329.876751 0.079493 0.077521\n", + "\t130 13204.256727 13274.552132 0.076387 0.074636\n", + "\t135 13181.230898 13253.288405 0.072152 0.070750\n", + "\t140 13140.558809 13181.363752 0.068672 0.068530\n", + "\t145 13100.506478 13181.575930 0.066137 0.064852\n", + "\t150 13031.203765 13082.422101 0.063556 0.062904\n", + "\t155 13001.184124 13062.551665 0.060949 0.060169\n", + "\t160 12942.508886 12987.751942 0.059074 0.058697\n", + "\t165 12912.024518 12945.550201 0.056502 0.056124\n", + "\t170 12891.299428 12940.320268 0.054352 0.053558\n", + "\t175 12829.756525 12896.507466 0.052885 0.051976\n", + "\t180 12724.080499 12797.825337 0.051865 0.050767\n", + "\t185 12693.543905 12735.338727 0.050000 0.049307\n", + "\t190 12633.066604 12652.175435 0.048720 0.048346\n", + "\t195 12597.395784 12589.600635 0.047184 0.047215\n", + "\t200 12524.995233 12547.627424 0.046111 0.045773\n", + "\t205 12411.617488 12449.180372 0.045838 0.045101\n", + "\t210 12367.171116 12369.571510 0.044366 0.044203\n", + "\t215 12267.188024 12289.644346 0.043741 0.043363\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 17/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.20561e+43 erg / s Luminosity absorbed = 8.28570e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13294.332464 13244.460538 0.639894 0.649708\n", + "\t5 13334.225375 13313.132071 0.545794 0.550698\n", + "\t10 13341.731754 13348.350349 0.489215 0.487904\n", + "\t15 13392.884162 13404.499117 0.433385 0.433249\n", + "\t20 13449.916507 13430.748739 0.384548 0.385039\n", + "\t25 13504.961340 13499.884145 0.345200 0.346484\n", + "\t30 13515.546456 13504.716737 0.311512 0.313246\n", + "\t35 13519.588682 13513.576476 0.284215 0.284651\n", + "\t40 13540.852305 13560.463567 0.259600 0.257079\n", + "\t45 13551.504707 13582.157063 0.237025 0.233521\n", + "\t50 13595.189072 13587.075628 0.214145 0.214888\n", + "\t55 13640.400297 13627.976702 0.194052 0.195310\n", + "\t60 13650.180777 13631.544645 0.178014 0.179027\n", + "\t65 13682.597008 13648.868190 0.163665 0.164849\n", + "\t70 13628.687357 13643.698829 0.153787 0.152570\n", + "\t75 13686.415370 13645.078826 0.140718 0.141862\n", + "\t80 13667.638767 13627.262692 0.131648 0.132550\n", + "\t85 13614.508146 13566.681998 0.123306 0.124692\n", + "\t90 13592.293834 13530.395581 0.115055 0.117295\n", + "\t95 13591.458834 13521.840968 0.106462 0.109117\n", + "\t100 13516.525837 13506.642129 0.102804 0.103689\n", + "\t105 13493.432984 13492.251957 0.097303 0.097032\n", + "\t110 13506.262863 13451.600943 0.089921 0.091607\n", + "\t115 13462.951218 13436.119810 0.084904 0.085909\n", + "\t120 13401.175756 13371.241690 0.081100 0.081890\n", + "\t125 13329.876751 13311.937285 0.077521 0.078795\n", + "\t130 13274.552132 13252.648553 0.074636 0.075308\n", + "\t135 13253.288405 13169.998818 0.070750 0.072800\n", + "\t140 13181.363752 13174.222148 0.068530 0.068996\n", + "\t145 13181.575930 13153.212857 0.064852 0.065518\n", + "\t150 13082.422101 13095.618126 0.062904 0.062932\n", + "\t155 13062.551665 13094.553925 0.060169 0.059735\n", + "\t160 12987.751942 13003.530962 0.058697 0.058432\n", + "\t165 12945.550201 12966.608407 0.056124 0.056056\n", + "\t170 12940.320268 12949.402109 0.053558 0.053719\n", + "\t175 12896.507466 12862.494358 0.051976 0.052531\n", + "\t180 12797.825337 12769.340474 0.050767 0.051329\n", + "\t185 12735.338727 12709.306732 0.049307 0.049652\n", + "\t190 12652.175435 12652.964504 0.048346 0.048401\n", + "\t195 12589.600635 12584.816637 0.047215 0.047114\n", + "\t200 12547.627424 12545.414915 0.045773 0.045822\n", + "\t205 12449.180372 12449.480425 0.045101 0.045031\n", + "\t210 12369.571510 12384.011229 0.044203 0.044008\n", + "\t215 12289.644346 12284.916605 0.043363 0.043442\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 18/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21325e+43 erg / s Luminosity absorbed = 8.21131e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13244.460538 13215.252735 0.649708 0.654200\n", + "\t5 13313.132071 13331.428875 0.550698 0.548599\n", + "\t10 13348.350349 13337.370576 0.487904 0.489688\n", + "\t15 13404.499117 13409.016011 0.433249 0.430614\n", + "\t20 13430.748739 13410.672754 0.385039 0.388442\n", + "\t25 13499.884145 13501.090419 0.346484 0.345055\n", + "\t30 13504.716737 13505.531068 0.313246 0.314231\n", + "\t35 13513.576476 13530.755819 0.284651 0.284280\n", + "\t40 13560.463567 13609.596856 0.257079 0.253643\n", + "\t45 13582.157063 13618.753022 0.233521 0.232414\n", + "\t50 13587.075628 13597.159113 0.214888 0.214130\n", + "\t55 13627.976702 13610.530973 0.195310 0.197545\n", + "\t60 13631.544645 13623.810035 0.179027 0.180485\n", + "\t65 13648.868190 13661.772054 0.164849 0.165192\n", + "\t70 13643.698829 13641.503930 0.152570 0.153796\n", + "\t75 13645.078826 13645.600774 0.141862 0.142685\n", + "\t80 13627.262692 13622.954322 0.132550 0.132334\n", + "\t85 13566.681998 13639.296322 0.124692 0.122790\n", + "\t90 13530.395581 13547.268529 0.117295 0.116846\n", + "\t95 13521.840968 13569.531356 0.109117 0.108375\n", + "\t100 13506.642129 13515.502271 0.103689 0.103755\n", + "\t105 13492.251957 13508.172150 0.097032 0.096633\n", + "\t110 13451.600943 13416.524575 0.091607 0.092537\n", + "\t115 13436.119810 13426.995476 0.085909 0.086135\n", + "\t120 13371.241690 13365.275241 0.081890 0.082669\n", + "\t125 13311.937285 13358.631205 0.078795 0.077631\n", + "\t130 13252.648553 13268.178553 0.075308 0.075116\n", + "\t135 13169.998818 13223.649500 0.072800 0.071572\n", + "\t140 13174.222148 13220.602955 0.068996 0.068243\n", + "\t145 13153.212857 13164.065526 0.065518 0.065354\n", + "\t150 13095.618126 13095.126934 0.062932 0.062998\n", + "\t155 13094.553925 13063.058963 0.059735 0.060406\n", + "\t160 13003.530962 13006.990221 0.058432 0.058529\n", + "\t165 12966.608407 12986.943190 0.056056 0.056053\n", + "\t170 12949.402109 12932.955739 0.053719 0.054015\n", + "\t175 12862.494358 12870.860913 0.052531 0.052642\n", + "\t180 12769.340474 12790.514749 0.051329 0.051470\n", + "\t185 12709.306732 12728.490336 0.049652 0.049786\n", + "\t190 12652.964504 12670.148291 0.048401 0.048367\n", + "\t195 12584.816637 12570.559127 0.047114 0.047467\n", + "\t200 12545.414915 12551.141851 0.045822 0.045916\n", + "\t205 12449.480425 12451.950510 0.045031 0.045138\n", + "\t210 12384.011229 12423.828932 0.044008 0.043706\n", + "\t215 12284.916605 12275.821494 0.043442 0.043662\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 19/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.20640e+43 erg / s Luminosity absorbed = 8.27703e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13215.252735 13275.593833 0.654200 0.641575\n", + "\t5 13331.428875 13320.464125 0.548599 0.550550\n", + "\t10 13337.370576 13370.286171 0.489688 0.486018\n", + "\t15 13409.016011 13423.974251 0.430614 0.431285\n", + "\t20 13410.672754 13405.029898 0.388442 0.388908\n", + "\t25 13501.090419 13480.899554 0.345055 0.347041\n", + "\t30 13505.531068 13493.738929 0.314231 0.313299\n", + "\t35 13530.755819 13526.947308 0.284280 0.285034\n", + "\t40 13609.596856 13573.130852 0.253643 0.256411\n", + "\t45 13618.753022 13564.868422 0.232414 0.235894\n", + "\t50 13597.159113 13596.162820 0.214130 0.215532\n", + "\t55 13610.530973 13606.961834 0.197545 0.197586\n", + "\t60 13623.810035 13613.799866 0.180485 0.180924\n", + "\t65 13661.772054 13670.974200 0.165192 0.165513\n", + "\t70 13641.503930 13609.754866 0.153796 0.155721\n", + "\t75 13645.600774 13620.099801 0.142685 0.142833\n", + "\t80 13622.954322 13622.149138 0.132334 0.131837\n", + "\t85 13639.296322 13622.629367 0.122790 0.123216\n", + "\t90 13547.268529 13550.206611 0.116846 0.116628\n", + "\t95 13569.531356 13576.442824 0.108375 0.108370\n", + "\t100 13515.502271 13550.661291 0.103755 0.101730\n", + "\t105 13508.172150 13553.415275 0.096633 0.095153\n", + "\t110 13416.524575 13440.968288 0.092537 0.092355\n", + "\t115 13426.995476 13449.937711 0.086135 0.085928\n", + "\t120 13365.275241 13406.827159 0.082669 0.081735\n", + "\t125 13358.631205 13387.097490 0.077631 0.076855\n", + "\t130 13268.178553 13293.561247 0.075116 0.074028\n", + "\t135 13223.649500 13246.846989 0.071572 0.071086\n", + "\t140 13220.602955 13249.313571 0.068243 0.067556\n", + "\t145 13164.065526 13155.060691 0.065354 0.065304\n", + "\t150 13095.126934 13129.179399 0.062998 0.062308\n", + "\t155 13063.058963 13105.788762 0.060406 0.059735\n", + "\t160 13006.990221 13052.993533 0.058529 0.057324\n", + "\t165 12986.943190 13049.537481 0.056053 0.054732\n", + "\t170 12932.955739 12964.141616 0.054015 0.053309\n", + "\t175 12870.860913 12873.365866 0.052642 0.052182\n", + "\t180 12790.514749 12837.826363 0.051470 0.050638\n", + "\t185 12728.490336 12704.945462 0.049786 0.050011\n", + "\t190 12670.148291 12654.429611 0.048367 0.048416\n", + "\t195 12570.559127 12591.048933 0.047467 0.047031\n", + "\t200 12551.141851 12538.868881 0.045916 0.045828\n", + "\t205 12451.950510 12478.435159 0.045138 0.044770\n", + "\t210 12423.828932 12404.402526 0.043706 0.043784\n", + "\t215 12275.821494 12299.938075 0.043662 0.043250\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 20/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.20984e+43 erg / s Luminosity absorbed = 8.24688e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13275.593833 13272.361297 0.641575 0.643131\n", + "\t5 13320.464125 13304.764532 0.550550 0.550883\n", + "\t10 13370.286171 13353.155833 0.486018 0.485906\n", + "\t15 13423.974251 13390.303881 0.431285 0.431612\n", + "\t20 13405.029898 13414.224321 0.388908 0.389333\n", + "\t25 13480.899554 13455.091925 0.347041 0.351015\n", + "\t30 13493.738929 13513.943048 0.313299 0.313074\n", + "\t35 13526.947308 13521.764429 0.285034 0.284974\n", + "\t40 13573.130852 13574.404278 0.256411 0.258177\n", + "\t45 13564.868422 13603.410987 0.235894 0.232270\n", + "\t50 13596.162820 13612.094041 0.215532 0.213702\n", + "\t55 13606.961834 13617.408271 0.197586 0.197195\n", + "\t60 13613.799866 13610.013488 0.180924 0.180863\n", + "\t65 13670.974200 13673.552054 0.165513 0.164200\n", + "\t70 13609.754866 13631.831420 0.155721 0.153669\n", + "\t75 13620.099801 13628.628207 0.142833 0.142200\n", + "\t80 13622.149138 13616.815750 0.131837 0.132482\n", + "\t85 13622.629367 13599.203244 0.123216 0.124053\n", + "\t90 13550.206611 13557.210091 0.116628 0.116238\n", + "\t95 13576.442824 13594.203419 0.108370 0.107687\n", + "\t100 13550.661291 13555.614817 0.101730 0.101546\n", + "\t105 13553.415275 13536.865702 0.095153 0.095418\n", + "\t110 13440.968288 13441.150339 0.092355 0.092386\n", + "\t115 13449.937711 13445.885394 0.085928 0.086177\n", + "\t120 13406.827159 13382.034029 0.081735 0.082523\n", + "\t125 13387.097490 13356.909386 0.076855 0.077689\n", + "\t130 13293.561247 13315.234449 0.074028 0.074039\n", + "\t135 13246.846989 13228.514457 0.071086 0.071362\n", + "\t140 13249.313571 13218.647649 0.067556 0.067926\n", + "\t145 13155.060691 13184.008520 0.065304 0.064832\n", + "\t150 13129.179399 13121.121278 0.062308 0.062483\n", + "\t155 13105.788762 13076.788591 0.059735 0.060584\n", + "\t160 13052.993533 13054.043317 0.057324 0.057635\n", + "\t165 13049.537481 13019.828137 0.054732 0.055219\n", + "\t170 12964.141616 12904.410073 0.053309 0.054136\n", + "\t175 12873.365866 12845.416495 0.052182 0.052463\n", + "\t180 12837.826363 12800.187003 0.050638 0.050817\n", + "\t185 12704.945462 12730.751566 0.050011 0.049486\n", + "\t190 12654.429611 12717.490013 0.048416 0.047431\n", + "\t195 12591.048933 12616.635162 0.047031 0.046664\n", + "\t200 12538.868881 12576.252299 0.045828 0.045482\n", + "\t205 12478.435159 12495.955151 0.044770 0.044613\n", + "\t210 12404.402526 12378.250834 0.043784 0.044170\n", + "\t215 12299.938075 12289.662367 0.043250 0.043437\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 21/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21090e+43 erg / s Luminosity absorbed = 8.23618e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13272.361297 13238.110707 0.643131 0.650874\n", + "\t5 13304.764532 13287.753370 0.550883 0.555871\n", + "\t10 13353.155833 13374.524718 0.485906 0.483023\n", + "\t15 13390.303881 13434.237643 0.431612 0.425922\n", + "\t20 13414.224321 13473.045124 0.389333 0.380603\n", + "\t25 13455.091925 13494.400242 0.351015 0.346455\n", + "\t30 13513.943048 13534.656312 0.313074 0.311422\n", + "\t35 13521.764429 13561.340870 0.284974 0.283608\n", + "\t40 13574.404278 13602.879384 0.258177 0.256415\n", + "\t45 13603.410987 13634.222625 0.232270 0.231800\n", + "\t50 13612.094041 13652.488522 0.213702 0.211430\n", + "\t55 13617.408271 13599.228313 0.197195 0.198330\n", + "\t60 13610.013488 13628.648194 0.180863 0.180805\n", + "\t65 13673.552054 13670.725712 0.164200 0.164859\n", + "\t70 13631.831420 13654.266746 0.153669 0.152613\n", + "\t75 13628.628207 13629.745042 0.142200 0.142187\n", + "\t80 13616.815750 13665.134288 0.132482 0.130750\n", + "\t85 13599.203244 13630.062407 0.124053 0.123503\n", + "\t90 13557.210091 13549.894759 0.116238 0.116418\n", + "\t95 13594.203419 13581.793685 0.107687 0.107581\n", + "\t100 13555.614817 13535.666076 0.101546 0.102364\n", + "\t105 13536.865702 13525.411134 0.095418 0.095427\n", + "\t110 13441.150339 13411.374333 0.092386 0.092997\n", + "\t115 13445.885394 13436.288355 0.086177 0.086271\n", + "\t120 13382.034029 13390.614836 0.082523 0.081941\n", + "\t125 13356.909386 13362.603164 0.077689 0.077358\n", + "\t130 13315.234449 13301.473749 0.074039 0.074348\n", + "\t135 13228.514457 13241.489428 0.071362 0.071113\n", + "\t140 13218.647649 13233.411060 0.067926 0.067449\n", + "\t145 13184.008520 13191.107749 0.064832 0.064329\n", + "\t150 13121.121278 13049.165675 0.062483 0.063520\n", + "\t155 13076.788591 13051.818930 0.060584 0.060553\n", + "\t160 13054.043317 12981.109020 0.057635 0.058468\n", + "\t165 13019.828137 12939.117333 0.055219 0.056205\n", + "\t170 12904.410073 12872.061882 0.054136 0.054645\n", + "\t175 12845.416495 12826.919565 0.052463 0.052669\n", + "\t180 12800.187003 12774.051094 0.050817 0.051284\n", + "\t185 12730.751566 12712.958930 0.049486 0.049636\n", + "\t190 12717.490013 12723.636192 0.047431 0.047333\n", + "\t195 12616.635162 12607.550913 0.046664 0.046856\n", + "\t200 12576.252299 12557.223085 0.045482 0.045825\n", + "\t205 12495.955151 12496.316576 0.044613 0.044607\n", + "\t210 12378.250834 12359.901379 0.044170 0.044432\n", + "\t215 12289.662367 12291.488426 0.043437 0.043549\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 22/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.20712e+43 erg / s Luminosity absorbed = 8.27261e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13238.110707 13242.589560 0.650874 0.650661\n", + "\t5 13287.753370 13302.236433 0.555871 0.553883\n", + "\t10 13374.524718 13384.295437 0.483023 0.482708\n", + "\t15 13434.237643 13461.126863 0.425922 0.422022\n", + "\t20 13473.045124 13470.826671 0.380603 0.381236\n", + "\t25 13494.400242 13484.247273 0.346455 0.347424\n", + "\t30 13534.656312 13508.470387 0.311422 0.313096\n", + "\t35 13561.340870 13543.424565 0.283608 0.283193\n", + "\t40 13602.879384 13580.176047 0.256415 0.255980\n", + "\t45 13634.222625 13614.362136 0.231800 0.233483\n", + "\t50 13652.488522 13644.537532 0.211430 0.212147\n", + "\t55 13599.228313 13614.396375 0.198330 0.196233\n", + "\t60 13628.648194 13645.330715 0.180805 0.179325\n", + "\t65 13670.725712 13671.254373 0.164859 0.164136\n", + "\t70 13654.266746 13646.107424 0.152613 0.154460\n", + "\t75 13629.745042 13616.420307 0.142187 0.143180\n", + "\t80 13665.134288 13666.532072 0.130750 0.130676\n", + "\t85 13630.062407 13575.991311 0.123503 0.125264\n", + "\t90 13549.894759 13537.593047 0.116418 0.116523\n", + "\t95 13581.793685 13519.324997 0.107581 0.109037\n", + "\t100 13535.666076 13533.503766 0.102364 0.102218\n", + "\t105 13525.411134 13524.296286 0.095427 0.095472\n", + "\t110 13411.374333 13408.597159 0.092997 0.093029\n", + "\t115 13436.288355 13433.139621 0.086271 0.086510\n", + "\t120 13390.614836 13346.445735 0.081941 0.082582\n", + "\t125 13362.603164 13314.078372 0.077358 0.078528\n", + "\t130 13301.473749 13292.870438 0.074348 0.074822\n", + "\t135 13241.489428 13212.983843 0.071113 0.071947\n", + "\t140 13233.411060 13187.883546 0.067449 0.068256\n", + "\t145 13191.107749 13154.962528 0.064329 0.064997\n", + "\t150 13049.165675 13089.139849 0.063520 0.063208\n", + "\t155 13051.818930 13028.832892 0.060553 0.061212\n", + "\t160 12981.109020 12955.890317 0.058468 0.058877\n", + "\t165 12939.117333 12919.786297 0.056205 0.056736\n", + "\t170 12872.061882 12831.233593 0.054645 0.055147\n", + "\t175 12826.919565 12821.111309 0.052669 0.052650\n", + "\t180 12774.051094 12755.985754 0.051284 0.051381\n", + "\t185 12712.958930 12699.330423 0.049636 0.049963\n", + "\t190 12723.636192 12658.081446 0.047333 0.048248\n", + "\t195 12607.550913 12616.588836 0.046856 0.046843\n", + "\t200 12557.223085 12520.004285 0.045825 0.046220\n", + "\t205 12496.316576 12469.494651 0.044607 0.044840\n", + "\t210 12359.901379 12336.831445 0.044432 0.044645\n", + "\t215 12291.488426 12292.263511 0.043549 0.043399\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 23/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21436e+43 erg / s Luminosity absorbed = 8.20267e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13242.589560 13258.427166 0.650661 0.646781\n", + "\t5 13302.236433 13299.514381 0.553883 0.554208\n", + "\t10 13384.295437 13387.580308 0.482708 0.480070\n", + "\t15 13461.126863 13417.789738 0.422022 0.427611\n", + "\t20 13470.826671 13478.947649 0.381236 0.380802\n", + "\t25 13484.247273 13483.853992 0.347424 0.346575\n", + "\t30 13508.470387 13481.971326 0.313096 0.315613\n", + "\t35 13543.424565 13537.341581 0.283193 0.283614\n", + "\t40 13580.176047 13545.349426 0.255980 0.258512\n", + "\t45 13614.362136 13606.415720 0.233483 0.234589\n", + "\t50 13644.537532 13633.155825 0.212147 0.212564\n", + "\t55 13614.396375 13620.012253 0.196233 0.196128\n", + "\t60 13645.330715 13634.573954 0.179325 0.179370\n", + "\t65 13671.254373 13674.376989 0.164136 0.163591\n", + "\t70 13646.107424 13696.720677 0.154460 0.151282\n", + "\t75 13616.420307 13628.110866 0.143180 0.142235\n", + "\t80 13666.532072 13664.116070 0.130676 0.130906\n", + "\t85 13575.991311 13593.787168 0.125264 0.123801\n", + "\t90 13537.593047 13574.594978 0.116523 0.115843\n", + "\t95 13519.324997 13545.986411 0.109037 0.109125\n", + "\t100 13533.503766 13541.716260 0.102218 0.102313\n", + "\t105 13524.296286 13557.042736 0.095472 0.095252\n", + "\t110 13408.597159 13528.911716 0.093029 0.090217\n", + "\t115 13433.139621 13477.779105 0.086510 0.085821\n", + "\t120 13346.445735 13422.628881 0.082582 0.081440\n", + "\t125 13314.078372 13348.244204 0.078528 0.078169\n", + "\t130 13292.870438 13320.460508 0.074822 0.074200\n", + "\t135 13212.983843 13253.512286 0.071947 0.071233\n", + "\t140 13187.883546 13192.120721 0.068256 0.068204\n", + "\t145 13154.962528 13155.656117 0.064997 0.064981\n", + "\t150 13089.139849 13077.428673 0.063208 0.063336\n", + "\t155 13028.832892 13036.653378 0.061212 0.060716\n", + "\t160 12955.890317 13014.311099 0.058877 0.057776\n", + "\t165 12919.786297 12900.924993 0.056736 0.056904\n", + "\t170 12831.233593 12849.453995 0.055147 0.054912\n", + "\t175 12821.111309 12859.855508 0.052650 0.052202\n", + "\t180 12755.985754 12760.002356 0.051381 0.051311\n", + "\t185 12699.330423 12700.724143 0.049963 0.050006\n", + "\t190 12658.081446 12650.853146 0.048248 0.048460\n", + "\t195 12616.588836 12569.990044 0.046843 0.047719\n", + "\t200 12520.004285 12505.132653 0.046220 0.046542\n", + "\t205 12469.494651 12440.163243 0.044840 0.045558\n", + "\t210 12336.831445 12347.349436 0.044645 0.044723\n", + "\t215 12292.263511 12297.391221 0.043399 0.043441\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 24/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21096e+43 erg / s Luminosity absorbed = 8.23273e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13258.427166 13278.010763 0.646781 0.641648\n", + "\t5 13299.514381 13348.641733 0.554208 0.549116\n", + "\t10 13387.580308 13389.995159 0.480070 0.480928\n", + "\t15 13417.789738 13420.336267 0.427611 0.429481\n", + "\t20 13478.947649 13465.269674 0.380802 0.385639\n", + "\t25 13483.853992 13432.804473 0.346575 0.352148\n", + "\t30 13481.971326 13478.152260 0.315613 0.315732\n", + "\t35 13537.341581 13546.768152 0.283614 0.281819\n", + "\t40 13545.349426 13627.278155 0.258512 0.251783\n", + "\t45 13606.415720 13604.882965 0.234589 0.232319\n", + "\t50 13633.155825 13625.845890 0.212564 0.211901\n", + "\t55 13620.012253 13602.102375 0.196128 0.197460\n", + "\t60 13634.573954 13666.469517 0.179370 0.178943\n", + "\t65 13674.376989 13631.745802 0.163591 0.166107\n", + "\t70 13696.720677 13632.637553 0.151282 0.154189\n", + "\t75 13628.110866 13585.427048 0.142235 0.143931\n", + "\t80 13664.116070 13609.745816 0.130906 0.133334\n", + "\t85 13593.787168 13624.246108 0.123801 0.122952\n", + "\t90 13574.594978 13590.553890 0.115843 0.115902\n", + "\t95 13545.986411 13594.756383 0.109125 0.107573\n", + "\t100 13541.716260 13549.898466 0.102313 0.102064\n", + "\t105 13557.042736 13581.169680 0.095252 0.094850\n", + "\t110 13528.911716 13545.934037 0.090217 0.090058\n", + "\t115 13477.779105 13474.183557 0.085821 0.085789\n", + "\t120 13422.628881 13394.318438 0.081440 0.081967\n", + "\t125 13348.244204 13350.571807 0.078169 0.078213\n", + "\t130 13320.460508 13348.176678 0.074200 0.073127\n", + "\t135 13253.512286 13248.295587 0.071233 0.071241\n", + "\t140 13192.120721 13223.932093 0.068204 0.067701\n", + "\t145 13155.656117 13196.467541 0.064981 0.064614\n", + "\t150 13077.428673 13086.152871 0.063336 0.063554\n", + "\t155 13036.653378 13042.460884 0.060716 0.060771\n", + "\t160 13014.311099 13046.647251 0.057776 0.057615\n", + "\t165 12900.924993 12940.265617 0.056904 0.056227\n", + "\t170 12849.453995 12889.311505 0.054912 0.054381\n", + "\t175 12859.855508 12849.617688 0.052202 0.052379\n", + "\t180 12760.002356 12778.700141 0.051311 0.050977\n", + "\t185 12700.724143 12705.841566 0.050006 0.050094\n", + "\t190 12650.853146 12649.616025 0.048460 0.048563\n", + "\t195 12569.990044 12616.346040 0.047719 0.047103\n", + "\t200 12505.132653 12531.869670 0.046542 0.046301\n", + "\t205 12440.163243 12465.741589 0.045558 0.045231\n", + "\t210 12347.349436 12411.682924 0.044723 0.043975\n", + "\t215 12297.391221 12321.315541 0.043441 0.043216\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 25/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21222e+43 erg / s Luminosity absorbed = 8.22471e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13278.010763 13270.696022 0.641648 0.642908\n", + "\t5 13348.641733 13322.383942 0.549116 0.552644\n", + "\t10 13389.995159 13347.242864 0.480928 0.485565\n", + "\t15 13420.336267 13374.912459 0.429481 0.433836\n", + "\t20 13465.269674 13432.753510 0.385639 0.389339\n", + "\t25 13432.804473 13457.984438 0.352148 0.351032\n", + "\t30 13478.152260 13472.554559 0.315732 0.317566\n", + "\t35 13546.768152 13504.185333 0.281819 0.286299\n", + "\t40 13627.278155 13558.293549 0.251783 0.258565\n", + "\t45 13604.882965 13517.497337 0.232319 0.238954\n", + "\t50 13625.845890 13556.174094 0.211901 0.217968\n", + "\t55 13602.102375 13576.181423 0.197460 0.197895\n", + "\t60 13666.469517 13645.945023 0.178943 0.178924\n", + "\t65 13631.745802 13615.151776 0.166107 0.167639\n", + "\t70 13632.637553 13628.125067 0.154189 0.154651\n", + "\t75 13585.427048 13622.936479 0.143931 0.143012\n", + "\t80 13609.745816 13578.172623 0.133334 0.134668\n", + "\t85 13624.246108 13627.023042 0.122952 0.122788\n", + "\t90 13590.553890 13561.024079 0.115902 0.116892\n", + "\t95 13594.756383 13579.888510 0.107573 0.107658\n", + "\t100 13549.898466 13534.426606 0.102064 0.102287\n", + "\t105 13581.169680 13502.800064 0.094850 0.096255\n", + "\t110 13545.934037 13492.870611 0.090058 0.091280\n", + "\t115 13474.183557 13490.398325 0.085789 0.085011\n", + "\t120 13394.318438 13395.721908 0.081967 0.081828\n", + "\t125 13350.571807 13350.634171 0.078213 0.078078\n", + "\t130 13348.176678 13327.064866 0.073127 0.073579\n", + "\t135 13248.295587 13255.401860 0.071241 0.070985\n", + "\t140 13223.932093 13256.925404 0.067701 0.066957\n", + "\t145 13196.467541 13201.182229 0.064614 0.064230\n", + "\t150 13086.152871 13093.443767 0.063554 0.063075\n", + "\t155 13042.460884 13030.185698 0.060771 0.060808\n", + "\t160 13046.647251 13023.151556 0.057615 0.057750\n", + "\t165 12940.265617 12962.788653 0.056227 0.055743\n", + "\t170 12889.311505 12850.336719 0.054381 0.055013\n", + "\t175 12849.617688 12862.277714 0.052379 0.052236\n", + "\t180 12778.700141 12754.450331 0.050977 0.051508\n", + "\t185 12705.841566 12704.562044 0.050094 0.050137\n", + "\t190 12649.616025 12631.226174 0.048563 0.048881\n", + "\t195 12616.346040 12617.823980 0.047103 0.047068\n", + "\t200 12531.869670 12502.116047 0.046301 0.046602\n", + "\t205 12465.741589 12448.945824 0.045231 0.045355\n", + "\t210 12411.682924 12392.114495 0.043975 0.044175\n", + "\t215 12321.315541 12313.877904 0.043216 0.043334\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 26/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21673e+43 erg / s Luminosity absorbed = 8.18157e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13270.696022 13271.530613 0.642908 0.644373\n", + "\t5 13322.383942 13336.456993 0.552644 0.550015\n", + "\t10 13347.242864 13390.353000 0.485565 0.479866\n", + "\t15 13374.912459 13423.306105 0.433836 0.428565\n", + "\t20 13432.753510 13450.367080 0.389339 0.385605\n", + "\t25 13457.984438 13487.378370 0.351032 0.345360\n", + "\t30 13472.554559 13493.026147 0.317566 0.313875\n", + "\t35 13504.185333 13474.848725 0.286299 0.288170\n", + "\t40 13558.293549 13541.466576 0.258565 0.259858\n", + "\t45 13517.497337 13545.038773 0.238954 0.237364\n", + "\t50 13556.174094 13572.532271 0.217968 0.216331\n", + "\t55 13576.181423 13580.512176 0.197895 0.198667\n", + "\t60 13645.945023 13612.481505 0.178924 0.180448\n", + "\t65 13615.151776 13634.562335 0.167639 0.166978\n", + "\t70 13628.125067 13656.348500 0.154651 0.153620\n", + "\t75 13622.936479 13648.706367 0.143012 0.142383\n", + "\t80 13578.172623 13614.335979 0.134668 0.132803\n", + "\t85 13627.023042 13624.661577 0.122788 0.123129\n", + "\t90 13561.024079 13601.743908 0.116892 0.116093\n", + "\t95 13579.888510 13596.847242 0.107658 0.107682\n", + "\t100 13534.426606 13537.100559 0.102287 0.101815\n", + "\t105 13502.800064 13528.220122 0.096255 0.095814\n", + "\t110 13492.870611 13494.146708 0.091280 0.091091\n", + "\t115 13490.398325 13473.791902 0.085011 0.084878\n", + "\t120 13395.721908 13387.298820 0.081828 0.082027\n", + "\t125 13350.634171 13385.907140 0.078078 0.077395\n", + "\t130 13327.064866 13331.066428 0.073579 0.073484\n", + "\t135 13255.401860 13303.936110 0.070985 0.069917\n", + "\t140 13256.925404 13224.287619 0.066957 0.067565\n", + "\t145 13201.182229 13189.834245 0.064230 0.064769\n", + "\t150 13093.443767 13130.766284 0.063075 0.062414\n", + "\t155 13030.185698 13079.978377 0.060808 0.060190\n", + "\t160 13023.151556 13006.542174 0.057750 0.058257\n", + "\t165 12962.788653 12994.794135 0.055743 0.055493\n", + "\t170 12850.336719 12917.381157 0.055013 0.054132\n", + "\t175 12862.277714 12871.990795 0.052236 0.052206\n", + "\t180 12754.450331 12787.392876 0.051508 0.051257\n", + "\t185 12704.562044 12739.362606 0.050137 0.049526\n", + "\t190 12631.226174 12654.117488 0.048881 0.048597\n", + "\t195 12617.823980 12621.624550 0.047068 0.046930\n", + "\t200 12502.116047 12517.729051 0.046602 0.046454\n", + "\t205 12448.945824 12475.287394 0.045355 0.045014\n", + "\t210 12392.114495 12421.244417 0.044175 0.043767\n", + "\t215 12313.877904 12305.788275 0.043334 0.043532\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 27/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.20760e+43 erg / s Luminosity absorbed = 8.26883e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13271.530613 13267.229895 0.644373 0.644613\n", + "\t5 13336.456993 13309.957848 0.550015 0.552350\n", + "\t10 13390.353000 13362.599180 0.479866 0.484647\n", + "\t15 13423.306105 13416.485413 0.428565 0.430518\n", + "\t20 13450.367080 13452.533705 0.385605 0.385570\n", + "\t25 13487.378370 13480.877666 0.345360 0.345231\n", + "\t30 13493.026147 13534.801383 0.313875 0.309576\n", + "\t35 13474.848725 13562.348057 0.288170 0.280234\n", + "\t40 13541.466576 13556.692056 0.259858 0.257486\n", + "\t45 13545.038773 13561.290843 0.237364 0.235385\n", + "\t50 13572.532271 13548.608594 0.216331 0.218351\n", + "\t55 13580.512176 13525.096450 0.198667 0.201994\n", + "\t60 13612.481505 13602.430049 0.180448 0.180981\n", + "\t65 13634.562335 13604.584566 0.166978 0.167370\n", + "\t70 13656.348500 13633.088350 0.153620 0.154084\n", + "\t75 13648.706367 13626.329192 0.142383 0.142836\n", + "\t80 13614.335979 13619.696143 0.132803 0.132119\n", + "\t85 13624.661577 13631.876027 0.123129 0.122202\n", + "\t90 13601.743908 13607.669741 0.116093 0.114933\n", + "\t95 13596.847242 13580.712969 0.107682 0.107904\n", + "\t100 13537.100559 13558.070068 0.101815 0.101611\n", + "\t105 13528.220122 13513.642542 0.095814 0.096155\n", + "\t110 13494.146708 13484.539851 0.091091 0.090864\n", + "\t115 13473.791902 13518.564728 0.084878 0.083821\n", + "\t120 13387.298820 13442.773703 0.082027 0.080492\n", + "\t125 13385.907140 13428.879039 0.077395 0.076252\n", + "\t130 13331.066428 13348.870947 0.073484 0.073431\n", + "\t135 13303.936110 13288.319864 0.069917 0.070795\n", + "\t140 13224.287619 13295.534754 0.067565 0.066347\n", + "\t145 13189.834245 13199.387702 0.064769 0.064713\n", + "\t150 13130.766284 13144.816061 0.062414 0.061923\n", + "\t155 13079.978377 13097.206394 0.060190 0.059794\n", + "\t160 13006.542174 13083.035223 0.058257 0.057139\n", + "\t165 12994.794135 13015.077160 0.055493 0.055308\n", + "\t170 12917.381157 12939.781727 0.054132 0.053766\n", + "\t175 12871.990795 12889.826041 0.052206 0.051931\n", + "\t180 12787.392876 12815.502907 0.051257 0.050555\n", + "\t185 12739.362606 12731.737958 0.049526 0.049591\n", + "\t190 12654.117488 12662.180603 0.048597 0.048226\n", + "\t195 12621.624550 12656.423731 0.046930 0.046356\n", + "\t200 12517.729051 12533.340286 0.046454 0.045948\n", + "\t205 12475.287394 12492.138421 0.045014 0.044581\n", + "\t210 12421.244417 12451.936111 0.043767 0.043184\n", + "\t215 12305.788275 12322.819248 0.043532 0.043061\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 28/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21615e+43 erg / s Luminosity absorbed = 8.17960e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13267.229895 13238.322676 0.644613 0.649373\n", + "\t5 13309.957848 13286.184894 0.552350 0.556429\n", + "\t10 13362.599180 13340.900794 0.484647 0.488148\n", + "\t15 13416.485413 13393.757553 0.430518 0.432761\n", + "\t20 13452.533705 13470.248408 0.385570 0.382713\n", + "\t25 13480.877666 13514.190373 0.345231 0.343651\n", + "\t30 13534.801383 13533.440429 0.309576 0.311547\n", + "\t35 13562.348057 13560.082221 0.280234 0.281544\n", + "\t40 13556.692056 13552.468579 0.257486 0.259401\n", + "\t45 13561.290843 13585.031702 0.235385 0.235486\n", + "\t50 13548.608594 13576.630341 0.218351 0.217249\n", + "\t55 13525.096450 13594.171150 0.201994 0.198141\n", + "\t60 13602.430049 13634.500032 0.180981 0.180635\n", + "\t65 13604.584566 13662.297529 0.167370 0.166281\n", + "\t70 13633.088350 13584.906890 0.154084 0.157331\n", + "\t75 13626.329192 13634.212626 0.142836 0.143265\n", + "\t80 13619.696143 13631.150560 0.132119 0.132716\n", + "\t85 13631.876027 13621.164051 0.122202 0.123512\n", + "\t90 13607.669741 13580.266266 0.114933 0.116573\n", + "\t95 13580.712969 13604.376545 0.107904 0.107431\n", + "\t100 13558.070068 13567.082796 0.101611 0.101743\n", + "\t105 13513.642542 13543.001457 0.096155 0.095716\n", + "\t110 13484.539851 13487.048499 0.090864 0.091156\n", + "\t115 13518.564728 13478.088798 0.083821 0.085159\n", + "\t120 13442.773703 13443.911065 0.080492 0.081045\n", + "\t125 13428.879039 13364.123108 0.076252 0.078233\n", + "\t130 13348.870947 13352.151229 0.073431 0.074123\n", + "\t135 13288.319864 13285.766577 0.070795 0.070851\n", + "\t140 13295.534754 13282.190208 0.066347 0.066687\n", + "\t145 13199.387702 13188.937557 0.064713 0.065227\n", + "\t150 13144.816061 13140.181851 0.061923 0.062637\n", + "\t155 13097.206394 13076.784967 0.059794 0.060149\n", + "\t160 13083.035223 13062.198936 0.057139 0.057810\n", + "\t165 13015.077160 12985.941742 0.055308 0.056064\n", + "\t170 12939.781727 12903.272755 0.053766 0.054406\n", + "\t175 12889.826041 12874.062323 0.051931 0.052334\n", + "\t180 12815.502907 12786.132777 0.050555 0.051414\n", + "\t185 12731.737958 12752.988933 0.049591 0.049518\n", + "\t190 12662.180603 12689.779763 0.048226 0.048065\n", + "\t195 12656.423731 12652.130967 0.046356 0.046490\n", + "\t200 12533.340286 12558.648843 0.045948 0.045720\n", + "\t205 12492.138421 12520.764854 0.044581 0.044458\n", + "\t210 12451.936111 12404.766604 0.043184 0.043961\n", + "\t215 12322.819248 12304.711932 0.043061 0.043405\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 29/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.20884e+43 erg / s Luminosity absorbed = 8.25550e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", + "\t t_rad next_t_rad w next_w\n", + "\tShell \n", + "\t0 13238.322676 13226.520939 0.649373 0.651729\n", + "\t5 13286.184894 13303.405286 0.556429 0.555673\n", + "\t10 13340.900794 13318.466941 0.488148 0.494104\n", + "\t15 13393.757553 13411.400215 0.432761 0.431949\n", + "\t20 13470.248408 13488.042544 0.382713 0.381133\n", + "\t25 13514.190373 13542.705706 0.343651 0.339527\n", + "\t30 13533.440429 13568.331422 0.311547 0.306915\n", + "\t35 13560.082221 13551.572869 0.281544 0.281839\n", + "\t40 13552.468579 13592.277138 0.259401 0.255392\n", + "\t45 13585.031702 13589.751369 0.235486 0.232627\n", + "\t50 13576.630341 13632.800307 0.217249 0.212804\n", + "\t55 13594.171150 13640.303834 0.198141 0.194787\n", + "\t60 13634.500032 13677.204511 0.180635 0.177897\n", + "\t65 13662.297529 13681.215583 0.166281 0.164354\n", + "\t70 13584.906890 13597.256415 0.157331 0.155568\n", + "\t75 13634.212626 13607.327783 0.143265 0.143908\n", + "\t80 13631.150560 13613.045200 0.132716 0.132719\n", + "\t85 13621.164051 13595.524194 0.123512 0.123664\n", + "\t90 13580.266266 13572.564003 0.116573 0.115506\n", + "\t95 13604.376545 13595.132272 0.107431 0.106680\n", + "\t100 13567.082796 13546.383010 0.101743 0.101658\n", + "\t105 13543.001457 13544.145309 0.095716 0.095267\n", + "\t110 13487.048499 13506.548914 0.091156 0.089785\n", + "\t115 13478.088798 13474.167961 0.085159 0.085114\n", + "\t120 13443.911065 13452.053813 0.081045 0.080430\n", + "\t125 13364.123108 13399.637044 0.078233 0.077133\n", + "\t130 13352.151229 13372.115285 0.074123 0.073564\n", + "\t135 13285.766577 13329.206300 0.070851 0.070258\n", + "\t140 13282.190208 13242.263561 0.066687 0.067495\n", + "\t145 13188.937557 13211.178475 0.065227 0.064506\n", + "\t150 13140.181851 13203.869308 0.062637 0.061557\n", + "\t155 13076.784967 13115.889609 0.060149 0.059769\n", + "\t160 13062.198936 13039.682821 0.057810 0.057924\n", + "\t165 12985.941742 13009.593050 0.056064 0.055611\n", + "\t170 12903.272755 12944.363874 0.054406 0.053647\n", + "\t175 12874.062323 12919.678296 0.052334 0.051682\n", + "\t180 12786.132777 12835.572161 0.051414 0.050506\n", + "\t185 12752.988933 12766.745125 0.049518 0.049311\n", + "\t190 12689.779763 12689.244008 0.048065 0.048149\n", + "\t195 12652.130967 12608.255163 0.046490 0.046950\n", + "\t200 12558.648843 12520.624767 0.045720 0.046167\n", + "\t205 12520.764854 12493.746296 0.044458 0.044625\n", + "\t210 12404.766604 12405.512267 0.043961 0.043934\n", + "\t215 12304.711932 12314.490441 0.043405 0.043273\n", + "\n", + " (\u001b[1mbase.py\u001b[0m:348)\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 30/30 (\u001b[1mbase.py\u001b[0m:266)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.20889e+43 erg / s Luminosity absorbed = 8.25019e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", + "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Simulation finished in 30 iterations and took 356.68 s (\u001b[1mbase.py\u001b[0m:306)\n" + ] + } + ], + "source": [ + "sim = run_tardis('blondin_model_compare.yml')" + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": {}, + "outputs": [], + "source": [ + "spec_artis_toy = pd.read_csv('spec_artis_toy.csv')" + ] + }, + { + "cell_type": "code", + "execution_count": 9, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/matplotlib/style/core.py:167: UserWarning: In /Users/wkerzend/.matplotlib/stylelib/talk_default.mplstyle: \n", + "The text.latex.unicode rcparam was deprecated in Matplotlib 2.2 and will be removed in 3.1.\n", + " styles = read_style_directory(stylelib_path)\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n", + "Populating the interactive namespace from numpy and matplotlib\n", + "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.\n", + " AstropyDeprecationWarning)\n", + " (\u001b[1mwarnings.py\u001b[0m:99)\n" + ] + }, + { + "data": { + "application/javascript": [ + "/* Put everything inside the global mpl namespace */\n", + "window.mpl = {};\n", + "\n", + "\n", + "mpl.get_websocket_type = function() {\n", + " if (typeof(WebSocket) !== 'undefined') {\n", + " return WebSocket;\n", + " } else if (typeof(MozWebSocket) !== 'undefined') {\n", + " return MozWebSocket;\n", + " } else {\n", + " alert('Your browser does not have WebSocket support.' +\n", + " 'Please try Chrome, Safari or Firefox ≥ 6. ' +\n", + " 'Firefox 4 and 5 are also supported but you ' +\n", + " 'have to enable WebSockets in about:config.');\n", + " };\n", + "}\n", + "\n", + "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n", + " this.id = figure_id;\n", + "\n", + " this.ws = websocket;\n", + "\n", + " this.supports_binary = (this.ws.binaryType != undefined);\n", + "\n", + " if (!this.supports_binary) {\n", + " var warnings = document.getElementById(\"mpl-warnings\");\n", + " if (warnings) {\n", + " warnings.style.display = 'block';\n", + " warnings.textContent = (\n", + " \"This browser does not support binary websocket messages. \" +\n", + " \"Performance may be slow.\");\n", + " }\n", + " }\n", + "\n", + " this.imageObj = new Image();\n", + "\n", + " this.context = undefined;\n", + " this.message = undefined;\n", + " this.canvas = undefined;\n", + " this.rubberband_canvas = undefined;\n", + " this.rubberband_context = undefined;\n", + " this.format_dropdown = undefined;\n", + "\n", + " this.image_mode = 'full';\n", + "\n", + " this.root = $('
');\n", + " this._root_extra_style(this.root)\n", + " this.root.attr('style', 'display: inline-block');\n", + "\n", + " $(parent_element).append(this.root);\n", + "\n", + " this._init_header(this);\n", + " this._init_canvas(this);\n", + " this._init_toolbar(this);\n", + "\n", + " var fig = this;\n", + "\n", + " this.waiting = false;\n", + "\n", + " this.ws.onopen = function () {\n", + " fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n", + " fig.send_message(\"send_image_mode\", {});\n", + " if (mpl.ratio != 1) {\n", + " fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n", + " }\n", + " fig.send_message(\"refresh\", {});\n", + " }\n", + "\n", + " this.imageObj.onload = function() {\n", + " if (fig.image_mode == 'full') {\n", + " // Full images could contain transparency (where diff images\n", + " // almost always do), so we need to clear the canvas so that\n", + " // there is no ghosting.\n", + " fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n", + " }\n", + " fig.context.drawImage(fig.imageObj, 0, 0);\n", + " };\n", + "\n", + " this.imageObj.onunload = function() {\n", + " fig.ws.close();\n", + " }\n", + "\n", + " this.ws.onmessage = this._make_on_message_function(this);\n", + "\n", + " this.ondownload = ondownload;\n", + "}\n", + "\n", + "mpl.figure.prototype._init_header = function() {\n", + " var titlebar = $(\n", + " '
');\n", + " var titletext = $(\n", + " '
');\n", + " titlebar.append(titletext)\n", + " this.root.append(titlebar);\n", + " this.header = titletext[0];\n", + "}\n", + "\n", + "\n", + "\n", + "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n", + "\n", + "}\n", + "\n", + "\n", + "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n", + "\n", + "}\n", + "\n", + "mpl.figure.prototype._init_canvas = function() {\n", + " var fig = this;\n", + "\n", + " var canvas_div = $('
');\n", + "\n", + " canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n", + "\n", + " function canvas_keyboard_event(event) {\n", + " return fig.key_event(event, event['data']);\n", + " }\n", + "\n", + " canvas_div.keydown('key_press', canvas_keyboard_event);\n", + " canvas_div.keyup('key_release', canvas_keyboard_event);\n", + " this.canvas_div = canvas_div\n", + " this._canvas_extra_style(canvas_div)\n", + " this.root.append(canvas_div);\n", + "\n", + " var canvas = $('');\n", + " canvas.addClass('mpl-canvas');\n", + " canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n", + "\n", + " this.canvas = canvas[0];\n", + " this.context = canvas[0].getContext(\"2d\");\n", + "\n", + " var backingStore = this.context.backingStorePixelRatio ||\n", + "\tthis.context.webkitBackingStorePixelRatio ||\n", + "\tthis.context.mozBackingStorePixelRatio ||\n", + "\tthis.context.msBackingStorePixelRatio ||\n", + "\tthis.context.oBackingStorePixelRatio ||\n", + "\tthis.context.backingStorePixelRatio || 1;\n", + "\n", + " mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n", + "\n", + " var rubberband = $('');\n", + " rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n", + "\n", + " var pass_mouse_events = true;\n", + "\n", + " canvas_div.resizable({\n", + " start: function(event, ui) {\n", + " pass_mouse_events = false;\n", + " },\n", + " resize: function(event, ui) {\n", + " fig.request_resize(ui.size.width, ui.size.height);\n", + " },\n", + " stop: function(event, ui) {\n", + " pass_mouse_events = true;\n", + " fig.request_resize(ui.size.width, ui.size.height);\n", + " },\n", + " });\n", + "\n", + " function mouse_event_fn(event) {\n", + " if (pass_mouse_events)\n", + " return fig.mouse_event(event, event['data']);\n", + " }\n", + "\n", + " rubberband.mousedown('button_press', mouse_event_fn);\n", + " rubberband.mouseup('button_release', mouse_event_fn);\n", + " // Throttle sequential mouse events to 1 every 20ms.\n", + " rubberband.mousemove('motion_notify', mouse_event_fn);\n", + "\n", + " rubberband.mouseenter('figure_enter', mouse_event_fn);\n", + " rubberband.mouseleave('figure_leave', mouse_event_fn);\n", + "\n", + " canvas_div.on(\"wheel\", function (event) {\n", + " event = event.originalEvent;\n", + " event['data'] = 'scroll'\n", + " if (event.deltaY < 0) {\n", + " event.step = 1;\n", + " } else {\n", + " event.step = -1;\n", + " }\n", + " mouse_event_fn(event);\n", + " });\n", + "\n", + " canvas_div.append(canvas);\n", + " canvas_div.append(rubberband);\n", + "\n", + " this.rubberband = rubberband;\n", + " this.rubberband_canvas = rubberband[0];\n", + " this.rubberband_context = rubberband[0].getContext(\"2d\");\n", + " this.rubberband_context.strokeStyle = \"#000000\";\n", + "\n", + " this._resize_canvas = function(width, height) {\n", + " // Keep the size of the canvas, canvas container, and rubber band\n", + " // canvas in synch.\n", + " canvas_div.css('width', width)\n", + " canvas_div.css('height', height)\n", + "\n", + " canvas.attr('width', width * mpl.ratio);\n", + " canvas.attr('height', height * mpl.ratio);\n", + " canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n", + "\n", + " rubberband.attr('width', width);\n", + " rubberband.attr('height', height);\n", + " }\n", + "\n", + " // Set the figure to an initial 600x600px, this will subsequently be updated\n", + " // upon first draw.\n", + " this._resize_canvas(600, 600);\n", + "\n", + " // Disable right mouse context menu.\n", + " $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n", + " return false;\n", + " });\n", + "\n", + " function set_focus () {\n", + " canvas.focus();\n", + " canvas_div.focus();\n", + " }\n", + "\n", + " window.setTimeout(set_focus, 100);\n", + "}\n", + "\n", + "mpl.figure.prototype._init_toolbar = function() {\n", + " var fig = this;\n", + "\n", + " var nav_element = $('
')\n", + " nav_element.attr('style', 'width: 100%');\n", + " this.root.append(nav_element);\n", + "\n", + " // Define a callback function for later on.\n", + " function toolbar_event(event) {\n", + " return fig.toolbar_button_onclick(event['data']);\n", + " }\n", + " function toolbar_mouse_event(event) {\n", + " return fig.toolbar_button_onmouseover(event['data']);\n", + " }\n", + "\n", + " for(var toolbar_ind in mpl.toolbar_items) {\n", + " var name = mpl.toolbar_items[toolbar_ind][0];\n", + " var tooltip = mpl.toolbar_items[toolbar_ind][1];\n", + " var image = mpl.toolbar_items[toolbar_ind][2];\n", + " var method_name = mpl.toolbar_items[toolbar_ind][3];\n", + "\n", + " if (!name) {\n", + " // put a spacer in here.\n", + " continue;\n", + " }\n", + " var button = $('