made the timo_0 in decay framework required

.mailmap

@@ -235,3 +235,11 @@ Le Truong <[email protected]>
 
 Kim Lingemann <[email protected]> kimsina <[email protected]>
 Kim Lingemann <[email protected]> kim <[email protected]>
+
+Le Truong <[email protected]>
+
+Kim Lingemann <[email protected]> kimsina <[email protected]>
+Kim Lingemann <[email protected]> kim <[email protected]>
+
+Soumya Tarafder <[email protected]>
+Soumya Tarafder <[email protected]>

docs/quickstart.ipynb

@@ -186,7 +186,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.8.13"
+   "version": "3.8.15"
   }
  },
  "nbformat": 4,

tardis/io/decay.py

@@ -9,11 +9,10 @@ class IsotopeAbundances(pd.DataFrame):
     _metadata = ["time_0"]
 
     def __init__(self, *args, **kwargs):
-        if "time_0" in kwargs:
-            time_0 = kwargs["time_0"]
-            kwargs.pop("time_0")
-        else:
-            time_0 = 0 * u.d
+        if "time_0" not in kwargs:
+            raise ValueError("The 'time_0' parameter must be specified when creating an instance of IsotopeAbundance.")
+        time_0 = kwargs["time_0"]
+        kwargs.pop("time_0")
         super(IsotopeAbundances, self).__init__(*args, **kwargs)
         self.time_0 = time_0
 

tardis/io/tests/test_decay.py

@@ -1,6 +1,6 @@
 import pytest
 import pandas as pd
-
+from astropy import units as u
 from tardis.io.decay import IsotopeAbundances
 from numpy.testing import assert_almost_equal
 
@@ -10,7 +10,7 @@ def simple_abundance_model():
     index = pd.MultiIndex.from_tuples(
         [(28, 56)], names=["atomic_number", "mass_number"]
     )
-    return IsotopeAbundances([[1.0, 1.0]], index=index)
+    return IsotopeAbundances([[1.0, 1.0]], index=index, time_0=0*u.s)
 
 
 def test_simple_decay(simple_abundance_model):

tardis/model/base.py

@@ -680,7 +680,7 @@ def from_config(cls, config, atom_data=None):
             abundance /= norm_factor
             isotope_abundance /= norm_factor
 
-        isotope_abundance = IsotopeAbundances(isotope_abundance)
+        isotope_abundance = IsotopeAbundances(isotope_abundance,time_0=0*u.s)
 
         elemental_mass = None
         if atom_data is not None:

tardis/model/tests/test_base.py

@@ -277,7 +277,7 @@ def simple_isotope_abundance():
         [(6, 14), (12, 28)], names=["atomic_number", "mass_number"]
     )
     abundance = [[0.2] * 20] * 2
-    return IsotopeAbundances(abundance, index=index)
+    return IsotopeAbundances(abundance, index=index, time_0=0*u.s)
 
 
 def test_model_decay(simple_isotope_abundance):