tAdd more grains and change length of consolidation and relaxation - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
  HTML git clone git://src.adamsgaard.dk/sphere
   DIR Log
   DIR Files
   DIR Refs
   DIR LICENSE
       ---
   DIR commit 382c95a6e4f9f333dc93531ec68fa9024a8d2bc9
   DIR parent 4aa7a4c98c7c77cd9d743df021b9d9f22878dd59
  HTML Author: Anders Damsgaard <andersd@riseup.net>
       Date:   Mon, 22 Jan 2018 06:57:01 -0800
       
       Add more grains and change length of consolidation and relaxation
       
       Diffstat:
         M python/shear-test-ocr.py            |      11 ++++++-----
       
       1 file changed, 6 insertions(+), 5 deletions(-)
       ---
   DIR diff --git a/python/shear-test-ocr.py b/python/shear-test-ocr.py
       t@@ -4,20 +4,21 @@
        import sphere
        
        ### EXPERIMENT SETUP ###
       -initialization = False
       +initialization = True
        consolidation  = True
       +relaxation     = True
        shearing       = True
        rendering      = False
        plots          = True
        
        # Number of particles
       -np = 1e4
       +np = 2e4
        
        # Common simulation id
        sim_id = "shear-test-ocr"
        
        # Effective normal stresses during consolidation [Pa]
       -Nlist = [5e3, 10e3, 25e3, 50e3, 100e3, 250e3, 500e3]
       +Nlist = [10e3, 25e3, 50e3, 100e3, 250e3, 500e3]
        
        # Effective normal stresses during relaxation and shear [Pa]
        Nshear = 10e3
       t@@ -94,7 +95,7 @@ for N in Nlist:
            cons.checkerboardColors(nx=cons.num[0]/2, ny=cons.num[1]/2, nz=cons.num[2]/2)
        
            # Set duration of simulation
       -    cons.initTemporal(total = 1.5)
       +    cons.initTemporal(total = 4.0)
        
            if (consolidation == True):
        
       t@@ -130,7 +131,7 @@ for N in Nlist:
                    nz=relax.num[2]/2)
        
            # Set duration of simulation
       -    relax.initTemporal(total = 1.0)
       +    relax.initTemporal(total = 3.0)
        
            if (relaxation == True):