tAdd stress BC to top wall in shear function - sphere - GPU-based 3D discrete element method algorithm with optional fluid coupling
 (HTM) git clone git://src.adamsgaard.dk/sphere
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) LICENSE
       ---
 (DIR) commit d3de1fea51e48156510b03c78bec3b8a6f56626c
 (DIR) parent cf01772825c92e36d827573770ad4ca04aebf450
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Mon, 31 Mar 2014 12:36:54 +0200
       
       Add stress BC to top wall in shear function
       
       Diffstat:
         M python/sphere.py                    |       3 ++-
       
       1 file changed, 2 insertions(+), 1 deletion(-)
       ---
 (DIR) diff --git a/python/sphere.py b/python/sphere.py
       t@@ -2087,7 +2087,7 @@ class sim:
                # Adjust grid and placement of upper wall
                self.adjustUpperWall()
        
       -        # Set the top wall BC to a value of deviatoric stress
       +        # Set the top wall BC to a value of normal stress
                self.wmode = numpy.array([1])
                self.w_devs = numpy.ones(1) * normal_stress
        
       t@@ -2168,6 +2168,7 @@ class sim:
        
                # Adjust grid and placement of upper wall
                self.adjustUpperWall()
       +        self.wmode = numpy.array([1])
        
                # Fix horizontal velocity to 0.0 of lowermost particles
                d_max_below = numpy.max(self.radius[numpy.nonzero(self.x[:,2] <