diff --git a/test/test_atom.cpp b/test/test_atom.cpp
index 7511143..da06393 100644
--- a/test/test_atom.cpp
+++ b/test/test_atom.cpp
@@ -128,7 +128,7 @@ TEST_CASE("Dirac Atom - transitions", "[DiracAtom]")
   REQUIRE(tmat.totalRate() * Physical::s == Approx(1.0775e+07).epsilon(1e-4));
 
   // Iron-like atom
-  DiracAtom daFe = DiracAtom(26, 1, 56, NuclearRadiusModel::SPHERE);
+  DiracAtom daFe = DiracAtom(26, 1, 56, -1 , NuclearRadiusModel::SPHERE);
 
   // 3d3/2 => 3p1/2
   tmat = daFe.getTransitionProbabilities(3, 2, false, 3, 1, false);
diff --git a/test/test_kappaa.cpp b/test/test_kappaa.cpp
index 1275fa8..711785d 100644
--- a/test/test_kappaa.cpp
+++ b/test/test_kappaa.cpp
@@ -35,7 +35,7 @@ int main(int argc, char *argv[])
     // Now, simulations
     for (int i = 0; i < dataZ.size(); ++i)
     {
-        DiracAtom da = DiracAtom(dataZ[i], Physical::m_mu, dataA[i], SPHERE);
+        DiracAtom da = DiracAtom(dataZ[i], Physical::m_mu, dataA[i], -1, SPHERE);
         DiracState p2 = da.getState(2, 1, false);
         DiracState s1 = da.getState(1, 0, false);