diff --git a/Scripts/Spin_Alignment/spinal_exp.py b/Scripts/Spin_Alignment/spinal_exp.py
index 59e9770..214222e 100644
--- a/Scripts/Spin_Alignment/spinal_exp.py
+++ b/Scripts/Spin_Alignment/spinal_exp.py
@@ -103,11 +103,11 @@ def spinal_experiment(pars, run):
                        
     pars['PROG'] = 'spinal_experiment'                       
                        
-   # 8-step phase cycle (1-14 modifided to deal with T1-recovery and extended for Re/Im imbalance)
-    pars['PH1'] = [0,  270,   0, 270,   180,  90, 180,  90] # 1st (90-degree) pulse
-    pars['PH3'] = [90, 180,  90, 180,    90, 180,  90, 180] # 2nd (90-degree) pulse
-    pars['PH4'] = [90,  90, 270, 270,     0,   0, 180, 180]  # 3rd (90-degree) pulse
-    pars['PH2'] = [0,  180, 180,   0,    90, 270, 270,  90]  # receiver
+   # 8-step phase cycle (Schaefer et al. J Magn Res A 115 (1995))
+    pars['PH1'] = [ 0,   0, 180, 180,  90, 90, 270, 270] # 1st (90-degree) pulse
+    pars['PH3'] = [90, 270,  90, 270, 180, 90, 180,  90] # 2nd (90-degree) pulse
+    pars['PH4'] = [90]  # 3rd (90-degree) pulse
+    pars['PH2'] = [90, 270]  # receiver
         
    # read in variables:
     P90 = pars['P90']
diff --git a/Scripts/Spin_Alignment_Four_Pulses/spinal4pulses_exp.py b/Scripts/Spin_Alignment_Four_Pulses/spinal4pulses_exp.py
index 0fc8dbf..076a717 100644
--- a/Scripts/Spin_Alignment_Four_Pulses/spinal4pulses_exp.py
+++ b/Scripts/Spin_Alignment_Four_Pulses/spinal4pulses_exp.py
@@ -104,12 +104,12 @@ def spinal4pulses_experiment(pars, run):
                        
     pars['PROG'] = 'spinal4pulses_experiment'                       
                        
-   # 8-step phase cycle (1-14 modifided to deal with T1-recovery and extended for Re/Im imbalance)
-    pars['PH1'] = [0,  270,   0, 270,   180,  90, 180,  90] # 1st (90-degree) pulse
-    pars['PH3'] = [90, 180,  90, 180,    90, 180,  90, 180] # 2nd (90-degree) pulse
-    pars['PH4'] = [90,  90, 270, 270,     0,   0, 180, 180]  # 3rd (90-degree) pulse
-    pars['PH5'] = [90,  90,  90,  90,   180, 180, 180, 180]  # 3rd (90-degree) pulse
-    pars['PH2'] = [0,  180, 180,   0,    90, 270, 270,  90]  # receiver
+   # 8-step phase cycle (Schaefer et al. J Magn Res A 115 (1995))
+    pars['PH1'] = [0,    0, 180, 180,  90, 90, 270, 270] # 1st (90-degree) pulse
+    pars['PH3'] = [90, 270,  90, 270, 180, 90, 180,  90] # 2nd (90-degree) pulse
+    pars['PH4'] = [90] # 3rd (90-degree) pulse
+    pars['PH5'] = [90] * 8 + [270] * 8
+    pars['PH2'] = [90, 270]  # receiver
         
    # read in variables:
     P90 = pars['P90']
@@ -128,9 +128,6 @@ def spinal4pulses_experiment(pars, run):
     PH2 = pars['PH2'][run%len(pars['PH2'])]
     PHA = pars['PHA']
 
-    if (run/8)%2 != 0:
-        PH5 += 180   
-
    # set sampling parameters:
     SI = pars['SI']
     SW = pars['SW']
diff --git a/Scripts/Zeeman_Order/zeeman_exp.py b/Scripts/Zeeman_Order/zeeman_exp.py
index 5646fec..5850711 100644
--- a/Scripts/Zeeman_Order/zeeman_exp.py
+++ b/Scripts/Zeeman_Order/zeeman_exp.py
@@ -101,11 +101,11 @@ def zeeman_experiment(pars, run):
         
     pars['PROG'] = 'zeeman_experiment'
 
-   # 8-step phase cycle (1-21 modifided to deal with T1-recovery and extended for Re/Im imbalance)
-    pars['PH1'] = [0, 270,   0, 270,   180,  90, 180, 90] # 1st (90-degree) pulse
-    pars['PH3'] = [0,  90,   0,  90,     0,  90,   0, 90] # 2nd (90-degree) pulse
-    pars['PH4'] = [0,   0, 180, 180,   270, 270,  90, 90]  # 3rd (90-degree) pulse
-    pars['PH2'] = [0, 180, 180,   0,    90, 270, 270, 90]  # receiver
+   # 8-step phase cycle (Schaefer et al. J Magn Res A 115 (1995))
+    pars['PH1'] = [0,   0, 180, 180, 90,  90, 270, 270] # 1st (90-degree) pulse
+    pars['PH3'] = [0, 180,   0, 180, 90, 180,  90, 180] # 2nd (90-degree) pulse
+    pars['PH4'] = [90]  # 3rd (90-degree) pulse
+    pars['PH2'] = [90, 270]  # receiver
 
    # read in variables:
     P90 = pars['P90']
diff --git a/Scripts/Zeeman_Order_Four_Pulses/zeeman4pulses_exp.py b/Scripts/Zeeman_Order_Four_Pulses/zeeman4pulses_exp.py
index 7f59f98..68cb0a2 100644
--- a/Scripts/Zeeman_Order_Four_Pulses/zeeman4pulses_exp.py
+++ b/Scripts/Zeeman_Order_Four_Pulses/zeeman4pulses_exp.py
@@ -102,12 +102,12 @@ def zeeman4pulses_experiment(pars, run):
         
     pars['PROG'] = 'zeeman4pulses_experiment'
     
-         # ok 8-step phase cycle (1-21 modifided to deal with T1-recovery and extended for Re/Im imbalance)
-    pars['PH1'] = [0, 270,   0, 270,   180,  90, 180,  90] # 1st (90-degree) pulse
-    pars['PH3'] = [0,  90,   0,  90,     0,  90,   0,  90] # 2nd (90-degree) pulse
-    pars['PH4'] = [0,   0, 180, 180,   270, 270,  90,  90]  # 3rd (90-degree) pulse
-    pars['PH5'] = [90, 90,  90,  90,   180, 180, 180, 180]  # 3rd (90-degree) pulse
-    pars['PH2'] = [0, 180, 180,   0,    90, 270, 270,  90]  # receiver
+    # ok 8-step phase cycle (Schaefer et al. J Magn Res A 115 (1995))
+    pars['PH1'] = [0,   0, 180, 180, 90,  90, 270, 270] # 1st (90-degree) pulse
+    pars['PH3'] = [0, 180,   0, 180, 90, 180,  90, 270] # 2nd (90-degree) pulse
+    pars['PH4'] = [0]  # 3rd (90-degree) pulse
+    pars['PH5'] = [90] * 8 + [270] * 8
+    pars['PH2'] = [90, 270]  # receiver
     
  
    # read in variables:
@@ -125,9 +125,6 @@ def zeeman4pulses_experiment(pars, run):
     PH2 = pars['PH2'][run%len(pars['PH2'])]
     PHA = pars['PHA']
     
-    if (run/8)%2 != 0:
-        PH5 += 180  
-
    # set sampling parameters:
     SI = pars['SI']
     SW = pars['SW']