diff --git a/Richards-Richards/multi-patch/RR-multi-patch-plus-gravity/run-simulation b/Richards-Richards/multi-patch/RR-multi-patch-plus-gravity/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Richards-Richards/multi-patch/RR-multi-patch-plus-gravity/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Richards-Richards/multi-patch/RR-multi-patch-with-inner-patch/run-simulation b/Richards-Richards/multi-patch/RR-multi-patch-with-inner-patch/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Richards-Richards/multi-patch/RR-multi-patch-with-inner-patch/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Richards-Richards/two-patch/RR-two-patch-symmetric-analytic-soltion/run-simulation b/Richards-Richards/two-patch/RR-two-patch-symmetric-analytic-soltion/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Richards-Richards/two-patch/RR-two-patch-symmetric-analytic-soltion/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Richards-Richards/two-patch/RR-two-patch-test-case/run-simulation b/Richards-Richards/two-patch/RR-two-patch-test-case/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Richards-Richards/two-patch/RR-two-patch-test-case/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Richards/multi-patch/TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting/run-simulation b/Two-phase-Richards/multi-patch/TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Richards/multi-patch/TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/run-simulation b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/run-simulation b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-const-solution/run-simulation b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-const-solution/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-const-solution/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch-constant-solution/run-simulation b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch-constant-solution/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch-constant-solution/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/mesh_study/run-simulation b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/mesh_study/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/mesh_study/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/run-simulation b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/run-simulation b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/run-simulation b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/multi-patch/TP-multi-patch-plus-gravity-with-same-wetting-phase-as-RR/run-simulation b/Two-phase-Two-phase/multi-patch/TP-multi-patch-plus-gravity-with-same-wetting-phase-as-RR/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/multi-patch/TP-multi-patch-plus-gravity-with-same-wetting-phase-as-RR/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/one-patch/TP-one-patch/mesh_study/run-simulation b/Two-phase-Two-phase/one-patch/TP-one-patch/mesh_study/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/one-patch/TP-one-patch/mesh_study/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/one-patch/TP-one-patch/mesh_study_for_fixed_timestep/run-simulation b/Two-phase-Two-phase/one-patch/TP-one-patch/mesh_study_for_fixed_timestep/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/one-patch/TP-one-patch/mesh_study_for_fixed_timestep/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/one-patch/TP-one-patch/run-simulation b/Two-phase-Two-phase/one-patch/TP-one-patch/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/one-patch/TP-one-patch/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/two-patch/injection/run-simulation b/Two-phase-Two-phase/two-patch/injection/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/two-patch/injection/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/two-patch/mesh_studies/run-simulation b/Two-phase-Two-phase/two-patch/mesh_studies/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/two-patch/mesh_studies/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE
diff --git a/Two-phase-Two-phase/two-patch/run-simulation b/Two-phase-Two-phase/two-patch/run-simulation
new file mode 100755
index 0000000000000000000000000000000000000000..0eb497502a082a0fec07a5449b1fe946d59c8cc7
--- /dev/null
+++ b/Two-phase-Two-phase/two-patch/run-simulation
@@ -0,0 +1,16 @@
+#!/bin/bash
+
+[ $# -eq 0 ] && { echo "Usage: $0 simulation_file [logfile_name]"; exit 1; }
+
+SIMULATION_FILE=$1
+SIMULATION=${SIMULATION_FILE%.py}
+LOGFILE_DEFAULT="$SIMULATION.log"
+
+DATE=$(date -I)
+LOGFILE=${2:-$DATE-$LOGFILE_DEFAULT}
+
+GREETING="Simulation $SIMULATION is run on $DATE by $USER"
+
+echo $GREETING
+echo "running $SIMULATION_FILE | tee $LOGFILE"
+./$SIMULATION_FILE | tee $LOGFILE