diff --git a/script/Co_xmcd_step_pshell.py b/script/Co_xmcd_step_pshell.py
new file mode 100644
index 0000000..ec846e3
--- /dev/null
+++ b/script/Co_xmcd_step_pshell.py
@@ -0,0 +1,40 @@
+
+set_pol_cplus(offset = 0)
+wait_beam()
+set_file('Co_plus')
+scan_e_mult([(760,774,1),(774.5,783,0.1),(783.25,790,0.25),(790.2,798,0.2),(798.5,840,1)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cminus(offset = 2)
+wait_beam()
+set_file('Co_minus')
+scan_e_mult([(760,774,1),(774.5,783,0.1),(783.25,790,0.25),(790.2,798,0.2),(798.5,840,1)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cminus(offset = 2)
+wait_beam()
+set_file('Co_minus')
+scan_e_mult([(760,774,1),(774.5,783,0.1),(783.25,790,0.25),(790.2,798,0.2),(798.5,840,1)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cplus(offset = 0)
+wait_beam()
+set_file('Co_plus')
+scan_e_mult([(760,774,1),(774.5,783,0.1),(783.25,790,0.25),(790.2,798,0.2),(798.5,840,1)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cplus(offset = 0)
+wait_beam()
+set_file('Co_plus')
+scan_e_mult([(760,774,1),(774.5,783,0.1),(783.25,790,0.25),(790.2,798,0.2),(798.5,840,1)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cminus(offset = 2)
+wait_beam()
+set_file('Co_minus')
+scan_e_mult([(760,774,1),(774.5,783,0.1),(783.25,790,0.25),(790.2,798,0.2),(798.5,840,1)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cminus(offset = 2)
+wait_beam()
+set_file('Co_minus')
+scan_e_mult([(760,774,1),(774.5,783,0.1),(783.25,790,0.25),(790.2,798,0.2),(798.5,840,1)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cplus(offset = 0)
+wait_beam()
+set_file('Co_plus')
+scan_e_mult([(760,774,1),(774.5,783,0.1),(783.25,790,0.25),(790.2,798,0.2),(798.5,840,1)], settling_time = 1, accumulation_time = 3)