diff --git a/parastell/cubit_io.py b/parastell/cubit_io.py
index 4a984ad..0e0a2bf 100644
--- a/parastell/cubit_io.py
+++ b/parastell/cubit_io.py
@@ -205,3 +205,14 @@ def export_dagmc_cubit_native(
     # exports
     if delete_upon_export:
         cubit.cmd(f"delete mesh volume all propagate")
+
+
+def create_new_cubit_instance():
+    """Creates new cubit instance checking if Cubit has been already
+    initialized.
+
+    """
+    if initialized:
+        cubit.cmd("new")
+    else:
+        init_cubit()