diff --git a/models/CLIpro/CLIpro_150k_15_1.py b/models/CLIpro/CLIpro_150k_15_1.py
index 066953944ba7d4f36e7be07f4594a4e3060b0e13..6a49f43cd44a1cfbf1f5a9d0cf933230290eab45 100644
--- a/models/CLIpro/CLIpro_150k_15_1.py
+++ b/models/CLIpro/CLIpro_150k_15_1.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/150k'
-    p['mandible'] = f'{path}/CLIpro_jaw_150k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_150k.stl'
     p['contact_pts'] = [[-13.654129, -90.974480, 4.162612],
                 [14.185822, -89.850540, 8.610174]]
     p['axis_pt1'] = [47.751984, 205.253143, 26.247160]
diff --git a/models/CLIpro/CLIpro_150k_15_2.py b/models/CLIpro/CLIpro_150k_15_2.py
index 6b37706afabddc9e7236092f353122fe2e54f2da..e21217fc0ef379dc4f2618f4583bf1598e0cee5d 100644
--- a/models/CLIpro/CLIpro_150k_15_2.py
+++ b/models/CLIpro/CLIpro_150k_15_2.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/150k'
-    p['mandible'] = f'{path}/CLIpro_jaw_150k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_150k.stl'
     p['contact_pts'] = [[-34.115112, 21.823414, 10.233394],
                 [32.383774, 21.721252, 11.580667]]
     p['axis_pt1'] = [47.751984, 205.253143, 26.247160]
diff --git a/models/CLIpro/CLIpro_150k_7_1.py b/models/CLIpro/CLIpro_150k_7_1.py
index a0a3fb50d70aa59cb16e2c6eb557ab80d52c6028..c5b6892863d02e2506a2814c67c5ee570242896e 100644
--- a/models/CLIpro/CLIpro_150k_7_1.py
+++ b/models/CLIpro/CLIpro_150k_7_1.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/150k'
-    p['mandible'] = f'{path}/CLIpro_jaw_150k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_150k.stl'
     p['contact_pts'] = [[-13.654129, -90.974480, 4.162612],
                 [14.185822, -89.850540, 8.610174]]
     p['axis_pt1'] = [47.751984, 205.253143, 26.247160]
diff --git a/models/CLIpro/CLIpro_150k_7_2.py b/models/CLIpro/CLIpro_150k_7_2.py
index 12e59fe517b541dc8055b33386097b774d2178bb..aa2ca36bcbc53e516a3071e637f14ec2c1e2591a 100644
--- a/models/CLIpro/CLIpro_150k_7_2.py
+++ b/models/CLIpro/CLIpro_150k_7_2.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/150k'
-    p['mandible'] = f'{path}/CLIpro_jaw_150k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_150k.stl'
     p['contact_pts'] = [[-34.115112, 21.823414, 10.233394],
                 [32.383774, 21.721252, 11.580667]]
     p['axis_pt1'] = [47.751984, 205.253143, 26.247160]
diff --git a/models/CLIpro/CLIpro_250k_15_1.py b/models/CLIpro/CLIpro_250k_15_1.py
index 2616637a12a30eeea222c20a4cb636eb524a1c2a..6a1c9e75f0adad9b4b1fdb946667d5e01930708e 100644
--- a/models/CLIpro/CLIpro_250k_15_1.py
+++ b/models/CLIpro/CLIpro_250k_15_1.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/250k'
-    p['mandible'] = f'{path}/CLIpro_jaw_250k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_250k.stl'
     p['contact_pts'] = [[-13.921872, -90.258591, 8.862718],
                 [13.158182, -91.377449, 7.355043]]
     p['axis_pt1'] = [47.686401, 205.691238, 26.261269]
diff --git a/models/CLIpro/CLIpro_250k_15_2.py b/models/CLIpro/CLIpro_250k_15_2.py
index 8e00f870673270527681863355ee0b0152cf8a20..68068656e78ecc1ba1dc03b1f9b2a4f876a908fc 100644
--- a/models/CLIpro/CLIpro_250k_15_2.py
+++ b/models/CLIpro/CLIpro_250k_15_2.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/250k'
-    p['mandible'] = f'{path}/CLIpro_jaw_250k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_250k.stl'
     p['contact_pts'] = [[-34.058887, 21.742476, 9.846808],
                 [33.013206, 24.755995, -0.564457]]
     p['axis_pt1'] = [47.686401, 205.691238, 26.261269]
diff --git a/models/CLIpro/CLIpro_250k_7_1.py b/models/CLIpro/CLIpro_250k_7_1.py
index cec986dc831f57e5939bad77caa42254f9970b2f..2d6612ab4bee6d1fa281ff0893d364bd0947fbbb 100644
--- a/models/CLIpro/CLIpro_250k_7_1.py
+++ b/models/CLIpro/CLIpro_250k_7_1.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/250k'
-    p['mandible'] = f'{path}/CLIpro_jaw_250k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_250k.stl'
     p['contact_pts'] = [[-13.921872, -90.258591, 8.862718],
                 [13.158182, -91.377449, 7.355043]]
     p['axis_pt1'] = [47.686401, 205.691238, 26.261269]
diff --git a/models/CLIpro/CLIpro_250k_7_2.py b/models/CLIpro/CLIpro_250k_7_2.py
index c9ad83bad855a54fe41bb1604c3fb4d91917d318..28f37180849bd93292d1166836a57e879a32e44d 100644
--- a/models/CLIpro/CLIpro_250k_7_2.py
+++ b/models/CLIpro/CLIpro_250k_7_2.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/250k'
-    p['mandible'] = f'{path}/CLIpro_jaw_250k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_250k.stl'
     p['contact_pts'] = [[-34.058887, 21.742476, 9.846808],
                 [33.013206, 24.755995, -0.564457]]
     p['axis_pt1'] = [47.686401, 205.691238, 26.261269]
diff --git a/models/CLIpro/CLIpro_500k_15_1.py b/models/CLIpro/CLIpro_500k_15_1.py
index b21de461bf3125b66de4eb4f8d196e68de2c465c..7e7d2cc9a0ac6725d33154ef53656522446ed573 100644
--- a/models/CLIpro/CLIpro_500k_15_1.py
+++ b/models/CLIpro/CLIpro_500k_15_1.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/500k'
-    p['mandible'] = f'{path}/CLIpro_jaw_500k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_500k.stl'
     p['contact_pts'] = [[-12.880455, -92.642342, -0.673375],
                 [12.990344, -91.013168, 7.432779]]
     p['axis_pt1'] = [45.043976, 199.899857, 33.671288]
diff --git a/models/CLIpro/CLIpro_500k_15_2.py b/models/CLIpro/CLIpro_500k_15_2.py
index e1d7ae5d377518b8dc158a4da852d5bd0a472877..4f4e75f77b9354107bfc18da0db376996faacb4b 100644
--- a/models/CLIpro/CLIpro_500k_15_2.py
+++ b/models/CLIpro/CLIpro_500k_15_2.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/500k'
-    p['mandible'] = f'{path}/CLIpro_jaw_500k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_500k.stl'
     p['contact_pts'] = [[-33.575844, 26.282110, -0.445806],
                 [32.137966, 22.818340, 7.904116]]
     p['axis_pt1'] = [45.043976, 199.899857, 33.671288]
diff --git a/models/CLIpro/CLIpro_500k_7_1.py b/models/CLIpro/CLIpro_500k_7_1.py
index cb26d6f0c97f717cc1481b67e8f260769f5c09f7..f2de1674cde70e9e4fe5a98799f1c27b5aea72c7 100644
--- a/models/CLIpro/CLIpro_500k_7_1.py
+++ b/models/CLIpro/CLIpro_500k_7_1.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/500k'
-    p['mandible'] = f'{path}/CLIpro_jaw_500k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_500k.stl'
     p['contact_pts'] = [[-12.880455, -92.642342, -0.673375],
                 [12.990344, -91.013168, 7.432779]]
     p['axis_pt1'] = [45.043976, 199.899857, 33.671288]
diff --git a/models/CLIpro/CLIpro_500k_7_2.py b/models/CLIpro/CLIpro_500k_7_2.py
index 6b1313d74f87dcc0221ed2feec6b725f9f179f64..ba7e03d7d7c6bf762889b21290e2b978962a0431 100644
--- a/models/CLIpro/CLIpro_500k_7_2.py
+++ b/models/CLIpro/CLIpro_500k_7_2.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/500k'
-    p['mandible'] = f'{path}/CLIpro_jaw_500k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_500k.stl'
     p['contact_pts'] = [[-33.575844, 26.282110, -0.445806],
                 [32.137966, 22.818340, 7.904116]]
     p['axis_pt1'] = [45.043976, 199.899857, 33.671288]
diff --git a/models/CLIpro/CLIpro_750k_15_1.py b/models/CLIpro/CLIpro_750k_15_1.py
index ff29c8fdc22a06f2ca0b525b4e6a3b23f8c9067d..3fe58238944b3cc36b16aa892bd62317f9ccf51b 100644
--- a/models/CLIpro/CLIpro_750k_15_1.py
+++ b/models/CLIpro/CLIpro_750k_15_1.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/750k'
-    p['mandible'] = f'{path}/CLIpro_jaw_750k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_750k.stl'
     p['contact_pts'] = [[-14.743377, -90.163612, 9.011174],
                 [13.049115, -90.677956, 7.640814]]
     p['axis_pt1'] = [46.317547, 203.956772, 35.086609]
diff --git a/models/CLIpro/CLIpro_750k_15_2.py b/models/CLIpro/CLIpro_750k_15_2.py
index cabb96f3da0dd8eb328523df93472cc3d150035f..530049e21004a2ed97635df1ec3a50a11326badb 100644
--- a/models/CLIpro/CLIpro_750k_15_2.py
+++ b/models/CLIpro/CLIpro_750k_15_2.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/750k'
-    p['mandible'] = f'{path}/CLIpro_jaw_750k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_750k.stl'
     p['contact_pts'] = [[-34.126896, 21.475969, 10.220996],
                 [32.139095, 22.748108, 8.446575]]
     p['axis_pt1'] = [46.317547, 203.956772, 35.086609]
