diff --git a/data/simulations/2018-05-08-ad-invivo-1-contrastrates-c025-s1.npz b/data/simulations/2018-05-08-ad-invivo-1-contrastrates-c025-s1.npz
new file mode 100644
index 0000000..93f8312
Binary files /dev/null and b/data/simulations/2018-05-08-ad-invivo-1-contrastrates-c025-s1.npz differ
diff --git a/data/simulations/2018-05-08-ad-invivo-1-contrastrates-c050-s1.npz b/data/simulations/2018-05-08-ad-invivo-1-contrastrates-c050-s1.npz
new file mode 100644
index 0000000..91e2451
Binary files /dev/null and b/data/simulations/2018-05-08-ad-invivo-1-contrastrates-c050-s1.npz differ
diff --git a/data/simulations/2018-05-08-ad-invivo-1-contrastspectrum-c025.npz b/data/simulations/2018-05-08-ad-invivo-1-contrastspectrum-c025.npz
new file mode 100644
index 0000000..5c3a056
Binary files /dev/null and b/data/simulations/2018-05-08-ad-invivo-1-contrastspectrum-c025.npz differ
diff --git a/data/simulations/2018-05-08-ad-invivo-1-contrastspectrum-c050.npz b/data/simulations/2018-05-08-ad-invivo-1-contrastspectrum-c050.npz
new file mode 100644
index 0000000..69a4603
Binary files /dev/null and b/data/simulations/2018-05-08-ad-invivo-1-contrastspectrum-c050.npz differ
diff --git a/regimes.py b/regimes.py
index 23b70b4..9ca37e5 100644
--- a/regimes.py
+++ b/regimes.py
@@ -10,11 +10,12 @@ from plotstyle import plot_style, lighter, darker
 
 
 model_cell = '2018-05-08-ad-invivo-1'      # 228Hz, CV=0.67
-alphas = [0.002, 0.01, 0.03, 0.06]
+#alphas = [0.002, 0.01, 0.03, 0.06]
+alphas = [0.002, 0.01, 0.025, 0.05]
 rmax = 500
 amax = 60
 cthresh1 = 1.2
-cthresh2 = 3.05
+cthresh2 = 2.6
 
 #model_cell = '2018-05-08-ab-invivo-1'      # 116, CV=0.68
 #alphas = [0.002, 0.008, 0.025, 0.05]
@@ -351,8 +352,7 @@ def amplitude_squarefit(contrast, power, max_contrast):
 
 
 def plot_peaks(ax, s, alphas, contrasts, powerf1, powerf2, powerfsum,
-               powerfdiff):
-    cmax = 10
+               powerfdiff, cmax=8):
     contrasts *= 100
     ax.plot(contrasts, amplitude_linearfit(contrasts, powerf1, cthresh1),
             **s.lsF01m)
@@ -391,7 +391,7 @@ def plot_peaks(ax, s, alphas, contrasts, powerf1, powerf2, powerfsum,
                 ha='center', va='bottom')
     ax.set_xlim(0, cmax)
     ax.set_ylim(0, amax)
-    ax.set_xticks_delta(2)
+    ax.set_xticks_delta(1)
     ax.set_yticks_delta(20)
     ax.set_xlabel('Contrast', r'\%')
     ax.set_ylabel('Amplitude', 'Hz')