diff --git a/lib/io/gmsh_io.h b/lib/io/gmsh_io.h
index 7827e6b..99cb155 100644
--- a/lib/io/gmsh_io.h
+++ b/lib/io/gmsh_io.h
@@ -101,6 +101,7 @@ namespace gctl
 		template <typename E, typename A> void save_mesh(const array<type_block<E>> &element, const array<vertex<point3dc, A>> &node);
 		template <typename E, typename A> void save_mesh(const array<type_tricone<E>> &element, const array<vertex<point3dc, A>> &node);
 		template <typename E, typename A> void save_mesh(const array<type_prism<E>> &element, const array<vertex<point3dc, A>> &node);
+		template <typename E, typename A> void save_mesh(const array<type_tesseroid<E>> &element, const array<vertex<point3dc, A>> &node);
 
 		template <typename E> void save_mesh(const array<type_tetrahedron<E>> &element);
 
@@ -410,6 +411,14 @@ namespace gctl
 		return;
 	}
 
+	template <typename E, typename A>
+	void gmshio::save_mesh(const array<type_tesseroid<E>> &element, const array<vertex<point3dc, A>> &node)
+	{
+		initialized(Output);
+		save2gmsh(gmsh_out, element, node, out_packed);
+		return;
+	}
+
 	template <typename E>
 	void gmshio::save_mesh(const array<type_tetrahedron<E>> &element)
 	{