diff --git a/models/CLIpro/CLIpro_750k_7_1.py b/models/CLIpro/CLIpro_750k_7_1.py
index a870262063d0d7c144ecdc4f6b69a3e3f5dbc4c9..0aa447cd10fc8829842a203e1a476e4b68574caa 100644
--- a/models/CLIpro/CLIpro_750k_7_1.py
+++ b/models/CLIpro/CLIpro_750k_7_1.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/750k'
-    p['mandible'] = f'{path}/CLIpro_jaw_750k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_750k.stl'
     p['contact_pts'] = [[-14.743377, -90.163612, 9.011174],
                 [13.049115, -90.677956, 7.640814]]
     p['axis_pt1'] = [46.317547, 203.956772, 35.086609]
diff --git a/models/CLIpro/CLIpro_750k_7_2.py b/models/CLIpro/CLIpro_750k_7_2.py
index d8e8d963909d6808c914a1d371f7962cc8c8d768..7adbfdc865a4bf0611a023637df11f9b54ae0286 100644
--- a/models/CLIpro/CLIpro_750k_7_2.py
+++ b/models/CLIpro/CLIpro_750k_7_2.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'CLIpro/750k'
-    p['mandible'] = f'{path}/CLIpro_jaw_750k.stl'
+    p['bone'] = f'{path}/CLIpro_jaw_750k.stl'
     p['contact_pts'] = [[-34.126896, 21.475969, 10.220996],
                 [32.139095, 22.748108, 8.446575]]
     p['axis_pt1'] = [46.317547, 203.956772, 35.086609]
diff --git a/models/Dolico_bon/Dolico_bon_150k_15_1.py b/models/Dolico_bon/Dolico_bon_150k_15_1.py
index d7ccc01a38cae33c4ad82ecb91c315bcf8267f78..0c0ef94c577fa1003ddae44cf3ab6bcf1320b86c 100644
--- a/models/Dolico_bon/Dolico_bon_150k_15_1.py
+++ b/models/Dolico_bon/Dolico_bon_150k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/150k'
-    p['mandible'] = f'{path}/Dolico_bon_150k.stl'
+    p['bone'] = f'{path}/Dolico_bon_150k.stl'
     p['contact_pts'] = [[4.212819, -267.574554, 28.792574],
                 [36.007748, -269.322662, 27.255693]]
     p['axis_pt1'] = [162.458069, 598.625000, -20.051981]
diff --git a/models/Dolico_bon/Dolico_bon_150k_15_2.py b/models/Dolico_bon/Dolico_bon_150k_15_2.py
index 13a4e459be1e7421e275bef4a7ab449f92c61f38..a03c73421eefb0149428c766c1d143adc7f7d760 100644
--- a/models/Dolico_bon/Dolico_bon_150k_15_2.py
+++ b/models/Dolico_bon/Dolico_bon_150k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/150k'
-    p['mandible'] = f'{path}/Dolico_bon_150k.stl'
+    p['bone'] = f'{path}/Dolico_bon_150k.stl'
     p['contact_pts'] = [[-14.645494, 62.509026, 23.022108],
                 [44.984268, 67.347466, -2.159538]]
     p['axis_pt1'] = [162.458069, 598.625000, -20.051981]
diff --git a/models/Dolico_bon/Dolico_bon_150k_7_1.py b/models/Dolico_bon/Dolico_bon_150k_7_1.py
index 0a9422e6eb6e22dbbf922a7e4d12dc040b2208db..160ae0c13b415eccc2f2e2e4eb91e2490164477a 100644
--- a/models/Dolico_bon/Dolico_bon_150k_7_1.py
+++ b/models/Dolico_bon/Dolico_bon_150k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/150k'
-    p['mandible'] = f'{path}/Dolico_bon_150k.stl'
+    p['bone'] = f'{path}/Dolico_bon_150k.stl'
     p['contact_pts'] = [[4.212819, -267.574554, 28.792574],
                 [36.007748, -269.322662, 27.255693]]
     p['axis_pt1'] = [162.458069, 598.625000, -20.051981]
diff --git a/models/Dolico_bon/Dolico_bon_150k_7_2.py b/models/Dolico_bon/Dolico_bon_150k_7_2.py
index ed997cc657f629e12d1daf28d5332a7cc5559770..45098943f35b6ae5afcfc0cc894daab2e1300f07 100644
--- a/models/Dolico_bon/Dolico_bon_150k_7_2.py
+++ b/models/Dolico_bon/Dolico_bon_150k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/150k'
-    p['mandible'] = f'{path}/Dolico_bon_150k.stl'
+    p['bone'] = f'{path}/Dolico_bon_150k.stl'
     p['contact_pts'] = [[-14.645494, 62.509026, 23.022108],
                 [44.984268, 67.347466, -2.159538]]
     p['axis_pt1'] = [162.458069, 598.625000, -20.051981]
diff --git a/models/Dolico_bon/Dolico_bon_250k_15_1.py b/models/Dolico_bon/Dolico_bon_250k_15_1.py
index d944e1f2ff9712aee6b85a35ede48879e5694ab3..59c6db2f1bb8c7d263f3bba4a38e2494c0103d71 100644
--- a/models/Dolico_bon/Dolico_bon_250k_15_1.py
+++ b/models/Dolico_bon/Dolico_bon_250k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/250k'
-    p['mandible'] = f'{path}/Dolico_bon_250k.stl'
+    p['bone'] = f'{path}/Dolico_bon_250k.stl'
     p['contact_pts'] = [[4.212819, -267.574554, 28.792574],
                 [36.279472, -269.853271, 26.685869]]
     p['axis_pt1'] = [168.908783, 603.534424, -17.702827]
diff --git a/models/Dolico_bon/Dolico_bon_250k_15_2.py b/models/Dolico_bon/Dolico_bon_250k_15_2.py
index 7636851498a43d62e6db1d5b56e647dab81b29c6..7954cd78e10f821cbcd7f578ffd89801691fb743 100644
--- a/models/Dolico_bon/Dolico_bon_250k_15_2.py
+++ b/models/Dolico_bon/Dolico_bon_250k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/250k'
-    p['mandible'] = f'{path}/Dolico_bon_250k.stl'
+    p['bone'] = f'{path}/Dolico_bon_250k.stl'
     p['contact_pts'] = [[-10.472356, 31.800030, -3.685302],
                 [50.958504, 38.902874, -5.352512]]
     p['axis_pt1'] = [168.908783, 603.534424, -17.702827]
diff --git a/models/Dolico_bon/Dolico_bon_250k_7_1.py b/models/Dolico_bon/Dolico_bon_250k_7_1.py
index 24e16755a82beb91177241e38bbaffb2070e801c..fc2f46ea1ddf48ed97104d974778bd1dc986597f 100644
--- a/models/Dolico_bon/Dolico_bon_250k_7_1.py
+++ b/models/Dolico_bon/Dolico_bon_250k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/250k'
-    p['mandible'] = f'{path}/Dolico_bon_250k.stl'
+    p['bone'] = f'{path}/Dolico_bon_250k.stl'
     p['contact_pts'] = [[4.212819, -267.574554, 28.792574],
                 [36.279472, -269.853271, 26.685869]]
     p['axis_pt1'] = [168.908783, 603.534424, -17.702827]
diff --git a/models/Dolico_bon/Dolico_bon_250k_7_2.py b/models/Dolico_bon/Dolico_bon_250k_7_2.py
index 909dc56de1447490d65a7be7a07ceb4fe7842c6e..1b5a728dbbbee84a124a46ab225ab83a6082d3b5 100644
--- a/models/Dolico_bon/Dolico_bon_250k_7_2.py
+++ b/models/Dolico_bon/Dolico_bon_250k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/250k'
-    p['mandible'] = f'{path}/Dolico_bon_250k.stl'
+    p['bone'] = f'{path}/Dolico_bon_250k.stl'
     p['contact_pts'] = [[-10.472356, 31.800030, -3.685302],
                 [50.958504, 38.902874, -5.352512]]
     p['axis_pt1'] = [168.908783, 603.534424, -17.702827]
diff --git a/models/Dolico_bon/Dolico_bon_500k_15_1.py b/models/Dolico_bon/Dolico_bon_500k_15_1.py
index 6af6f7520be408461c700e8ff9c70127c3b02d22..c85448023d33d1fd41cb2761102983c4b15667bd 100644
--- a/models/Dolico_bon/Dolico_bon_500k_15_1.py
+++ b/models/Dolico_bon/Dolico_bon_500k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/500k'
-    p['mandible'] = f'{path}/Dolico_bon_500k.stl'
+    p['bone'] = f'{path}/Dolico_bon_500k.stl'
     p['contact_pts'] = [[4.858893, -266.720367, 28.316463],
                 [36.279472, -269.853271, 26.685869]]
     p['axis_pt1'] = [164.323242, 605.580627, -20.302845]
diff --git a/models/Dolico_bon/Dolico_bon_500k_15_2.py b/models/Dolico_bon/Dolico_bon_500k_15_2.py
index 99cf92f6025b410c458aeb18c58858e86f33b432..f10f8eedcb21313a5062b1bdc63cf36014e901a2 100644
--- a/models/Dolico_bon/Dolico_bon_500k_15_2.py
+++ b/models/Dolico_bon/Dolico_bon_500k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/500k'
-    p['mandible'] = f'{path}/Dolico_bon_500k.stl'
+    p['bone'] = f'{path}/Dolico_bon_500k.stl'
     p['contact_pts'] = [[-16.827028, 42.577381, -1.938522],
                 [44.807407, 48.136395, -0.112404]]
     p['axis_pt1'] = [164.323242, 605.580627, -20.302845]
