diff --git a/user_models/Mont_Blanc.py b/user_models/Mont_Blanc.py
new file mode 100644
index 00000000..2b87f621
--- /dev/null
+++ b/user_models/Mont_Blanc.py
@@ -0,0 +1,34 @@
+from pathlib import Path
+
+import gprMax
+
+fn = Path(__file__)
+
+title = gprMax.Title(name=fn.with_suffix('').name)
+domain = gprMax.Domain(p1=(0.4,0.4 ,0.12))
+dxdydz = gprMax.Discretisation(p1=(0.002, 0.002, 0.002))
+time_window = gprMax.TimeWindow(time=10e-9)
+
+g_stl = gprMax.GeometryObjectsReadSTL(stl_file='data_LasoMax/Mont_Blanc.stl', mat_index=2, discretization=(0.002,0.002,0.002))
+
+g_read = gprMax.GeometryObjectsRead(p1=(0.04,0.04,0.04), geofile=fn.with_suffix('.h5') , matfile='user_models/materials.txt')
+
+gv = gprMax.GeometryView(p1=(0, 0, 0),
+                         p2=(0.4, 0.4 , 0.12),
+                         dl=(0.002, 0.002, 0.002),
+                         filename=fn.with_suffix('').name,
+                         output_type='n')
+                         
+# create a scene
+scene = gprMax.Scene()
+# add the simulation objects to the scene
+scene.add(title)
+scene.add(domain)
+scene.add(dxdydz)
+scene.add(time_window)
+scene.add(g_stl)
+scene.add(g_read)
+scene.add(gv)
+
+# run the simulation
+gprMax.run(scenes=[scene], geometry_only=True, n=1, outputfile=fn)
diff --git a/user_models/amphora.py b/user_models/amphora.py
new file mode 100644
index 00000000..ef07ac2a
--- /dev/null
+++ b/user_models/amphora.py
@@ -0,0 +1,34 @@
+from pathlib import Path
+
+import gprMax
+
+fn = Path(__file__)
+
+title = gprMax.Title(name=fn.with_suffix('').name)
+domain = gprMax.Domain(p1=(0.08,0.07,0.11))
+dxdydz = gprMax.Discretisation(p1=(0.001, 0.001, 0.001))
+time_window = gprMax.TimeWindow(time=10e-9)
+
+g_stl = gprMax.GeometryObjectsReadSTL(stl_file='data_LasoMax/amphora.stl', mat_index=2, discretization =(0.001,0.001,0.001))
+
+g_read = gprMax.GeometryObjectsRead(p1=(0.001,0.001,0.0), geofile=fn.with_suffix('.h5') , matfile= 'user_models/materials.txt')
+
+gv = gprMax.GeometryView(p1=(0, 0, 0),
+                         p2=(0.08, 0.07,0.11),
+                         dl=(0.001, 0.001, 0.001),
+                         filename=fn.with_suffix('').name,
+                         output_type='n')
+                         
+# create a scene
+scene = gprMax.Scene()
+# add the simulation objects to the scene
+scene.add(title)
+scene.add(domain)
+scene.add(dxdydz)
+scene.add(time_window)
+scene.add(g_stl)
+scene.add(g_read)
+scene.add(gv)
+
+# run the simulation
+gprMax.run(scenes=[scene], geometry_only=True, n=1, outputfile=fn)
diff --git a/user_models/bunny.py b/user_models/bunny.py
new file mode 100644
index 00000000..1eb824cb
--- /dev/null
+++ b/user_models/bunny.py
@@ -0,0 +1,34 @@
+from pathlib import Path
+
+import gprMax
+
+fn = Path(__file__)
+
+title = gprMax.Title(name=fn.with_suffix('').name)
+domain = gprMax.Domain(p1=(0.12,0.1,0.12))
+dxdydz = gprMax.Discretisation(p1=(0.001, 0.001, 0.001))
+time_window = gprMax.TimeWindow(time=10e-9)
+
+g_stl = gprMax.GeometryObjectsReadSTL(stl_file='data_LasoMax/Stanford_Bunny.stl', mat_index=2, discretization=(0.001,0.001,0.001) )
+
+g_read = gprMax.GeometryObjectsRead(p1=(0 ,0 ,0), geofile=fn.with_suffix('.h5') , matfile= 'user_models/materials.txt')
+
+gv = gprMax.GeometryView(p1=(0, 0, 0),
+                         p2=(0.12, 0.1,0.12),
+                         dl=(0.001, 0.001, 0.001),
+                         filename=fn.with_suffix('').name,
+                         output_type='n')
+                         
+# create a scene
+scene = gprMax.Scene()
+# add the simulation objects to the scene
+scene.add(title)
+scene.add(domain)
+scene.add(dxdydz)
+scene.add(time_window)
+scene.add(g_stl)
+scene.add(g_read)
+scene.add(gv)
+
+# run the simulation
+gprMax.run(scenes=[scene], geometry_only=True, n=1, outputfile=fn)
diff --git a/user_models/dubai.py b/user_models/dubai.py
new file mode 100644
index 00000000..4b7127f0
--- /dev/null
+++ b/user_models/dubai.py
@@ -0,0 +1,34 @@
+from pathlib import Path
+
+import gprMax
+
+fn = Path(__file__)
+
+title = gprMax.Title(name=fn.with_suffix('').name)
+domain = gprMax.Domain(p1=(0.14,0.12,0.12))
+dxdydz = gprMax.Discretisation(p1=(0.0008, 0.0008, 0.0008))
+time_window = gprMax.TimeWindow(time=10e-9)
+
+g_stl = gprMax.GeometryObjectsReadSTL(stl_file='data_LasoMax/Dubai.stl', mat_index=1, discretization=(0.0008,0.0008,0.0008) )
+
+g_read = gprMax.GeometryObjectsRead(p1=(0.01 ,0.01,0.01), geofile=fn.with_suffix('.h5') , matfile= 'user_models/materials.txt')
+
+gv = gprMax.GeometryView(p1=(0, 0, 0),
+                         p2=(0.14, 0.12,0.12),
+                         dl=(0.0008, 0.0008, 0.0008),
+                         filename=fn.with_suffix('').name,
+                         output_type='n')
+                         
+# create a scene
+scene = gprMax.Scene()
+# add the simulation objects to the scene
+scene.add(title)
+scene.add(domain)
+scene.add(dxdydz)
+scene.add(time_window)
+scene.add(g_stl)
+scene.add(g_read)
+scene.add(gv)
+
+# run the simulation
+gprMax.run(scenes=[scene], geometry_only=True, n=1, outputfile=fn)
diff --git a/user_models/house.py b/user_models/house.py
new file mode 100644
index 00000000..6f0693b4
--- /dev/null
+++ b/user_models/house.py
@@ -0,0 +1,34 @@
+from pathlib import Path
+
+import gprMax
+
+fn = Path(__file__)
+
+title = gprMax.Title(name=fn.with_suffix('').name)
+domain = gprMax.Domain(p1=(0.112,0.104,0.082))
+dxdydz = gprMax.Discretisation(p1=(0.001, 0.001, 0.001))
+time_window = gprMax.TimeWindow(time=10e-9)
+
+g_stl = gprMax.GeometryObjectsReadSTL(stl_file='data_LasoMax/house.stl', mat_index=2,discretization=(0.001,0.001,0.001))
+
+g_read = gprMax.GeometryObjectsRead(p1=(0,0,0), geofile=fn.with_suffix('.h5') , matfile= 'user_models/materials.txt')
+
+gv = gprMax.GeometryView(p1=(0, 0, 0),
+                         p2=(0.112, 0.104,0.082),
+                         dl=(0.001, 0.001, 0.001),
+                         filename=fn.with_suffix('').name,
+                         output_type='n')
+                         
+# create a scene
+scene = gprMax.Scene()
+# add the simulation objects to the scene
+scene.add(title)
+scene.add(domain)
+scene.add(dxdydz)
+scene.add(time_window)
+scene.add(g_stl)
+scene.add(g_read)
+scene.add(gv)
+
+# run the simulation
+gprMax.run(scenes=[scene], geometry_only=True, n=1, outputfile=fn)
diff --git a/user_models/materials.txt b/user_models/materials.txt
new file mode 100644
index 00000000..53964258
--- /dev/null
+++ b/user_models/materials.txt
@@ -0,0 +1,4 @@
+#material: 8 0 1 0 pseudo_mat
+#material: 9 0 1 0 mat1
+#material: 8.5 0 1 0 mat2
+#material: 7 0 1 0 mat3
\ No newline at end of file
diff --git a/user_models/newyorkcity.py b/user_models/newyorkcity.py
new file mode 100644
index 00000000..b03d3013
--- /dev/null
+++ b/user_models/newyorkcity.py
@@ -0,0 +1,34 @@
+from pathlib import Path
+
+import gprMax
+
+fn = Path(__file__)
+
+title = gprMax.Title(name=fn.with_suffix('').name)
+domain = gprMax.Domain(p1=(0.11,0.115,0.07))
+dxdydz = gprMax.Discretisation(p1=(0.001, 0.001, 0.001))
+time_window = gprMax.TimeWindow(time=10e-9)
+
+g_stl = gprMax.GeometryObjectsReadSTL(stl_file='data_LasoMax/newyorkcity.stl', mat_index=2, discretization=(0.001,0.001,0.001))
+
+g_read = gprMax.GeometryObjectsRead(p1=(0.001,0.001,0.0), geofile=fn.with_suffix('.h5') , matfile= 'user_models/materials.txt')
+
+gv = gprMax.GeometryView(p1=(0, 0, 0),
+                         p2=(0.11, 0.115,0.07),
+                         dl=(0.001, 0.001, 0.001),
+                         filename=fn.with_suffix('').name,
+                         output_type='n')
+
+# create a scene
+scene = gprMax.Scene()
+# add the simulation objects to the scene
+scene.add(title)
+scene.add(domain)
+scene.add(dxdydz)
+scene.add(time_window)
+scene.add(g_stl)
+scene.add(g_read)
+scene.add(gv)
+
+# run the simulation
+gprMax.run(scenes=[scene], geometry_only=True, n=1, outputfile=fn)
\ No newline at end of file