diff --git a/script/Cr_xmcd_step_pshell.py b/script/Cr_xmcd_step_pshell.py
new file mode 100644
index 0000000..4515882
--- /dev/null
+++ b/script/Cr_xmcd_step_pshell.py
@@ -0,0 +1,41 @@
+
+set_pol_cplus(offset = -1)
+wait_beam()
+set_file('Cr_plus')
+scan_e_mult([(560,573,1),(573.1,589,0.1),(589.5,625,0.5)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cminus(offset = 0.5)
+wait_beam()
+set_file('Cr_minus')
+scan_e_mult([(560,573,1),(573.1,589,0.1),(589.5,625,0.5)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cminus(offset = 0.5)
+wait_beam()
+set_file('Cr_minus')
+scan_e_mult([(560,573,1),(573.1,589,0.1),(589.5,625,0.5)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cplus(offset = -1)
+wait_beam()
+set_file('Cr_plus')
+scan_e_mult([(560,573,1),(573.1,589,0.1),(589.5,625,0.5)], settling_time = 1, accumulation_time = 3)
+
+
+set_pol_cplus(offset = -1)
+wait_beam()
+set_file('Cr_plus')
+scan_e_mult([(560,573,1),(573.1,589,0.1),(589.5,625,0.5)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cminus(offset = 0.5)
+wait_beam()
+set_file('Cr_minus')
+scan_e_mult([(560,573,1),(573.1,589,0.1),(589.5,625,0.5)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cminus(offset = 0.5)
+wait_beam()
+set_file('Cr_minus')
+scan_e_mult([(560,573,1),(573.1,589,0.1),(589.5,625,0.5)], settling_time = 1, accumulation_time = 3)
+
+set_pol_cplus(offset = -1)
+wait_beam()
+set_file('Cr_plus')
+scan_e_mult([(560,573,1),(573.1,589,0.1),(589.5,625,0.5)], settling_time = 1, accumulation_time = 3)