diff --git a/include/teqp/models/multifluid.hpp b/include/teqp/models/multifluid.hpp
index 611ccd1822529bea03b7d9a3a73cbe4c2ebf8abf..8c3db016bc82275a0affb9315c5473a4ab912582 100644
--- a/include/teqp/models/multifluid.hpp
+++ b/include/teqp/models/multifluid.hpp
@@ -772,16 +772,16 @@ auto build_multifluid_mutant(Model& model, const nlohmann::json& jj) {
 
                 auto BIP = entry["BIP"];
                 // Set the reducing function parameters in the copy
-                betaT(i, j) = BIP["betaT"];
+                betaT(i, j) = BIP.at("betaT");
                 betaT(j, i) = 1 / red.betaT(i, j);
-                betaV(i, j) = BIP["betaV"];
+                betaV(i, j) = BIP.at("betaV");
                 betaV(j, i) = 1 / red.betaV(i, j);
-                gammaT(i, j) = BIP["gammaT"]; gammaT(j, i) = gammaT(i, j);
-                gammaV(i, j) = BIP["gammaV"]; gammaV(j, i) = gammaV(i, j);
+                gammaT(i, j) = BIP.at("gammaT"); gammaT(j, i) = gammaT(i, j);
+                gammaV(i, j) = BIP.at("gammaV"); gammaV(j, i) = gammaV(i, j);
 
                 // Build the matrix of F and departure functions
                 auto dep = entry["departure"];
-                F(i, j) = BIP["Fij"];
+                F(i, j) = BIP.at("Fij");
                 F(j, i) = F(i, j);
                 funcs[i][j] = build_departure_function(dep);
                 funcs[j][i] = build_departure_function(dep);
@@ -827,19 +827,19 @@ auto build_multifluid_mutant_invariant(Model& model, const nlohmann::json& jj) {
 
                 auto BIP = entry["BIP"];
                 // Set the reducing function parameters in the copy
-                phiT(i, j) = BIP["phiT"];
+                phiT(i, j) = BIP.at("phiT");
                 phiT(j, i) = phiT(i, j);
-                lambdaT(i, j) = BIP["lambdaT"]; 
+                lambdaT(i, j) = BIP.at("lambdaT"); 
                 lambdaT(j, i) = -lambdaT(i, j);
 
-                phiV(i, j) = BIP["phiV"];
+                phiV(i, j) = BIP.at("phiV");
                 phiV(j, i) = phiV(i, j);
-                lambdaV(i, j) = BIP["lambdaV"]; 
+                lambdaV(i, j) = BIP.at("lambdaV"); 
                 lambdaV(j, i) = -lambdaV(i, j);
 
                 // Build the matrix of F and departure functions
                 auto dep = entry["departure"];
-                F(i, j) = BIP["Fij"];
+                F(i, j) = BIP.at("Fij");
                 F(j, i) = F(i, j);
                 funcs[i][j] = build_departure_function(dep);
                 funcs[j][i] = build_departure_function(dep);
diff --git a/setup.py b/setup.py
index fed119e872dc44c07c8723fbaceefad1ef717480..b4c8ca2118f44f528de9fd759e3fa78bd60e5c3c 100644
--- a/setup.py
+++ b/setup.py
@@ -10,7 +10,7 @@ from setuptools import setup, Extension
 from setuptools.command.build_ext import build_ext
 from distutils.version import LooseVersion
 
-VERSION = '0.1.0'
+VERSION = '0.1.1'
 with open('interface/teqpversion.hpp','w') as fpver:
     fpver.write(f'#include <string>\nconst std::string TEQPVERSION = "{VERSION}";')
 
@@ -72,4 +72,4 @@ setup(
     ext_modules=[CMakeExtension('teqp')],
     cmdclass=dict(build_ext=CMakeBuild),
     zip_safe=False,
-)
+)
\ No newline at end of file