diff --git a/wk/test_ExB_shear.m b/wk/test_ExB_shear.m
new file mode 100644
index 0000000000000000000000000000000000000000..0e203aa6d9b657c65649da6973035fd8745da576
--- /dev/null
+++ b/wk/test_ExB_shear.m
@@ -0,0 +1,17 @@
+resdir = '/Users/ahoffmann/gyacomo/results/dev/ExB_SF/shearing_picture/';
+MVIN  = ['cd ',resdir];
+% RUN   = './gyacomo23_dp 0';
+% Picture, const. zonal mode in phi in nonlin term
+% RUN   = ['/opt/homebrew/bin/mpirun -np 6 ./gyacomo23_dp 1 6 1 0'];
+% Picture, bckg ExB shear
+% RUN   = '/opt/homebrew/bin/mpirun -np 6 ./gyacomo23_dp 1 6 1 1';
+% RUN   = '/opt/homebrew/bin/mpirun -np 1 ./gyacomo23_dp 1 1 1 1';
+% Test with parallel modes initialization (should not trigger nonlinear
+% coupling)
+% RUN   = '/opt/homebrew/bin/mpirun -np 6 ./gyacomo23_dp 1 6 1 3';
+RUN   = '/opt/homebrew/bin/mpirun -np 1 ./gyacomo23_dp 1 1 1 3';
+% RUN   = '/opt/homebrew/bin/mpirun -np 1 ./gyacomo23_debug 1 1 1 3';
+MVOUT = 'cd /Users/ahoffmann/gyacomo/wk';
+system([MVIN,'; ',RUN,'; ',MVOUT]);
+%%
+fast_analysis