diff --git a/xbanalysis/model/xba_reaction.py b/xbanalysis/model/xba_reaction.py
index e483fee3fccab608621b74d3038b5c5433a58f4d..31335edfe64ffc0e9eb9dc25f3a54e98b87602da 100644
--- a/xbanalysis/model/xba_reaction.py
+++ b/xbanalysis/model/xba_reaction.py
@@ -34,7 +34,7 @@ class XbaReaction:
             self.kinetic_law = s_reaction['kineticLaw'] if type(s_reaction['kineticLaw']) == str else ''
             # convert to numpy, replace local parameters with numerical values, inline lambda functions:
             expanded_kl = sbmlxdf.misc.mathml2numpy(self.kinetic_law)
-            for pid, val in self.local_params.items():
+            for pid, [val, _] in self.local_params.items():
                 expanded_kl = re.sub(r'\b' + pid + r'\b', str(val), expanded_kl)
             self.expanded_kl = expand_kineticlaw(expanded_kl, functions)
             for cid in compartments.keys():
diff --git a/xbanalysis/utils/utils.py b/xbanalysis/utils/utils.py
index 348bf00fd6356654b45b3ffcfdbf2d7fbc12b67b..b91343f0690f3bc342baa8a0e1bc456a4659db7f 100644
--- a/xbanalysis/utils/utils.py
+++ b/xbanalysis/utils/utils.py
@@ -92,7 +92,7 @@ def get_local_params(reaction):
     if 'localParams' in reaction:
         for lp in record_generator(reaction['localParams']):
             params = extract_params(lp)
-            local_params[params['id']] = float(params['value'])
+            local_params[params['id']] = [float(params['value']), params.get('units', 'dimensionless')]
     return local_params