Implement the extraction of the discovery year

Author: php1ic

include/inch/nuclide.hpp

@@ -267,7 +267,13 @@ class Nuclide
    *
    * \return Nothing
    */
-  inline void setYear() const { year = Converter::StringToInt(full_data, NUBASE_START_YEAR, NUBASE_END_YEAR); }
+  inline void setYear() const
+  {
+    // Some isotopes have no value for the year so we need to watch for that.
+    // Leave it as the default if no year is given
+    const auto value = full_data.substr(NUBASE_START_YEAR, NUBASE_END_YEAR - NUBASE_START_YEAR);
+    year             = std::isspace(value.front()) != 0 ? 1900 : Converter::StringToInt(value);
+  }
 
   /**
    * Extract the half life value from the data file

src/massTable.cpp

@@ -196,6 +196,12 @@ bool MassTable::readNUBASE(const std::filesystem::path& nubaseTable)
 
       isotope.setHalfLife();
 
+      // Discovery year was added after 2003
+      if (table_year != 3)
+        {
+          isotope.setYear();
+        }
+
       isotope.setDecayMode(pnSide, table_year);
 
       isotope.setNeutronOrProtonRich(pnSide);