a narch (nb of archives of metafor) as a parameter

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

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'])