diff --git a/examples/cpr_sns.py b/examples/cpr_sns.py
index 4845b45108e55441ce6967a853375da2245ce1d1..5f5b0fda4fbc03ce555168b4396c57ccc36436d5 100644
--- a/examples/cpr_sns.py
+++ b/examples/cpr_sns.py
@@ -95,7 +95,7 @@ def Gamma_to_alpha(Gamma_DP, Gamma_ST):
     return alpha_soc, eta
 
 
-def do(W_xi=6, multin=False):
+def do(W_xi=6, multin=True):
     T = 0.1
     Gamma_DP = 10
     Gamma_ST = 1
@@ -154,7 +154,7 @@ def do(W_xi=6, multin=False):
     ax = axs[0]
     ax.plot(phi / pi, Jx_mean.T / abs(Jx_mean).max(axis=1))
     if multin:
-        ax.plot(mphi / pi, Jx0_mean.T / abs(Jx0_mean).max(axis=1), "k:", alpha=0.25)
+        # ax.plot(mphi / pi, Jx0_mean.T / abs(Jx0_mean).max(axis=1), "k:", alpha=0.25)
         ax.plot(mphi / pi, Jx1_mean.T / abs(Jx1_mean).max(axis=1), "k:")
     ax.set_xlabel(r"$\varphi / \pi$")
     ax.set_ylabel(r"$I / I_{\mathrm{max}}$")
@@ -168,13 +168,13 @@ def do(W_xi=6, multin=False):
     ax = axs[1]
     ax.plot(L / xi, 100 * eff(Jx_mean))
     if multin:
-        ax.plot(L / xi, 100 * eff(Jx0_mean), "k:", alpha=0.25)
+        # ax.plot(L / xi, 100 * eff(Jx0_mean), "k:", alpha=0.25)
         ax.plot(L / xi, 100 * eff(Jx1_mean), "k:")
     ax.set_xlabel(r"$L / \xi$")
     ax.set_ylabel(r"$\eta$  [%]")
 
     fig.suptitle(
-        rf"$W = {W/xi} \xi_0$ $\Gamma_{{DP}} = {Gamma_DP} \Delta_0$, $\Gamma_{{ST}} = {Gamma_ST} \Delta_0$   ($\tilde{{\eta}} = {eta:.3g}$, $\tilde{{\alpha}} = {alpha:.3g}$)"
+        rf"$h = {h} \Delta_0$, $T = {T} \Delta_0$, $W = {W/xi} \xi_0$ $\Gamma_{{DP}} = {Gamma_DP} \Delta_0$, $\Gamma_{{ST}} = {Gamma_ST} \Delta_0$   ($\tilde{{\eta}} = {eta:.3g}$, $\tilde{{\alpha}} = {alpha:.3g}$)"
     )
 
     fig.savefig("cpr_sns.pdf")