tfix vtk files - 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 a46ed4ed9cb15cae01a5897578a35aeace71c035
DIR parent db378368459b79a351971445588510f7479c9455
HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Tue, 14 Oct 2014 13:49:48 +0200
fix vtk files
Diffstat:
M python/sphere.py | 13 ++++++++-----
1 file changed, 8 insertions(+), 5 deletions(-)
---
DIR diff --git a/python/sphere.py b/python/sphere.py
t@@ -5091,19 +5091,19 @@ class sim:
# First iteration: Allocate arrays and find constant values
if (i == 0):
# Shear displacement
- self.xdisp = numpy.zeros(lastfile+1, dtype=numpy.float64)
+ self.xdisp = numpy.zeros(lastfile+1, dtype=numpy.float64)
# Normal stress
- self.sigma_eff = numpy.zeros(lastfile+1, dtype=numpy.float64)
+ self.sigma_eff= numpy.zeros(lastfile+1, dtype=numpy.float64)
# Normal stress
- self.sigma_def = numpy.zeros(lastfile+1, dtype=numpy.float64)
+ self.sigma_def= numpy.zeros(lastfile+1, dtype=numpy.float64)
# Shear stress
- self.tau = numpy.zeros(lastfile+1, dtype=numpy.float64)
+ self.tau = numpy.zeros(lastfile+1, dtype=numpy.float64)
# Upper wall position
- self.dilation = numpy.zeros(lastfile+1, dtype=numpy.float64)
+ self.dilation = numpy.zeros(lastfile+1, dtype=numpy.float64)
# Upper wall position
self.tau_p = 0.0 # Peak shear stress
t@@ -5116,6 +5116,9 @@ class sim:
w_x0 = sb.w_x[0] # Original height
A = sb.L[0] * sb.L[1] # Upper surface area
+ if (i == 1):
+ w_x0 = sb.w_x[0] # Original height
+
# Summation of shear stress contributions
for j in fixvel[0]:
if (sb.vel[j,0] > 0.0):