diff --git a/models/Dolico_bon/Dolico_bon_500k_7_1.py b/models/Dolico_bon/Dolico_bon_500k_7_1.py
index 2746e93387f87ef8cc186f5102ae42b102a359ff..cfce80e0d90da38a645f7e1face767651a32f660 100644
--- a/models/Dolico_bon/Dolico_bon_500k_7_1.py
+++ b/models/Dolico_bon/Dolico_bon_500k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/500k'
-    p['mandible'] = f'{path}/Dolico_bon_500k.stl'
+    p['bone'] = f'{path}/Dolico_bon_500k.stl'
     p['contact_pts'] = [[4.858893, -266.720367, 28.316463],
                 [36.279472, -269.853271, 26.685869]]
     p['axis_pt1'] = [164.323242, 605.580627, -20.302845]
diff --git a/models/Dolico_bon/Dolico_bon_500k_7_2.py b/models/Dolico_bon/Dolico_bon_500k_7_2.py
index 3fbe9b08ace2767443d31b96e4d84b1ed42068b6..458cb378240ba61a716c5f3097d3b3a41fbb83c8 100644
--- a/models/Dolico_bon/Dolico_bon_500k_7_2.py
+++ b/models/Dolico_bon/Dolico_bon_500k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/500k'
-    p['mandible'] = f'{path}/Dolico_bon_500k.stl'
+    p['bone'] = f'{path}/Dolico_bon_500k.stl'
     p['contact_pts'] = [[-16.827028, 42.577381, -1.938522],
                 [44.807407, 48.136395, -0.112404]]
     p['axis_pt1'] = [164.323242, 605.580627, -20.302845]
diff --git a/models/Dolico_bon/Dolico_bon_750k_15_1.py b/models/Dolico_bon/Dolico_bon_750k_15_1.py
index 8f8f4be71ce37cc316b05aa2c4f88577e42ccb37..9511f78cf812fa5eed4441cc9bed2ac994e7cae8 100644
--- a/models/Dolico_bon/Dolico_bon_750k_15_1.py
+++ b/models/Dolico_bon/Dolico_bon_750k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/750k'
-    p['mandible'] = f'{path}/Dolico_bon_750k.stl'
+    p['bone'] = f'{path}/Dolico_bon_750k.stl'
     p['contact_pts'] = [[5.189603, -266.977203, 28.278206],
                 [36.279472, -269.853271, 26.685869]]
     p['axis_pt1'] = [163.730804, 595.499817, -20.042530]
diff --git a/models/Dolico_bon/Dolico_bon_750k_15_2.py b/models/Dolico_bon/Dolico_bon_750k_15_2.py
index 5ed97a7bfe1abc814e334e3939bdff9cf48df2fd..04f50c7964e495e3c6a280a1ccff48bca45f5ba9 100644
--- a/models/Dolico_bon/Dolico_bon_750k_15_2.py
+++ b/models/Dolico_bon/Dolico_bon_750k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/750k'
-    p['mandible'] = f'{path}/Dolico_bon_750k.stl'
+    p['bone'] = f'{path}/Dolico_bon_750k.stl'
     p['contact_pts'] = [[-18.989965, 40.841747, -1.935307],
                 [34.655277, 41.826595, -7.642198]]
     p['axis_pt1'] = [163.730804, 595.499817, -20.042530]
diff --git a/models/Dolico_bon/Dolico_bon_750k_7_1.py b/models/Dolico_bon/Dolico_bon_750k_7_1.py
index 575253d410065cb1b12ff07e0294ae68b1f8e27e..708b93115567b667384d9e040c4bb6e3587c4b6a 100644
--- a/models/Dolico_bon/Dolico_bon_750k_7_1.py
+++ b/models/Dolico_bon/Dolico_bon_750k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/750k'
-    p['mandible'] = f'{path}/Dolico_bon_750k.stl'
+    p['bone'] = f'{path}/Dolico_bon_750k.stl'
     p['contact_pts'] = [[5.189603, -266.977203, 28.278206],
                 [36.279472, -269.853271, 26.685869]]
     p['axis_pt1'] = [163.730804, 595.499817, -20.042530]
diff --git a/models/Dolico_bon/Dolico_bon_750k_7_2.py b/models/Dolico_bon/Dolico_bon_750k_7_2.py
index d7cbd8609bb122bfefe582ad8c98abf2c8c66ad2..1fd3809e45d750e8a48b21e3c549a3799a5513bb 100644
--- a/models/Dolico_bon/Dolico_bon_750k_7_2.py
+++ b/models/Dolico_bon/Dolico_bon_750k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_bon/750k'
-    p['mandible'] = f'{path}/Dolico_bon_750k.stl'
+    p['bone'] = f'{path}/Dolico_bon_750k.stl'
     p['contact_pts'] = [[-18.989965, 40.841747, -1.935307],
                 [34.655277, 41.826595, -7.642198]]
     p['axis_pt1'] = [163.730804, 595.499817, -20.042530]
diff --git a/models/Dolico_osb/Dolico_osb_150k_15_1.py b/models/Dolico_osb/Dolico_osb_150k_15_1.py
index 77d8a0db4aa369bc8bf1f5a38c1f1d47a6618aa4..a7e112c9a9c4cc8e6b6545fdad655b0d6cf3f78a 100644
--- a/models/Dolico_osb/Dolico_osb_150k_15_1.py
+++ b/models/Dolico_osb/Dolico_osb_150k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/150k'
-    p['mandible'] = f'{path}/Dolico_osb_150k.stl'
+    p['bone'] = f'{path}/Dolico_osb_150k.stl'
     p['contact_pts'] = [[70.939697, -369.730804, 305.997803],
                 [72.499474, -367.223450, 278.139832]]
     p['axis_pt1'] = [39.231628, 126.756012, 209.439102]
diff --git a/models/Dolico_osb/Dolico_osb_150k_15_2.py b/models/Dolico_osb/Dolico_osb_150k_15_2.py
index 4378765fec4d837b2fbaf80b42a3a4196b674303..77dd876d5c8beebb2d7881f236cc3c10216e1d8f 100644
--- a/models/Dolico_osb/Dolico_osb_150k_15_2.py
+++ b/models/Dolico_osb/Dolico_osb_150k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/150k'
-    p['mandible'] = f'{path}/Dolico_osb_150k.stl'
+    p['bone'] = f'{path}/Dolico_osb_150k.stl'
     p['contact_pts'] = [[52.009193, -201.873154, 307.677521],
                 [64.623055, -199.101791, 282.517853]]
     p['axis_pt1'] = [39.231628, 126.756012, 209.439102]
diff --git a/models/Dolico_osb/Dolico_osb_150k_7_1.py b/models/Dolico_osb/Dolico_osb_150k_7_1.py
index d795576e48ce289e30d1fccd0961f71e14599c32..874bec86f19f8aebc65dee5645632b8f2fe694e5 100644
--- a/models/Dolico_osb/Dolico_osb_150k_7_1.py
+++ b/models/Dolico_osb/Dolico_osb_150k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/150k'
-    p['mandible'] = f'{path}/Dolico_osb_150k.stl'
+    p['bone'] = f'{path}/Dolico_osb_150k.stl'
     p['contact_pts'] = [[70.939697, -369.730804, 305.997803],
                 [72.499474, -367.223450, 278.139832]]
     p['axis_pt1'] = [39.231628, 126.756012, 209.439102]
diff --git a/models/Dolico_osb/Dolico_osb_150k_7_2.py b/models/Dolico_osb/Dolico_osb_150k_7_2.py
index 3d8cf339832a465736f438fb4d801473416f6c7f..ecd205e752928e8d673b754bc2aae576a4073343 100644
--- a/models/Dolico_osb/Dolico_osb_150k_7_2.py
+++ b/models/Dolico_osb/Dolico_osb_150k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/150k'
-    p['mandible'] = f'{path}/Dolico_osb_150k.stl'
+    p['bone'] = f'{path}/Dolico_osb_150k.stl'
     p['contact_pts'] = [[52.009193, -201.873154, 307.677521],
                 [64.623055, -199.101791, 282.517853]]
     p['axis_pt1'] = [39.231628, 126.756012, 209.439102]
diff --git a/models/Dolico_osb/Dolico_osb_250k_15_1.py b/models/Dolico_osb/Dolico_osb_250k_15_1.py
index 15e269f027e73d0d71f2b6a12f6a2a52c85b371c..eff7c305c670dd4831f90b94d8be1fabe98f9222 100644
--- a/models/Dolico_osb/Dolico_osb_250k_15_1.py
+++ b/models/Dolico_osb/Dolico_osb_250k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/250k'
-    p['mandible'] = f'{path}/Dolico_osb_250k.stl'
+    p['bone'] = f'{path}/Dolico_osb_250k.stl'
     p['contact_pts'] = [[71.647064, -367.307281, 303.918365],
                 [71.428085, -366.684906, 278.119751]]
     p['axis_pt1'] = [25.142067, 117.287262, 198.309219]
diff --git a/models/Dolico_osb/Dolico_osb_250k_15_2.py b/models/Dolico_osb/Dolico_osb_250k_15_2.py
index 02fc9715a354b6bc0e4a2c0911c8fb126fe01291..a678573c65eeaf642d7e7e4ec58383a9d4dea717 100644
--- a/models/Dolico_osb/Dolico_osb_250k_15_2.py
+++ b/models/Dolico_osb/Dolico_osb_250k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/250k'
-    p['mandible'] = f'{path}/Dolico_osb_250k.stl'
+    p['bone'] = f'{path}/Dolico_osb_250k.stl'
     p['contact_pts'] = [[58.887959, -197.486603, 305.862701],
                 [64.476349, -199.363113, 282.813385]]
     p['axis_pt1'] = [25.142067, 117.287262, 198.309219]
