trefer to correct velocity array - 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 18c620d69ca54e3b0de0faf91314ea164f0f0d2c
 (DIR) parent cb034e35510563d64958f7a01f84ea6c5f932e90
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Tue, 10 Feb 2015 20:52:08 +0100
       
       refer to correct velocity array
       
       Diffstat:
         M python/sphere.py                    |       4 ++--
       
       1 file changed, 2 insertions(+), 2 deletions(-)
       ---
 (DIR) diff --git a/python/sphere.py b/python/sphere.py
       t@@ -5945,7 +5945,7 @@ class sim:
        
                        if (i > 0):
                            self.xdisp[i] = sb.xyzsum[fixvel,0].max()
       -                    self.v[i] = sb.vel[fixvel,0].max()
       +                    self.v[i]     = sb.vel[fixvel,0].max()
        
                        self.sigma_eff[i] = sb.w_force[0]/A
                        self.sigma_def[i] = sb.currentNormalStress()
       t@@ -6000,7 +6000,7 @@ class sim:
        
                        # Middle plot
                        ax5 = plt.subplot(3, 1, 2, sharex=ax1)
       -                ax5.semilogy(time[1:], self.vel[1:], label='Shear velocity')
       +                ax5.semilogy(time[1:], self.v[1:], label='Shear velocity')
                        ax5.set_ylabel('Shear velocity [ms$^{-1}$]')
        
                        # Lower plot