diff --git a/models/Panthera_pardus/Panthera_pardus_RC_30.py b/models/Panthera_pardus/Panthera_pardus_RC_30.py
index 0726ed9729c2a19642609800304d50de956795c8..ede6e27972eb78429be8271d5e363a2857b2a583 100644
--- a/models/Panthera_pardus/Panthera_pardus_RC_30.py
+++ b/models/Panthera_pardus/Panthera_pardus_RC_30.py
@@ -70,6 +70,7 @@ def parms(d={}):
     p['Poisson'] = 0.3       # [-]
 
     # p['use_gmshOld'] = True
+    p['narch'] = 2          # number of archives (Metafor)
 
     p.update(d)
     return p
diff --git a/models/bonemodel.py b/models/bonemodel.py
index 6193e0bfdd2e1c53e3fafbff30f5b8817cbc6174..b95a1bba57386717f76a696be5871fc450c8b2e9 100644
--- a/models/bonemodel.py
+++ b/models/bonemodel.py
@@ -81,7 +81,8 @@ def parms(d={}):
     # numerical parameters
     p['tolNR'] = 1e-6        # [-] equilibrium tolerance
     p['dt0'] = 1.0           # [s] time step size
-
+    p['narch'] = 1           # number of archives
+    
     # gmsh toolbox
     p['use_gmshOld'] = False  # use old gmsh interface
 
@@ -232,7 +233,7 @@ def getMetafor(p={}):
 
     tsm = metafor.getTimeStepManager()
     tsm.setInitialTime(0.0, p['dt0'])
-    tsm.setNextTime(1.0, 1, 1.0)
+    tsm.setNextTime(1.0, p['narch'], 1.0)
 
     mim = metafor.getMechanicalIterationManager()
     mim.setResidualTolerance(p['tolNR'])