diff --git a/models/Dolico_osb/Dolico_osb_250k_7_1.py b/models/Dolico_osb/Dolico_osb_250k_7_1.py
index 167ecf5ad2f7b3265ce9ffd079951b998c4e3aca..a0aa2f704bd672c39503e499b86b67130cd23353 100644
--- a/models/Dolico_osb/Dolico_osb_250k_7_1.py
+++ b/models/Dolico_osb/Dolico_osb_250k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/250k'
-    p['mandible'] = f'{path}/Dolico_osb_250k.stl'
+    p['bone'] = f'{path}/Dolico_osb_250k.stl'
     p['contact_pts'] = [[71.647064, -367.307281, 303.918365],
                 [71.428085, -366.684906, 278.119751]]
     p['axis_pt1'] = [25.142067, 117.287262, 198.309219]
diff --git a/models/Dolico_osb/Dolico_osb_250k_7_2.py b/models/Dolico_osb/Dolico_osb_250k_7_2.py
index 914085bcae7612580f43386f71011c5502550b5f..dda249c890b25eed41d7bcdd18b3f89e90ef79f2 100644
--- a/models/Dolico_osb/Dolico_osb_250k_7_2.py
+++ b/models/Dolico_osb/Dolico_osb_250k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/250k'
-    p['mandible'] = f'{path}/Dolico_osb_250k.stl'
+    p['bone'] = f'{path}/Dolico_osb_250k.stl'
     p['contact_pts'] = [[58.887959, -197.486603, 305.862701],
                 [64.476349, -199.363113, 282.813385]]
     p['axis_pt1'] = [25.142067, 117.287262, 198.309219]
diff --git a/models/Dolico_osb/Dolico_osb_500k_15_1.py b/models/Dolico_osb/Dolico_osb_500k_15_1.py
index 64368c79ebd6ed375b2602fea826b7ad1f4dfa31..0cbe7b49c59a7e17412768e240a40f33d0e23b4d 100644
--- a/models/Dolico_osb/Dolico_osb_500k_15_1.py
+++ b/models/Dolico_osb/Dolico_osb_500k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/500k'
-    p['mandible'] = f'{path}/Dolico_osb_500k.stl'
+    p['bone'] = f'{path}/Dolico_osb_500k.stl'
     p['contact_pts'] = [[72.492805, -367.768738, 304.363922],
                 [72.060631, -366.996246, 278.112061]]
     p['axis_pt1'] = [27.536142, 119.829124, 198.032684]
diff --git a/models/Dolico_osb/Dolico_osb_500k_15_2.py b/models/Dolico_osb/Dolico_osb_500k_15_2.py
index a9038e7113b8077a3af547643a6237daf2279bff..d4323937aa595eaca502f74a00ea99cf26474cc5 100644
--- a/models/Dolico_osb/Dolico_osb_500k_15_2.py
+++ b/models/Dolico_osb/Dolico_osb_500k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/500k'
-    p['mandible'] = f'{path}/Dolico_osb_500k.stl'
+    p['bone'] = f'{path}/Dolico_osb_500k.stl'
     p['contact_pts'] = [[58.887959, -197.486603, 305.862701],
                 [58.952919, -197.350830, 280.833832]]
     p['axis_pt1'] = [27.536142, 119.829124, 198.032684]
diff --git a/models/Dolico_osb/Dolico_osb_500k_7_1.py b/models/Dolico_osb/Dolico_osb_500k_7_1.py
index eb6e1e5c9641239be7ce592fd546704c61a9a677..b71a2207dd5740e72493ca75bff1c3d69b1911db 100644
--- a/models/Dolico_osb/Dolico_osb_500k_7_1.py
+++ b/models/Dolico_osb/Dolico_osb_500k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/500k'
-    p['mandible'] = f'{path}/Dolico_osb_500k.stl'
+    p['bone'] = f'{path}/Dolico_osb_500k.stl'
     p['contact_pts'] = [[72.492805, -367.768738, 304.363922],
                 [72.060631, -366.996246, 278.112061]]
     p['axis_pt1'] = [27.536142, 119.829124, 198.032684]
diff --git a/models/Dolico_osb/Dolico_osb_500k_7_2.py b/models/Dolico_osb/Dolico_osb_500k_7_2.py
index b1bcd8542c8575ce58de1e049f67cfe1dfe88425..4176d789e9a277822690c641c20c51721d741dff 100644
--- a/models/Dolico_osb/Dolico_osb_500k_7_2.py
+++ b/models/Dolico_osb/Dolico_osb_500k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/500k'
-    p['mandible'] = f'{path}/Dolico_osb_500k.stl'
+    p['bone'] = f'{path}/Dolico_osb_500k.stl'
     p['contact_pts'] = [[58.887959, -197.486603, 305.862701],
                 [58.952919, -197.350830, 280.833832]]
     p['axis_pt1'] = [27.536142, 119.829124, 198.032684]
diff --git a/models/Dolico_osb/Dolico_osb_750k_15_1.py b/models/Dolico_osb/Dolico_osb_750k_15_1.py
index 8b631a4310549ae40d20dd2c88e262163b34a0e4..3ce53de519edd96cb2cceb4ec70b38739adfa51e 100644
--- a/models/Dolico_osb/Dolico_osb_750k_15_1.py
+++ b/models/Dolico_osb/Dolico_osb_750k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/750k'
-    p['mandible'] = f'{path}/Dolico_osb_750k.stl'
+    p['bone'] = f'{path}/Dolico_osb_750k.stl'
     p['contact_pts'] = [[72.896973, -368.044189, 304.529999],
                 [72.543152, -367.867920, 278.888367]]
     p['axis_pt1'] = [28.769444, 127.225159, 198.837112]
diff --git a/models/Dolico_osb/Dolico_osb_750k_15_2.py b/models/Dolico_osb/Dolico_osb_750k_15_2.py
index 38e46ff67399daa771a6ebbd83c7920e3468ec4e..6c95f1269eb4b7d15dbb735e6f89d56631060d98 100644
--- a/models/Dolico_osb/Dolico_osb_750k_15_2.py
+++ b/models/Dolico_osb/Dolico_osb_750k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/750k'
-    p['mandible'] = f'{path}/Dolico_osb_750k.stl'
+    p['bone'] = f'{path}/Dolico_osb_750k.stl'
     p['contact_pts'] = [[58.887959, -197.486603, 305.862701],
                 [58.952919, -197.350830, 280.833832]]
     p['axis_pt1'] = [28.769444, 127.225159, 198.837112]
diff --git a/models/Dolico_osb/Dolico_osb_750k_7_1.py b/models/Dolico_osb/Dolico_osb_750k_7_1.py
index 0dca00c14f1ea3940566ef984748180c15c59456..ddbcc99735b6b352cd543fcacf22aa79e2054f6f 100644
--- a/models/Dolico_osb/Dolico_osb_750k_7_1.py
+++ b/models/Dolico_osb/Dolico_osb_750k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/750k'
-    p['mandible'] = f'{path}/Dolico_osb_750k.stl'
+    p['bone'] = f'{path}/Dolico_osb_750k.stl'
     p['contact_pts'] = [[72.896973, -368.044189, 304.529999],
                 [72.543152, -367.867920, 278.888367]]
     p['axis_pt1'] = [28.769444, 127.225159, 198.837112]
diff --git a/models/Dolico_osb/Dolico_osb_750k_7_2.py b/models/Dolico_osb/Dolico_osb_750k_7_2.py
index e5e6eb8e4c56aa516ae4c175517c0db339121401..e62167716e081a423d936de6e76ed34f70c536a9 100644
--- a/models/Dolico_osb/Dolico_osb_750k_7_2.py
+++ b/models/Dolico_osb/Dolico_osb_750k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Dolico_osb/750k'
-    p['mandible'] = f'{path}/Dolico_osb_750k.stl'
+    p['bone'] = f'{path}/Dolico_osb_750k.stl'
     p['contact_pts'] = [[58.887959, -197.486603, 305.862701],
                 [58.952919, -197.350830, 280.833832]]
     p['axis_pt1'] = [28.769444, 127.225159, 198.837112]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_150k_15_1.py b/models/PLTtymFrankencarpus/PLTtympFranken_150k_15_1.py
index 5d81f9160b7790e3dfe626058236979426c60bb4..3c8cfdea1f075421dd81b92caab8c07f7f9c43b2 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_150k_15_1.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_150k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/150k'
-    p['mandible'] = f'{path}/PLTtym_150k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_150k_symph.stl'
     p['contact_pts'] = [[-140.484802, 32.752762, 185.435608],
                 [-141.219009, -31.436172, 178.064987]]
     p['axis_pt1'] = [207.172852, -79.520950, 187.358459]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_150k_15_2.py b/models/PLTtymFrankencarpus/PLTtympFranken_150k_15_2.py
index a66b9a9ae9adcf455a46b996f03c92091a0c63e4..c68301c1a4a7feab385f1d90b83a10581bd7f038 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_150k_15_2.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_150k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/150k'
-    p['mandible'] = f'{path}/PLTtym_150k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_150k_symph.stl'
     p['contact_pts'] = [[-47.591175, 59.634678, 182.979019],
                 [-47.909733, -48.301151, 176.236832]]
     p['axis_pt1'] = [207.172852, -79.520950, 187.358459]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_150k_7_1.py b/models/PLTtymFrankencarpus/PLTtympFranken_150k_7_1.py
index afeae9b396e4ad5ad04460aeb0ef62a03ed4afd2..db8496528f84760860fb31f074b5e0115ffbdde5 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_150k_7_1.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_150k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/150k'
-    p['mandible'] = f'{path}/PLTtym_150k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_150k_symph.stl'
     p['contact_pts'] = [[-140.484802, 32.752762, 185.435608],
                 [-141.219009, -31.436172, 178.064987]]
     p['axis_pt1'] = [207.172852, -79.520950, 187.358459]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_150k_7_2.py b/models/PLTtymFrankencarpus/PLTtympFranken_150k_7_2.py
