tdilation in mean grain diameters as default - 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 d4606c13b8efe33dfd5d4dc972f9221374aaef5e
 (DIR) parent 565130ebc1f63424aff2f72cf769fa7ed75aa7f3
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Mon, 31 Mar 2014 13:29:51 +0200
       
       dilation in mean grain diameters as default
       
       Diffstat:
         M python/sphere.py                    |      12 ++++++++++--
       
       1 file changed, 10 insertions(+), 2 deletions(-)
       ---
 (DIR) diff --git a/python/sphere.py b/python/sphere.py
       t@@ -4182,8 +4182,16 @@ class sim:
                            xdisp[i]    = xdisp[i-1] + sb.time_file_dt[0] * shearvel
                        sigma_eff[i] = sb.w_force[0] / A
                        sigma_def[i] = sb.w_devs[0]
       -                #dilation[i] = sb.w_x[0] - w_x0   # dilation in meters
       -                dilation[i] = (sb.w_x[0] - w_x0)/w_x0 * 100.0 # dilation in percent
       +
       +                # dilation in meters
       +                #dilation[i] = sb.w_x[0] - w_x0
       +
       +                # dilation in percent
       +                #dilation[i] = (sb.w_x[0] - w_x0)/w_x0 * 100.0 # dilation in percent
       +
       +                # dilation in number of mean particle diameters
       +                d_bar = numpy.mean(self.radius)*2.0
       +                dilation[i] = (sb.w_x[0] - w_x0)/d_var
        
                        # Test if this was the max. shear stress
                        if (tau[i] > tau_u):