diff --git a/src/cmd/sim.cpp b/src/cmd/sim.cpp
index ec32c5b9..b3c37aa8 100644
--- a/src/cmd/sim.cpp
+++ b/src/cmd/sim.cpp
@@ -82,7 +82,7 @@ int main_sim(args::Subparser &parser)
   Sim::Result result;
   if (ng) {
     Sim::Parameter const gamma{ng.Get(), gLo.Get(), gHi.Get(), false};
-    result = Sim::Eddy(T1, beta, gamma, B1, seq, randomSamp.Get());
+    result = Sim::Eddy(T1, beta, gamma, seq, randomSamp.Get());
   } else if (mupa) {
     Sim::Parameter const T2{65, 0.02, 0.2, true};
     result = Sim::MUPA(T1, T2, B1, seq, randomSamp.Get());
diff --git a/src/sim/eddy.cpp b/src/sim/eddy.cpp
index 0a87aabf..50de2943 100644
--- a/src/sim/eddy.cpp
+++ b/src/sim/eddy.cpp
@@ -9,7 +9,6 @@ Result Eddy(
   Parameter const T1p,
   Parameter const betap,
   Parameter const gammap,
-  Parameter const B1p,
   Sequence const seq,
   Index const nRand)
 {
@@ -25,12 +24,11 @@ Result Eddy(
   Log::Print(FMT_STRING("{} values of T1 from {} to {}s"), T1p.N, T1p.lo, T1p.hi);
   Log::Print(FMT_STRING("{} values of β from {} to {}"), betap.N, betap.lo, betap.hi);
   Log::Print(FMT_STRING("{} values of ɣ from {} to {}"), gammap.N, gammap.lo, gammap.hi);
-  Log::Print(FMT_STRING("{} values of B1 from {} to {}"), B1p.N, B1p.lo, B1p.hi);
-  ParameterGenerator<4> gen(T1p, betap, gammap, B1p);
+  ParameterGenerator<3> gen(T1p, betap, gammap);
   Index totalN = (nRand > 0) ? nRand : gen.totalN();
   Result result;
   result.dynamics.resize(totalN, 4 * seq.sps);
-  result.parameters.resize(totalN, 4);
+  result.parameters.resize(totalN, 3);
   result.Mz_ss.resize(totalN);
 
   auto task = [&](Index const lo, Index const hi, Index const ti) {
diff --git a/src/sim/eddy.h b/src/sim/eddy.h
index 9fde40cb..09d21656 100644
--- a/src/sim/eddy.h
+++ b/src/sim/eddy.h
@@ -10,7 +10,6 @@ Result Eddy(
   Parameter const T1,
   Parameter const beta,
   Parameter const gamma,
-  Parameter const B1,
   Sequence const seq,
   Index const nRand);