tadd output when sim isn't found, increase height of fluid grid in diffusivity experiments - 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 20a203aa918a4e84b93380acc6f5effeb7475f16
   DIR parent 31815a035aca67067bbcf99b7f56fbf28546935e
  HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Fri, 29 Aug 2014 13:47:40 +0200
       
       add output when sim isn't found, increase height of fluid grid in diffusivity experiments
       
       Diffstat:
         M python/diffusivity-starter.py       |       6 +++++-
         M python/permeability-results.py      |       2 ++
       
       2 files changed, 7 insertions(+), 1 deletion(-)
       ---
   DIR diff --git a/python/diffusivity-starter.py b/python/diffusivity-starter.py
       t@@ -17,7 +17,7 @@ for sigma0_str in sys.argv[4:]:
            sim.readlast()
        
            sim.sid = 'diffusivity-sigma0=' + str(sigma0) + '-c_phi=' + str(c_phi) + \
       -            '-c_grad_p=' + str(c_grad_p)
       +            '-c_grad_p=' + str(c_grad_p) + '-tall'
            print(sim.sid)
        
            sim.checkerboardColors()
       t@@ -28,6 +28,10 @@ for sigma0_str in sys.argv[4:]:
        
            sim.consolidate(normal_stress = sigma0)
        
       +    # Increase height of grid
       +    sim.L[2] *= 2.0
       +    sim.num[2] *= 2
       +
            sim.initFluid(mu = 17.87e-4, p = 1.0e5, hydrostatic = True)
            sim.setFluidBottomNoFlow()
            sim.setFluidTopFixedPressure()
   DIR diff --git a/python/permeability-results.py b/python/permeability-results.py
       t@@ -47,6 +47,8 @@ for c_grad_p in cvals:
                    Q[c][i] = pc.Q[2]
                    pc.findMeanPorosity()
                    phi_bar[c][i] = pc.phi_bar
       +        else:
       +            print(sid + ' not found')
        
                i += 1