set up layered soil simulations

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")
 

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")
 

Two-phase-Richards/multi-patch/layered_soil/mesh_study/TP-R-layered_soil.py → 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