diff --git a/.gitignore b/.gitignore
index 81c7539..27de75d 100644
--- a/.gitignore
+++ b/.gitignore
@@ -29,3 +29,4 @@ figs/
 figs/*
 archive/
 archive/*
+driver.sh
diff --git a/run/Ising_2024-11-12.jl b/run/Ising_2024-11-12.jl
index c1ac621..508bce2 100644
--- a/run/Ising_2024-11-12.jl
+++ b/run/Ising_2024-11-12.jl
@@ -18,10 +18,10 @@ end
 cases = Any[];
 κs = [0.0];
 #Ks = zip(1e-2*[0.0; 1.0; 0.0; 10.0], 1e-2*[1.0; 0.0; 10.0; 0.0]);
-Ks = zip(1e-1*[0.0; 1.0], 1e-1*[1.0; 0.0]);
+Ks = zip(1e-1*[0.0; 1.0; 0.0; 4.0], 1e-1*[1.0; 0.0; 4.0; 0.0]);
 #E0s = vcat(0.0, 0.005:0.005:1.0);
 E0s = [0.1; 1.0; 10.0];
-ns = Int[25];
+ns = Int[25; 100];
 bs = [1.0];
 kTs = [0.1; 1.0];
 Fzs = [0.0; 0.1; 0.2; 0.3; 0.4; 0.5; 1.0; 2.0; 3.0; 4.0; 5.0; 7.5; 10.0; 12.5; 15.0; 20.0; 25.0; 30.0; 
diff --git a/run/Ising_2024-11-15.jl b/run/Ising_2024-11-15.jl
new file mode 100644
index 0000000..bebf0d4
--- /dev/null
+++ b/run/Ising_2024-11-15.jl
@@ -0,0 +1,54 @@
+println("Hello world!");
+
+@everywhere using Printf;
+
+@show workdir = if length(ARGS) > 0
+  ARGS[1];
+else
+  ".";
+end
+
+@everywhere fmt(x) = @sprintf("%07d", round(Int, 1e3*x));
+@everywhere fmt_int(x) = @sprintf("%03d", x);
+
+@everywhere function prefix(case)
+    "E0-$(fmt(case[:E0]))_K1-$(fmt(case[:K1]))_K2-$(fmt(case[:K2]))_kT-$(fmt(case[:kT]))_Fz-$(fmt(case[:Fz]))_Fx-$(fmt(case[:Fx]))_n-$(fmt(case[:n]))_b-$(fmt(case[:b]))_kappa-$(fmt(case[:kappa]))_run-$(case[:run])";
+end
+
+cases = Any[];
+κs = [0.0];
+#Ks = zip(1e-2*[0.0; 1.0; 0.0; 10.0], 1e-2*[1.0; 0.0; 10.0; 0.0]);
+Ks = zip(1e-1*[0.0; 1.0; 0.0; 4.0], 1e-1*[1.0; 0.0; 4.0; 0.0]);
+#E0s = vcat(0.0, 0.005:0.005:1.0);
+E0s = [0.1; 1.0; 10.0];
+ns = Int[25; 100];
+bs = [1.0];
+kTs = [0.1; 1.0];
+Fxs = [0.0; 0.1; 0.2; 0.3; 0.4; 0.5; 1.0; 2.0; 3.0; 4.0; 5.0; 7.5; 10.0; 12.5; 15.0; 20.0; 25.0; 30.0; 
+       35.0; 40.0; 45.0; 50.0];
+Fzs = [0.0];
+runs = 1:10
+for b in bs, κ in κs, n in ns, Kvec in Ks, E0 in E0s, kT in kTs, Fx in Fxs, Fz in Fzs, run in runs
+  K1, K2 = Kvec;
+  push!(cases, Dict(:E0 => E0, :K1 => K1, :K2 => K2,
+                    :kT => kT, :Fz => Fz, :kappa => κ,
+                    :Fx => Fx, :n => n, :b => b, :run => run));
+end
+
+@info "total number of cases to run: $(length(cases))";
+
+mkpath(workdir);
+
+pmap(case -> begin;
+
+  outfile = joinpath(workdir, "$(prefix(case)).out");
+  if !isfile(outfile)
+    println("Running case: $case.");
+    command = `julia -t 1 -O 3 mcmc_clustering_eap_chain.jl --chain-type dielectric --energy-type Ising -b $(case[:b]) --bend-mod $(case[:kappa]) --E0 $(case[:E0]) --K1 $(case[:K1]) --K2 $(case[:K2]) --kT $(case[:kT]) --Fz $(case[:Fz]) --Fx $(case[:Fx]) -n $(case[:n]) --num-steps 2500000 --burn-in 100000 -v 2 --prefix $(joinpath(workdir, prefix(case))) --stepout 250`;
+    output = read(command, String);
+    write(outfile, output); 
+  else
+    println("Case: $case has already been run.");
+  end
+
+end, cases);
diff --git a/run/Ising_2024-11-19.jl b/run/Ising_2024-11-19.jl
new file mode 100644
index 0000000..83ebbd8
--- /dev/null
+++ b/run/Ising_2024-11-19.jl
@@ -0,0 +1,54 @@
+println("Hello world!");
+
+@everywhere using Printf;
+
+@show workdir = if length(ARGS) > 0
+  ARGS[1];
+else
+  ".";
+end
+
+@everywhere fmt(x) = @sprintf("%07d", round(Int, 1e3*x));
+@everywhere fmt_int(x) = @sprintf("%03d", x);
+
+@everywhere function prefix(case)
+    "E0-$(fmt(case[:E0]))_K1-$(fmt(case[:K1]))_K2-$(fmt(case[:K2]))_kT-$(fmt(case[:kT]))_Fz-$(fmt(case[:Fz]))_Fx-$(fmt(case[:Fx]))_n-$(fmt(case[:n]))_b-$(fmt(case[:b]))_kappa-$(fmt(case[:kappa]))_run-$(case[:run])";
+end
+
+cases = Any[];
+κs = [0.0];
+#Ks = zip(1e-2*[0.0; 1.0; 0.0; 10.0], 1e-2*[1.0; 0.0; 10.0; 0.0]);
+Ks = zip(1e-1*[0.0; 1.0; 0.0; 4.0], 1e-1*[1.0; 0.0; 4.0; 0.0]);
+#E0s = vcat(0.0, 0.005:0.005:1.0);
+E0s = [0.1; 1.0; 10.0];
+ns = Int[25; 100];
+bs = [1.0];
+kTs = [0.1; 1.0];
+Fxs = [0.0; 0.1; 0.2; 0.3; 0.4; 0.5; 1.0; 2.0; 3.0; 4.0; 5.0; 7.5; 10.0; 12.5; 15.0; 20.0; 25.0; 30.0; 
+       35.0; 40.0; 45.0; 50.0];
+Fzs = [0.0];
+runs = 1:10
+for b in bs, κ in κs, n in ns, Kvec in Ks, E0 in E0s, kT in kTs, Fx in Fxs, Fz in Fzs, run in runs
+  K1, K2 = Kvec;
+  push!(cases, Dict(:E0 => E0, :K1 => K1, :K2 => K2,
+                    :kT => kT, :Fz => Fx / sqrt(2), :kappa => κ,
+                    :Fx => Fx / sqrt(2), :n => n, :b => b, :run => run));
+end
+
+@info "total number of cases to run: $(length(cases))";
+
+mkpath(workdir);
+
+pmap(case -> begin;
+
+  outfile = joinpath(workdir, "$(prefix(case)).out");
+  if !isfile(outfile)
+    println("Running case: $case.");
+    command = `julia -t 1 -O 3 mcmc_clustering_eap_chain.jl --chain-type dielectric --energy-type Ising -b $(case[:b]) --bend-mod $(case[:kappa]) --E0 $(case[:E0]) --K1 $(case[:K1]) --K2 $(case[:K2]) --kT $(case[:kT]) --Fz $(case[:Fz]) --Fx $(case[:Fx]) -n $(case[:n]) --num-steps 2500000 --burn-in 100000 -v 2 --prefix $(joinpath(workdir, prefix(case))) --stepout 250`;
+    output = read(command, String);
+    write(outfile, output); 
+  else
+    println("Case: $case has already been run.");
+  end
+
+end, cases);