diff --git a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-all-params-one.py b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-all-params-one.py
index cb44482a0b31d6f7a4414d4badc6de0689e6c603..60f94ca84baa2b10b6bf89720047e30776f88f8b 100755
--- a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-all-params-one.py
+++ b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-all-params-one.py
@@ -17,6 +17,15 @@ import datetime
 import os
 import pandas as pd
 
+# check if output directory exists
+if not os.path.exists('./output'):
+    os.mkdir('./output')
+    print("Directory ", './output',  " created ")
+else:
+    print("Directory ", './output',  " already exists. Will use as output \
+    directory")
+
+
 date = datetime.datetime.now()
 datestr = date.strftime("%Y-%m-%d")
 
diff --git a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
index c05a2447c7c35a0e7255be2d3c7e13a2266c8408..8f3f087cf35cdc105021421c9bb247ee76122342 100755
--- a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
+++ b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
@@ -17,6 +17,15 @@ import datetime
 import os
 import pandas as pd
 
+# check if output directory exists
+if not os.path.exists('./output'):
+    os.mkdir('./output')
+    print("Directory ", './output',  " created ")
+else:
+    print("Directory ", './output',  " already exists. Will use as output \
+    directory")
+
+
 date = datetime.datetime.now()
 datestr = date.strftime("%Y-%m-%d")
 
diff --git a/Two-phase-Richards/multi-patch/layered_soil/mesh_study/TP-R-layered_soil.py b/Two-phase-Richards/multi-patch/layered_soil/mesh_study/TP-R-layered_soil-mesh-study.py
similarity index 98%
rename from Two-phase-Richards/multi-patch/layered_soil/mesh_study/TP-R-layered_soil.py
rename to Two-phase-Richards/multi-patch/layered_soil/mesh_study/TP-R-layered_soil-mesh-study.py
index 1d56fc1f6dabd48340891ee315ddd614306df6a4..0e9c966abbb7d08ade59ca7acbb95799587aa990 100755
--- a/Two-phase-Richards/multi-patch/layered_soil/mesh_study/TP-R-layered_soil.py
+++ b/Two-phase-Richards/multi-patch/layered_soil/mesh_study/TP-R-layered_soil-mesh-study.py
@@ -17,6 +17,14 @@ import datetime
 import os
 import pandas as pd
 
+# check if output directory exists
+if not os.path.exists('./output'):
+    os.mkdir('./output')
+    print("Directory ", './output',  " created ")
+else:
+    print("Directory ", './output',  " already exists. Will use as output \
+    directory")
+
 date = datetime.datetime.now()
 datestr = date.strftime("%Y-%m-%d")
 
@@ -25,7 +33,8 @@ sym.init_printing()
 # solver_tol = 6E-7
 use_case = "TP-R-layered-soil-realistic"
 # name of this very file. Needed for log output.
-thisfile = "TP-R-layered_soil.py"
+thisfile = "TP-R-layered_soil-mesh-study.py"
+
 max_iter_num = 300
 FEM_Lagrange_degree = 1
 mesh_study = True
@@ -352,8 +361,8 @@ lambda_param = {
 intrinsic_permeability = {
     1: 0.1,  # sand
     2: 0.1,  # sand, there is a range
-    3: 0.001,  #10e-2,  # clay has a range
-    4: 0.001,  #10e-3
+    3: 0.001,  # 10e-2,  # clay has a range
+    4: 0.001,  # 10e-3
 }