tread noshear file as input when c_grad_p = 0.1 and sigma0=20e3 - 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 53ebe460d5a88ab9eb12115197c0ca687da5c0f2
   DIR parent 66b4f94d6377edc0dbd6c2722bc3a7281e8a08e5
  HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Wed, 24 Sep 2014 10:09:12 +0200
       
       read noshear file as input when c_grad_p = 0.1 and sigma0=20e3
       
       Diffstat:
         M python/shear-starter.py             |      10 +++++++---
         M python/sphere.py                    |       4 ++--
       
       2 files changed, 9 insertions(+), 5 deletions(-)
       ---
   DIR diff --git a/python/shear-starter.py b/python/shear-starter.py
       t@@ -25,9 +25,13 @@ sim = sphere.sim('diffusivity-sigma0=' + str(sigma0),
                fluid=True)
        sim.readlast()
        
       +if sigma0 == 20.0e3 and c_phi == 1.0 and c_grad_p == 0.1:
       +    sim.sid = 'shear-sigma0=20000.0-c_phi=1.0-c_grad_p=0.1-hi_mu-lo_visc-hw-noshear'
       +    sim.readlast()
       +
        if fluid:
            sim.sid = 'shear-sigma0=' + str(sigma0) + '-c_phi=' + str(c_phi) + \
       -            '-c_grad_p=' + str(c_grad_p) + '-hi_mu-lo_visc-hw-noshear'
       +            '-c_grad_p=' + str(c_grad_p) + '-hi_mu-lo_visc-hw'
        else:
            sim.sid = 'shear-sigma0=' + str(sigma0) + '-hw'
        
       t@@ -39,8 +43,8 @@ sim.cleanup()
        sim.adjustUpperWall()
        sim.zeroKinematics()
        
       -#sim.shear(1.0/20.0)
       -sim.shear(0.0)
       +sim.shear(1.0/20.0)
       +#sim.shear(0.0)
        
        if fluid:
            #sim.num[2] *= 2
   DIR diff --git a/python/sphere.py b/python/sphere.py
       t@@ -3659,7 +3659,7 @@ class sim:
                :type disp: str
                '''
        
       -        self.writebin()
       +        self.writebin(verbose=False)
        
                nd = ''
                if (disp == '2d'):
       t@@ -3668,7 +3668,7 @@ class sim:
                subprocess.call("cd .. && ./forcechains " + nd + "-f " + outformat \
                        + " -lc " + str(lc) + " -uc " + str(uc) + " input/" + self.sid \
                        + ".bin > python/tmp.gp", shell=True)
       -        subprocess.call("gnuplot tmp.gp && rm tmp.bin && rm tmp.gp", shell=True)
       +        subprocess.call("gnuplot tmp.gp && rm tmp.gp", shell=True)
        
        
            def forcechainsRose(self, lower_limit=0.25, graphics_format='pdf'):