diff --git a/python/otto_motor/subprojects/cycle_shift/cycle_shift.org b/python/otto_motor/subprojects/cycle_shift/cycle_shift.org
index 9b36ccb..1d5255e 100644
--- a/python/otto_motor/subprojects/cycle_shift/cycle_shift.org
+++ b/python/otto_motor/subprojects/cycle_shift/cycle_shift.org
@@ -988,7 +988,7 @@ hamiltonian.
 
 #+begin_src jupyter-python :tangle tangle/rot.py
   rot_models = []
-  weights = [.3]
+  weights = [.5]
   for weight in weights:
       off_ax = sc.make_model(0, 0)
       off_ax.H_bias = ConstantMatrix(weight / 2 * qt.sigmax().full())
@@ -1001,17 +1001,17 @@ hamiltonian.
 #+begin_src jupyter-python
   τs = rot_models[0].t
   #plt.plot(τs, np.einsum('tij,ij->t', rot_models[0].H(τs), qt.sigmay().full()).real)
-  plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 0]))
-  plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 1]))
-  plt.plot(τs, abs(rot_models[0].H.operator_norm(τs)))
+  # plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 0]))
+  # plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 1]))
+  # plt.plot(τs, abs(rot_models[0].H.operator_norm(τs)))
   H = rot_models[0].H
   plt.plot(τs, list(map(lambda t: get_energy_gap(H(t)), τs)), color="black")
 #+end_src
 
 #+RESULTS:
 :RESULTS:
-| <matplotlib.lines.Line2D | at | 0x7f36244f3340> |
-[[file:./.ob-jupyter/e2600fe61c160903c60b7159b388b82b5ce667d0.svg]]
+| <matplotlib.lines.Line2D | at | 0x7f36470cdd00> |
+[[file:./.ob-jupyter/d5d7c229833877657c2abc95a30a085919c05d40.svg]]
 :END:
 
 ** Integration
diff --git a/python/otto_motor/subprojects/cycle_shift/tangle/cycle_shift.py b/python/otto_motor/subprojects/cycle_shift/tangle/cycle_shift.py
index 45cb7a7..8f29f1e 100644
--- a/python/otto_motor/subprojects/cycle_shift/tangle/cycle_shift.py
+++ b/python/otto_motor/subprojects/cycle_shift/tangle/cycle_shift.py
@@ -373,9 +373,9 @@ fs.export_fig(f"energy_change_off_axis", x_scaling=2, y_scaling=0.7)
 
 τs = rot_models[0].t
 #plt.plot(τs, np.einsum('tij,ij->t', rot_models[0].H(τs), qt.sigmay().full()).real)
-plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 0]))
-plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 1]))
-plt.plot(τs, abs(rot_models[0].H.operator_norm(τs)))
+# plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 0]))
+# plt.plot(τs, abs(rot_models[0].H(τs)[:, 0, 1]))
+# plt.plot(τs, abs(rot_models[0].H.operator_norm(τs)))
 H = rot_models[0].H
 plt.plot(τs, list(map(lambda t: get_energy_gap(H(t)), τs)), color="black")
 
diff --git a/python/otto_motor/subprojects/cycle_shift/tangle/rot.py b/python/otto_motor/subprojects/cycle_shift/tangle/rot.py
index 1b8ae2c..ae83931 100644
--- a/python/otto_motor/subprojects/cycle_shift/tangle/rot.py
+++ b/python/otto_motor/subprojects/cycle_shift/tangle/rot.py
@@ -22,7 +22,7 @@ import itertools
 from hops.util.dynamic_matrix import ConstantMatrix
 
 rot_models = []
-weights = [.3]
+weights = [.5]
 for weight in weights:
     off_ax = sc.make_model(0, 0)
     off_ax.H_bias = ConstantMatrix(weight / 2 * qt.sigmax().full())