diff --git a/models/bonemodel.py b/models/bonemodel.py
index dadbd6f5406c5c2e5165aa9bff1c73e845dc9718..8222c07fb8883980007975d65074707499a0d5de 100644
--- a/models/bonemodel.py
+++ b/models/bonemodel.py
@@ -68,6 +68,9 @@ def parms(d={}):
     p['tolNR'] = 1e-6        # [-] equilibrium tolerance
     p['dt0'] = 1.0           # [s] time step size
 
+    # gmsh toolbox
+    p['use_gmshOld'] = True  # use old gmsh interface
+
     p.update(d)
     return p
 
@@ -283,8 +286,10 @@ def import_mesh(domain, filename):
         raise Exception(f'Unknown extension: {ext}, please use .ply, .stl, .msh or .geo')
 
     print(f'importing {mandible_gmsh}')
-    # from toolbox.gmsh import GmshImport # FIXME: does not work!
-    from toolbox.gmshOld import GmshImport
+    if p['use_gmshOld']:
+        from toolbox.gmshOld import GmshImport
+    else:
+        from toolbox.gmsh import GmshImport # FIXME: does not work!
     f = os.path.join(os.path.dirname(__file__), mandible_gmsh)
     mesher = GmshImport(f, domain)
     mesher.execute()