************************************************************
Model run description: testrun

#########################################################
###### cdem2D output parameters text file ############
#########################################################
 
This model run was started, and the parameters file was created, on Mon Mar  9 11:14:01 2026

 
###### cdem2D  -  Code Origin, Version and Configuration ######
Stuart Hardy's 2D Discrete Element Code
Copyright: Stuart Hardy 2005-2025.
email: stuart.hardy@icrea.cat
 
 
2D Code Version =   2.0.7
 
OpenMP parallelized 2D Code using 32 threads in addition to master thread
 
2D Code compiled using Microsoft Visual Studio 2022
 
 
###### Basic Geometric/Boundary Information ######
Number of (initial) assembly elements read in from input file <med.txt> = 9253
Number of (initial) assembly elements including base, walls, faults etc = 10891
Walls start at element 9254
Right wall starts at element 9404
Walls finish at element 9553
Base starts at element 9554
Base finishes at element 10253
Faults start at element 10254
Second fault starts at element 10414
Faults end at element 10891
 
Number of elements in each wall 150
Number of base elements 700
Number of fault elements 160
 
Initial model width  (m) = 5000.000000
Initial model height (m) = 2769.739250
 
Minimum element radius (m) = 12.500000
Maximum element radius (m) = 31.250000
Average element radius (m) = 19.463282
Boundary/wall element radius (m) = 19.463282
 
Minimum element mass (kg) = 1227184.630309
Maximum element mass (kg) = 7669903.939428
Average element mass (kg) = 2975240.242389
 
 
###### Physical-Material Constants, Calculation/Assembly Parameters and Constraints ######
Spring constant normal interaction (N/m) = 1.544870E+09
Spring constant frictional shear interaction (N/m) = 1.544870E+08
Gravitational constant (m/s2) = 9.810000
Viscosity (damping) constant (N.s.m-1) = 1.500000E+07
Element Density (kg/m3) = 2500.000000
Element Density - Salt (kg/m3) = 2500.000000
Unit scaling length (m) = 250.000000
Max number of neighbours considered when detecting contacts = 10
Max distance (m) searched for neighbours in x and y directions = 93.750000
Max number of possible assembly/array elements = 500000
Sorting element index = 500000
Neighbour recalculation increment (no. of timesteps) = 200
Ratio of neighbour search distance/max element radius  = 3
Number of (pre-growth) marker layers = 24
 
 
 
###### Time-step, Runtime and any Pre-Deformation Stages etc ######
Time-step (s) = 0.010000
Using a timestep that is close to stable limit, for speedier runs - be careful with sedimentation steps
Critical Stable Time-step (s) = 0.028184
Total number of time steps = 100
Total Run Time (s) = 1.000000
Total Run Time (hours) = 0.000278
Total Run Time (days) = 0.000012
 
Equilibrating the assembly for 20000 timesteps before deformation starts
Equilibrating without friction & cohesion
Equilibrating without any bonding
Equilibrating without any cohesion
 
 
 
###### Boundary Displacements, Deformation Modes and File Saving ######
###### Boundary Displacement is -ve for down and towards right (including horizontal)  ######
###### Boundary Displacement is +ve for up and towards left (including horizontal)  ######
Boundary displacement (m/s) = -0.250000
Boundary displacement (m/timestep) = -0.002500
Ratio of boundary displacement/min element radius = 0.000200
Extension mode
Display/Recording increment (m) of boundary displacement = 2.000000
Total (nominal) displacement (m) = 49.750000
 
Display/Recording increment in model seconds = 8.000000
Total number of geometry output files (including 0 time) = 1
First 26 output files (including 0) are solely equilibration
 
###### Fault Configuration ######
Basement fault dip (degrees) = 60
Fault is simulated using a pulled mylar sheet
Fault location (m) = 2235.991030
2nd Fault offset (m) = 1000.000000
 
###### Other Simulation Flags ######
Initial packed elements allowed to overlap
but equilibrium separation to calculate forces still equal to sum of both radii
 
 
Final model geometry saved (newpack.txt) as a possible starting geometry
 
 
###### Material Properties - Frictional: Element Interaction Flags ######
Material is frictional (touching elements with no elastic unbroken bond)
Material is frictional after equilibration (touching elements with no elastic unbroken bond)
    Material is frictional (touching elements) even when initial bond not yet broken
Internal (element-element) Coeff of friction = 0.250000
C0 - Inter-element shear-cohesive force (N) = 5.000000E+06
Cohesion lost after elements separate (i.e. lose initial contact)
Base is frictional
Frictional fault plane
Element-Base Coeff of friction = 0.250000
Walls are frictional
Element-Wall Coeff of friction = 0.250000
Element-Fault Coeff of friction = 0.000000

Testing weak faults option, using strain softening fraction 1.000000



#################################

###### Material Properties - Elastic: Element Interaction Flags ######

Material is initially elastically bonded = tensile strength
Homogeneous internal breaking strain = 0
Average force to break a symmetric bond (N) = 0
Minimum force to break a symmetric bond (N) = 0
Average symmetric bond strength (MPa) = 0
Walls are initially bonded
Wall breaking strain = 0
Base is initially bonded
Basal breaking strain = 0
Faults are initially bonded
Fault breaking strain = 0
 
#################################
Growth sedimentation parameters and flags
No growth sedimentation
 
 
###### cdem2D  -  end of model parameters file ######

###### cdem2D  -  model run statistics ######
Model run started: Mon Mar  9 11:14:01 2026
Model run finished: Mon Mar  9 11:14:01 2026
Model run took 0 seconds