index 4322ed435bb6bebfd7ead83dae5c64b12568fab7..dba9446c53053c76606396f221b80329892eb9ac 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_150k_7_2.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_150k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/150k'
-    p['mandible'] = f'{path}/PLTtym_150k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_150k_symph.stl'
     p['contact_pts'] = [[-47.591175, 59.634678, 182.979019],
                 [-47.909733, -48.301151, 176.236832]]
     p['axis_pt1'] = [207.172852, -79.520950, 187.358459]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_20k.py b/models/PLTtymFrankencarpus/PLTtympFranken_20k.py
index c815c9a845b269a98e5733b71ef3ed2a8eeab0ff..8deed3c9536ba6b684d4e128ce96ea8a31959a1e 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_20k.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_20k.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/20k'
-    p['mandible'] = f'{path}/PLTtym_20k_symph_lfs.stl'
+    p['bone'] = f'{path}/PLTtym_20k_symph_lfs.stl'
     p['contact_pts'] = [[-139.862946, 34.411610, 186.322037],
                   [-140.781250, -29.009081, 185.878799]]
     p['axis_pt1'] = [214.665802, -70.867393, 186.018097]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_250k.py b/models/PLTtymFrankencarpus/PLTtympFranken_250k.py
index 42dd60a6627527f86581c7a50eac8823454a68df..c9e61bf3cccffd7551ffaaee3baeb196bf91a495 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_250k.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_250k.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/250k'
-    p['mandible'] = f'{path}/PLTtym_250k_symph_lfs.stl'
+    p['bone'] = f'{path}/PLTtym_250k_symph_lfs.stl'
     p['contact_pts'] = [[-140.616379, 34.634388, 178.193710],
                 [-141.087341, -31.283962, 177.457901]]
     p['axis_pt1'] = [206.616394, -80.626038, 187.317291]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_250k_15_1.py b/models/PLTtymFrankencarpus/PLTtympFranken_250k_15_1.py
index 225be6af95590f15c9fde6adaf9a13235581dc8c..b878d713a22f65ada4f2ed08442fdad1b2615785 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_250k_15_1.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_250k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/250k'
-    p['mandible'] = f'{path}/PLTtym_250k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_250k_symph.stl'
     p['contact_pts'] = [[-140.616379, 34.634388, 178.193710],
                 [-141.087341, -31.283962, 177.457901]]
     p['axis_pt1'] = [206.616394, -80.626038, 187.317291]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_250k_15_2.py b/models/PLTtymFrankencarpus/PLTtympFranken_250k_15_2.py
index e0ff55c505c30661175110634c2fa50f8a954bc4..6818487a5f1fb538516b0a96ac6267bcd029e0f0 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_250k_15_2.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_250k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/250k'
-    p['mandible'] = f'{path}/PLTtym_250k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_250k_symph.stl'
     p['contact_pts'] = [[-47.591175, 59.634678, 182.979019],
                 [-47.556026, -48.933891, 180.280045]]
     p['axis_pt1'] = [206.616394, -80.626038, 187.317291]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_250k_7_1.py b/models/PLTtymFrankencarpus/PLTtympFranken_250k_7_1.py
index 5984686751335534c48b7e98af6a71c1e4bd94a5..abdf8eafd8630e981073bbe4f86b764e27678d06 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_250k_7_1.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_250k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/250k'
-    p['mandible'] = f'{path}/PLTtym_250k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_250k_symph.stl'
     p['contact_pts'] = [[-140.616379, 34.634388, 178.193710],
                 [-141.087341, -31.283962, 177.457901]]
     p['axis_pt1'] = [206.616394, -80.626038, 187.317291]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_250k_7_2.py b/models/PLTtymFrankencarpus/PLTtympFranken_250k_7_2.py
index 0008e0f8fd8a6d0d1a8fdc4cf5fee07e3b7def3c..100bdba0a053bd585e422b7c5289b6fb90a9320f 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_250k_7_2.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_250k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/250k'
-    p['mandible'] = f'{path}/PLTtym_250k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_250k_symph.stl'
     p['contact_pts'] = [[-47.591175, 59.634678, 182.979019],
                 [-47.556026, -48.933891, 180.280045]]
     p['axis_pt1'] = [206.616394, -80.626038, 187.317291]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_500k_15_1.py b/models/PLTtymFrankencarpus/PLTtympFranken_500k_15_1.py
index 36313dcfe6a466756786c5cfa70d77796b823314..5aa1bdc326296585908bc6169c9cabb7e606ae83 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_500k_15_1.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_500k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/500k'
-    p['mandible'] = f'{path}/PLTtym_500k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_500k_symph.stl'
     p['contact_pts'] = [[-140.616379, 34.634388, 178.193710],
                 [-141.554642, -31.049093, 176.841354]]
     p['axis_pt1'] = [209.550949, -84.342560, 175.265411]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_500k_15_2.py b/models/PLTtymFrankencarpus/PLTtympFranken_500k_15_2.py
index 74fd2bbd15f9bbdf9b862ffb1f1bd590cf7107f6..0005e5079bd2f77a2623440d8372da436a70e498 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_500k_15_2.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_500k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/500k'
-    p['mandible'] = f'{path}/PLTtym_500k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_500k_symph.stl'
     p['contact_pts'] = [[-48.250050, 60.041912, 182.285355],
                 [-47.561031, -48.934483, 180.753647]]
     p['axis_pt1'] = [209.550949, -84.342560, 175.265411]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_500k_7_1.py b/models/PLTtymFrankencarpus/PLTtympFranken_500k_7_1.py
index 38e8bccc9233a57ec26b13144020c3a4f5267ba0..f20663dc4d762ae0664aa3febc64669be3fbcd3d 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_500k_7_1.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_500k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/500k'
-    p['mandible'] = f'{path}/PLTtym_500k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_500k_symph.stl'
     p['contact_pts'] = [[-140.616379, 34.634388, 178.193710],
                 [-141.554642, -31.049093, 176.841354]]
     p['axis_pt1'] = [209.550949, -84.342560, 175.265411]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_500k_7_2.py b/models/PLTtymFrankencarpus/PLTtympFranken_500k_7_2.py
index f0c33bb6e277520c56decc4596c6a5f8e30bf359..d38057a575d41c0851f62b3ecfb61b9e91be93c3 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_500k_7_2.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_500k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/500k'
-    p['mandible'] = f'{path}/PLTtym_500k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_500k_symph.stl'
     p['contact_pts'] = [[-48.250050, 60.041912, 182.285355],
                 [-47.561031, -48.934483, 180.753647]]
     p['axis_pt1'] = [209.550949, -84.342560, 175.265411]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_750k.py b/models/PLTtymFrankencarpus/PLTtympFranken_750k.py
index 8f664dcbbe690940fc432da92026001029ae4d29..c54850201ebbfbb756729021380c09be3ff54d74 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_750k.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_750k.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/750k'
-    p['mandible'] = f'{path}/PLTtym_750k_symph_lfs.stl'
+    p['bone'] = f'{path}/PLTtym_750k_symph_lfs.stl'
     p['contact_pts'] = [[-140.392929, 34.519073, 177.620209],
                 [-141.554642, -31.049093, 176.841354]]
     p['axis_pt1'] = [207.311279, -80.298553, 187.387451]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_750k_15_1.py b/models/PLTtymFrankencarpus/PLTtympFranken_750k_15_1.py
index 3549e64575d954b607cb4c3657e65f7f692cb78f..51f224e0d25f9c7f50737b80c6a936ccbdc9bbd1 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_750k_15_1.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_750k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/750k'
-    p['mandible'] = f'{path}/PLTtym_750k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_750k_symph.stl'
     p['contact_pts'] = [[-140.392929, 34.519073, 177.620209],
                 [-141.554642, -31.049093, 176.841354]]
     p['axis_pt1'] = [207.311279, -80.298553, 187.387451]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_750k_15_2.py b/models/PLTtymFrankencarpus/PLTtympFranken_750k_15_2.py
index 0aa2251ee042dcd04a3fbe9e3190b3333e23efe0..864ce1b6f9977378c487357df483a702bf03365f 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_750k_15_2.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_750k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/750k'
-    p['mandible'] = f'{path}/PLTtym_750k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_750k_symph.stl'
     p['contact_pts'] = [[-48.250050, 60.041912, 182.285355],
                 [-47.400860, -48.840836, 180.895721]]
     p['axis_pt1'] = [207.311279, -80.298553, 187.387451]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_750k_7_1.py b/models/PLTtymFrankencarpus/PLTtympFranken_750k_7_1.py
index 0170107da7063638d42437327dc70e181d3bac82..9d33889af1063ed74d18064a696494d4770d30ce 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_750k_7_1.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_750k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/750k'
-    p['mandible'] = f'{path}/PLTtym_750k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_750k_symph.stl'
     p['contact_pts'] = [[-140.392929, 34.519073, 177.620209],
                 [-141.554642, -31.049093, 176.841354]]
     p['axis_pt1'] = [207.311279, -80.298553, 187.387451]
diff --git a/models/PLTtymFrankencarpus/PLTtympFranken_750k_7_2.py b/models/PLTtymFrankencarpus/PLTtympFranken_750k_7_2.py
index 9f6435fa0356e5deeb8cc393d711b1d240096dc4..1f98bc77e3ae8bc4a44eb20fbcca609bee35e690 100644
--- a/models/PLTtymFrankencarpus/PLTtympFranken_750k_7_2.py
+++ b/models/PLTtymFrankencarpus/PLTtympFranken_750k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'PLTtymFrankencarpus/750k'
-    p['mandible'] = f'{path}/PLTtym_750k_symph.stl'
+    p['bone'] = f'{path}/PLTtym_750k_symph.stl'
     p['contact_pts'] = [[-48.250050, 60.041912, 182.285355],
                 [-47.400860, -48.840836, 180.895721]]
     p['axis_pt1'] = [207.311279, -80.298553, 187.387451]
