diff --git a/models/Smilodon_fatalis_juvenile/500k/Smilodon_fatalis_juvenile_LACMHC2002LR3_mandible_500k_lfs.ply b/models/Smilodon_fatalis_juvenile/500k/Smilodon_fatalis_juvenile_LACMHC2002LR3_mandible_500k_lfs.ply
new file mode 100644
index 0000000000000000000000000000000000000000..704b629538cd2877895724727494eb74c8bd7ce8
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/Smilodon_fatalis_juvenile_LACMHC2002LR3_mandible_500k_lfs.ply
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:81f9e34fce122fb420ccc282b69be73238a179eb52706bc5e1aa5ea5d849b1f9
+size 9499786
diff --git a/models/Smilodon_fatalis_juvenile/500k/m. masseter left.stl b/models/Smilodon_fatalis_juvenile/500k/m. masseter left.stl
new file mode 100644
index 0000000000000000000000000000000000000000..3d1b61b4d1be3cae40a8bc85660a0fa21b9d2c35
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/m. masseter left.stl	
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:766f2501ed240819add0958920c104c4f7627140348aeed3a0ace4b7bccad8ae
+size 2099797
diff --git a/models/Smilodon_fatalis_juvenile/500k/m. masseter right.stl b/models/Smilodon_fatalis_juvenile/500k/m. masseter right.stl
new file mode 100644
index 0000000000000000000000000000000000000000..e030bd4c42de64c992f6e0fece2be9e199561621
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/m. masseter right.stl	
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:96adcd27e29d766622f54f4a9cda724be8db454120b99461f9c73b4c50d4fe03
+size 2442162
diff --git a/models/Smilodon_fatalis_juvenile/500k/m. pterygoideus left.stl b/models/Smilodon_fatalis_juvenile/500k/m. pterygoideus left.stl
new file mode 100644
index 0000000000000000000000000000000000000000..96341ff1c39bc9540a89f6fddfbfd113a00a42af
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/m. pterygoideus left.stl	
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:0fd6e556989938bc7fb061246119894b5f09451444135b8ebdf472f17618e906
+size 425077
diff --git a/models/Smilodon_fatalis_juvenile/500k/m. pterygoideus right.stl b/models/Smilodon_fatalis_juvenile/500k/m. pterygoideus right.stl
new file mode 100644
index 0000000000000000000000000000000000000000..652b4f6937120921398e60f9f4d0605dd5ce4ccb
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/m. pterygoideus right.stl	
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:4dab0b4312d219619e69abd75085e7f55ce73d907397fbfb0e0c9760a1ff343d
+size 351754
diff --git a/models/Smilodon_fatalis_juvenile/500k/m. temporalis left.stl b/models/Smilodon_fatalis_juvenile/500k/m. temporalis left.stl
new file mode 100644
index 0000000000000000000000000000000000000000..617c9e30b4498d9fa3d840926df175943189a7f2
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/m. temporalis left.stl	
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:199b75986ecae6a42fe432fd213323ea8b7956913fb20a678ddad1037776d915
+size 2761401
diff --git a/models/Smilodon_fatalis_juvenile/500k/m. temporalis right.stl b/models/Smilodon_fatalis_juvenile/500k/m. temporalis right.stl
new file mode 100644
index 0000000000000000000000000000000000000000..7e1e68d96a177e53e7849a7cec9bbb5afa5b8d7e
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/500k/m. temporalis right.stl	
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:5c189fc79791e9bd2cd7aac6e8580e948e587be313c882ddd31a72b27e33f20b
+size 3375261
diff --git a/models/Smilodon_fatalis_juvenile/LACMHC_2002_L_R_3_30.py b/models/Smilodon_fatalis_juvenile/LACMHC_2002_L_R_3_30.py
new file mode 100644
index 0000000000000000000000000000000000000000..2105643f1a9300c871259b956f0ca8e832dac662
--- /dev/null
+++ b/models/Smilodon_fatalis_juvenile/LACMHC_2002_L_R_3_30.py
@@ -0,0 +1,80 @@
+#! \usr\bin\env python3
+# -*- coding: utf-8 -*-
+# Smilodon fatalis juvenile L-R-3
+
+
+def parms(d={}):
+    p = {}
+    import os
+    path = os.path.join(os.path.dirname(__file__),'500k')
+    p['bone'] = f'{path}\Smilodon_fatalis_juvenile_LACMHC2002LR3_mandible_500k_lfs.ply'
+    p['contact_pts'] = [ 
+                   [-3.128700, 44.665501, -9.424970] ,
+                   [3.667650, 45.900398, -8.284160],
+                   [9.765930, 44.752102, -6.164700]]
+    p['axis_pt1'] = [-27.858801, -35.784000, -1.227750]
+    p['axis_pt2'] = [22.619600, -43.196201, 4.224330] 
+    p['muscles'] = [
+        {
+            'file': f'{path}\m. masseter left.stl', 
+            'force': 60.7224,
+            'focalpt': [-41.7800, -33.8660, 14.6611],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}\m. masseter right.stl',
+            'force': 75.40299,
+            'focalpt': [41.2037, -39.5222, 14.6911],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}\m. pterygoideus left.stl',
+            'force': 12.51531,
+            'focalpt': [-12.7555, -31.5549, 12.3858],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}\m. pterygoideus right.stl',
+            'force': 9.85383,
+            'focalpt': [10.5836, -32.2269, 12.1060],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}\m. temporalis left.stl',
+            'force': 86.78298,
+            'focalpt': [-22.6034, -79.5285, 22.4788],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        },
+        {
+            'file': f'{path}\m. temporalis right.stl',
+            'force': 104.6819,
+            'focalpt': [15.2406, -81.9674, 20.9924],
+            'method': 'T+N'                              # 'U', 'T', 'T+N'
+        }
+    ]
+    p['fixations'] = {
+        'contact_pts': ['z'],
+        'axis_pt1': ['z','y'],
+        'axis_pt2': ['z','y']
+    }
+
+    # material properties
+    p['density'] = 1.662e-9  # [T\mm³]
+    p['Young'] = 18000.      # [MPa]
+    p['Poisson'] = 0.3       # [-]
+
+    p['dt0'] = 0.3       # [s]
+    p['tolNR'] = 1e-3    # [-]
+
+    p.update(d)
+    return p
+
+
+def getMetafor(p={}):
+    import bonemodel as model
+    return model.getMetafor(parms(p))
+
+
+if __name__ == "__main__":
+    import models.bonemodel2 as model
+    model.solve(parms())