tplot force chains in a time series - 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 47e1a28dd303b7c5f4cf7e2a796e903c4e7d05b7
 (DIR) parent a8dfa931c6944978116c1c0a446925e9e09821e8
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Wed, 24 Sep 2014 10:55:37 +0200
       
       plot force chains in a time series
       
       Diffstat:
         A python/shear-results-fc.py          |      22 ++++++++++++++++++++++
       
       1 file changed, 22 insertions(+), 0 deletions(-)
       ---
 (DIR) diff --git a/python/shear-results-fc.py b/python/shear-results-fc.py
       t@@ -0,0 +1,22 @@
       +#!/usr/bin/env python
       +import sphere
       +import numpy
       +import subprocess
       +
       +sid = 'shear-sigma0=20000.0-hw'
       +imgformat = 'png'
       +sim = sphere.sim(sid, fluid=False)
       +subprocess.call('mkdir -p ' + sid + '-fc', shell=True)
       +
       +d = 0
       +for i in numpy.arange(210,300):
       +#for i in numpy.arange(240,260):
       +    print("File: %d, output: %s-fc/%05d.png" % (i, sid, d))
       +    sim.readstep(i, verbose=False)
       +    #sim.forcechains(lc = 3.0e1, uc=1.0e2)
       +    sim.forcechains(lc = 30.0, uc=1.0e2)
       +    subprocess.call('mv ' + sid + '-fc.' + imgformat \
       +            + ' ' + sid + '-fc/%05d.png' % (d), shell=True)
       +    d += 1
       +
       +subprocess.call('mv ' + sid + '-fc && sh ./make_video.sh', shell=True)