diff --git a/models/Panthera_pardus/Panthera_pardus_RC_30.py b/models/Panthera_pardus/Panthera_pardus_RC_30.py
index 891dca0cd4c6a7f4fc97b098c9346df97f1d6f89..8c15ebdc80017de8cee6eebd98e5079158a75977 100644
--- a/models/Panthera_pardus/Panthera_pardus_RC_30.py
+++ b/models/Panthera_pardus/Panthera_pardus_RC_30.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Panthera_pardus/150k'
-    p['mandible'] = f'{path}/Panthera_pardus_AMNH_113745_mandible.stl'
+    p['bone'] = f'{path}/Panthera_pardus_AMNH_113745_mandible.stl'
     p['contact_pts'] = [16.4418, 17.3567, -53.2842]
     p['axis_pt1'] = [-37.433147, -7.689619, 48.622185]
     p['axis_pt2'] = [37.795021,	-7.337307, 49.882557]
diff --git a/models/README.md b/models/README.md
index f8c32bda369678a4322d0f1c13a236164771c6b1..ae2c3f5f0bf0a736c21b53171005584659829952 100644
--- a/models/README.md
+++ b/models/README.md
@@ -30,7 +30,7 @@ The mesh files are specified in the input file:
 
 ```
 path = 'dolicorhynchops/10k'
-p['mandible'] = f'{path}/mandible.ply'
+p['bone'] = f'{path}/mandible.ply'
 p['contact_pts'] = f'{path}/teeth.off'
 p['axis_pt1'] = f'{path}/LTMJ.off'
 p['axis_pt2'] = f'{path}/RTMJ.off'
diff --git a/models/TYLber/TYLber_150k_15_1.py b/models/TYLber/TYLber_150k_15_1.py
index 318378fc55e6441c16d292331230eddbb097c18b..a443d46cf8089d624561d92ce3532abc2d10833f 100644
--- a/models/TYLber/TYLber_150k_15_1.py
+++ b/models/TYLber/TYLber_150k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/150k'
-    p['mandible'] = f'{path}/TYLber_jaw_150k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_150k.stl'
     p['contact_pts'] = [[74.470337, 531.928772, 38.113079],
                 [-71.465904, 531.618591, 37.830475]]
     p['axis_pt1'] = [-232.980133, -518.113647, 14.254566]
diff --git a/models/TYLber/TYLber_150k_15_2.py b/models/TYLber/TYLber_150k_15_2.py
index fcb3baab0acccb932cec6d9f0376d710f52a1935..30bcedf32965354bd8437c110dcb175707437a6c 100644
--- a/models/TYLber/TYLber_150k_15_2.py
+++ b/models/TYLber/TYLber_150k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/150k'
-    p['mandible'] = f'{path}/TYLber_jaw_150k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_150k.stl'
     p['contact_pts'] = [[159.952377, 145.793655, 10.386530],
                 [-147.823532, 152.052460, 16.779087]]
     p['axis_pt1'] = [-232.980133, -518.113647, 14.254566]
diff --git a/models/TYLber/TYLber_150k_7_1.py b/models/TYLber/TYLber_150k_7_1.py
index 45ba836d9124e1b805dc5fe23409c1aa118d6860..4fc0e07596e53293b17edbcbf31733ca7f48a2a3 100644
--- a/models/TYLber/TYLber_150k_7_1.py
+++ b/models/TYLber/TYLber_150k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/150k'
-    p['mandible'] = f'{path}/TYLber_jaw_150k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_150k.stl'
     p['contact_pts'] = [[74.470337, 531.928772, 38.113079],
                 [-71.465904, 531.618591, 37.830475]]
     p['axis_pt1'] = [-232.980133, -518.113647, 14.254566]
diff --git a/models/TYLber/TYLber_150k_7_2.py b/models/TYLber/TYLber_150k_7_2.py
index 27d6ec9e4696cf99b837b5b02f377ce196ed7826..8aa55048cb5868b2d1796ff897d97cbfd7a16207 100644
--- a/models/TYLber/TYLber_150k_7_2.py
+++ b/models/TYLber/TYLber_150k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/150k'
-    p['mandible'] = f'{path}/TYLber_jaw_150k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_150k.stl'
     p['contact_pts'] = [[159.952377, 145.793655, 10.386530],
                 [-147.823532, 152.052460, 16.779087]]
     p['axis_pt1'] = [-232.980133, -518.113647, 14.254566]
diff --git a/models/TYLber/TYLber_250k_15_1.py b/models/TYLber/TYLber_250k_15_1.py
index aaac1bf38ff02bf4a529f81396fa528e493b3adb..20fb0b5924069b12389965c7865a6ba17840a13b 100644
--- a/models/TYLber/TYLber_250k_15_1.py
+++ b/models/TYLber/TYLber_250k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/250k'
-    p['mandible'] = f'{path}/TYLber_jaw_250k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_250k.stl'
     p['contact_pts'] = [[76.895927, 530.584045, 32.612164],
                 [-73.001038, 536.202087, 37.195301]]
     p['axis_pt1'] = [-249.095215, -528.940735, -12.830746]
diff --git a/models/TYLber/TYLber_250k_15_2_bad.py b/models/TYLber/TYLber_250k_15_2_bad.py
index fe8961d858f4b60ac7cff450509a7e0dba87ec49..cc337963c98e56c9ae5118ac90ddb234241d7577 100644
--- a/models/TYLber/TYLber_250k_15_2_bad.py
+++ b/models/TYLber/TYLber_250k_15_2_bad.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/250k'
-    p['mandible'] = f'{path}/TYLber_jaw_250k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_250k.stl'
     p['contact_pts'] = [[159.370773, 150.709351, 18.192839],
                 [-152.182220, 150.355057, 22.975760]]
     p['axis_pt1'] = [-249.095215, -528.940735, -12.830746]
diff --git a/models/TYLber/TYLber_250k_7_1.py b/models/TYLber/TYLber_250k_7_1.py
index f175bd116393cb155a0231502eb6140c1ec14b67..91e960daaed7f492f424450a146bf79d2e8f19a0 100644
--- a/models/TYLber/TYLber_250k_7_1.py
+++ b/models/TYLber/TYLber_250k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/250k'
-    p['mandible'] = f'{path}/TYLber_jaw_250k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_250k.stl'
     p['contact_pts'] = [[76.895927, 530.584045, 32.612164],
                 [-73.001038, 536.202087, 37.195301]]
     p['axis_pt1'] = [-249.095215, -528.940735, -12.830746]
diff --git a/models/TYLber/TYLber_250k_7_2_bad.py b/models/TYLber/TYLber_250k_7_2_bad.py
index 5a0257aa87d970c6b9c0f1619b8347b0650f2ab0..fb5cefeed2875238983052e1c5c71ab776446b05 100644
--- a/models/TYLber/TYLber_250k_7_2_bad.py
+++ b/models/TYLber/TYLber_250k_7_2_bad.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/250k'
-    p['mandible'] = f'{path}/TYLber_jaw_250k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_250k.stl'
     p['contact_pts'] = [[159.370773, 150.709351, 18.192839],
                 [-152.182220, 150.355057, 22.975760]]
     p['axis_pt1'] = [-249.095215, -528.940735, -12.830746]
diff --git a/models/TYLber/TYLber_500k_15_1.py b/models/TYLber/TYLber_500k_15_1.py
index 415a5d6bc4c183ed2d21b45ac0d48043b62c5974..57daab1c1a4d9e3b00bc3075570c4c7f12426acb 100644
--- a/models/TYLber/TYLber_500k_15_1.py
+++ b/models/TYLber/TYLber_500k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/500k'
-    p['mandible'] = f'{path}/TYLber_jaw_500k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_500k.stl'
     p['contact_pts'] = [[77.305153, 544.409424, 30.836790],
                 [-73.187263, 536.783264, 36.419006]]
     p['axis_pt1'] = [-228.034103, -519.674866, 16.911722]
diff --git a/models/TYLber/TYLber_500k_15_2.py b/models/TYLber/TYLber_500k_15_2.py
index b69d369598e5d4e745fc5d73207edb5f1f1110be..3f7878a5c4ed5b476bceeb61dc35e48819350e01 100644
--- a/models/TYLber/TYLber_500k_15_2.py
+++ b/models/TYLber/TYLber_500k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/500k'
-    p['mandible'] = f'{path}/TYLber_jaw_500k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_500k.stl'
     p['contact_pts'] = [[153.858017, 155.593140, 17.697323],
                 [-150.512238, 145.369690, 17.743414]]
     p['axis_pt1'] = [-228.034103, -519.674866, 16.911722]
diff --git a/models/TYLber/TYLber_500k_7_1.py b/models/TYLber/TYLber_500k_7_1.py
index 38d7836472c642124ebc3144495d09d65ddfc9a5..9d6b2bcf5c4742c121b7e8ed0a51a67932cd3aa7 100644
--- a/models/TYLber/TYLber_500k_7_1.py
+++ b/models/TYLber/TYLber_500k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/500k'
-    p['mandible'] = f'{path}/TYLber_jaw_500k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_500k.stl'
     p['contact_pts'] = [[77.305153, 544.409424, 30.836790],
                 [-73.187263, 536.783264, 36.419006]]
     p['axis_pt1'] = [-228.034103, -519.674866, 16.911722]
diff --git a/models/TYLber/TYLber_500k_7_2.py b/models/TYLber/TYLber_500k_7_2.py
index 0fdf18193d278df29541beeb3831420beba6e496..a9e54960f46f8a6704da45c544e45e10d9467e5c 100644
--- a/models/TYLber/TYLber_500k_7_2.py
+++ b/models/TYLber/TYLber_500k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/500k'
-    p['mandible'] = f'{path}/TYLber_jaw_500k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_500k.stl'
     p['contact_pts'] = [[153.858017, 155.593140, 17.697323],
                 [-150.512238, 145.369690, 17.743414]]
     p['axis_pt1'] = [-228.034103, -519.674866, 16.911722]
diff --git a/models/TYLber/TYLber_750k_15_1.py b/models/TYLber/TYLber_750k_15_1.py
index 94140911da8578cf838962bb1218015ad7bb47cc..c0b2e16fb390fadf6f7c1f917fa538bcc9b58a95 100644
--- a/models/TYLber/TYLber_750k_15_1.py
+++ b/models/TYLber/TYLber_750k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/750k'
-    p['mandible'] = f'{path}/TYLber_jaw_750k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_750k.stl'
     p['contact_pts'] = [[78.026253, 541.746887, 33.878143],
                 [-73.820549, 540.698853, 33.361473]]
     p['axis_pt1'] = [-212.512512, -567.425232, -59.523800]
diff --git a/models/TYLber/TYLber_750k_15_2.py b/models/TYLber/TYLber_750k_15_2.py
index 6b47b4c2f111a9398a77748cd56453cb6a379c25..5534a4018019a5f1ff107a19f009f3a0d8551cfa 100644
--- a/models/TYLber/TYLber_750k_15_2.py
+++ b/models/TYLber/TYLber_750k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/750k'
-    p['mandible'] = f'{path}/TYLber_jaw_750k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_750k.stl'
     p['contact_pts'] = [[196.424347, 90.526627, -41.083321],
                 [-159.784760, 98.047890, 30.349947]]
     p['axis_pt1'] = [-212.512512, -567.425232, -59.523800]
diff --git a/models/TYLber/TYLber_750k_7_1.py b/models/TYLber/TYLber_750k_7_1.py
index d2f162ef01eab1b8fc50555d1c76e17c81222928..7a802a52b0cdb403072ba5ff41619e316d30a29f 100644
--- a/models/TYLber/TYLber_750k_7_1.py
+++ b/models/TYLber/TYLber_750k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/750k'
-    p['mandible'] = f'{path}/TYLber_jaw_750k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_750k.stl'
     p['contact_pts'] = [[78.026253, 541.746887, 33.878143],
                 [-73.820549, 540.698853, 33.361473]]
     p['axis_pt1'] = [-212.512512, -567.425232, -59.523800]
diff --git a/models/TYLber/TYLber_750k_7_2.py b/models/TYLber/TYLber_750k_7_2.py
index bcee311f23a5d161523d49a38a4eb9c85b37ee93..a3fec53bfbe8ce301002df40091987b32526a283 100644
--- a/models/TYLber/TYLber_750k_7_2.py
+++ b/models/TYLber/TYLber_750k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLber/750k'
-    p['mandible'] = f'{path}/TYLber_jaw_750k.stl'
+    p['bone'] = f'{path}/TYLber_jaw_750k.stl'
     p['contact_pts'] = [[196.424347, 90.526627, -41.083321],
                 [-159.784760, 98.047890, 30.349947]]
     p['axis_pt1'] = [-212.512512, -567.425232, -59.523800]
diff --git a/models/TYLnep/TYLnep_150k_15_1.py b/models/TYLnep/TYLnep_150k_15_1.py
index 5441cfc0eff85e58951992f3ec6fad9f0ec89b04..c041f562bf312644a4e93cf48eac2a2fdce985a7 100644
--- a/models/TYLnep/TYLnep_150k_15_1.py
+++ b/models/TYLnep/TYLnep_150k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/150k'
-    p['mandible'] = f'{path}/TYLnep_jaw_150k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_150k.stl'
     p['contact_pts'] = [[-411.039917, 227.579102, 43.466621],
                 [-426.836792, 165.685944, 36.533264]]
     p['axis_pt1'] = [380.324341, 44.402424, 89.259888]
diff --git a/models/TYLnep/TYLnep_150k_15_2.py b/models/TYLnep/TYLnep_150k_15_2.py
index 8caf1cdaf4eb413b5ae97cc80060e71e95cc2c7f..95288ff2c8766adc64d1c2243449dedee7e313ea 100644
--- a/models/TYLnep/TYLnep_150k_15_2.py
+++ b/models/TYLnep/TYLnep_150k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/150k'
-    p['mandible'] = f'{path}/TYLnep_jaw_150k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_150k.stl'
     p['contact_pts'] = [[-87.998215, 318.365540, 63.300953],
                 [-99.972519, 80.537491, 43.511547]]
     p['axis_pt1'] = [380.324341, 44.402424, 89.259888]
diff --git a/models/TYLnep/TYLnep_150k_7_1.py b/models/TYLnep/TYLnep_150k_7_1.py
index 984bf270820a9a6f82e1d3283e242661e52dc526..48fa28713315736235e594ca7ae619479dcc286b 100644
--- a/models/TYLnep/TYLnep_150k_7_1.py
+++ b/models/TYLnep/TYLnep_150k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/150k'
-    p['mandible'] = f'{path}/TYLnep_jaw_150k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_150k.stl'
     p['contact_pts'] = [[-411.039917, 227.579102, 43.466621],
                 [-426.836792, 165.685944, 36.533264]]
     p['axis_pt1'] = [380.324341, 44.402424, 89.259888]
diff --git a/models/TYLnep/TYLnep_150k_7_2.py b/models/TYLnep/TYLnep_150k_7_2.py
index 41d36296d0bce4aec4208541e08b497935d02008..a4b5b33fb1e903045d7a9d773d34052d8814cd91 100644
--- a/models/TYLnep/TYLnep_150k_7_2.py
+++ b/models/TYLnep/TYLnep_150k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/150k'
-    p['mandible'] = f'{path}/TYLnep_jaw_150k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_150k.stl'
     p['contact_pts'] = [[-87.998215, 318.365540, 63.300953],
                 [-99.972519, 80.537491, 43.511547]]
     p['axis_pt1'] = [380.324341, 44.402424, 89.259888]
diff --git a/models/TYLnep/TYLnep_250k_15_1.py b/models/TYLnep/TYLnep_250k_15_1.py
index 9bd78aecec8ebefccdb1364b28b5fe3cf8467b64..d32983436e0137f238f1ce164aaad4d39ee24a54 100644
--- a/models/TYLnep/TYLnep_250k_15_1.py
+++ b/models/TYLnep/TYLnep_250k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/250k'
-    p['mandible'] = f'{path}/TYLnep_jaw_250k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_250k.stl'
     p['contact_pts'] = [[-420.489166, 226.380966, 16.269716],
                 [-426.962891, 159.505264, 35.790138]]
     p['axis_pt1'] = [384.577301, 39.049088, 84.218117]
diff --git a/models/TYLnep/TYLnep_250k_15_2.py b/models/TYLnep/TYLnep_250k_15_2.py
index bdea96351629e916156e6a80bd52e1936536859e..7d8837f38a828da32c37ce9cd286f11be6a7e08f 100644
--- a/models/TYLnep/TYLnep_250k_15_2.py
+++ b/models/TYLnep/TYLnep_250k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/250k'
-    p['mandible'] = f'{path}/TYLnep_jaw_250k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_250k.stl'
     p['contact_pts'] = [[-86.279778, 320.183533, 52.348541],
                 [-96.408920, 78.982399, 52.918877]]
     p['axis_pt1'] = [384.577301, 39.049088, 84.218117]
diff --git a/models/TYLnep/TYLnep_250k_7_1.py b/models/TYLnep/TYLnep_250k_7_1.py
index d11f1bb2f7bab851d69c71a11f0f1eeb859fbb0d..42a97dd8858d08dc34e3433611f90f564f1f8ecc 100644
--- a/models/TYLnep/TYLnep_250k_7_1.py
+++ b/models/TYLnep/TYLnep_250k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/250k'
-    p['mandible'] = f'{path}/TYLnep_jaw_250k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_250k.stl'
     p['contact_pts'] = [[-420.489166, 226.380966, 16.269716],
                 [-426.962891, 159.505264, 35.790138]]
     p['axis_pt1'] = [384.577301, 39.049088, 84.218117]
diff --git a/models/TYLnep/TYLnep_250k_7_2.py b/models/TYLnep/TYLnep_250k_7_2.py
index 24d02f16a7f3320ebf07529988ad94cd6198ecee..ccd415c3aef2a56e72b3e94c4cf04b52b20c2d97 100644
--- a/models/TYLnep/TYLnep_250k_7_2.py
+++ b/models/TYLnep/TYLnep_250k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/250k'
-    p['mandible'] = f'{path}/TYLnep_jaw_250k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_250k.stl'
     p['contact_pts'] = [[-86.279778, 320.183533, 52.348541],
                 [-96.408920, 78.982399, 52.918877]]
     p['axis_pt1'] = [384.577301, 39.049088, 84.218117]
diff --git a/models/TYLnep/TYLnep_500k_15_1.py b/models/TYLnep/TYLnep_500k_15_1.py
index 13380d92893199c092a514bbf89dc0ff68318045..33ab448940b64856769f0dae06fbae32fd314887 100644
--- a/models/TYLnep/TYLnep_500k_15_1.py
+++ b/models/TYLnep/TYLnep_500k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/500k'
-    p['mandible'] = f'{path}/TYLnep_jaw_500k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_500k.stl'
     p['contact_pts'] = [[-417.999359, 222.242569, 31.007950],
                 [-428.224182, 162.031265, 36.896519]]
     p['axis_pt1'] = [381.261627, 39.182289, 82.731766]
diff --git a/models/TYLnep/TYLnep_500k_15_2.py b/models/TYLnep/TYLnep_500k_15_2.py
index 6f4e0d763077400c0f4aefb53f2714e960a49a2a..d383b91b2b8c0034361d31a3891f29406b44de73 100644
--- a/models/TYLnep/TYLnep_500k_15_2.py
+++ b/models/TYLnep/TYLnep_500k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/500k'
-    p['mandible'] = f'{path}/TYLnep_jaw_500k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_500k.stl'
     p['contact_pts'] = [[-88.511429, 320.280304, 40.494282],
                 [-95.066902, 77.562195, 41.572414]]
     p['axis_pt1'] = [381.261627, 39.182289, 82.731766]
diff --git a/models/TYLnep/TYLnep_500k_7_1.py b/models/TYLnep/TYLnep_500k_7_1.py
index cb182b505113a8e45921bee8faaaedaf13978e3d..fdd255eea85899bbe751ce4f198d9a59bd0ed23b 100644
--- a/models/TYLnep/TYLnep_500k_7_1.py
+++ b/models/TYLnep/TYLnep_500k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/500k'
-    p['mandible'] = f'{path}/TYLnep_jaw_500k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_500k.stl'
     p['contact_pts'] = [[-417.999359, 222.242569, 31.007950],
                 [-428.224182, 162.031265, 36.896519]]
     p['axis_pt1'] = [381.261627, 39.182289, 82.731766]
diff --git a/models/TYLnep/TYLnep_500k_7_2.py b/models/TYLnep/TYLnep_500k_7_2.py
index 7d2db4cdf4b99646011419f0ec98af22bcc23587..dbb35614cd75d00992a6eb7cafe02183efb495ff 100644
--- a/models/TYLnep/TYLnep_500k_7_2.py
+++ b/models/TYLnep/TYLnep_500k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/500k'
-    p['mandible'] = f'{path}/TYLnep_jaw_500k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_500k.stl'
     p['contact_pts'] = [[-88.511429, 320.280304, 40.494282],
                 [-95.066902, 77.562195, 41.572414]]
     p['axis_pt1'] = [381.261627, 39.182289, 82.731766]
diff --git a/models/TYLnep/TYLnep_750k_15_1.py b/models/TYLnep/TYLnep_750k_15_1.py
index 87e27f513b085720337b6f8c4d95787ee5da6c73..b319cd9439ea9d566605cddeca94de8b47d70e48 100644
--- a/models/TYLnep/TYLnep_750k_15_1.py
+++ b/models/TYLnep/TYLnep_750k_15_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/750k'
-    p['mandible'] = f'{path}/TYLnep_jaw_750k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_750k.stl'
     p['contact_pts'] = [[-440.282074, 226.035706, 8.676176],
                 [-453.118317, 156.467865, 11.607960]]
     p['axis_pt1'] = [392.782898, 35.557762, 84.161575]
diff --git a/models/TYLnep/TYLnep_750k_15_2.py b/models/TYLnep/TYLnep_750k_15_2.py
index d0c0aaca68f51284964aa2c74e3a1e66980ad099..8b9e32db8f32491975e08f8fa491c03ad67fed48 100644
--- a/models/TYLnep/TYLnep_750k_15_2.py
+++ b/models/TYLnep/TYLnep_750k_15_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/750k'
-    p['mandible'] = f'{path}/TYLnep_jaw_750k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_750k.stl'
     p['contact_pts'] = [[-91.923141, 317.501709, 54.304459],
                 [-96.786438, 79.207016, 55.289452]]
     p['axis_pt1'] = [392.782898, 35.557762, 84.161575]
diff --git a/models/TYLnep/TYLnep_750k_7_1.py b/models/TYLnep/TYLnep_750k_7_1.py
index f6f88fc4203634cf38b6b6080fdae0d04a1a7a21..8e618b310f325ef4987ee22f62113a00e28400d2 100644
--- a/models/TYLnep/TYLnep_750k_7_1.py
+++ b/models/TYLnep/TYLnep_750k_7_1.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/750k'
-    p['mandible'] = f'{path}/TYLnep_jaw_750k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_750k.stl'
     p['contact_pts'] = [[-440.282074, 226.035706, 8.676176],
                 [-453.118317, 156.467865, 11.607960]]
     p['axis_pt1'] = [392.782898, 35.557762, 84.161575]
diff --git a/models/TYLnep/TYLnep_750k_7_2.py b/models/TYLnep/TYLnep_750k_7_2.py
index 3614d05c874d601909872333bee669a1656f6669..a9458c5a56f08f4cb5f36eaaa38b32547a8caf37 100644
--- a/models/TYLnep/TYLnep_750k_7_2.py
+++ b/models/TYLnep/TYLnep_750k_7_2.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'TYLnep/750k'
-    p['mandible'] = f'{path}/TYLnep_jaw_750k.stl'
+    p['bone'] = f'{path}/TYLnep_jaw_750k.stl'
     p['contact_pts'] = [[-91.923141, 317.501709, 54.304459],
                 [-96.786438, 79.207016, 55.289452]]
     p['axis_pt1'] = [392.782898, 35.557762, 84.161575]
diff --git a/models/Tapirus_indicus/T_indicus_shoulder50k.py b/models/Tapirus_indicus/T_indicus_shoulder50k.py
index 9f3313916dc325fc0ef08d5134121dd8125d2b50..c52e31c90c8362d8d8150f6dddb6e5d95ce9a8ed 100644
--- a/models/Tapirus_indicus/T_indicus_shoulder50k.py
+++ b/models/Tapirus_indicus/T_indicus_shoulder50k.py
@@ -6,7 +6,7 @@
 def parms(d={}):
     p = {}
     path = 'Tapirus_indicus/50k'
-    p['mandible'] = f'{path}/TAPind901 humerus.ply'
+    p['bone'] = f'{path}/TAPind901 humerus.ply'
     p['contact_pts'] = [43.003590, 88.870224, -205.887726] # Distal humerus
     p['axis_pt1'] = [18.711651, -26.742256, -93.027802] #Humeral head
     p['axis_pt2'] = [11.976404, -9.600277, -102.244888] #Humeral head
diff --git a/models/bonemodel.py b/models/bonemodel.py
index 1652645bc263f4c967b6a938c3d19e2d137f774e..b266fd9bc8855761179e8f844f6b0cfee573cb74 100644
--- a/models/bonemodel.py
+++ b/models/bonemodel.py
@@ -21,7 +21,7 @@ def parms(d={}):
     # Geomagic/MeshLab input files
     path = 'dolicorhynchops/10k'
     # mandible mesh (.ply/.geo/.msh)
-    p['mandible'] = f'{path}/mandible.ply' 
+    p['bone'] = f'{path}/mandible.ply' 
     # one or several vertices (.ply/.off) (used if p['food']=='fixteeth')
     p['contact_pts'] = f'{path}/teeth.off'
     # left temporomandibular joint - 1 vertex (.ply/.off/coordinates)
@@ -88,7 +88,7 @@ def getMetafor(p={}):
     geometry.setDim3D()
 
     # load main mandible mesh (volume or surface - it will be meshed in 3D)
-    import_mesh(domain, p['mandible'], p['use_gmshOld'])
+    import_mesh(domain, p['bone'], p['use_gmshOld'])
 
     # define groups
 
diff --git a/models/dolicorhynchops/dolicorhynchops_10k.py b/models/dolicorhynchops/dolicorhynchops_10k.py
index 496f75df9360aed996340d0e4b140b05ae814214..6087c273a6ee9a90004106a5cb195006791338af 100644
--- a/models/dolicorhynchops/dolicorhynchops_10k.py
+++ b/models/dolicorhynchops/dolicorhynchops_10k.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'dolicorhynchops/10k'
-    p['mandible'] = f'{path}/mandible.stl'
+    p['bone'] = f'{path}/mandible.stl'
     p['contact_pts'] = [[-10.20362, -17.46838, -229.9061],
                   [-11.92466, 26.3042, -229.5354]]
     p['axis_pt1'] = f'{path}/LTMJ.off'
diff --git a/models/dolicorhynchops/dolicorhynchops_10k_contact.py b/models/dolicorhynchops/dolicorhynchops_10k_contact.py
index 92354193ed28d9b77f4aaedd088c781e2b677e33..f8b7c6c97a31805a43bbb724d22afaf9cf6d9406 100644
--- a/models/dolicorhynchops/dolicorhynchops_10k_contact.py
+++ b/models/dolicorhynchops/dolicorhynchops_10k_contact.py
@@ -8,7 +8,7 @@
 def parms(d={}):
     p = {}
     path = 'dolicorhynchops/10k'
-    p['mandible'] = f'{path}/mandible.ply'
+    p['bone'] = f'{path}/mandible.ply'
     # p['contact_pts'] = f'{path}/teeth.off'
     p['axis_pt1'] = [0.1458893, -73.13895, 227.3909]
     p['axis_pt2'] = f'{path}/RTMJ.off'
diff --git a/models/dolicorhynchops/dolicorhynchops_149k.py b/models/dolicorhynchops/dolicorhynchops_149k.py
index ee24fa77b8e00f966c3cbe134564219839e48586..cb399138d9673494b28606187b130abbe645800b 100644
--- a/models/dolicorhynchops/dolicorhynchops_149k.py
+++ b/models/dolicorhynchops/dolicorhynchops_149k.py
@@ -7,7 +7,7 @@
 def parms(d={}):
     p = {}
     path = 'dolicorhynchops/149k'
-    p['mandible'] = f'{path}/mandible_lfs.ply'
+    p['bone'] = f'{path}/mandible_lfs.ply'
     p['contact_pts'] = f'{path}/teeth.ply'
     p['axis_pt1'] = f'{path}/LTMJ.ply'
     p['axis_pt2'] = f'{path}/RTMJ.ply'
diff --git a/models/view_model.py b/models/view_model.py
index ad093cf8778f911110e8a65f4446b9132dac0555..ffeb8a0a5849e512c78da399bb24ef54db22f111 100644
--- a/models/view_model.py
+++ b/models/view_model.py
@@ -12,7 +12,7 @@ if __name__=="__main__":
     p = model.parms()
 
     # load meshes
-    surface, _, _ = SurfMesh.load(p['mandible'], vertices=False)
+    surface, _, _ = SurfMesh.load(p['bone'], vertices=False)
     ltmj, _, _ = SurfMesh.load(p['axis_pt1'])
     rtmj, _, _ = SurfMesh.load(p['axis_pt2'])
     teeth, _, _ = SurfMesh.load(p['contact_pts'])