diff --git a/gprMax/input_cmd_funcs.py b/gprMax/input_cmd_funcs.py
index 6f9a7793..5f13d923 100644
--- a/gprMax/input_cmd_funcs.py
+++ b/gprMax/input_cmd_funcs.py
@@ -695,3 +695,20 @@ def rx_steps(dx=0, dy=0, dz=0):
     c = Coordinate(dx, dy, dz)
     command('rx_steps', str(c))
     return c
+
+def geometry_objects_read(x, y, z, file1, file2):
+    """Prints the #geometry_objects_read command.
+
+    Args:
+        x y z are the lower left (x,y,z) coordinates in the domain where the lower left corner of the geometry array should be placed.
+		file1 is the path to and filename of the HDF5 file that contains an integer array which defines the geometry.
+		file2 is the path to and filename of the text file that contains #material commands.
+		not used: c1 is an optional parameter which can be y or n, used to switch on and off dielectric smoothing. Dielectric smoothing can only be turned on if the geometry objects that are being read were originally generated by gprMax, i.e. via the #geometry_objects_write command.
+
+    Returns:
+        coordinates (tuple): namedtuple Coordinate in the domain where the lower left corner of the geometry array is placed.
+    """
+
+    c = Coordinate(x, y, z)
+    command('geometry_objects_read', str(c), file1, file2)
+    return c