granular-basin

commit bfb1778678ea1d379dbd47473c69de546c8d5a3f
parent 9dc3378f0fe774ce724fa45814a1f3db335c6aff
Author: esbenpalmstrom <esbenpalmstroem@gmail.com>
Date:   Thu, 25 Nov 2021 12:08:26 +0100

minor changes

Diffstat:
M
compact_basin.jl
 | 
2
+-
M
init_basin.jl
 | 
2
+-
M
layer_basin.jl
 | 
4
++--

3 files changed, 4 insertions(+), 4 deletions(-)
diff --git a/compact_basin.jl b/compact_basin.jl
@@ -8,7 +8,7 @@ t_start = Dates.now() # Save the start time, print the end time later.
 # lav en lille test? se om dit appendede carpet stadig er forbundet til hoved-
 # simulationsobjektet
 
-id = "simulation1000"    # id of simulation to load
+id = "simulation500"    # id of simulation to load
 N = 20e3                # amount of stress to be applied
 t_comp = 3.0            # compaction max duration [s]
 
diff --git a/init_basin.jl b/init_basin.jl
@@ -12,7 +12,7 @@ t_stack = 0.8                   # duration for each stack to settle [s]
 
 g = [0.,-9.8]                   # vector for direction and magnitude of gravitational acceleration of grains
 
-ngrains = 1000                 # total number of grains
+ngrains = 500                  # total number of grains
 aspect_ratio = 4                # should be x times as wide as it is tall
 
 mkpath("simulation$(ngrains)")
diff --git a/layer_basin.jl b/layer_basin.jl
@@ -3,7 +3,7 @@ import JLD2
 import PyPlot
 import Dates
 
-id = "simulation1000"    # id of simulation to load, just write the folder
+id = "simulation500"    # id of simulation to load, just write the folder
                         # name here
 
 # Layer interface positions
@@ -103,7 +103,7 @@ for grain in sim.grains
     end
 end
 
-Granular.findContactsAllToAll!(carpet) # find the grain contacts
+Granular.findContactsAllToAll!(sim) # find the grain contacts
 
 #reduce the contact radius again
 for i = 1:size(sim.grains,1)