diff --git a/Applications/Utils/Tests.cmake b/Applications/Utils/Tests.cmake
index b7dcc2304814f001d1b2731451d3ef8025faa6c5..c0e8e0db9c082f0a85d3d92cf9763f23d72eee0e 100644
--- a/Applications/Utils/Tests.cmake
+++ b/Applications/Utils/Tests.cmake
@@ -602,3 +602,15 @@ AddTest(
     cube_1x1x1_hex_1e3_1.vtu cube_1x1x1_hex_1e3_1.vtu 1e-16
     cube_1x1x1_hex_1e3_2.vtu cube_1x1x1_hex_1e3_2.vtu 1e-16
 )
+
+AddTest(
+    NAME geometryToGmshGeo
+    PATH MeshGeoToolsLib/geometryToGmshGeo/
+    WORKING_DIRECTORY ${Data_SOURCE_DIR}/MeshGeoToolsLib/geometryToGmshGeo
+    EXECUTABLE geometryToGmshGeo
+    EXECUTABLE_ARGS -i square_1x1.gml -o ${Data_BINARY_DIR}/MeshGeoToolsLib/geometryToGmshGeo/square_1x1.geo
+    REQUIREMENTS NOT OGS_USE_MPI
+    TESTER diff
+    DIFF_DATA
+    square_1x1.geo
+)
diff --git a/Tests/Data/MeshGeoToolsLib/geometryToGmshGeo/square_1x1.geo b/Tests/Data/MeshGeoToolsLib/geometryToGmshGeo/square_1x1.geo
new file mode 100644
index 0000000000000000000000000000000000000000..b46deb659b94c95c51660157cd25060be65fa8ba
--- /dev/null
+++ b/Tests/Data/MeshGeoToolsLib/geometryToGmshGeo/square_1x1.geo
@@ -0,0 +1,12 @@
+// GMSH input file created by OpenGeoSys 6.4.0-427-gdca679359.dirty
+
+Point(0) = {0, 0, 0, 0.2};
+Point(1) = {0, 1, 0, 0.2};
+Point(2) = {1, 0, 0, 0.2};
+Point(3) = {1, 1, 0, 0.2};
+Line(0) = {0,1};
+Line(1) = {1,3};
+Line(2) = {3,2};
+Line(3) = {2,0};
+Line Loop(4) = {0,1,2,3};
+Plane Surface(0) = {4};
diff --git a/Tests/Data/MeshGeoToolsLib/geometryToGmshGeo/square_1x1.gml b/Tests/Data/MeshGeoToolsLib/geometryToGmshGeo/square_1x1.gml
new file mode 100644
index 0000000000000000000000000000000000000000..76f0ac21eba0702028ff3d1b01ae249c86d5e4e4
--- /dev/null
+++ b/Tests/Data/MeshGeoToolsLib/geometryToGmshGeo/square_1x1.gml
@@ -0,0 +1,38 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml-stylesheet type="text/xsl" href="OpenGeoSysGLI.xsl"?>
+
+<OpenGeoSysGLI xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ogs="http://www.opengeosys.org">
+    <name>square_1x1_geometry</name>
+    <points>
+        <point id="0" x="0" y="0" z="0"/>
+        <point id="1" x="0" y="1" z="0"/>
+        <point id="2" x="1" y="0" z="0"/>
+        <point id="3" x="1" y="1" z="0"/>
+    </points>
+
+    <polylines>
+        <polyline id="0" name="left">
+            <pnt>0</pnt>
+            <pnt>1</pnt>
+        </polyline>
+        <polyline id="1" name="right">
+            <pnt>2</pnt>
+            <pnt>3</pnt>
+        </polyline>
+        <polyline id="2" name="bottom">
+            <pnt>0</pnt>
+            <pnt>2</pnt>
+        </polyline>
+        <polyline id="3" name="top">
+            <pnt>1</pnt>
+            <pnt>3</pnt>
+        </polyline>
+        <polyline id="4" name="boundary">
+            <pnt>0</pnt>
+            <pnt>1</pnt>
+            <pnt>3</pnt>
+            <pnt>2</pnt>
+            <pnt>0</pnt>
+        </polyline>
+    </polylines>
+</OpenGeoSysGLI>