diff --git a/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys1_mc_moab.sh b/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys1_mc_moab.sh
old mode 100644
new mode 100755
diff --git a/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys2_bromide_mc_moab.sh b/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys2_bromide_mc_moab.sh
new file mode 100755
index 00000000..7600458a
--- /dev/null
+++ b/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys2_bromide_mc_moab.sh
@@ -0,0 +1,20 @@
+#!/bin/bash
+#PBS -l nodes=1:ppn=1
+#PBS -l walltime=24:00:00
+#PBS -l pmem=8000mb
+#PBS -N svat_lys2_bromide_mc
+#PBS -m bea
+#PBS -M robin.schwemmle@hydrology.uni-freiburg.de
+ 
+eval "$(conda shell.bash hook)"
+conda activate roger
+cd /home/fr/fr_fr/fr_rs1092/roger/examples/plot_scale/reckenholz/svat_monte_carlo
+ 
+# adapt command to your available scheduler / MPI implementation
+python svat.py -b numpy -d cpu -lys lys2_bromide -td "${TMPDIR}"
+# Write output to temporary SSD of computing node
+echo "Write output to $TMPDIR"
+# Move output from temporary SSD to workspace
+echo "Move output to /beegfs/work/workspace/ws/fr_rs1092-workspace-0/reckenholz/svat_monte_carlo"
+mkdir -p /beegfs/work/workspace/ws/fr_rs1092-workspace-0/reckenholz/svat_monte_carlo
+mv "${TMPDIR}"/*.nc /beegfs/work/workspace/ws/fr_rs1092-workspace-0/reckenholz/svat_monte_carlo
diff --git a/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys2_mc_moab.sh b/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys2_mc_moab.sh
old mode 100644
new mode 100755
diff --git a/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys3_mc_moab.sh b/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys3_mc_moab.sh
old mode 100644
new mode 100755
diff --git a/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys4_mc_moab.sh b/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys4_mc_moab.sh
old mode 100644
new mode 100755
diff --git a/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys8_bromide_mc_moab.sh b/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys8_bromide_mc_moab.sh
new file mode 100755
index 00000000..cb1e1a4b
--- /dev/null
+++ b/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys8_bromide_mc_moab.sh
@@ -0,0 +1,20 @@
+#!/bin/bash
+#PBS -l nodes=1:ppn=1
+#PBS -l walltime=24:00:00
+#PBS -l pmem=8000mb
+#PBS -N svat_lys8_bromide_mc
+#PBS -m bea
+#PBS -M robin.schwemmle@hydrology.uni-freiburg.de
+ 
+eval "$(conda shell.bash hook)"
+conda activate roger
+cd /home/fr/fr_fr/fr_rs1092/roger/examples/plot_scale/reckenholz/svat_monte_carlo
+ 
+# adapt command to your available scheduler / MPI implementation
+python svat.py -b numpy -d cpu -lys lys8_bromide -td "${TMPDIR}"
+# Write output to temporary SSD of computing node
+echo "Write output to $TMPDIR"
+# Move output from temporary SSD to workspace
+echo "Move output to /beegfs/work/workspace/ws/fr_rs1092-workspace-0/reckenholz/svat_monte_carlo"
+mkdir -p /beegfs/work/workspace/ws/fr_rs1092-workspace-0/reckenholz/svat_monte_carlo
+mv "${TMPDIR}"/*.nc /beegfs/work/workspace/ws/fr_rs1092-workspace-0/reckenholz/svat_monte_carlo
diff --git a/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys8_mc_moab.sh b/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys8_mc_moab.sh
old mode 100644
new mode 100755
diff --git a/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys9_bromide_mc_moab.sh b/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys9_bromide_mc_moab.sh
new file mode 100755
index 00000000..5b49923b
--- /dev/null
+++ b/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys9_bromide_mc_moab.sh
@@ -0,0 +1,20 @@
+#!/bin/bash
+#PBS -l nodes=1:ppn=1
+#PBS -l walltime=24:00:00
+#PBS -l pmem=8000mb
+#PBS -N svat_lys9_bromide_mc
+#PBS -m bea
+#PBS -M robin.schwemmle@hydrology.uni-freiburg.de
+ 
+eval "$(conda shell.bash hook)"
+conda activate roger
+cd /home/fr/fr_fr/fr_rs1092/roger/examples/plot_scale/reckenholz/svat_monte_carlo
+ 
+# adapt command to your available scheduler / MPI implementation
+python svat.py -b numpy -d cpu -lys lys9_bromide -td "${TMPDIR}"
+# Write output to temporary SSD of computing node
+echo "Write output to $TMPDIR"
+# Move output from temporary SSD to workspace
+echo "Move output to /beegfs/work/workspace/ws/fr_rs1092-workspace-0/reckenholz/svat_monte_carlo"
+mkdir -p /beegfs/work/workspace/ws/fr_rs1092-workspace-0/reckenholz/svat_monte_carlo
+mv "${TMPDIR}"/*.nc /beegfs/work/workspace/ws/fr_rs1092-workspace-0/reckenholz/svat_monte_carlo
diff --git a/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys9_mc_moab.sh b/examples/plot_scale/reckenholz/svat_monte_carlo/svat_lys9_mc_moab.sh
old mode 100644
new mode 100755
diff --git a/examples/plot_scale/reckenholz/svat_monte_carlo/write_moab_jobs.py b/examples/plot_scale/reckenholz/svat_monte_carlo/write_moab_jobs.py
index a252137e..294da63a 100644
--- a/examples/plot_scale/reckenholz/svat_monte_carlo/write_moab_jobs.py
+++ b/examples/plot_scale/reckenholz/svat_monte_carlo/write_moab_jobs.py
@@ -9,7 +9,8 @@ def main(job_type):
     base_path = Path(__file__).parent
     base_path_binac = '/home/fr/fr_fr/fr_rs1092/roger/examples/plot_scale/reckenholz/svat_monte_carlo'
     base_path_ws = Path('/beegfs/work/workspace/ws/fr_rs1092-workspace-0')
-    lysimeters = ['lys1', 'lys2', 'lys3', 'lys4', 'lys8', 'lys9']
+    lysimeters = ['lys1', 'lys2', 'lys3', 'lys4', 'lys8', 'lys9', 'lys2_bromide',
+                  'lys8_bromide', 'lys9_bromide']
     for lys in lysimeters:
         script_name = f'svat_{lys}_mc'
         output_path_ws = base_path_ws / 'reckenholz' / 'svat_monte_carlo'
@@ -38,7 +39,7 @@ def main(job_type):
         file = open(file_path, "w")
         file.writelines(lines)
         file.close()
-        subprocess.Popen(f"chmod +x {file_path}.sh", shell=True)
+        subprocess.Popen(f"chmod +x {file_path}", shell=True)
 
     return